Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Navier Stokes”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Wall-Modeled Lattice Boltzmann and Navier-Stokes Approaches for Separated Flows

Lattice Boltzmann (LB) and hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) methods within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are applied to NASA's Revolutionary Computational Aerosciences (RCA) standard test cases for separated flows. A detailed comparison between the performance and accuracy of the two emerging numerical methodologies for turbulence resolving simulations, i.e. the LB and hybrid RANS/LES methods will be presented. This contribution addresses the RCA technical challenge to identify and down-select critical turbulence, transition, and numerical method technologies for 40% reduction in predictive error for standard turbulence separated flow test cases. Results for the 2D NASA wall-mounted hump and the axisymmetric transonic bump including time-averaged pressure coefficient, skin friction, and velocity pro les, as well as resolved and modeled Reynolds stresses for both numerical approaches will be presented and differences between LB and hybrid RANS/LES will be discussed.

Lattice Boltzmann↗
NASA NTRS1991

Transonic Navier-Stokes solutions about a generic hypersonic configuration

Three-dimensional transonic viscous flow computations are presented for a generic high-speed accelerator model that includes wing, body, fillets, and a no-flow-through engine nacelle. Solutions are obtained from an algorithm for the compressible Navier-Stokes equations that incorporates an upwind-biased, flux-vector-splitting approach along with longitudinally patched grids. Results are presented for fully turbulent flow assumptions and include correlations with wind-tunnel data. A good quantitative agreement for the forebody surface pressure distribution is achieved between computations and the available wind-tunnel measurements at M∞ = 0.9. Furthermore, it is demonstrated that the flow is stagnating around the boattail region due to separation from the aft-engine cowl lip.

Hypersonic flows↗
NASA NTRS1987

Evaluation of Navier-Stokes and Euler solutions for leading-edge separation vortices

Extensive study on the numerical simulation of the vortical flow over a double delta wing is carried out using the thin layer Navier-Stokes and Euler equations. Two important flow characteristics, vortex interaction and vortex breakdown, are successfully simulated. Grid resolution is one of the most important factors associated with the vortex problem. Computations were performed on a series of grids with various levels of refinement, coarse, medium, and fine. Computations using either the coarse or medium grids fail to capture the proper physical phenomena. The computed result using a fine grid shows flow unsteadiness once the vortex breakdown takes place. The C sub L - alpha characteristics are well predicted up to the breakdown angle of attack for all the grid distributions. The Euler solutions show fairly good agreement with the experiment on the C sub L - alpha characteristics. However, other aspects of the solution at each angle of attack, such as the locus of the leading edge separation vortex, are not consistent with the experiment. Even for the fine grid Navier-Stokes computations, further grid resolution is required to obtain good quantitative agreement with the experiment.

Fujii, K.↗
NASA NTRS1987

Evaluation of Navier-Stokes and Euler solutions for leading-edge separation vortices

Extensive study on the numerical simulation of the vortical flow over a double delta wing is carried out using the thin layer Navier-Stokes and Euler equations. Two important flow characteristics, vortex interaction and vortex breakdown, are successfully simulated. Grid resolution is one of the most important factors associated with the vortex problem. Computations were performed on a series of grids with various levels of refinement, coarse, medium, and fine. Computations using either the coarse or medium grids fail to capture the proper physical phenomena. The computed result using a fine grid shows flow unsteadiness once the vortex breakdown takes place. The C sub L - alpha characteristics are well predicted up to the breakdown angle of attack for all the grid distributions. The Euler solutions show fairly good agreement with the experiment on the C sub L - alpha characteristics. However, other aspects of the solution at each angle of attack, such as the locus of the leading edge separation vortex, are not consistent with the experiment. Even for the fine grid Navier-Stokes computations, further grid resolution is required to obtain good quantitative agreement with the experiment.

Fujii, K.↗
NASA NTRS1988

A diagonally inverted LU implicit multigrid scheme for the 3-D Navier-Stokes equations and a two equation model of turbulence

An LU implicit multigrid algorithm is developed to calculate 3-D compressible viscous flows. This scheme solves the full 3-D Reynolds-Averaged Navier-Stokes equation with a two-equation kappa-epsilon model of turbulence. The flow equations are integrated by an efficient, diagonally inverted, LU implicit multigrid scheme while the kappa-epsilon equations are solved, uncoupled from the flow equations, by a block LU implicit algorithm. The flow equations are solved within the framework of the multigrid method using a four-grid level W-cycle, while the kappa-epsilon equations are iterated only on the finest grid. This treatment of the Reynolds-Averaged Navier-Stokes equations proves to be an efficient method for calculating 3-D compressible viscous flows.

