Search NASA⌕ Search

SEARCH · Search NASA

Results for “Navier Stokes simulation”

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 307 records · Page 17

Nonlinear Aeroelastic Analysis Using a Time-Accurate Navier-Stokes Equations Solver

A method to simulate limit cycle oscillation (LCO) due to control surface freeplay using a modified CFL3D, a time-accurate Navier-Stokes computational fluid dynamics (CFD) analysis code with structural modeling capability, is presented. This approach can be used to analyze aeroelastic response of aircraft with structural behavior characterized by nonlinearity in the force verses displacement curve. A limited validation of the method, using very low Mach number experimental data for a three-degrees-of-freedom (pitch/plunge/flap deflection) airfoil model with flap freeplay, is also presented.

Kuruvila, Geojoe↗

Wake Management Strategies for Reduction of Turbomachinery Fan Noise

The primary objective of our work was to evaluate and test several wake management schemes for the reduction of turbomachinery fan noise. Throughout the course of this work we relied on several tools. These include 1) Two-dimensional steady boundary-layer and wake analyses using MISES (a thin-shear layer Navier-Stokes code), 2) Two-dimensional unsteady wake-stator interaction simulations using UNSFLO, 3) Three-dimensional, steady Navier-Stokes rotor simulations using NEWT, 4) Internal blade passage design using quasi-one-dimensional passage flow models developed at MIT, 5) Acoustic modeling using LINSUB, 6) Acoustic modeling using VO72, 7) Experiments in a low-speed cascade wind-tunnel, and 8) ADP fan rig tests in the MIT Blowdown Compressor.

Waitz, Ian A.↗

Tip vortices of wings in subsonic and transonic flow: A numerical simulation

Thin layer Navier-Stokes and Euler equations are numerically solved using a multi-block zonal approach to simulate the formation and roll up of tip vortices of wings in subsonic and transonic flows. Several wing planforms were considered to examine the influence of tip-cap shape, planform geometry and free stream Mach number on the formation process. A good definition of the formation and qualitative roll up of tip vortices was achieved.

Srinivasan, G. R.↗

Effects of Chemistry on Blunt-Body Wake Structure

Results of a numerical study are presented for hypersonic low-density flow about a 70-deg blunt cone using direct simulation Monte Carlo (DSMC) and Navier-Stokes calculations. Particular emphasis is given to the effects of chemistry on the near-wake structure and on the surface quantities and the comparison of the DSMC results with the Navier-Stokes calculations. The flow conditions simulated are those experienced by a space vehicle at an altitude of 85 km and a velocity of 7 km/s during Earth entry. A steady vortex forms in the near wake for these freestream conditions for both chemically reactive and nonreactive air gas models. The size (axial length) of the vortex for the reactive air calculations is 25% larger than that of the nonreactive air calculations. The forebody surface quantities are less sensitive to the chemistry than the base surface quantities. The presence of the afterbody has no effect on the forebody flow structure or the surface quantities. The comparisons of DSMC and Navier-Stokes calculations show good agreement for the wake structure and the forebody surface quantities.

Dogra, Virendra K.↗

Time-Accurate Simulations and Acoustic Analysis of Slat Free-Shear Layer

A detailed computational aeroacoustic analysis of a high-lift flow field is performed. Time-accurate Reynolds Averaged Navier-Stokes (RANS) computations simulate the free shear layer that originates from the slat cusp. Both unforced and forced cases are studied. Preliminary results show that the shear layer is a good amplifier of disturbances in the low to mid-frequency range. The Ffowcs-Williams and Hawkings equation is solved to determine the acoustic field using the unsteady flow data from the RANS calculations. The noise radiated from the excited shear layer has a spectral shape qualitatively similar to that obtained from measurements in a corresponding experimental study of the high-lift system.

Khorrami, Mehdi R.↗

Study of Convective Flow Effects in Endwall Casing Treatments in Transonic Compressor Rotors

