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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS2000

A gas-kinetic BGK scheme for the compressible Navier-Stokes equations

This paper presents an improved gas-kinetic scheme based on the Bhatnagar-Gross-Krook (BGK) model for the compressible Navier-Stokes equations. The current method extends the previous gas-kinetic Navier-Stokes solver developed by Xu and Prendergast by implementing a general nonequilibrium state to represent the gas distribution function at the beginning of each time step. As a result, the requirement in the previous scheme, such as the particle collision time being less than the time step for the validity of the BGK Navier-Stokes solution, is removed. Therefore, the applicable regime of the current method is much enlarged and the Navier-Stokes solution can be obtained accurately regardless of the ratio between the collision time and the time step. The gas-kinetic Navier-Stokes solver developed by Chou and Baganoff is the limiting case of the current method, and it is valid only under such a limiting condition. Also, in this paper, the appropriate implementation of boundary condition for the kinetic scheme, different kinetic limiting cases, and the Prandtl number fix are presented. The connection among artificial dissipative central schemes, Godunov-type schemes, and the gas-kinetic BGK method is discussed. Many numerical tests are included to validate the current method.

Xu, Kun↗

DOE OSTI2021

Modeling and simulation of transitional Rayleigh–Taylor flow with partially averaged Navier–Stokes equations

In this work, the partially averaged Navier–Stokes (PANS) equations are used to predict the variable-density Rayleigh–Taylor (RT) flow at Atwood number 0.5 and maximum Reynolds number 500. This is a prototypical problem of material mixing, featuring laminar, transitional, and turbulent flow, instabilities and coherent structures, density fluctuations, and production of turbulence kinetic energy by both shear and buoyancy mechanisms. These features pose numerous challenges to modeling and simulation, making the RT flow ideal to develop the validation space of the recently proposed PANS Besnard–Harlow–Rauenzahn-linear eddy viscosity model closure. The numerical simulations are conducted at different levels of physical resolution and test three approaches to set the parameters $f_\phi$ defining the range of physically resolved scales. The computations demonstrate the efficiency (accuracy vs cost) of the PANS model predicting the spatiotemporal development of the RT flow. Results comparable to large-eddy simulations and direct numerical simulations are obtained at significantly lower physical resolution without the limitations of the Reynolds-averaged Navier–Stokes equations in these transitional flows. The data also illustrate the importance of appropriate selection of the physical resolution and the resolved fraction of each dependent quantity $\phi$ of the turbulent closure, $f_\phi$. These two aspects determine the ability of the model to resolve the flow phenomena not amenable to modeling by the closure and, as such, the computations’ fidelity.

Navier Stokes equations↗

NASA NTRS1986

Numerical solutions of the complete Navier-Stokes equations

Using ideas from the kinetic theory, the Navier-Stokes equations are modified in such a way that they can be cast as a set of first order hyperbolic equations. This is achieved by incorporating time dependent terms into the definition of the stress tensor and the heat flux vectors. The boundary conditions are then determined from the theory of characteristics. Because the resulting equations reduce to the traditional Navier-Stokes equations when the steady state is reached, the present approach provides a straightforward scheme for the determination of inflow and outflow boundary conditions. The method is validated by comparing its predictions with known exact solutions of the steady Navier-Stokes equations.

Hassan, H. A.↗

NASA NTRS1986

Navier-Stokes computations for circulation controlled airfoils

Navier-Stokes computations of subsonic to transonic flow past airfoils with augmented lift due to rearward jet blowing over a curved trailing edge are presented. The approach uses a spiral grid topology. Solutions are obtained using a Navier-Stokes code which employs an implicit finite difference method, an algebraic turbulence model, and developments which improve stability, convergence, and accuracy. Results are compared against experiments for no jet blowing and moderate jet pressures and demonstrate the capability to compute these complicated flows.

Pulliam, T. H.↗

NASA NTRS1987

Navier-Stokes computations for circulation control airfoils

Pulliam, Thomas H.↗

NASA NTRS1986

