SEARCH · Engineering Papers

Results for “Stokes equations”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1988

A streamwise upwind algorithm for the Euler and Navier-Stokes equations applied to transonic flows

A new algorithm was developed for the Euler and Navier-Stokes equations that uses upwind differencing based on the streamwise direction. This algorithm is time accurate and can be used in codes for calculating unsteady transonic flows over wings. Such codes can be used for the flutter analysis of wings. In this algorithm, the coordinate system is locally rotated to align with the streamwise direction. For differencing the convective terms in the streamwise direction, a new form of flux splitting is employed, in which the biasing depends on the local Mach number. In the plane perpendicular to the stream direction, the new flux splitting uses the condition of no flow in that local plane. By using a locally rotated coordinate system, the convective flux vector biasing depends on the total Mach number. Hence, the switching of the flux vector biasing occurs across shock waves and the proper domain of dependence is used in supersonic regions. For comparison, many other upwind methods switch differencing based on Mach number of shock wave in multidimensional flows. The formulas for the convective flux vector differencing do not contain any user specified parameters. So, the amount of numerical dissipation is automatically determined.

Goorjian, P. M.↗

NASA NTRS1989

Spectral solution of the incompressible Navier-Stokes equations on the Connection Machine2

The issue of solving the time-dependent incompressible Navier-Stokes equations on the Connection Machine 2 is addressed, for the problem of transition to turbulence of the steady flow in a channel. The spectral algorithm used serially requires O(N4) operations when solving the equations on an N x N x N grid; using the massive parallelism of the CM, it becomes an O(N2) problem. Preliminary timings of the code, written in LISP, are included and compared with a corresponding code optimized for the Cray-2 for a 128 x 128 x 101 grid.

Tomboulian, Sherryl↗

NASA NTRS1989

Unstructured mesh solution of the Euler and Navier-Stokes equations

Mesh generation procedures as well as solution algorithms for solving the Euler and Navier-Stokes equations on unstructured meshes are presented. The solution algorithms discussed utilize approximate Riemann solver, upwind differencing to achieve high spatial accuracy. Numerical results for Euler flow over single and multi-element airfoils are presented.

Barth, Timothy J.↗

NASA NTRS1991

Spectral methods for the Navier-Stokes equations with one infinite and two periodic directions

The time-dependent, 3D incompressible Navier-Stokes equations in (1) boundary layers, the semiinfinite domain, and (2) mixing layers or wakes, the fully infinite domain, are respectively solved by two numerical methods which employ rapidly decaying spectral basis functions to approximate the vertical dependence of the solutions. These functions are then combined with one, for method (1), and two, for method (2), slowly decaying 'extra functions' for each wave vector. Each extra function can exactly represent the solution's irrotational component at large distances. The two methods have been applied to extensive direct-numerical simulation of transition and turbulence.

Spalart, Philippe R.↗

NASA NTRS1992

Investigation of vortex breakdown on delta wings using Navier-Stokes equations

An efficient finite-difference scheme solving for the three-dimensional incompressible Navier-Stokes equations is described. Numerical simulations of vortex breakdown are then carried out for a sharp-edged delta wing and a round-edged double-delta wing at high Reynolds numbers. Computed results show that several major features of vortex breakdown are qualitatively in agreement with observations made in experiments.

Hsu, C.-H.↗

NASA NTRS1995

Fractional Step and Pseudocompressibility Methods for the Solution of Incompressible Navier-Stokes Equations

The fractional step and the pseudocompressibility methods for the solution of the incompressible Navier-Stokes equations are outlined. The fractional step method is based on finite-volume formulation and uses the pressure and the volume fluxes across the faces of each cell as dependent variables. The momentum equations are solved implicitly and the Poisson equation for the pressure is solved by using the multigrid method. The pseudocompressibility approach uses an implicit-higher-order-upwind differencing scheme for the convective terms together with the Gauss-Seidel line relaxation method. The dependent variables in the pseudocompressibility approach are the pressure and the cartesian velocity components in unstaggered mesh orientation. The 90-degree square duct flow, the wing-tip vortex wake flow and unsteady turbulent flows over an oscillating NACA 0015 airfoil are computed using both the fractional step and the pseudocompressibility methods. The results obtained from two different schemes are compared against experimental measurements.

Kiris, Cetin↗

NASA NTRS1991

A multiple-block multigrid method for the solution of the three-dimensional Euler and Navier-Stokes equations