The unsteady convective flow effects in a transonic compressor rotor with a circumferential-groove casing treatment are investigated in this paper. Experimental results show that the circumferential-groove casing treatment increases the compressor stall margin by almost 50% for the current transonic compressor rotor. Steady flow simulation of the current casing treatment, however, yields only a 15% gain in stall margin. The flow field at near-stall operation is highly unsteady due to several self-induced flow phenomena. These include shock oscillation, vortex shedding at the trailing edge, and interaction between the passage shock and the tip clearance vortex. The primary focus of the current investigation is to assess the effects of flow unsteadiness and unsteady flow convection on the circumferential-groove casing treatment. Unsteady Reynolds-averaged Navier-Stokes (URANS) and Large Eddy Simulation (LES) techniques were applied in addition to steady Reynolds-averaged Navier-Stokes (RANS) to simulate the flow field at near-stall operation and to determine changes in stall margin. The current investigation reveals that unsteady flow effects are as important as steady flow effects on the performance of the circumferential grooves casing treatment in extending the stall margin of the current transonic compressor rotor. The primary unsteady flow mechanism is unsteady flow injection from the grooves into the main flow near the casing. Flows moving into and out of the grooves are caused due to local pressure difference near the grooves. As the pressure field becomes transient due to self-induced flow oscillation, flow injection from the grooves also becomes unsteady. The unsteady flow simulation shows that this unsteady flow injection from the grooves is substantial and contributes significantly to extending the compressor stall margin. Unsteady flows into and out of the grooves have as large a role as steady flows in the circumferential grooves. While the circumferential-groove casing treatment seems to be a steady flow device, unsteady flow effects should be included to accurately assess its performance as the flow is transient at near-stall operation.

Hah, Chunill↗

Calculation of a simulated 3-D high speed inlet using the Navier-Stokes equations

A hybrid numerical algorithm, developed to solve the full three-dimensional Navier-Stokes equations, is applied to the computation of the flowfield in a simulated three-dimensional high speed aircraft inlet at a Mach number of 2.5 and Reynolds number of 1.4 x 10 to the 7th based on inlet length. The numerical algorithm incorporates a coordinate transformation in order to handle general flow geometries, and utilizes the algebraic turbulent eddy viscosity model of Baldwin and Lomax. The hybrid algorithm has been vectorized on the CDC CYBER 203 computer using the SL/1 vector programming language developed at NASA Langley. The computed results are compared with experimental measurements of the ramp and cowl static pressures, and boundary layer pitot profiles. The results are also compared with a previous two-dimensional Navier-Stokes computation of the same configuration. The agreement with the experimental data is generally good; however, additional improvements in turbulence modeling are needed.

Knight, D. D.↗

Reciprocal interactions of hairpin-shaped vortices and a boundary layer

The present effort simulates the flow dynamics resulting from hairpin-shaped vortices in a boundary layer through the solution of the time-dependent, three-dimensional, compressible Navier-Stokes equations. Two simulations have been carried out. In the first case, the initial condition contains only one imposed hairpin-shaped vortex. In the second case, there are two incipient vortices separated by a short distance in the streamwise direction. The results add to the detailed understanding of the dynamic role of hairpin vortices in the entrainment, production and continual replenishment of turbulent structure, as well as the process of the interaction, intertwining and eventual conglomeration of vortices in boundary layers.

Liu, N.-S.↗

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗

Application of Lattice Boltzmann and Navier-Stokes Methods to NASA's Wall Mounted Hump

Lattice Boltzmann (LB) based Large Eddy Simulation (LES), Reynolds-averaged Navier-Stokes (RANS) as well as hybrid RANS/LES methods within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are applied to NASA's wall-mounted hump. Computational results are compared with experiments performed by Greenblatt et al. A detailed comparison between the accuracy and resolution requirements of the two approaches for turbulence resolving simulations, as well as the suitability of different grid paradigms (body-fitted curvilinear and block structured Cartesian) are presented. This test case is part of NASA's Revolutionary Computational Aerosciences (RCA) sub-project which addresses the technical challenge of predicting flow separation and reattachment accurately. Improvements in predictive accuracy by as much as 90% are demonstrated using LB as well as hybrid RANS/LES approaches compared to state-of-the-art steady state RANS simulations.

Lattice Boltzmann↗

Numerical simulation of particle-wave interaction in boundary layers