Yokota, Jeffrey W.↗
NASA NTRS1989

A diagonally inverted LU implicit multigrid scheme for the 3-D Navier-Stokes equations and a two equation model of turbulence

An LU implicit multigrid algorithm is developed to calculate 3-D compressible viscous flows. This scheme solves the full 3-D Reynolds-Averaged Navier-Stokes equation with a two-equation kappa-epsilon model of turbulence. The flow equations are integrated by an efficient, diagonally inverted, LU implicit multigrid scheme while the kappa-epsilon equations are solved, uncoupled from the flow equations, by a block LU implicit algorithm. The flow equations are solved within the framework of the multigrid method using a four-grid level W-cycle, while the kappa-epsilon equations are iterated only on the finest grid. This treatment of the Reynolds-Averaged Navier-Stokes equations proves to be an efficient method for calculating 3-D compressible viscous flows.

Yokota, Jeffrey W.↗
NASA NTRS1995

Parametrics on 2D Navier-Stokes analysis of a Mach 2.68 bifurcated rectangular mixed-compression inlet

The supersonic diffuser of a Mach 2.68 bifurcated, rectangular, mixed-compression inlet was analyzed using a two-dimensional (2D) Navier-Stokes flow solver. Parametric studies were performed on turbulence models, computational grids and bleed models. The computer flowfield was substantially different from the original inviscid design, due to interactions of shocks, boundary layers, and bleed. Good agreement with experimental data was obtained in many aspects. Many of the discrepancies were thought to originate primarily from 3D effects. Therefore, a balance should be struck between expending resources on a high fidelity 2D simulation, and the inherent limitations of 2D analysis. The solutions were fairly insensitive to turbulence models, grids and bleed models. Overall, the k-e turbulence model, and the bleed models based on unchoked bleed hole discharge coefficients or uniform velocity are recommended. The 2D Navier-Stokes methods appear to be a useful tool for the design and analysis of supersonic inlets, by providing a higher fidelity simulation of the inlet flowfield than inviscid methods, in a reasonable turnaround time.

Mizukami, M.↗
NASA NTRS1989

Computation of aircraft component flow fields at transonic Mach numbers using a three-dimensional Navier-Stokes algorithm

A computer analysis was developed for calculating steady (or unsteady) three-dimensional aircraft component flow fields. This algorithm, called ENS3D, can compute the flow field for the following configurations: diffuser duct/thrust nozzle, isolated wing, isolated fuselage, wing/fuselage with or without integrated inlet and exhaust, nacelle/inlet, nacelle (fuselage) afterbody/exhaust jet, complete transport engine installation, and multicomponent configurations using zonal grid generation technique. Solutions can be obtained for subsonic, transonic, or hypersonic freestream speeds. The algorithm can solve either the Euler equations for inviscid flow, the thin shear layer Navier-Stokes equations for viscous flow, or the full Navier-Stokes equations for viscous flow. The flow field solution is determined on a body-fitted computational grid. A fully-implicit alternating direction implicit method is employed for the solution of the finite difference equations. For viscous computations, either a two layer eddy-viscosity turbulence model or the k-epsilon two equation transport model can be used to achieve mathematical closure.

Shrewsbury, George D.↗
NASA NTRS2003

Evaluation of a Second-Order Accurate Navier-Stokes Code for Detached Eddy Simulation Past a Circular Cylinder

We evaluate the applicability of a production computational fluid dynamics code for conducting detached eddy simulation for unsteady flows. A second-order accurate Navier-Stokes code developed at NASA Langley Research Center, known as TLNS3D, is used for these simulations. We focus our attention on high Reynolds number flow (Re = 5 x 10(sup 4) - 1.4 x 10(sup 5)) past a circular cylinder to simulate flows with large-scale separations. We consider two types of flow situations: one in which the flow at the separation point is laminar, and the other in which the flow is already turbulent when it detaches from the surface of the cylinder. Solutions are presented for two- and three-dimensional calculations using both the unsteady Reynolds-averaged Navier-Stokes paradigm and the detached eddy simulation treatment. All calculations use the standard Spalart-Allmaras turbulence model as the base model.

Vatsa, Veer N.↗
NASA NTRS1986

Transonic Navier-Stokes wing solution using a zonal approach. Part 1: Solution methodology and code validation