A multiple block multigrid method for the solution of the three dimensional Euler and Navier-Stokes equations is presented. The basic flow solver is a cell vertex method which employs central difference spatial approximations and Runge-Kutta time stepping. The use of local time stepping, implicit residual smoothing, multigrid techniques and variable coefficient numerical dissipation results in an efficient and robust scheme is discussed. The multiblock strategy places the block loop within the Runge-Kutta Loop such that accuracy and convergence are not affected by block boundaries. This has been verified by comparing the results of one and two block calculations in which the two block grid is generated by splitting the one block grid. Results are presented for both Euler and Navier-Stokes computations of wing/fuselage combinations.

Atkins, Harold↗

NASA NTRS1975

Use of numerically generated body-fitted coordinate systems for solution of the Navier-Stokes equations

A procedure for numerical solution of the time-dependent, two-dimensional incompressible Navier-Stokes equations that can treat the unsteady laminar flow about bodies of arbitrary shape, such as two-dimensional airfoils, multiple airfoils, and submerged hydrofoils, as naturally as it can deal with the flow about simple bodies. The solution is based on a method of automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multiconnected region containing any number of arbitrarily shaped bodies. The curvilinear coordinates are generated as the solution of two elliptical partial differential equations with Dirichlet boundary conditions, one coordinate being specified to be constant on each of the boundaries, and a distribution of the other being specified along the boundaries. The solution compares excellently with the Blasius boundary layer solution for the flow past a semiinfinite flat plate.

Thompson, J. F.↗

NASA NTRS1977

An implicit factored scheme for the compressible Navier-Stokes equations

An implicit finite difference scheme is developed for the numerical solution of the compressible Navier-Stokes equations in conservation-law form. The algorithm is second-order-time accurate, noniterative, and spatially factored. In order to obtain an efficient factored algorithm, the spatial cross-derivatives are evaluated explicitly. However, the algorithm is unconditionally stable and, although a three-time-level scheme, requires only two-time-levels of data storage. The algorithm is constructed in a 'delta' form (i.e., increments of the conserved variables and fluxes) that provides a direct derivation of the scheme and leads to an efficient computational algorithm. In addition, the delta form has the advantageous property of a steady-state (if one exists) independent of the size of the time step. Numerical results are presented for a two-dimensional shock boundary-layer interaction problem.

Beam, R. M.↗

NASA NTRS1978

Numerical solution of the Navier-Stokes equations for super-sonic flows with strong shocks

The numerical solution of the full Navier-Stokes Equations for viscous flows with high Mach numbers and a strong detached bow shock was obtained. Two dimensional flows around a circular cylinder, and a circular cylinder with an aft-body in the form of a fairing, were considered. The solution of the compressible N.S. equations was accomplished by the method of finite differences. An implicit scheme of solution, the S.O.R., was used with the optimum acceleration parameters determined by trial and error. The tensor notation was used in writing the N-S Equations transformed into general curvilinear coordinates. The equations for the generation of the coordinate system were solved, followed by the solution of the N.S. equations, at the end of a set of given number of time steps. "Wiggles", constituted the one major problem that needed to be overcome. These oscillations give rise to quantities such as negative temperatures, which ultimately caused the computational program to break down. Certain dissipative finite-difference schemes damped these oscillations.

Devarayalu, K.↗

NASA NTRS1978

Numerical solution of the Navier-Stokes equations for arbitrary blunt bodies in supersonic flows

A time-dependent, two-dimensional Navier-Stokes code employing the body-fitted coordinate technique has been developed for supersonic flows past blunt bodies of arbitrary shape. The computer program is based on the finite-difference approximation of the compressible Navier-Stokes equations transformed to nonorthogonal curvilinear coordinates with the contravariant components of the velocity vector as dependent variables. The bow shock ahead of the body is obtained as part of the solution, by 'shock capturing'. Numerical solutions of the complete equations are presented in detail for free-stream Mach number 4.6, Reynolds number 10,000, and an isothermal wall temperature of 556 K for a circular cylinder with the free-stream outer boundaries forming a hyperbola in the front and a circular arc in the back.

Warsi, Z. U. A.↗

NASA NTRS1984

Evaluation of a velocity-split solution of the Navier-Stokes equations for fighter forebodies

A velocity-splitting method of solving the three-dimensional Navier-Stokes equations, originally designed for afterbody flows, is examined for its applicability for predicting fighter forebody flows. Results from the AFTEND Code are compared with wind tunnel data for two fighter configurations at a Mach number of 0.9 and angles-of-attack from 0 deg to 20 deg. Results compare well with data, and in areas where data do not exist, the viscous AFTEND results show realistic effects of viscosity compared to inviscid predictions. The inviscid results themselves are in general superior to results obtained from a small-disturbance transonic potential code.