The effects of wall injection and particle motion on the spatial stability of two-dimensional plane channel flow are investigated. For this purpose, an accurate Navier-Stokes solver to simulate the space-time evolution of disturbances in three-dimensional flows has been developed. The code is operational on the NASA Langley CRAY2 and can be ported to any other supercomputer. The code has been tested extensively in tracking the spatial evolution of two-dimensional disturbances in plane channel flow and provided excellent agreement with the linear theory including at the inflow/outflow boundaries. Preliminary calculations have been performed to investigate the effects of stationary and moving sources of vortical disturbances simulating a particle traveling in the flow field. Results suggest that even at very low amplitudes, vortical disturbances act as amplifiers on the Tollmien-Schlichting waves promoting rapid instability. It is also found that slow moving particles are more dangerous than both stationary and fast moving particles for the same disturbance levels.

Biringen, S.↗

Large Eddy Simulation of High-Speed, Premixed Ethylene Combustion

A large-eddy simulation / Reynolds-averaged Navier-Stokes (LES/RANS) methodology is used to simulate premixed ethylene-air combustion in a model scramjet designed for dual mode operation and equipped with a cavity for flameholding. A 22-species reduced mechanism for ethylene-air combustion is employed, and the calculations are performed on a mesh containing 93 million cells. Fuel plumes injected at the isolator entrance are processed by the isolator shock train, yielding a premixed fuel-air mixture at an equivalence ratio of 0.42 at the cavity entrance plane. A premixed flame is anchored within the cavity and propagates toward the opposite wall. Near complete combustion of ethylene is obtained. The combustor is highly dynamic, exhibiting a large-scale oscillation in global heat release and mass flow rate with a period of about 2.8 ms. Maximum heat release occurs when the flame front reaches its most downstream extent, as the flame surface area is larger. Minimum heat release is associated with flame propagation toward the cavity and occurs through a reduction in core flow velocity that is correlated with an upstream movement of the shock train. Reasonable agreement between simulation results and available wall pressure, particle image velocimetry, and OH-PLIF data is obtained, but it is not yet clear whether the system-level oscillations seen in the calculations are actually present in the experiment.

Ramesh, Kiran↗

Aeroacoustic Simulations of a Nose Landing Gear Using FUN3D on Pointwise Unstructured Grids

Numerical simulations have been performed for a partially-dressed, cavity-closed (PDCC) nose landing gear configuration that was tested in the University of Florida's open-jet acoustic facility known as the UFAFF. The unstructured-grid flow solver FUN3D is used to compute the unsteady flow field for this configuration. Mixed-element grids generated using the Pointwise(TradeMark) grid generation software are used for these simulations. Particular care is taken to ensure quality cells and proper resolution in critical areas of interest in an effort to minimize errors introduced by numerical artifacts. A hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) turbulence model is used for these simulations. Solutions are also presented for a wall function model coupled to the standard turbulence model. Time-averaged and instantaneous solutions obtained on these Pointwise grids are compared with the measured data and previous numerical solutions. The resulting CFD solutions are used as input to a Ffowcs Williams-Hawkings noise propagation code to compute the farfield noise levels in the flyover and sideline directions. The computed noise levels compare well with previous CFD solutions and experimental data.

Vatsa, Veer N.↗

Proper Orthogonal Decomposition in Optimal Control of Fluids

In this article, we present a reduced order modeling approach suitable for active control of fluid dynamical systems based on proper orthogonal decomposition (POD). The rationale behind the reduced order modeling is that numerical simulation of Navier-Stokes equations is still too costly for the purpose of optimization and control of unsteady flows. We examine the possibility of obtaining reduced order models that reduce computational complexity associated with the Navier-Stokes equations while capturing the essential dynamics by using the POD. The POD allows extraction of certain optimal set of basis functions, perhaps few, from a computational or experimental data-base through an eigenvalue analysis. The solution is then obtained as a linear combination of these optimal set of basis functions by means of Galerkin projection. This makes it attractive for optimal control and estimation of systems governed by partial differential equations. We here use it in active control of fluid flows governed by the Navier-Stokes equations. We show that the resulting reduced order model can be very efficient for the computations of optimization and control problems in unsteady flows. Finally, implementational issues and numerical experiments are presented for simulations and optimal control of fluid flow through channels.

