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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1995

The 3D 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 three-dimensional (3D) Navier-Stokes flow solver. A two-equation turbulence model, and a porous bleed model based on unchoked bleed hole discharge coefficients were used. Comparisons were made with experimental data, inviscid theory, and two-dimensional Navier-Stokes analyses. The main objective was to gain insight into the inlet fluid dynamics. Examination of the computational results along with the experimental data suggest that the cowl shock-sidewall boundary layer interaction near the leading edge caused a substantial separation in the wind tunnel inlet model. As a result, the inlet performance may have been compromised by increased spillage and higher bleed mass flow requirements. The internal flow contained substantial waves that were not in the original inviscid design. 3D effects were fairly minor for this inlet at on-design conditions. Navier-Stokes analysis appears to be an useful tool for gaining insight into the inlet fluid dynamics. It provides a higher fidelity simulation of the flowfield than the original inviscid design, by taking into account boundary layers, porous bleed, and their interactions with shock waves.

Mizukami, M.↗

NASA NTRS1993

Hypersonic blunt body wake computations using DSMC and Navier-Stokes solvers

Numerical results obtained with direct simulation Monte Carlo (DSMC) and Navier-Stokes methods are presented for Mach 20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how does the wake structure change as a function of rarefaction, what are the afterbody levels of heating, and to what limits are continuum models realistic as rarefaction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the afterbody sting are emphasized in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Moss, James N.↗

NASA NTRS2001

Prediction of Business Jet Airloads Using The Overflow Navier-Stokes Code

The objective of this work is to evaluate the application of Navier-Stokes computational fluid dynamics technology, for the purpose of predicting off-design condition airloads on a business jet configuration in the transonic regime. The NASA Navier-Stokes flow solver OVERFLOW with Chimera overset grid capability, availability of several numerical schemes and convergence acceleration techniques was selected for this work. A set of scripts which have been compiled to reduce the time required for the grid generation process are described. Several turbulence models are evaluated in the presence of separated flow regions on the wing. Computed results are compared to available wind tunnel data for two Mach numbers and a range of angles-of-attack. Comparisons of wing surface pressure from numerical simulation and wind tunnel measurements show good agreement up to fairly high angles-of-attack.

Bounajem, Elias↗

NASA NTRS1994

On Bi-Grid Local Mode Analysis of Solution Techniques for 3-D Euler and Navier-Stokes Equations

A procedure is presented for utilizing a bi-grid stability analysis as a practical tool for predicting multigrid performance in a range of numerical methods for solving Euler and Navier-Stokes equations. Model problems based on the convection, diffusion and Burger's equation are used to illustrate the superiority of the bi-grid analysis as a predictive tool for multigrid performance in comparison to the smoothing factor derived from conventional von Neumann analysis. For the Euler equations, bi-grid analysis is presented for three upwind difference based factorizations, namely Spatial, Eigenvalue and Combination splits, and two central difference based factorizations, namely LU and ADI methods. In the former, both the Steger-Warming and van Leer flux-vector splitting methods are considered. For the Navier-Stokes equations, only the Beam-Warming (ADI) central difference scheme is considered. In each case, estimates of multigrid convergence rates from the bi-grid analysis are compared to smoothing factors obtained from single-grid stability analysis. Effects of grid aspect ratio and flow skewness are examined. Both predictions are compared with practical multigrid convergence rates for 2-D Euler and Navier-Stokes solutions based on the Beam-Warming central scheme.

Ibraheem, S. O.↗

NASA NTRS1995

Convergence Acceleration of a Navier-Stokes Solver for Efficient Static Aeroelastic Computations

New capabilities have been developed for a Navier-Stokes solver to perform steady-state simulations more efficiently. The flow solver for solving the Navier-Stokes equations is based on a combination of the lower-upper factored symmetric Gauss-Seidel implicit method and the modified Harten-Lax-van Leer-Einfeldt upwind scheme. A numerically stable and efficient pseudo-time-marching method is also developed for computing steady flows over flexible wings. Results are demonstrated for transonic flows over rigid and flexible wings.

Obayashi, Shigeru↗

NASA NTRS1987

A comparison of numerical flux formulas for the Euler and Navier-Stokes equations

Numerical flux formulas for the convection terms in the Euler or Navier-Stokes equations are analyzed with regard to their accuracy in representing steady nonlinear and linear waves (shocks and entropy/shear waves, respectively). Numerical results are obtained for a one-dimensional conical Navier-Stokes flow including both a shock and a boundary layer. Analysis and experiments indicate that for an accurate representation of both layers the flux formula must include information about all different waves by which neighboring cells interact, as in Roe's flux-difference splitting. In comparison, Van Leer's flux-vector splitting, which ignores the linear waves, badly diffuses the boundary layer. The results of MacCormack's scheme, if properly tuned, are significantly better. The use of a sufficiently detailed flux formula appears to reduce the number of cells required to resolve a boundary layer by a factor 1/2 to 1/4 and thus pays off.

