Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stokes parameters”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Advanced Low-Noise Research Fan Stage Design

This report describes the design of the Advanced Low-Noise Research Fan stage. The fan is a variable pitch design, which is designed at the cruise pitch condition. Relative to the cruise setting, the blade is closed at takeoff and opened for reverse thrust operation. The fan stage is a split flow design with fan exit guide vanes (FEGVs) and core stators. The fan stage design is combined with a nacelle and engine core duct to form a powered fan/nacelle subscale model. This model is intended for use in combined aerodynamic, acoustic, and structural testing in a wind tunnel. The fan has an outer diameter of 22 in. and a hub-to-tip of 0.426 in., which allows the use of existing NASA fan and cowl force balance and rig drive systems. The design parameters were selected to permit valid acoustic and aerodynamic comparisons with the Pratt & Whitney (P&W) 17- and 22-in. rigs previously tested under NASA contract. The fan stage design is described in detail. The results of the design axisymmetric and Navier-Stokes aerodynamic analysis are presented at the critical design conditions. The structural analysis of the fan rotor and attachment is included. The blade and attachment are predicted to have adequate low-cycle fatigue life and an acceptable operating range without resonant stress or flutter. The stage was acoustically designed with airfoil counts in the FEGV and core stator to minimize noise. A fan/FEGV tone analysis developed separately under NASA contract was used to determine the optimum airfoil counts. The fan stage was matched to the existing nacelle, designed under the previous P&W low-noise contract, to form a fan/nacelle model for wind tunnel testing. It is an axisymmetric nacelle for convenience in testing and analysis. Previous testing confirmed that the nacelle performed as required at various aircraft operating conditions.

Neubert, Robert↗
DOE OSTI2023

Two self-similar Reynolds-stress transport models with anisotropic eddy viscosity

Two Reynolds-averaged Navier-Stokes models with full Reynolds-stress transport (RST) and tensor eddy viscosity are presented. These new models represent RST extensions of the $k−2L−a−\mathscr{C}$ and $k−ϕ−L−a−\mathscr{C}$ models by Morgan. Self-similarity analysis is used to derive constraints on model coefficients required to reproduce expected growth parameters for a variety of canonical flows, including Rayleigh-Taylor (RT) and Kelvin-Helmholtz (KH) mixing layers. Both models are then applied in one-dimensional simulation of RT and KH mixing layers, and the expected self-similar growth rates and anisotropy are obtained. Next, models are applied in two-dimensional simulation of the so-called “tilted rocket rig” inclined RT experiment and in simulation of a shock-accelerated localized patch of turbulence. Here it is found that RST is required to capture the qualitative growth of the shock-accelerated patch, and an anisotropic eddy viscosity provides substantial improvement over a Boussinesq treatment for the tilted rocket rig problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Soot-particle core-shell and fractal structures from small-angle X-ray scattering measurements in a flame

We have characterized soot particles measured in situ in a laminar co-flow ethylene-air diffusion flame using small-angle X-ray scattering (SAXS). The analysis includes temperature measurements made with coherent anti-Stokes Raman spectroscopy (CARS) and complements soot volume-fraction and maturity measurements made with laser-induced incandescence (LII). We compared the results of fits to the SAXS measurements using a unified model and a fractal core-shell model. Power-law parameters yielded by the unified model indicate that aggregates of primary particles are in the mass-fractal regime, whereas the primary particles are in the surface-fractal regime in the middle of the flame. Higher and lower in the flame, the primary-particle power-law parameter approaches 4, suggesting smooth primary particles. Furthermore, these trends are consistent with fits using the fractal core-shell model, which indicate that particles have an established core-shell structure in the middle of the flame and are internally homogeneous at higher and lower heights in the flame. Primary-particle size distributions derived using the fractal core-shell model demonstrate excellent agreement with distributions inferred from transmission electron microscopy (TEM) images in the middle of the flame. Higher in the flame, a second small mode appears in the size distributions, suggesting particle fragmentation during oxidation. Surface oxidation would explain (1) aggregate fragmentation and (2) loss of core-shell structure leading to smoother primary-particle surfaces by removal of carbon overlayers. SAXS measurements are much more sensitive to incipient and young soot particles than LII and demonstrate significant volume fraction from particles low in the flame where the LII signal is negligible.

42 ENGINEERING↗
NASA NTRS1980

A Navier-Stokes fast solver for turbulence modeling applications