Ravindran, S. S.↗

Down-to-earth supercomputing at NASA Ames

The requirements for computers for aerospace vehicle design problems, in particular fluid dynamics, are examined. The levels of approximation to Navier-Stokes equations are: (1) linearized inviscid form, (2) nonlinear inviscid equation, (3) Reynolds-averaged Navier-Stokes equation, (4) large eddy simulation with models for subgrid scale turbulence, and (5) full Navier-Stokes equations to calculate laminar/turbulent transition and turbulence dissipation. The Cray 2 with 256 megaword memory and 4.1 nsec clock and the Cray X-MP/48 with a 128 megaword memory and a 9.5 nsec clock of the Ames facility are described, and an example revealing their applicability is presented. The use of zonal methods to simulate aerodynamic flow fields numerically is discussed. Consideration is given to the development of a program for simulation of aerodynamic flow fields, the compilers, symbolic debuggers, the scheduling of multiple CPUs, and the analysis and storage of the data.

Stevens, Kenneth G., Jr.↗

Isolating Curvature Effects in Computing Wall-Bounded Turbulent Flows

The flow over the zero-pressure-gradient So-Mellor convex curved wall is simulated using the Navier-Stokes equations. An inviscid effective outer wall shape, undocumented in the experiment, is obtained by using an adjoint optimization method with the desired pressure distribution on the inner wall as the cost function. Using this wall shape with a Navier-Stokes method, the abilities of various turbulence models to simulate the effects of curvature without the complicating factor of streamwise pressure gradient can be evaluated. The one-equation Spalart-Allmaras turbulence model overpredicts eddy viscosity, and its boundary layer profiles are too full. A curvature-corrected version of this model improves results, which are sensitive to the choice of a particular constant. An explicit algebraic stress model does a reasonable job predicting this flow field. However, results can be slightly improved by modifying the assumption on anisotropy equilibrium in the model's derivation. The resulting curvature-corrected explicit algebraic stress model possesses no heuristic functions or additional constants. It lowers slightly the computed skin friction coefficient and the turbulent stress levels for this case (in better agreement with experiment), but the effect on computed velocity profiles is very small.

Rumsey, Christopher L.↗

Real gas effects on the numerical simulation of a hypersonic inlet

The present use of a parabolized Navier-Stokes solver to accurately simulate the flowfield in a supersonic inlet yields good agreement between numerical analysis and experiment for a Mach 7.4 inlet under cruise conditions, with an internal compression ratio of 8. The significance of real gas effects on the performance calculation of a hypersonic inlet is demonstrated, with small changes in the ratio of specific heats resulting in a substantial change in the calculated pitot pressure ratio.

Ng, W.-F.↗

Aerothermodynamic Analyses of Towed Ballutes

A ballute (balloon-parachute) is an inflatable, aerodynamic drag device for application to planetary entry vehicles. Two challenging aspects of aerothermal simulation of towed ballutes are considered. The first challenge, simulation of a complete system including inflatable tethers and a trailing toroidal ballute, is addressed using the unstructured-grid, Navier-Stokes solver FUN3D. Auxiliary simulations of a semi-infinite cylinder using the rarefied flow, Direct Simulation Monte Carlo solver, DSV2, provide additional insight into limiting behavior of the aerothermal environment around tethers directly exposed to the free stream. Simulations reveal pressures higher than stagnation and corresponding large heating rates on the tether as it emerges from the spacecraft base flow and passes through the spacecraft bow shock. The footprint of the tether shock on the toroidal ballute is also subject to heating amplification. Design options to accommodate or reduce these environments are discussed. The second challenge addresses time-accurate simulation to detect the onset of unsteady flow interactions as a function of geometry and Reynolds number. Video of unsteady interactions measured in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel and CFD simulations using the structured grid, Navier-Stokes solver LAURA are compared for flow over a rigid spacecraft-sting-toroid system. The experimental data provides qualitative information on the amplitude and onset of unsteady motion which is captured in the numerical simulations. The presence of severe unsteady fluid - structure interactions is undesirable and numerical simulation must be able to predict the onset of such motion.

Gnoffo, Peter A.↗