Van Leer, Bram↗

NASA NTRS1986

Upwind Navier-Stokes solutions for separated periodic flows

The application of an upwind implicit approximate factorization Navier-Stokes algorithm to highly separated flow is described. Using both the thin-layer and complete forms of the Navier-Stokes equations, the low Reynolds number laminar flow around a circular cylinder with periodic shedding is solved. The effect of grid density, grid extent, and time step on the Strouhal number is shown. Results from both sets of equations agree within the experimental data band. Unsteady, laminar flow computations around inclined plates and separated airfoils are also described. Strouhal numbers agree to within 5 percent of experiments for inclined plates. Differences between the complete equations and the thin-layer approximation for separated periodic flows are discussed. Computations of an impulsively started circular cylinder and airfoil yield time-accurate flowfield shapes in good agreement with experimental flow visualizations. The turbulent computation of an airfoil at a high angle-of-attack is massively separated, but shows no evidence of periodicity.

Rumsey, C. L.↗

NASA NTRS1990

Navier-Stokes computations on swept-tapered wings, including flexibility

A procedure to couple the Navier-Stokes solutions with modal structural equations of motion is presented for computing aeroelastic responses of flexible fighter wings. The Navier-Stokes flow equations are solved by a finite-difference scheme with dynamic grids. The coupled aeroelastic equations of motion are solved using the linear-acceleration method. The configuration-adaptive dynamic grids are time-accurately generated using the aeroelastically deformed shape of the wing. The coupled calculations are compared with experiments when available. Effects of flexibility and pitch rate are demonstrated for flows with vortices. Turbulent flow computations are also compared with laminar flow computations.

Guruswamy, Guru P.↗

NASA NTRS1993

Three-dimensional flow analysis inside turbomachinery stages with steady and unsteady Navier-Stokes method

This study presents a numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed in the present study. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow field of a subsonic turbine stage and the unsteady flow field inside a transonic, high-through-flow, axial compressor stage. The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid.

Copenhaver, W. W.↗

NASA NTRS1991

Navier-Stokes analysis of a liquid rocket engine disk cavity

This paper presents a Navier-Stokes analysis of hydrodynamic phenomena occurring in the aft disk cavity of a liquid rocket engine turbine. The cavity analyzed in the Space Shuttle Main Engine Alternate Turbopump currently being developed by NASA and Pratt and Whitney. Comparison of results obtained from the Navier-Stokes code for two rotating disk datasets available in the literature are presented as benchmark validations. The benchmark results obtained using the code show good agreement relative to experimental data, and the turbine disk cavity was analyzed with comparable grid resolution, dissipation levels, and turbulence models. Predicted temperatures in the cavity show that little mixing of hot and cold fluid occurs in the cavity and the flow is dominated by swirl and pumping up the rotating disk.

Benjamin, Theodore G.↗

NASA NTRS1974

Numerical solution of the steady-state Navier-Stokes equations for hypersonic flow about blunt axisymmetric bodies

The steady-state Navier-Stokes equations are solved for hypersonic flow about blunt axisymmetric bodies. The equations of motion are solved by successive approximations using an implicit finite-difference scheme. The results are compared with viscous shock-layer theory, experimental data, and time-dependent solutions of the Navier-Stokes equations. It is demonstrated that viscous shock-layer theory is sufficiently accurate for the range of flight conditions normally encountered by entry vehicles.

Anderson, E. C.↗

NASA NTRS1976

Numerical solution of the steady-state Navier-Stokes equations for hypersonic flow about blunt axisymmetric bodies

Anderson, E. C.↗

DOE OSTI2020

Critique on “A Revisit of Navier-Stokes Equation”

In a recently published work (Sheng, 2020), Sheng identified inconsistencies in the concepts and principles leading to the Navier–Stokes equations and proposed a new asymmetric “friction tensor” to resolve these inconsistencies. We show that this friction tensor can be interpreted as the non-isotropic part of a vorticity-influenced viscous stress tensor following the analysis in Berdahl and Strang (1986). The asymmetric stress tensor however violates the principles of mechanical equilibrium in fluid flows that are free of external or body couples and is therefore non-physical in classical continuum hydrodynamics. We refute the arguments in favour of the asymmetric stress tensor in Sheng (2020) whilst also demonstrating that the friction tensor leads to the incorrect form of the Navier–Stokes equations in non-Cartesian co-ordinates and its ability to recover the correct governing equations in Cartesian co-ordinates is fortuitous.

asymmetric↗

NASA NTRS1983

A fully vectorized numerical solution of the incompressible Navier-Stokes equations