The calculation of supersonic viscous flows using the parabolized Navier-Stokes equations

Solution of the parabolic Navier-Stokes (PNS) equations for supersonic flows is discussed, and compatibility of the PNS method with the triple-deck theory of Stewartson (1974) is demonstrated. Characteristic and stability analyses show that use of an appropriate filter on the pressure term in the x-momentum equation can suppress the departure solutions, giving the usual desired weak interaction solution. An Alternating Direction Explicit procedure, with minimal computer storage requirements compared to the full Navier-Stokes solvers, is proposed to calculate strongly interacting flows using a global iteration procedure for the PNS equations. The PNS equations are used to solve the hypersonic viscous interaction problem, and good agreement is found with experimental results.

Davis, R. T.↗

NASA NTRS1987

Transonic Navier-Stokes solutions for a fighter-like configuration

The Transonic Navier-Stokes wing code is extended to a 16-zone TNS wing-fuselage code and used to solve the transonic viscous flow over a modified F-16A. The computer code, called Transonic Navier-Stokes Wing/Fuselage uses a zonal approach to solve the three-dimensional Euler and Navier-Stokes equations. With the zonal implementation, clustering suitable for viscous calculations is achieved on all solid surfaces. The transonic case has flow conditions of free-stream M = 0.9, alpha = 4.12 deg, and a Reynolds number based on root chord of 4.5 million. This case required about 3,000 iterations to reduce the L2-norm of the residual by three orders, which takes about 15 hr of cpu time on the Cray X-MP/48 processor. Pressure distributions, as well as separaton patterns, compare favorably with experiment for this transonic case.

Flores, Jolen↗

NASA NTRS1993

Calculation of AGARD Wing 445.6 Flutter Using Navier-Stokes Aerodynamics

The flutter characteristics of the first AGARD standard aeroelastic configuration for dynamic response, Wing 445.6, are studied using an unsteady Navier-Stokes algorithm in order to investigate a previously noted discrepancy between Euler flutter characteristics and the experimental data. The algorithm, which is a three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1), was previously modified for the time-marching, aeroelastic analysis of wings using the unsteady Euler equations. These modifications include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time integration with the governing flow equations. In this paper, the aeroelastic method is extended and evaluated for applications that use the Navier- Stokes aerodynamics. The paper presents a brief description of the aeroelastic method and presents unsteady calculations which verify this method for Navier-Stokes calculations. A linear stability analysis and a time-marching aeroelastic analysis are used to determine the flutter characteristics of the isolated 45 deg. swept-back wing. Effects of fluid viscosity, structural damping, and number of modes in the structural model are investigated. For the linear stability analysis, the unsteady generalized aerodynamic forces of the wing are computed for a range of reduced frequencies using the pulse transfer-function approach. The flutter characteristics of the wing are determined using these unsteady generalized aerodynamic forces in a traditional V-g analysis. This stability analysis is used to determine the flutter characteristics of the wing at free-stream Mach numbers of 0.96 and 1.141 using the generalized aerodynamic forces generated by solving the Euler equations and the Navier-Stokes equations. Time-marching aeroelastic calculations are performed at a free-stream Mach number of 1.141 using the Euler and Navier-Stokes equations to compare with the linear V-g flutter analysis method. The V-g analysis, which is used in conjunction with the time-marching analysis, indicates that the fluid viscosity has a significant effect on the supersonic flutter boundary for this wing while the structural damping and number of modes in the structural model have a lesser effect.

Lee-Rausch, Elizabeth M.↗

NASA NTRS1995

Progress in incompressible Navier-Stokes computations for propulsion flows and its dual-use applications

