Search NASA⌕ Search

SEARCH · Search NASA

Results for “Navier Stokes algorithm”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

A numerical simulation of the flow in the diffuser of the NASA Lewis Icing Research Tunnel

The flow in the diffuser section of the Icing Research Tunnel at the NASA Lewis Research Center is numerically investigated. To accomplish this, an existing computer code is utilized. The code, known as PARC3D, is based on the Beam-Warming algorithm applied to the strong conservation law form of the complete Navier-Stokes equations. The first portion of the paper consists of a brief description of the diffuser and its current flow characteristics. A brief discussion of the code work follows. Predicted velocity patterns are then compared with the measured values.

Addy, Harold E., Jr.↗
NASA NTRS1991

Choice of velocity variables for complex flow computation

The issue of adopting the velocity components as dependent velocity variables for the Navier-Stokes flow computations is investigated. The viewpoint advocated is that a numerical algorithm should preferably honor both the physical conservation law in differential form and the geometric conservation law in discrete form. With the use of Cartesian velocity vector, the momentum equations in curvilinear coordinates can retain the full conservation-law form and satisfy the physical conservation laws. With the curvilinear velocity components, source terms appear in differential equations and hence the full conservation law form can not be retained. In discrete expressions, algorithms based on the Cartesian components can satisfy the geometric conservation-law form for convection terms but not for viscous terms; those based on the curvilinear components, on the other hand, cannot satisfy the geometric conservation-law form for either convection or viscous terms. Several flow solutions for domain with 90 and 360 degree turnings are presented to illustrate the issues of using the Cartesian velocity components and the staggered grid arrangement.

Shyy, W.↗
NASA NTRS1992

Parabolized Navier-Stokes methods for hypersonic flows

The development and numerical solution of the Parabolized Navier-Stokes (PNS) equations is discussed. Special emphasis is placed on numerical integration algorithms for PNS solvers, most notably recent applications of upwind schemes. Application of this technology to high speed, high Reynolds number flows was the focus, and modeling techniques for real gas effects and turbulence were discussed. Experiences in validating the UPS upwind PNS flow solver for relatively simple body shapes are presented. The test cases range in complexity from two dimensional perfect gas cases to three dimensional finite rate chemistry cases.

Lawrence, S. L.↗
NASA NTRS1993

Accuracy of least-squares methods for the Navier-Stokes equations

Recently there has been substantial interest in least-squares finite element methods for velocity-vorticity-pressure formulations of the incompressible Navier-Stokes equations. The main cause for this interest is the fact that algorithms for the resulting discrete equations can be devised which require the solution of only symmetric, positive definite systems of algebraic equations. On the other hand, it is well-documented that methods using the vorticity as a primary variable often yield very poor approximations. Thus, here we study the accuracy of these methods through a series of computational experiments, and also comment on theoretical error estimates. It is found, despite the failure of standard methods for deriving error estimates, that computational evidence suggests that these methods are, at the least, nearly optimally accurate. Thus, in addition to the desirable matrix properties yielded by least-squares methods, one also obtains accurate approximations.

Bochev, Pavel B.↗
NASA NTRS1993

Numerical study of shock-wave/boundary-layer interactions with bleed

A numerical study was conducted to investigate how bleed through a two-dimensional slot affects shock-wave induced, boundary-layer separation on a flat plate. This study is based on the ensemble-averaged, compressible, Navier-Stokes equations closed by the Baldwin-Lomax, algebraic turbulence model. The algorithm used to obtain solutions was the implicit, partially split, two-factored scheme of Steger. This study examined the effects of the following parameters in controlling shock-wave induced flow separation: location of slot in relation to where the incident shock wave impinged on the boundary layer, size of slot in relation to the boundary-layer thickness, number of slots, spacings between slots, and strength of the incident shock wave. This study also showed the nature of the very complex flowfield about the slot or slots and how the plenum affects the bleed process. The results of this study are relevant to problems where bleed is used to control shock-wave induced, boundary-layer separation (e.g., inside jet engine inlets and wind tunnels).

Hahn, T. O.↗
NASA NTRS1993