A fast diagonalized Beam-Warming algorithm is coupled with a zonal approach to solve the three dimensional Euler/Navier-Stokes equations. The computer code, called Transonic Navier-Stokes (TNS), uses a total of four zones for wing configurations (or can be extended to complete aircraft configurations by adding zones). In the inner blocks near the wing surface, the thin-layer Navier-Stokes equations are solved, while in the outer two blocks the Euler equations are solved. The diagonal algorithm yields a speedup of as much as a factor of 40 over the original algorithm/zonal method code. The TNS code, in addition, has the capability to model wind tunnel walls. Transonic viscous solutions are obtained on a 150,000-point mesh for a NACA 0012 wing. A three-order-of-magnitude drop in the L2-norm of the residual requires approximately 500 iterations, which takes about 45 min of CPU time on a Cray-XMP processor. Simulations are also conducted for a different geometrical wing called WING C. All cases show good agreement with experimental data.

Flores, J.↗
NASA NTRS1996

Three-dimensional Navier-Stokes analysis and redesign of an imbedded bellmouth nozzle in a turbine cascade inlet section

Experimental measurements in the inlet of a transonic turbine blade cascade showed unacceptable pitchwise flow non-uniformity. A three-dimensional, Navier-Stokes computational fluid dynamics (CFD) analysis of the imbedded bellmouth inlet in the facility was performed to identify and eliminate the source of the flow non-uniformity. The blockage and acceleration effects of the blades were accounted for by specifying a periodic static pressure exit condition interpolated from a separate three-dimensional Navier-Stokes CFD solution of flow around a single blade in an infinite cascade. Calculations of the original inlet geometry showed total pressure loss regions consistent in strength and location to experimental measurements. The results indicate that the distortions were caused by a pair of streamwise vortices that originated as a result of the interaction of the flow with the imbedded bellmouth. Computations were performed for an inlet geometry which eliminated the imbedded bellmouth by bridging the region between it and the upstream wall. This analysis indicated that eliminating the imbedded bellmouth nozzle also eliminates the pair of vortices, resulting in a flow with much greater pitchwise uniformity. Measurements taken with an installed redesigned inlet verify that the flow non-uniformity has indeed been eliminated.

VanFossen, G. J.↗
NASA NTRS1992

A study of hypersonic swept shock wave/turbulent boundary layer interactions using a conical Navier-Stokes code

A computational study has been performed of sharp fin-induced swept shock wave/turbulent boundary layer interactions at low hypersonic Mach numbers. The objective was to determine if results obtained using a conical Navier-Stokes code, particularly the peak heating and pressure, are adequate for engineering predictions. The advantage of the conical approach is that the problem becomes two-dimensional and requires much less computational effort than a fully three-dimensional calculation. In this code, the standard Baldwin-Lomax model is used for turbulent closure and its performance is studied in some detail. To assess the approach interactive flowfields generated by unswept sharp fins at two angles of attack at each of three hyprsonic freestream Mach numbers (5, 6, 11) have been calculated and the results compared with experimental wall pressure and heat transfer data. Although the conical Navier-Stokes/Baldwin-Lomax approach is reasonably successful at Mach numbers up to 5, the performance deteriorates as the Mach number is increased. Nevertheless, the approach could be a valuable tool in preliminary parametric design studies.

Rodi, Patrick E.↗
NASA NTRS1986

Transonic Navier-Stokes wing solution using a zonal approach. Part 1: Solution methodology and code validation

A fast diagonalized Beam-Warming algorithm is coupled with a zonal approach to solve the three-dimensional Euler/Navier-Stokes equations. The computer code, called Transonic Navier-Stokes (TNS), uses a total of four zones for wing configurations (or can be extended to complete aircraft configurations by adding zones). In the inner blocks near the wing surface, the thin-layer Navier-Stokes equations are solved, while in the outer two blocks the Euler equations are solved. The diagonal algorithm yields a speedup of as much as a factor of 40 over the original algorithm/zonal method code. The TNS code, in addition, has the capability to model wind tunnel walls. Transonic viscous solutions are obtained on a 150,000-point mesh for a NACA 0012 wing. A three-order-of-magnitude drop in the L2-norm of the residual requires approximately 500 iterations, which takes about 45 min of CPU time on a Cray-XMP processor. Simulations are also conducted for a different geometrical wing called WING C. All cases show good agreement with experimental data.

Flores, J.↗
NASA NTRS1993