Development of an incompressible Navier-Stokes solution procedure was performed for the analysis of a liquid rocket engine pump components and for the mechanical heart assist devices. The solution procedure for the propulsion systems is applicable to incompressible Navier-Stokes flows in a steadily rotating frame of reference for any general complex configurations. The computer codes were tested on different complex configurations such as liquid rocket engine inducer and impellers. As a spin-off technology from the turbopump component simulations, the flow analysis for an axial heart pump was conducted. The baseline Left Ventricular Assist Device (LVAD) design was improved by adding an inducer geometry by adapting from the liquid rocket engine pump. The time-accurate mode of the incompressible Navier-Stokes code was validated with flapping foil experiment by using different domain decomposition methods. In the flapping foil experiment, two upstream NACA 0025 foils perform high-frequency synchronized motion and generate unsteady flow conditions for a downstream larger stationary foil. Fairly good agreement was obtained between unsteady experimental data and numerical results from two different moving boundary procedures. Incompressible Navier-Stokes code (INS3D) has been extended for heat transfer applications. The temperature equation was written for both forced and natural convection phenomena. Flow in a square duct case was used for the validation of the code in both natural and forced convection.

Kiris, Cetin↗

NASA NTRS1990

Navier-Stokes and Euler solutions for lee-side flows over supersonic delta wings. A correlation with experiment

An Euler flow solver and a thin layer Navier-Stokes flow solver were used to numerically simulate the supersonic leeside flow fields over delta wings which were observed experimentally. Three delta wings with 75, 67.5, and 60 deg leading edge sweeps were computed over an angle-of-attack range of 4 to 20 deg at a Mach number 2.8. The Euler code and Navier-Stokes code predict equally well the primary flow structure where the flow is expected to be separated or attached at the leading edge based on the Stanbrook-Squire boundary. The Navier-Stokes code is capable of predicting both the primary and the secondary flow features for the parameter range investigated. For those flow conditions where the Euler code did not predict the correct type of primary flow structure, the Navier-Stokes code illustrated that the flow structure is sensitive to boundary layer model. In general, the laminar Navier-Stokes solutions agreed better with the experimental data, especially for the lower sweep delta wings. The computational results and a detailed re-examination of the experimental data resulted in a refinement of the flow classifications. This refinement in the flow classification results in the separation bubble with the shock flow type as the intermediate flow pattern between separated and attached flows.

Mcmillin, S. Naomi↗

NASA NTRS1993

A new flux conserving Newton's method scheme for the two-dimensional, steady Navier-Stokes equations

A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.

Scott, James R.↗

NASA NTRS1989

Vortical flow computations on swept flexible wings using Navier-Stokes equations

A procedure to couple the Navier-Stokes solutions with modal structural equations of motion is presented for computing aeroelastic responses of swept flexible 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 calculations are compared with the experiment when available. Effects of flexibility, sweep angle, and pitch rate are demonstrated for flows with vortices.

Guruswamy, Guru P.↗

NASA NTRS1985

Viscous-inviscid interaction computations using a pseudo Navier-Stokes approach

A new method is presented for the computation of viscous-inviscid interaction. The idea is to treat rotational inviscid flow (of which flows are almost entirely composed) in a thorough manner, and accept an approximation treatment of vorticity as introduced by viscous effects. The approach is to numerically solve the Navier-Stokes equations with the viscous terms determined from an inverse boundary-layer solution. The method falls somewhere between a Navier-Stokes approach and an Euler and boundary-layer equation coupling approach; consequently, it is referred to as a pseudo Navier-Stokes approach. Results from both the Navier-Stokes equations and the pseudo Navier-Stokes approach are presented.

Whitfield, D. L.↗

NASA NTRS1988

On multigrid methods for the Navier-Stokes Computer

The overall architecture of the multipurpose parallel-processing Navier-Stokes Computer (NSC) being developed by Princeton and NASA Langley (Nosenchuck et al., 1986) is described and illustrated with extensive diagrams, and the NSC implementation of an elementary multigrid algorithm for simulating isotropic turbulence (based on solution of the incompressible time-dependent Navier-Stokes equations with constant viscosity) is characterized in detail. The present NSC design concept calls for 64 nodes, each with the performance of a class VI supercomputer, linked together by a fiber-optic hypercube network and joined to a front-end computer by a global bus. In this configuration, the NSC would have a storage capacity of over 32 Gword and a peak speed of over 40 Gflops. The multigrid Navier-Stokes code discussed would give sustained operation rates of about 25 Gflops.