Accuracy of least-squares methods for the Navier-Stokes equations

Recently there has been substantial interest in least-squares finite element methods for velocity-vorticity-pressure formulations of the incompressible Navier-Stokes equations. The main cause for this interest is the fact that algorithms for the resulting discrete equations can be devised which require the solution of only symmetric, positive definite systems of algebraic equations. On the other hand, it is well-documented that methods using the vorticity as a primary variable often yield very poor approximations. Thus, here we study the accuracy of these methods through a series of computational experiments, and also comment on theoretical error estimates. It is found, despite the failure of standard methods for deriving error estimates, that computational evidence suggests that these methods are, at the least, nearly optimally accurate. Thus, in addition to the desirable matrix properties yielded by least-squares methods, one also obtains accurate approximations.

Bochev, Pavel B.↗
NASA NTRS1984

A vectorized solution for incompressible flow

An algorithm is developed to obtain solutions to the unsteady Reynolds-averaged incompressible Navier-Stokes equations in general curvilinear coordinates on a vector processor. The governing equations are in nonconservative form with the velocity and pressure as dependent variables. Two momentum equations and the Poisson equation for pressure form a set of three governing equations for three flow field unknowns: u, v, and p. The governing equations and boundary conditions are expressed in terms of boundary-conforming curvilinear coordinates, and a checkerboard SOR iteration is used to solve the governing equations. Several possible sequences for a checkerboard SOR iteration are investigated for finding the best overall convergence rate. The efficiency and capability of the present algorithm was assessed using the example of an 18 percent thick NACA 66(3)018 airfoil at zero degree angle of attack for chord Reynolds number range 1000-40,000.

Patel, N. R.↗
NASA NTRS1989

Computational fluid dynamics research in three-dimensional zonal techniques

Patched-grid algorithms for the analysis of complex configurations with an implicit, upwind-biased Navier-Stokes solver were investigated. Conservative and non-conservative approaches for performing zonal interpolations were implemented. The latter approach yields the most flexible technique in that it can handle both patched and overlaid grids. Results for a two-dimensional blunt body problem show that either approach yield accurate steady-state shock locations and jump conditions. In addition, calculations of the turbulent flow through a hypersonic inlet on a three-zone grid show that the numerical prediction is in good agreement with the experimental results. Through the use of a generalized coordinate transformation at the zonal interface between two or more blocks, the algorithm can be applied to highly stretched viscous grids and to arbitrarily-shaped zonal boundaries. Applications were made to the F-18 aircraft at subsonic, high-alpha conditions, in support of the NASA High-Alpha Research Program. The calculations were compared to ground-based and flight test experiments and were used as a guide to understanding the ground-based tests, which are laminar and transitional, and their relationship to flight. Calculations about a complete reconnaissance aircraft were also performed in order to further demonstrate the capability of the patched-grid algorithm.

Walters, Robert W.↗
NASA NTRS1981

A numerical simulation of hypersonic viscous flow over arbitrary geometries at angle of attack

An implicit conservative, noniterative, finite-difference algorithm that predicts the supersonic, laminar or turbulent viscous flow about arbitrary geometries at large angles of attack is presented. The three-dimensional parabolized form of the thin-layer Navier-Stokes equations are written in generalized coordinates. These equations are solved using the delta form of the Beam-Warming implicit algorithm. Flow field simulations have been obtained for a blunt biconic with windward and leeward cuts and an X-24C lifting body for both laminar and turbulent flow at various Mach numbers and angles of attack. When compared with experiment or with previous theories, these computational predictions show good agreement.

Chaussee, D. S.↗
NASA NTRS1977

The finite element method in computational fluid mechanics

A finite element solution algorithm is established for a general statement of the time-averaged Navier-Stokes equations governing multi-dimensional turbulent flows. Numerical results are presented which evaluate factors affecting solution accuracy for a broad spectrum of linear and non-linear problem classes.

Baker, A. J.↗
NASA NTRS1986

High speed viscous flow calculations about complex configurations

