Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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.

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

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Numerical solution of the compressible Navier-Stokes equations using density gradients as additional dependent variables

Numerical solution of two dimensional, time dependent, compressible viscous Navier-Stokes equations about arbitrary bodies was treated using density gradients as additional dependent variables. Thus, six dependent variables were computed with the SOR iteration method. Besides formulation for pressure gradient terms, a formulation for computing the body density was presented. To approximate the governing equations, an implicit finite difference method was employed. In computing the solution for the flow about a circular cylinder, a problem arose near the wall at both stagnation points. Thus, computations with various conditions were tried to examine the problem. Also, computations with and without formulations are compared. The flow variables were computed on 37 by 40 field first, then on an 81 by 40 field.

Kwon, J. H.↗
NASA NTRS1978

Forebody and afterbody solutions of the Navier-Stokes equations for supersonic flow over blunt bodies in a generalized orthogonal coordinate system

A coordinate transformation, which can approximate many different two-dimensional and axisymmetric body shapes with an analytic function, is used as a basis for solving the Navier-Stokes equations for the purpose of predicting 0 deg angle of attack supersonic flow fields. The transformation defines a curvilinear, orthogonal coordinate system in which coordinate lines are perpendicular to the body and the body is defined by one coordinate line. This system is mapped in to a rectangular computational domain in which the governing flow field equations are solved numerically. Advantages of this technique are that the specification of boundary conditions are simplified and, most importantly, the entire flow field can be obtained, including flow in the wake. Good agreement has been obtained with experimental data for pressure distributions, density distributions, and heat transfer over spheres and cylinders in supersonic flow. Approximations to the Viking aeroshell and to a candidate Jupiter probe are presented and flow fields over these shapes are calculated.

Gnoffo, P. A.↗
NASA NTRS1978

A numerical solution of the Navier-Stokes equations for chemically nonequilibrium, merged stagnation shock layers on spheres and two-dimensional cylinders in air

Results of solving the Navier-Stokes equations for chemically nonequilibrium, merged stagnation shock layers on spheres and two-dimensional cylinders are presented. The effects of wall catalysis and slip are also examined. The thin shock layer assumption is not made, and the thick viscous shock is allowed to develop within the computational domain. The results show good comparison with existing data. Due to the more pronounced merging of shock layer and boundary layer for the sphere, the heating rates for spheres become higher than those for cylinders as the altitude is increased.

Johnston, K. D.↗
NASA NTRS1978

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

A time dependent, two dimensional Navier-Stokes code employing the method of body fitted coordinate technique was developed for supersonic flows past blunt bodies of arbitrary shapes. The bow shock ahead of the body is obtained as part of the solution, viz., by shock capturing. A first attempt at mesh refinement in the shock region was made by using the forcing function in the coordinate generating equations as a linear function of the density gradients. The technique displaces a few lines from the neighboring region into the shock region. Numerical calculations for Mach numbers 2 and 4.6 and Reynolds numbers from 320 to 10,000 were performed for a circular cylinder with and without a fairing. Results of Mach number 4.6 and Reynolds number 10,000 for an isothermal wall temperature of 556 K are presented in detail.

Warsi, Z. U. A.↗
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 NTRS1977

A Navier-Stokes solution of the three-dimensional viscous compressible flow in a centrifugal compressor impeller

A two-dimensional time-dependent computer code was utilized to calculate the three-dimensional steady flow within the impeller blading. The numerical method is an explicit time marching scheme in two spatial dimensions. Initially, an inviscid solution is generated on the hub blade-to-blade surface by the method of Katsanis and McNally (1973). Starting with the known inviscid solution, the viscous effects are calculated through iteration. The approach makes it possible to take into account principal impeller fluid-mechanical effects. It is pointed out that the second iterate provides a complete solution to the three-dimensional, compressible, Navier-Stokes equations for flow in a centrifugal impeller. The problems investigated are related to the study of a radial impeller and a backswept impeller.

Harp, J. L., Jr.↗
NASA NTRS1977

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of the nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The continuous cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and symptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Cooke, C. H.↗
NASA NTRS1978

Navier-Stokes calculations for laminar and turbulent hypersonic flow over indented nosetips

A time-accurate finite-difference Navier-Stokes code has been used to calculate the viscous flow over a severely indented blunt body in a supersonic stream. An algebraic turbulence model is used and the results are compared with experimental data from wind-tunnel tests. Qualitative agreement is obtained for the surface pressure distribution and flow-field structure, including the separated bubble in the indented region. However, uncertainties still exist in the heating calculations, which are attributed to the turbulence model. For both laminar and turbulent calculations, the flow exhibits a fundamental unsteady character at a frequency of about 50 kHz.

Rakich, J. V.↗
NASA NTRS1978

Approximate solutions to the Navier-Stokes initial value problem

A Galerkin-Ritz procedure for any arbitrary system of field equations is shown to follow generically from 'two-functional' variational conditions. The initial-value problem for boundary-free incompressible Navier-Stokes flow is solved analytically in the two-parameter-function approximation.

Rosen, G.↗
NASA NTRS1977

Navier-Stokes solution for the laminar and turbulent supersonic flow near the trailing edge of a flat plate