Yaros, S. F.↗

NASA NTRS1985

On the accuracy of the pseudocompressibility method in solving the incompressible Navier-Stokes equations

The method of pseudocompressibility is being tested for its accuracy in solving the incompressible Navier-Stokes equations. An implicit, finite-difference computer code is used to solve the equations in a three-dimensional, curvilinear coordinate system. The code employs artificial compressibility for solving the pressure field, coupled with an implicit, approximate-factorization scheme. This coupling is known as the pseudocompressibility method. The pseudocompressibility method introduces pressure waves of finite speed into the fluid medium that would otherwise have an infinite sound speed. The waves die out as the solution converges, and the steady state solution approaches a divergence-free condition. However, these waves limit the time-accuracy of the computations. The effects of these waves are analyzed and criteria are set for choosing the pseudocompressibility parameters that govern the pressure wave speed. Test cases are presented that verify these criteria. The code is tested by computing laminar flow over a two-dimensional, backward-facing step and over a two-dimensional, impulsively started circular cylinder.

Rogers, S. E.↗

NASA NTRS1984

Multi-grid technique for the solution of incompressible Navier-Stokes equations

A computationally efficient method is proposed for obtaining fine-mesh solutions of the Navier-Stokes equations for high-Re flow, including separation. The method involves implementing a multi-grid solution procedure with suitably chosen elements, such that the resulting overall solution technique is efficient as well as robust. The robustness and efficiency of the solution technique are demonstrated by applying it to three model problems: flow in a driven cavity, downstream asymptotic flow in curved ducts of square and polar sections, and Newmann boundary-value problem in clustered curvilinear orthogonal coordinates.

Ghia, K. N.↗

NASA NTRS1986

About the coupling of turbulence closure models with averaged Navier-Stokes equations

The MacCormack implicit predictor-corrector model (1981) for numerical solution of the coupled Navier-Stokes equations for turbulent flows is extended to nonconservative multiequation turbulence models, as well as the inclusion of second-order Reynolds stress turbulence closure. A scalar effective pressure turbulent contribution to the pressure field is defined to approximate the effects of the Reynolds stress in strongly sheared flows. The Jacobian matrices of the transport equations are diagonalized to reduce the required computer memory and run time. Techniques are defined for including turbulence in the diagonalization. Application of the method is demonstrated with solutions generated for transonic nozzle flow and for the interaction between a supersonic flat plate boundary layer and a 12 deg compression-expansion ramp.

Vandromme, D.↗

NASA NTRS1988

On the smallest scale for the incompressible Navier-Stokes equations

It is proven that for solutions to the two- and three-dimensional incompressible Navier-Stokes equations the minimum scale is inversely proportional to the square root of the Reynolds number based on the kinematic viscosity and the maximum of the velocity gradients. The bounds on the velocity gradients can be obtained for two-dimensional flows, but have to be assumed to be three-dimensional. Numerical results in two dimensions are given which illustrate and substantiate the features of the proof. Implications of the minimum scale result to the decay rate of the energy spectrum are discussed.

Henshaw, W. D.↗

NASA NTRS1989

Spectral solution of the incompressible Navier-Stokes equations on the Connection Machine 2

The issue of solving the time-dependent incompressible Navier-Stokes equations on the Connection Machine 2 is addressed, for the problem of transition to turbulence of the steady flow in a channel. The spectral algorithm used serially requires O(N(4)) operations when solving the equations on an N x N x N grid; using the massive parallelism of the CM, it becomes an O(N(2)) problem. Preliminary timings of the code, written in LISP, are included and compared with a corresponding code optimized for the Cray-2 for a 128 x 128 x 101 grid.

Tomboulian, Sherryl↗

NASA NTRS1992

Direct Coupling Method for Time-Accurate Solution of Incompressible Navier-Stokes Equations

A noniterative finite difference numerical method is presented for the solution of the incompressible Navier-Stokes equations with second order accuracy in time and space. Explicit treatment of convection and diffusion terms and implicit treatment of the pressure gradient give a single pressure Poisson equation when the discretized momentum and continuity equations are combined. A pressure boundary condition is not needed on solid boundaries in the staggered mesh system. The solution of the pressure Poisson equation is obtained directly by Gaussian elimination. This method is tested on flow problems in a driven cavity and a curved duct.

Soh, Woo Y.↗