A vectorizable algorithm is presented for the implicit finite difference solution of the incompressible Navier-Stokes equations in general curvilinear coordinates. The unsteady Reynolds averaged Navier-Stokes equations solved are in two dimension and non-conservative primitive variable form. A two-layer algebraic eddy viscosity turbulence model is used to incorporate the effects of turbulence. Two momentum equations and a Poisson pressure equation, which is obtained by taking the divergence of the momentum equations and satisfying the continuity equation, are solved simultaneously at each time step. An elliptic grid generation approach is used to generate a boundary conforming coordinate system about an airfoil. The governing equations are expressed in terms of the curvilinear coordinates and are solved on a uniform rectangular computational domain. A checkerboard SOR, which can effectively utilize the computer architectural concept of vector processing, is used for iterative solution of the governing equations.

Patel, N.↗

NASA NTRS2014

A Modular Approach to Model Oscillating Control Surfaces Using Navier Stokes Equations

The use of active controls for rotorcraft is becoming more important for modern aerospace configurations. Efforts to reduce the vibrations of helicopter blades with use of active-controls are in progress. Modeling oscillating control surfaces using the linear aerodynamics theory is well established. However, higher-fidelity methods are needed to account for nonlinear effects, such as those that occur in transonic flow. The aeroelastic responses of a wing with an oscillating control surface, computed using the transonic small perturbation (TSP) theory, have been shown to cause important transonic flow effects such as a reversal of control surface effectiveness that occurs as the shock wave crosses the hinge line. In order to account for flow complexities such as blade-vortex interactions of rotor blades higher-fidelity methods based on the Navier-Stokes equations are used. Reference 6 presents a procedure that uses the Navier-Stokes equations with moving-sheared grids and demonstrates up to 8 degrees of control-surface amplitude, using a single grid. Later, this procedure was extended to accommodate larger amplitudes, based on sliding grid zones. The sheared grid method implemented in EulerlNavier-Stokes-based aeroelastic code ENS AERO was successfully applied to active control design by industry. Recently there are several papers that present results for oscillating control surface using Reynolds Averaged Navier-Stokes (RANS) equations. References 9 and 10 report 2-D cases by filling gaps with overset grids. Reference 9 compares integrated forces with the experiment at low oscillating frequencies whereas Ref. 10 reports parametric studies but with no validation. Reference II reports results for a 3D case by modeling the gap region with a deformed grid and compares force results with the experiment only at the mid-span of flap. In Ref. II grid is deformed to match the control surface deflections at the section where the measurements are made. However, there is no indication in Ref. II that the gaps are explicitly modeled as in Ref. 6. Computations using overset grids are reported in Ref. 12 for a case by adding moving control surface to an existing blade but with no validation either with an experiment or another computation.

Navier↗

NASA NTRS1986

Pseudo-time algorithms for the Navier-Stokes equations

A pseudo-time method is introduced to integrate the compressible Navier-Stokes equations to a steady state. This method is a generalization of a method used by Crocco and also by Allen and Cheng. We show that for a simple heat equation that this is just a renormalization of the time. For a convection-diffusion equation the renormalization is dependent only on the viscous terms. We implement the method for the Navier-Stokes equations using a Runge-Kutta type algorithm. This permits the time step to be chosen based on the inviscid model only. We also discuss the use of residual smoothing when viscous terms are present.

Swanson, R. C.↗

NASA NTRS1992

Transonic airfoil and wing design using Navier-Stokes codes

An iterative design method has been implemented into 2D and 3D Navier-Stokes codes for the design of airfoils or wings with given target pressure distributions. The method begins with the analysis of an initial geometry, and obtains the analysis pressure distributions of that geometry. The differences between analysis pressures and target pressures are used to drive geometry changes through the use of a streamline curvature method. This paper describes the procedure that makes the iterative design method work for Navier-Stokes codes. Examples of 2D airfoil design, and 3D wing design are included. It is demonstrated that the method is highly effective for airfoil or wing design at flow conditions where no substantial separation occurs. Problems encountered in the airfoil design with shock induced flow separations are discussed.

Yu, N. J.↗

NASA NTRS1988

Navier-Stokes simulation of transonic wing flow fields using a zonal grid approach

The transonic Navier-Stokes code was used to simulate flow fields about isolated wings for workshop wind-tunnel and free-air cases using the thin-layer Reynolds-averaged Navier-Stokes equations. An implicit finite-difference scheme based on a diagonal version of the Beam-Warming algorithm was used to integrate the governing equations. A zonal grid approach was used to allow efficient grid refinement near the wing surface. The flow field was sensitive to the turbulent transition model, and flow unsteadiness was observed for a wind-tunnel case but not for the corresponding free-air case. The specification of experimental pressure at the wind-tunnel exit plane is the primary reason for the difference of these two numerical solutions.

Chaderjian, Neal M.↗