Applications of the NASA Ames Parabolized Navier-Stokes (PNS) code to a variety of complex generic configurations is presented. The algorithm, boundary conditions, initial conditions, and grid generators are discussed as applied to these configurations. The PNS code was used as the mainline procedure to numerically simulate the viscous supersonic flow over these generic configurations. The turbulence model that was used in this study is the Baldwin-Lomax model. The boundary conditions are the usual viscous no slip at the wall, and a characteristic procedure is used to fit the bow shock wave which is the outermost boundary. An elliptic grid generator is employed to discretize the flow domain. In addition, an equilibrium air capability has been incorporated into the code. It uses the curve fits of Tannehill, et al. The flow regimes vary from a Mach number of 2 up to 25. Both laminar and turbulent flow are considered. Varying angles of attack have also been computed. Configurations vary from simple cone-type bodies to lifting winged bodies, such as the space shuttle or the generic supersonic cruise fighter.

Chaussee, D. S.↗
NASA NTRS1989

Non-linear evolution of a second mode wave in supersonic boundary layers

The nonlinear time evolution of a second mode instability in a Mach 4.5 wall-bounded flow is computed by solving the full compressible, time-dependent Navier-Stokes equations. High accuracy is achieved by using a Fourier-Chebyshev collocation algorithm. Primarily inviscid in nature, second modes are characterized by high frequency and high growth rates compared to first modes. Time evolution of growth rate as a function of distance from the plate suggests this problem is amenable to the Stuart-Watson perturbation theory as generalized by Herbert.

Erlebacher, Gordon↗
NASA NTRS1989

Progress on a Taylor weak statement finite element algorithm for high-speed aerodynamic flows

A new finite element numerical Computational Fluid Dynamics (CFD) algorithm has matured to the point of efficiently solving two-dimensional high speed real-gas compressible flow problems in generalized coordinates on modern vector computer systems. The algorithm employs a Taylor Weak Statement classical Galerkin formulation, a variably implicit Newton iteration, and a tensor matrix product factorization of the linear algebra Jacobian under a generalized coordinate transformation. Allowing for a general two-dimensional conservation law system, the algorithm has been exercised on the Euler and laminar forms of the Navier-Stokes equations. Real-gas fluid properties are admitted, and numerical results verify solution accuracy, efficiency, and stability over a range of test problem parameters.

Baker, A. J.↗
NASA NTRS1991

Reynolds number effects on supersonic asymmetrical flows over a cone at high angle of attack

The supersonic viscous flow over a 5-degree half-angle cone at an angle of attack of four times the cone half-angle is studied computationally using both the conical and the three-dimensional Navier-Stokes equations. The numerical solutions were obtained with an implicit, upwind-biased algorithm. Asymmetrical flowfields of the absolute-instability type are found using the conical-flow equations which agree with published results. However, the absolute instabilities of the originally symmetric flow found with the conical equations do not occur in the three-dimensional simulations, although spurious asymmetric three-dimensional flows for symmetric bodies arise if the grid resolution is insufficient in the nose region. The asymmetric flows computed with the three-dimensional equations are convective instabilities and are possible if the local Reynolds number exceeds a critical value and a fixed geometric asymmetry is imposed. A continuous range of asymmetries can be developed, depending on the size of the disturbance and the Reynolds number. As the Reynolds number is increased, the asymmetries demonstrate a bistable behavior at levels of side force consistent with those predicted using the conical equations. Below a certain critical Reynolds number, any flow asymmetries arising from a geometrical asymmetry are damped with increasing distance downstream from the geometrical asymmetry.

Thomas, J. L.↗
NASA NTRS1991

A new finite element formulation for computational fluid dynamics. X - The compressible Euler and Navier-Stokes equations

A space-time element method is presented for solving the compressible Euler and Navier-Stokes equations. The proposed formulation includes the variational equation, predictor multi-corrector algorithms and boundary conditions. The variational equation is based on the time-discontinuous Galerkin method, in which the physical entropy variables are employed. A least-squares operator and a discontinuity-capturing operator are added, resulting in a high-order accurate and unconditionally stable method. Implicit/explicit predictor multi-corrector algorithms, applicable to steady as well as unsteady problems, are presented; techniques are developed to enhance their efficiency. Implementation of boundary conditions is addressed; in particular, a technique is introduced to satisfy nonlinear essential boundary conditions, and a consistent method is presented to calculate boundary fluxes. Numerical results are presented to demonstrate the performance of the method.