A computer code for the evaluation and/or optimization of the predicative potential of second order turbulent closure models in simple two dimensional flow configurations is discussed. A procedure for the numerical solution of the steady constant property Navier-Stokes equations are described together with algebraic, one dimensional and two dimensional equations of turbulence closure models. Four turbulence models are compared with several sets of experimental data. The effects of initial conditions and boundary conditions are also described. The effects of purely numerical parameters, such as mesh size, boundary locations, and convergence criteria are presented.

Murphy, J. D.↗
NASA NTRS1981

Basic lubrication equations

Lubricants, usually Newtonian fluids, are assumed to experience laminar flow. The basic equations used to describe the flow are the Navier-Stokes equation of motion. The study of hydrodynamic lubrication is, from a mathematical standpoint, the application of a reduced form of these Navier-Stokes equations in association with the continuity equation. The Reynolds equation can also be derived from first principles, provided of course that the same basic assumptions are adopted in each case. Both methods are used in deriving the Reynolds equation, and the assumptions inherent in reducing the Navier-Stokes equations are specified. Because the Reynolds equation contains viscosity and density terms and these properties depend on temperature and pressure, it is often necessary to couple the Reynolds with energy equation. The lubricant properties and the energy equation are presented. Film thickness, a parameter of the Reynolds equation, is a function of the elastic behavior of the bearing surface. The governing elasticity equation is therefore presented.

Hamrock, B. J.↗
NASA NTRS1984

Comparison of experimental and computational compressible flow in a S-duct

This paper describes experimental measurements of secondary flow in a constant area, circular cross-section 30-30 deg S-duct, and compares the results obtained with the computations performed using the PEPSIG code, a parabolized Navier-Stokes code. The flow entering the duct was turbulent, with entrance Mach number of 0.6, and the boundary layer thickness at the duct entrance was 10 percent of the duct diameter. The duct mean radius of curvature to the duct diameter was 5.077. Flow parameters were measured at six stations along the length of the duct. These measurements were made using a five-port cone probe. At least ten radial traverses were made at each station on both sides of the symmetry plane. Wall static pressures along three azimuth angles of zero, 90, and 180 deg along the duct were measured. Plots presenting the secondary velocity field as well as contour plots of the total and static-pressure fields have been obtained. Strong secondary flows were observed in the first bend, and these continued into the second bend with the formation of new vorticity in the opposite sense in the second bend. The flow exiting the duct contained two pairs of counter-rotating vortices. The computational results are in general agreement with the experiments. However, it appears that the computations underestimate the extent of the pressure distortion, due to simplifications made in the pressure field calculations.

Vakili, A.↗
NASA NTRS2016

Coaxial Rotor Flow Phenomena in Forward Flight

Coaxial rotors are nding use in advanced rotorcraft concepts. Combined with lift offset rotor technology, they offer a solution to the problems of dynamic stall and re-verse ow that often limit single rotor edgewise forward ight speeds. Lower tip speed means reduced high-speed impulsive noise. The need for an anti-torque tail rotor is eliminated, a major boon during operation in conned areas. However, the operation of two counter-rotating ro-tors in close proximity generates many possibilities for aerodynamic interactions between rotor blades, blades and vortices, and between vortices. The parameter de-sign space is very large, and requires efcient computations as well as basic experiments to explore important physics to determine performance, loads, and acoustics. Computations are done on the classic HarringtonDingeldein rotor test case from the 1950s using the ROTUNS Navier Stokes code. Two regimes are explored: very low advance ratio as a perturbation from hover, and high advance ratio. Flow eld properties from RotUNS are used with 2-D OVERFLOW computations to capture blade crossing effects including those of higher subsonic Mach numbers. Bladeblade and bladevortex intersection events are captured using a MatLab-based predictor.

Flow Phenomena↗
DOE OSTI2023

Scaled ILU Smoothers for Navier-Stokes Pressure Projection

Incomplete LU (ILU) smoothers are effective in the algebraic multigrid (AMG) V-cycle for reducing high-frequency components of the error. However, the requisite direct triangular solves are comparatively slow on GPUs. Previous work has demonstrated the advantages of Jacobi iteration as an alternative to direct solution of these systems. Depending on the threshold and fill-level parameters chosen, the factors can be highly nonnormal and Jacobi is unlikely to converge in a low number of iterations. We demonstrate that row scaling can reduce the departure from normality, allowing us to replace the inherently sequential solve with a rapidly converging Richardson iteration. There are several advantages beyond the lower compute time. Scaling is performed locally for a diagonal block of the global matrix because it is applied directly to the factor. Further, an ILUT Schur complement smoother maintains a constant GMRES iteration count as the number of MPI ranks increases, and thus parallel strong-scaling is improved. Our algorithms have been incorporated into hypre, and we demonstrate improved time to solution for linear systems arising in the Nalu-Wind and PeleLM pressure solvers. For large problem sizes, GMRES+AMG executes at least five times faster when using iterative triangular solves compared with direct solves on massively parallel GPUs.

