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 217 records · Page 12

All sources NASA NTRS DOE OSTI

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↗

DOE OSTI2024

Note on two formulations of Crank-Nicolson method for Navier-Stokes equations

Here, we consider two formulations of the Crank-Nicolson (CN) method for the Navier-Stokes equations (NSE). The “natural” way of implementing CN for NSE is formally second order accurate in time for both velocity and pressure, whereas another formulation approximates pressure with only first order accuracy in time. Both versions of the method are applied to the benchmark problem of computing drag and lift in the flow around a cylinder. We show that the presumably more accurate version of the CN can create a solution with nonphysical oscillations and give incorrect predictions for the maximal drag coefficient, whereas the other formulation of the method predicts the drag and lift coefficients more accurately and does not introduce nonphysical oscillations. We locate the source of the issue and suggest several remedies.

Crank-Nicolson↗

DOE OSTI2021

An embedded variable step IMEX scheme for the incompressible Navier–Stokes equations

This report presents a series of implicit–explicit (IMEX) variable stepsize algorithms for the incompressible Navier–Stokes equations (NSE). Here, IMEX means the nonlinear term is treated fully explicitly while the remaining terms are treated implicitly. With the advent of new computer architectures there has been growing demand for low memory solvers of this type. The addition of time adaptivity improves the accuracy and greatly enhances the efficiency of the algorithm. We prove energy stability of an embedded first–second order IMEX pair. For the first order member of the pair, we prove stability for variable stepsizes, and analyze convergence. We believe this to be the first proof of this type for a variable stepsize IMEX scheme for the incompressible NSE. We then define and test a variable stepsize, variable order IMEX scheme using these methods. Our work contributes several firsts for IMEX NSE schemes, including an energy argument and error analysis of a two-step, variable stepsize method, and embedded error estimation for an IMEX multi-step method.

42 ENGINEERING↗

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.↗

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 NTRS2001

Time Integration Schemes for the Unsteady Navier-stokes Equations

The efficiency and accuracy of several time integration schemes are investigated for the unsteady Navier-Stokes equations. This study focuses on the efficiency of higher-order Runge-Kutta schemes in comparison with the popular Backward Differencing Formulations. For this comparison an unsteady two-dimensional laminar flow problem is chosen, i.e., flow around a circular cylinder at Re = 1200. It is concluded that for realistic error tolerances (smaller than 10(exp -1)) fourth-and fifth-order Runge-Kutta schemes are the most efficient. For reasons of robustness and computer storage, the fourth-order Runge-Kutta method is recommended. The efficiency of the fourth-order Runge-Kutta scheme exceeds that of second-order Backward Difference Formula by a factor of 2.5 at engineering error tolerance levels (10(exp -1) to 10(exp -2)). Efficiency gains are more dramatic at smaller tolerances.

Bijl, Hester↗