Nosenchuck, D. M.↗

NASA NTRS1981

Two-boundary grid generation for the solution of the three dimensional compressible Navier-Stokes equations

A grid generation technique called the two boundary technique is developed and applied for the solution of the three dimensional Navier-Stokes equations. The Navier-Stokes equations are transformed from a cartesian coordinate system to a computational coordinate system, and the grid generation technique provides the Jacobian matrix describing the transformation. The two boundary technique is based on algebraically defining two distinct boundaries of a flow domain and the distribution of the grid is achieved by applying functions to the uniform computational grid which redistribute the computational independent variables and consequently concentrate or disperse the grid points in the physical domain. The Navier-Stokes equations are solved using a MacCormack time-split technique. Grids and supersonic laminar flow solutions are obtained for a family of three dimensional corners and two spike-nosed bodies.

Smith, R. E.↗

NASA NTRS1999

A Comparison of Three Navier-Stokes Solvers for Exhaust Nozzle Flowfields

A comparison of the NPARC, PAB, and WIND (previously known as NASTD) Navier-Stokes solvers is made for two flow cases with turbulent mixing as the dominant flow characteristic, a two-dimensional ejector nozzle and a Mach 1.5 elliptic jet. The objective of the work is to determine if comparable predictions of nozzle flows can be obtained from different Navier-Stokes codes employed in a multiple site research program. A single computational grid was constructed for each of the two flows and used for all of the Navier-Stokes solvers. In addition, similar k-e based turbulence models were employed in each code, and boundary conditions were specified as similarly as possible across the codes. Comparisons of mass flow rates, velocity profiles, and turbulence model quantities are made between the computations and experimental data. The computational cost of obtaining converged solutions with each of the codes is also documented. Results indicate that all of the codes provided similar predictions for the two nozzle flows. Agreement of the Navier-Stokes calculations with experimental data was good for the ejector nozzle. However, for the Mach 1.5 elliptic jet, the calculations were unable to accurately capture the development of the three dimensional elliptic mixing layer.

Georgiadis, Nicholas J.↗

NASA NTRS1986

A comparison of interactive boundary-layer and thin-layer Navier-Stokes procedures

Numerical procedures for solving the thin-shear-layer Navier-Stokes equations and for the interaction of solutions to inviscid and boundary-layer equations are described and evaluated. To allow appraisal of the numerical and fluid dynamic abilities of the two schemes, they have been applied to one airfoil as a function of angle of attack at two slightly different Reynolds numbers. The NACA 0012 airfoil has been chosen because it allows comparison with measured lift, drag, and moment and with surface-pressure distributions. Calculations have been performed with algebraic eddy-viscosity formulations, and they include consideration of transition. The results are presented in a form that allows easy appraisal of the accuracy of both procedures and of the relative costs. The interactive procedure is computationally efficient but restrictive relative to the thin-layer Navier-Stokes procedure. The latter procedure does a better job of predicting drag than does the former. In both procedures, the location of transition is crucial for accurate or detailed computations, particularly at high angles of attack. When the upstream influence of pressure field through the shear layer is important, the thin-layer Navier-Stokes procedure has an edge over the interactive procedure.

Mehta, U.↗

NASA NTRS1990

Prediction of vortical flows on wings using incompressible Navier-Stokes equations

Numerical solutions to the incompressible Navier-Stokes equations are in greater demand than ever before as the field of computational fluid dynamics (CFD) increases its impact as an engineering tool. An implicit upwind-relation finite-difference scheme is developed for solving unsteady incompressible Navier-Stokes equations incorporating an artificial time derivative of the pressure to the continuity equation. The essential features observed in experiments for complicated vortical flows past three round-edged low-aspect-ratio wings at high angles of attack are reproduced both qualitatively and quantitatively.

Hsu, C.-H.↗