algebraic multigrid↗
NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗
DOE OSTI2023

The myth of URANS

Since the 1990s, RANS practitioners have observed spontaneous unsteadiness in RANS simulations. Some have suggested deliberately using this as a method of resolving large turbulent structures. However, to date, no one has produced a theoretical justification for this unsteady RANS (URANS) approach. Here, we extend the dynamical system fixed point analysis to create a theoretical model for URANS dynamics. The results are compared to URANS simulations for homogeneous isotropic decaying turbulence. The model shows that URANS can predict incorrect decay rates and that the solution tends towards steady RANS over time. Similar analysis for forced turbulence shows a fixed modelled energy of about 30% of total energy, regardless of the model parameters. The same analysis can be used to show how hybrid type models can begin to address these issues.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

Details of the computed flowfield over a circular cylinder at Reynolds number 1200

The application of an upwind-biased implicit approximate factorization Navier-Stokes algorithm to the unsteady impulsive start-up flow over a circular cylinder at Reynolds number 1200 is described. The complete form of the compressible Navier-Stokes equations is used, and the algorithm is second-order accurate in both space and time. The development with time of the shape and size of the separated vortical flow region is computed, as well as the time-variation of several boundary layer parameters and profile shapes. Computations, in general, show excellent agreement with experiment, although the present method predicts a more rapid onset of reversed flow on the cylinder than evidenced in experiment. The changes that the vortical region behind the cylinder undergoes as the symmetric flow transitions to periodic vortex shedding are discussed. The flow becomes periodic with a Strouhal frequency of 0.222, which compares well with the experimental value of approximately 0.21. The effect of grid density on the development of the unsteady flow is also shown.

Rumsey, C. L.↗
NASA NTRS2018

Aerodynamics and Aeroacoustic Sources of a Coaxial Rotor

Vehicles with coaxial, contra-rotating rotor systems (CACR) are being considered for a range of applications. Design parameters are identified that effect the fundamental aerodynamics and fluid dynamic features of a CACR in hover, vertical, and edgewise flight. Understanding the fluid dynamic features is a precursor to studying the aeroacoustics of a coaxial rotor. Rotor performance was computed initially using Navier-Stokes solver with prescribed blade section aerodynamic properties, the results validated against generic experimental test cases. Using a A 2D potential flow code and 2D OVERFLOW compressible-flow Navier-Stokes solver, the fluid dynamics of blade interactions was simplified and broken into a 2D blade crossing problem, with crossing locations and velocity fields from the rotor results to understand circulation, thickness, compressibility, shed vorticity, downwash, and viscosity effects. A calculation tool has been developed to identify time and location of blade overlap and BVI time location for CACR. Specific aerodynamic phenomena that occur for each noise source relevant to CACR are presented, along with computational tools to predict these occurrences.

Aeroacoustic↗
DOE OSTI2023

Multifidelity Approach to Sensitivity Estimation in Large-Eddy Simulation

A novel approach to compute affordable approximate sensitivities in a large-eddy simulation (LES) is proposed and assessed. The approach is based on solving a Reynolds-averaged Navier–Stokes (RANS) problem that has been linearized around the mean LES solution, with closure modeling required for the linearized changes in turbulent Reynolds stresses. In the present study, the closure modeling is based on a linearized algebraic turbulence model. The method is assessed for the flow over a NACA0012 airfoil at a fixed angle of attack, with the Reynolds number as the variable parameter. The results show that, provided an accurate linearized closure model, the method predicts the correct sensitivity of the skin friction coefficient and the mean velocity field at the cost of a linearized RANS, which provides an important proof-of-concept for this approach. Here, the linearized algebraic turbulence model with standard model coefficients produces a reasonably accurate sensitivity, but the results also suggest that accuracy could be gained from recalibrating the model coefficients for this new use of the model.

42 ENGINEERING↗
DOE OSTI2022

Accurate and Efficient Parametric Model-Order Reduction for Turbulent Thermal Transport