CNSFV code development, virtual zone Navier-Stokes computations of oscillating control surfaces and computational support of the laminar flow supersonic wind tunnel

The work performed during the past year on this cooperative agreement covered two major areas and two lesser ones. The two major items included further development and validation of the Compressible Navier-Stokes Finite Volume (CNSFV) code and providing computational support for the Laminar Flow Supersonic Wind Tunnel (LFSWT). The two lesser items involve a Navier-Stokes simulation of an oscillating control surface at transonic speeds and improving the basic algorithm used in the CNSFV code for faster convergence rates and more robustness. The work done in all four areas is in support of the High Speed Research Program at NASA Ames Research Center.

Klopfer, Goetz H.↗
NASA NTRS1993

Chemical nonequilibrium Navier-Stokes solutions for hypersonic flow over an ablating graphite nosetip

The general boundary conditions including mass and energy balances of chemically equilibrated or nonequilibrated gas adjacent to ablating surfaces have been derived. A computer procedure based on these conditions was developed and interfaced with the Navier-Stokes solver for predictions of the flow field, surface temperature, and surface ablation rates over re-entry space vehicles with ablating Thermal Protection Systems (TPS). The Navier-Stokes solver with general surface thermochemistry boundary conditions can predict more realistic solutions and provide useful information for the design of TPS. A test case with a proposed hypersonic test vehicle configuration and associated free stream conditions was developed. Solutions with various surface boundary conditions were obtained, and the effect of nonequilibrium gas as well as surface chemistry on surface heating and ablation rate were examined. The solutions of the GASP code with complete ablating surface conditions were compared with those of the ASC code. The direction of future work is also discussed.

Chen, Y. K.↗
NASA NTRS2003

Algorithmic Enhancements to the VULCAN Navier-Stokes Solver

VULCAN (Viscous Upwind aLgorithm for Complex flow ANalysis) is a cell centered, finite volume code used to solve high speed flows related to hypersonic vehicles. Two algorithms are presented for expanding the range of applications of the current Navier-Stokes solver implemented in VULCAN. The first addition is a highly implicit approach that uses subiterations to enhance block to block connectivity between adjacent subdomains. The addition of this scheme allows more efficient solution of viscous flows on highly-stretched meshes. The second algorithm addresses the shortcomings associated with density-based schemes by the addition of a time-derivative preconditioning strategy. High speed, compressible flows are typically solved with density based schemes, which show a high level of degradation in accuracy and convergence at low Mach numbers (M less than or equal to 0.1). With the addition of preconditioning and associated modifications to the numerical discretization scheme, the eigenvalues will scale with the local velocity, and the above problems will be eliminated. With these additions, VULCAN now has improved convergence behavior for multi-block, highly-stretched meshes and also can solve the Navier-Stokes equations for very low Mach numbers.

Litton, D. K.↗
NASA NTRS1994

A hybrid Navier-Stokes/full-potential method for the prediction of iced wing aerodynamics

A hybrid method for computing compressible viscous flows is presented. This method divides the computational domain into two zones. In the outer zone, the unsteady full-potential equation (FPE) is solved. In the inner zone, the Navier-Stokes equations are solved. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. The resulting CPU times are less than 50 percent of the required for a full-blown Navier-Stokes analysis. Sample applications of the method to an unswept iced wing at 4 deg and 8 deg angle of attack are presented. Surface pressures are in good agreement with the measurements obtained by Bragg et al. at the University of Illinois.

Mello, O. A. F.↗
NASA NTRS1972

Numerical solutions of the Navier-Stokes equations for the supersonic laminar flow over a two-dimensional compression corner

Numerical solutions have been obtained for the supersonic, laminar flow over a two-dimensional compression corner. These solutions were obtained as steady-state solutions to the unsteady Navier-Stokes equations using the finite difference method of Brailovskaya, which has second-order accuracy in the spatial coordinates. Good agreement was obtained between the computed results and wall pressure distributions measured experimentally for Mach numbers of 4 and 6.06, and respective Reynolds numbers, based on free-stream conditions and the distance from the leading edge to the corner. In those calculations, as well as in others, sufficient resolution was obtained to show the streamline pattern in the separation bubble. Upstream boundary conditions to the compression corner flow were provided by numerically solving the unsteady Navier-Stokes equations for the flat plate flow field, beginning at the leading edge. The compression corner flow field was enclosed by a computational boundary with the unknown boundary conditions supplied by extrapolation from internally computed points.

Carter, J. E.↗