Shakib, Farzin↗
NASA NTRS1993

Hypersonic flow past open cavities

The hypersonic flow over a cavity is investigated. The time-dependent compressible Navier-Stokes equations, in terms of mass averaged variables, are numerically solved. An implicit algorithm, with a subiteration procedure to recover time-accuracy, is used to perform the time-accurate computations. The objective of the study is to investigate the effects of Reynolds number and cavity dimensions. The comparison of the computations with available experimental data, in terms of time mean static pressure, heat transfer, and Mach number show good agreement. In the computations large vortex structures, which adversely affect the cavity flow characteristics, are observed at the rear of the cavity. A self-sustained oscillatory motion occurs within the cavity over a range of Reynolds number and cavity dimensions. The frequency spectra of the oscillations show good agreement with a modified semi-empirical relation.

Morgenstern, Alagacyr, Jr.↗
NASA NTRS1995

Computational analysis of forebody tangential slot blowing on the high alpha research vehicle

A numerical analysis of forebody tangential slot blowing as a means of generating side force and yawing moment is conducted using an aircraft geometry. The Reynolds-averaged, thin-layer, Navier-Stokes equations are solved using a partially flux-split, approximately-factored algorithm. An algebraic turbulence model is used to determine the turbulent eddy viscosity values. Solutions are obtained using both patched and overset grid systems. In the patched grid model, and actuator plane is used to introduce jet variables into the flow field. The overset grid model is used to model the physical slot geometry and facilitate modeling of the full aircraft configuration. A slot optimization study indicates that a short slot located close to the nose of the aircraft provided the most side force and yawing moment per unit blowing coefficient. Comparison of computed surface pressure with that obtained in full-scale wind tunnel tests produce good agreement, indicating the numerical method and grid system used in the study are valid. Full aircraft computations resolve the changes in vortex burst point due to blowing. A time-accurate full-aircraft solution shows the effect of blowing on the changes in the frequency of the aerodynamic loads over the vertical tails. A study of the effects of freestream Mach number and various jet parameters indicates blowing remains effective through the transonic Mach range. An investigation of the force onset time lag associated with forebody blowing shows the lag to be minimal. The knowledge obtained in this study may be applied to the design of a forebody tangential slot blowing system for use on flight aircraft.

Gee, Ken↗
NASA NTRS1993

Optimization of contoured hypersonic scramjet inlets with a least-squares parabolized Navier-Stokes procedure

A new optimization procedure, in which a parabolized Navier-Stokes solver is coupled with a non-linear least-squares optimization algorithm, is applied to the design of a Mach 14, laminar two-dimensional hypersonic subscale flight inlet with an internal contraction ratio of 15:1 and a length-to-throat half-height ratio of 150:1. An automated numerical search of multiple geometric wall contours, which are defined by polynomical splines, results in an optimal geometry that yields the maximum total-pressure recovery for the compression process. Optimal inlet geometry is obtained for both inviscid and viscous flows, with the assumption that the gas is either calorically or thermally perfect. The analysis with a calorically perfect gas results in an optimized inviscid inlet design that is defined by two cubic splines and yields a mass-weighted total-pressure recovery of 0.787, which is a 23% improvement compared with the optimized shock-canceled two-ramp inlet design. Similarly, the design procedure obtains the optimized contour for a viscous calorically perfect gas to yield a mass-weighted total-pressure recovery value of 0.749. Additionally, an optimized contour for a viscous thermally perfect gas is obtained to yield a mass-weighted total-pressure recovery value of 0.768. The design methodology incorporates both complex fluid dynamic physics and optimal search techniques without an excessive compromise of computational speed; hence, this methodology is a practical technique that is applicable to optimal inlet design procedures.

Korte, J. J.↗