This project produced new algorithms and software for low-cost exploration of turbulent thermal-fluids behavior under parametric variation by using reduced-order models (ROMs). The ROMs numerically solve the governing equations for fluid motion by using a small set of basis functions (typically, N=20-200 modes) to represent the solution. The base modes are computed as optimal combinations of solutions from expensive high-fidelity "anchor-point" solutions involving millions of unknowns, which are typically generated by solving the full Navier-Stokes equations on a supercomputer. The ROM solution is itself a combination of the base modes, where the basis coefficients are determined by evolving an NxN system of nonlinear equations. The overarching idea is to use the inexpensive ROM to predict solutions under conditions where the parameters differ from the anchor-point conditions. Several ingredients are required to make ROMs useful for thermal hydraulics analysis. These include: a stable and accurate ROM that is capable of reproducing the large-scale dynamics of turbulent flow, error indicators than can guide the choice of anchor points, and low-cost mechanisms for evaluating nonlinear terms in the reduced equations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1974

Computer extension and analytic continuation of Stokes' expansion for gravity waves

Stokes' infinitesimal-wave expansion for steady progressive free-surface waves has been extended to high order using a computer to perform the coefficient arithmetic. Stokes' expansion has been found to be incapable of yielding the highest wave for any value of the water depth since convergence is limited by a square-root branch-point some distance short of the maximum. By reformulating the problem using a different independent parameter, the highest waves are obtained correctly. Series summation and analytic continuation are facilitated by the use of Pade approximants. The method is valid in principle for any finite value of the wavelength and solutions of high accuracy can be obtained for most values of the wave height and water depth.

Schwartz, L. W.↗
NASA NTRS1993

2D axisymmetric analysis of SRM ignition transient

To analyze ignition transient of Space Shuttle solid rocket motor, a transient two-dimensional numerical model based on turbulent compressible Navier-Stokes equations in a generalized coordinate system was developed. One-dimensional numerical models (Peretz et al., 1973; Han, 1992; Pardue and Han, 1992) with empirical correlations data obtained from steady turbulent boundary layer flows agrees reasonably well with test rocket data by adjusting a few parameters. However, a 1D model can not provide a physical insight into the complex multidimensional thermal fields and flowfields in the chamber and the converging-diverging rocket nozzle. As an interim step, a 2D model was developed and compared with test data. A modified version of SIMPLE algorithm was used for the numerical model, and the standard k-epsilon model with a wall function was used for turbulence closure. Transient flowfields and thermal fields in the combustion chamber and the attached nozzle were obtained for a selected rocket geometry and propellant. Transient behaviors of the flow and thermal fields were analyzed, and were found to be in good agreement with physical expectations.

Bai, S. D.↗
NASA NTRS1990

Comparisons of a Three-Dimensional, Full Navier Stokes Computer Model with High Mach Number Combuster Test Data

Comparisons between scramjet combustor data and a three-dimensional full Navier-Stokes calculation have been made to verify and substantiate computational fluid dynamics (CFD) codes and application procedures. High Mach number scramjet combustor development will rely heavily on CFD applications to provide wind tunnel-equivalent data of quality sufficient to design, build and fly hypersonic aircraft. Therefore. detailed comparisons between CFD results and test data are imperative. An experimental case is presented, for which combustor wall static pressures were measured and flow-fieid interferograms were obtained. A computer model was done of the experiment, and counterpart parameters are compared with experiment. The experiment involved a subscale combustor designed and fabricated for the National Aero-Space Plane Program, and tested in the Calspan Corporation 96" hypersonic shock tunnel. The combustor inlet ramp was inclined at a 20 angle to the shock tunnel nozzle axis, and resulting combustor entrance flow conditions simulated freestream M=10. The combustor body and cowl walls were instrumented with static pressure transducers, and the combustor lateral walls contained windows through which flowfield holographic interferograms were obtained. The CFD calculation involved a three-dimensional time-averaged full Navier-Stokes code applied to the axial flow segment containing fuel injection and combustion. The full Navier-Stokes approach allowed for mixed supersonic and subsonic flow, downstream-upstream communication in subsonic flow regions, and effects of adverse pressure gradients. The code included hydrogen-air chemistry in the combustor segment which begins near fuel injection and continues through combustor exhaust. Combustor ramp and inlet segments on the combustor lateral centerline were modelled as two dimensional. Comparisons to be shown include calculated versus measured wall static pressures as functions of axial flow coordinate, and calculated path-averaged density contours versus an holographic Interferogram.

Watkins, William B.↗
NASA NTRS1989

"Smart" Electromechanical Shock Absorber

Shock-absorbing apparatus includes electromechanical actuator and digital feedback control circuitry rather than springs and hydraulic damping as in conventional shock absorbers. Device not subject to leakage and requires little or no maintenance. Attenuator parameters adjusted in response to sensory feedback and predictive algorithms to obtain desired damping characteristic. Device programmed to decelerate slowly approaching vehicle or other large object according to prescribed damping characteristic.

Stokes, Lebarian↗