The present paper deals with problems of two-dimensional laminar and turbulent flow of a viscous incompressible fluid near the trailing edge of a thin flat plate. The complete system of Navier-Stokes equations is solved by means of MacCormack's finite-difference scheme. In the laminar case, the results are compared with the triple deck solution of Daniels (1973). The comparisons indicate that the asymptotic triple deck theory for supersonic trailing-edge flow is accurate within five percent for Reynolds numbers greater than 1000. In the turbulent case, the Prandtl-Van Driest-Clausner algebraic eddy viscosity model is used, showing that the region of upstream influence is roughly of the order of the boundary layer thickness.

Hussaini, M. Y.↗
NASA NTRS1978

Navier-Stokes solutions using Stetter's method

Compressible Navier-Stokes equations for quasi-one-dimensional flow in a converging-diverging nozzle have been solved using Stetter's three-step predictor-corrector technique. Particular emphasis is given to the minimum iterative step feature for steady-state solutions. It is found that for the nonviscous-dominated case, Stetter's method attains the steady-state solution in the fewest steps when compared to four other currently used techniques.

Graves, R. A., Jr.↗
NASA NTRS1978

Development of a three-dimensional Navier-Stokes code on CDC star-100 computer

A three-dimensional code in body-fitted coordinates was developed using MacCormack's algorithm. The code is structured to be compatible with any general configuration, provided that the metric coefficients for the transformation are available. The governing equations are developed in primitive variables in order to facilitate the incorporation of physical boundary conditions and turbulence-closure models. MacCormack's two-step, unsplit, time-marching algorithm is used to solve the unsteady Navier-Stokes equations until steady-state solution is achieved. Cases discussed include (1) flat plate in supersonic free stream; (2) supersonic flow along an axial corner; (3) subsonic flow in an axial corner at M infinity = 0.95; and (4) supersonic flow in an axial corner at M infinity 1.5.

Vatsa, V. N.↗
NASA NTRS1978

Navier-Stokes calculations with a coupled strongly implicit method. Part 2: Spline solutions

A coupled strongly implicit method is combined with a deferred-corrector spline solver for the vorticity-stream function form of the Navier-Stokes equation. Solutions for cavity, channel and cylinder flows are obtained with the fourth-order spline 4 procedure. The strongly coupled spline corrector method converges as rapidly as the finite difference calculations and also allows for arbitrary large time increments for the Reynolds numbers considered. In some cases fourth-order smoothing or filtering is required in order to suppress high frequency oscillations.

Rubin, S. G.↗
NASA NTRS1979

A split finite element algorithm for the compressible Navier-Stokes equations

An accurate and efficient numerical solution algorithm is established for solution of the high Reynolds number limit of the Navier-Stokes equations governing the multidimensional flow of a compressible essentially inviscid fluid. Finite element interpolation theory is used within a dissipative formulation established using Galerkin criteria within the Method of Weighted Residuals. An implicit iterative solution algorithm is developed, employing tensor product bases within a fractional steps integration procedure, that significantly enhances solution economy concurrent with sharply reduced computer hardware demands. The algorithm is evaluated for resolution of steep field gradients and coarse grid accuracy using both linear and quadratic tensor product interpolation bases. Numerical solutions for linear and nonlinear, one, two and three dimensional examples confirm and extend the linearized theoretical analyses, and results are compared to competitive finite difference derived algorithms.

Baker, A. J.↗
NASA NTRS1979

Numerical solution of the Navier-Stokes equations for laminar, transonic flows

An implicit finite difference solution of the Navier-Stokes equations yielded time histories of the transonic laminar flow development about a circular cylinder and NACA-0018 airfoil. Reynolds numbers ranged from those corresponding to purely laminar flow to those corresponding to significant turbulence in the boundary layer. Body thermal conditions of an adiabatic wall and a specified body temperature were considered. Versatility in treating arbitrary bodies was incorporated by using numerically generated, body-fitted coordinate transformations. Solution of the simultaneous difference equations for the dependent variables was obtained using an accelerated Gauss-Seidel iterative scheme. Computational results are presented in the form of velocity vector fields, Mach number contours, aerodynamic coefficients, heat transfer rates at the body surface, and body temperature distributions. Truncation analyses of first and second derivative difference approximations resulted in general criteria for numerically generated coordinates so that flow near a body is more accurately represented.

Turner, L., III↗
NASA NTRS1979

A compressible solution of the Navier-Stokes equations for turbulent flow about an airfoil

A compressible time dependent solution of the Navier-Stokes equations including a transition turbulence model is obtained for the isolated airfoil flow field problem. The equations are solved by a consistently split linearized block implicit scheme. A nonorthogonal body-fitted coordinate system is used which has maximum resolution near the airfoil surface and in the region of the airfoil leading edge. The transition turbulence model is based upon the turbulence kinetic energy equation and predicts regions of laminar, transitional, and turbulent flow. Mean flow field and turbulence field results are presented for an NACA 0012 airfoil at zero and nonzero incidence angles of Reynolds number up to one million and low subsonic Mach numbers.

Shamroth, S. J.↗
NASA NTRS1979

Numerical solution of the azimuthal-invariant thin-layer Navier-Stokes equations

The paper reports a numerical procedure developed for a two-dimensional azimuthal (or planar) invariant form of the thin-layer Navier-Stokes equations. Generalization of the governing equations is described, and the equations are solved with an implicit approximate factorization finite-difference scheme. Inviscid and viscous results are presented for both external and internal flows and for spinning and nonspinning bodies.

Nietubicz, C. J.↗