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 595 records · Page 33

Numerical solutions of Navier-Stokes equations for a Butler wing

The flow field is simulated on the surface of a given delta wing (Butler wing) at zero incident in a uniform stream. The simulation is done by integrating a set of flow field equations. This set of equations governs the unsteady, viscous, compressible, heat conducting flow of an ideal gas. The equations are written in curvilinear coordinates so that the wing surface is represented accurately. These equations are solved by the finite difference method, and results obtained for high-speed freestream conditions are compared with theoretical and experimental results. In this study, the Navier-Stokes equations are solved numerically. These equations are unsteady, compressible, viscous, and three-dimensional without neglecting any terms. The time dependency of the governing equations allows the solution to progress naturally for an arbitrary initial initial guess to an asymptotic steady state, if one exists. The equations are transformed from physical coordinates to the computational coordinates, allowing the solution of the governing equations in a rectangular parallel-piped domain. The equations are solved by the MacCormack time-split technique which is vectorized and programmed to run on the CDC VPS 32 computer.

Abolhassani, J. S.↗

Computational aspects of zonal algorithms for solving the compressible Navier-Stokes equations in three dimensions

Transonic flow fields about wing geometries are computed using an Euler/Navier-Stokes approach in which the flow field is divided into several zones. The flow field immediately adjacent to the wing surface is resolved with fine grid zones and solved using a Navier-Stokes algorithm. Flow field regions removed from the wing are resolved with less finely clustered grid zones and are solved with an Euler algorithm. Computational issues associated with this zonal approach, including data base management aspects, are discussed. Solutions are obtained that are in good agreement with experiment, including cases with significant wind tunnel wall effects. Additional cases with significant shock induced separation on the upper wing surface are also presented.

Holst, T. L.↗

A fourth-order scheme for the unsteady compressible Navier-Stokes equations

A computational scheme is described which is second-order accurate in time and fourth-order accurate in space (2-4). This method is applied to study the stability of compressible boundary layers. The laminar compressible Navier-Stokes equations are solved with a time harmonic inflow superimposed on the steady state solution. This results in spatially unstable modes. It is shown that the second-order methods are inefficient for calculating the growth rates and phases of the unstable modes. In contrast the fourth-order method yields accurate results on relatively course meshes.

Bayliss, A.↗

Marching iterative methods for the parabolized and thin layer Navier-Stokes equations

Downstream marching iterative schemes for the solution of the Parabolized or Thin Layer (PNS or TL) Navier-Stokes equations are described. Modifications of the primitive equation global relaxation sweep procedure result in efficient second-order marching schemes. These schemes take full account of the reduced order of the approximate equations as they behave like the SLOR for a single elliptic equation. The improved smoothing properties permit the introduction of Multi-Grid acceleration. The proposed algorithm is essentially Reynolds number independent and therefore can be applied to the solution of the subsonic Euler equations. The convergence rates are similar to those obtained by the Multi-Grid solution of a single elliptic equation; the storage is also comparable as only the pressure has to be stored on all levels. Extensions to three-dimensional and compressible subsonic flows are discussed. Numerical results are presented.

Israeli, M.↗

Steady supersonic Navier-Stokes solutions of a 75 deg delta wing

Steady solutions about a slender sharp edged delta wing in a supersonic freestream for moderate and high angles of attack are obtained numerically by time integration of the unsteady compressible three dimensional laminar Navier-Stokes equations. The main features of the flow, including primary and secondary separation, and vortex position and strength, are adequately simulated in the numerical solutions. Improved resolution of the computational grid in the leading edge region from a previous solution had considerable effect on the accuracy of the solutions. Good agreement between numerical solutions and experimental data was obtained for two cases. A local timestepping procedure is used to speed convergence by approximately a factor of two.

Buter, T. A.↗

Efficient solution methods for the Navier-Stokes equations

Implicit finite difference schemes for solving two-dimensional and three-dimensional Euler and thin layer Navier-Stokes equations are addressed. The methods are demonstrated in fully vectorized codes for a Cray type architecture. The Beam and Warming implicit approximate factorization algorithm in generalized coordinates is used. The methods are either time accurate or accelerated non-time accurate steady state schemes. Acceleration and efficiency modifications such as matrix reduction, diagonalization, and flux split schemes are presented. Two dimensional inviscid and viscous calculations (e.g., airfoils with a deflected spoiler, circulation control airfoils, and unsteady buffeting) and of three dimensional viscous elliptical bodies, exhausting boattails, and generic oblique wing computations are discussed.

Pulliam, T. H.↗

Navier-Stokes computations for exotic airfoils

A hyperbolic grid generator with improvements for handling sharp corners and concave surfaces combined with a Navier-Stokes flow solver is applied to complex two dimensional airfoil configurations. Steady separated flow about an iced leading edge of an airfoil is presented. Unsteady viscous separated flows past an airfoil at two angles of attack with a spoiler deployed at 60 deg are compared with experiment. The spoiler computations are performed with two different topological maps of the physical domain to the computational domain. Graphical techniques for the static and unsteady display of the flow fields are discussed.

Barth, T. J.↗

Patched-grid calculations with the Euler and Navier-Stokes equations: Theory and application

The Rai (1984,85) patch-boundary scheme for the Euler equations is described. The integration methods used to update the interior grid points are are discussed. Stability of patch-boundary schemes and the use of these schemes in Navier-Stokes calculations are mentioned. Results for inviscid, supersonic flow over a cylinder, blast wave diffraction by ramp, and the motion of a vortex in a freestream are presented. These test cases demonstrate the quality of solutions possible with the scheme.

Rai, M. M.↗

A spectral collocation method for the Navier-Stokes equations

A Fourier-Chebyshev spectral method for the incompressible Navier-Stokes equations is described. It is applicable to a variety of problems including some with fluid properties which vary strongly both in the normal direction and in time. In this fully spectral algorithm, a preconditioned iterative technique is used for solving the implicit equations arising from semi-implicit treatment of pressure, mean advection and vertical diffusion terms. The algorithm is tested by applying it to hydrodynamic stability problems in channel flow and in external boundary layers with both constant and variable viscosity.

Malik, M. R.↗

A comparative study of the parabolized Navier-Stokes code using various grid-generation techniques

The parabolized Navier-Stokes (PNS) equations are used to calculate the flow-field characteristics about the hypersonic research aircraft X-24C. A comparison of the results obtained using elliptic, hyperbolic and algebraic grid generators is presented. The outer bow shock is treated as a sharp discontinuity, and the discontinuities within the shock layer are captured. Surface pressures and heat-transfer results at angles of attack of 6 deg and 20 deg, obtained using the three grid generators, are compared. The PNS equations are marched downstream over the body in both Cartesian and cylindrical base coordinate systems, and the results are compared. A robust marching procedure is demonstrated by successfully using large marching-step sizes with the implicit shock fitting procedure. A correlation is found between the marching-step size, Reynolds number and the angle of attack at fixed values of smoothing and stability coefficients for the marching scheme.

Kaul, U. K.↗

Parabolized Navier-Stokes analysis of ducted turbulent mixing problems with finite-rate chemistry

A hybrid implicit/explicit approach for analyzing 2D ducted turbulent mixing problems with finite-rate chemistry is presented. The approach combines a fully-implicit parabolic mixing algorithm, an explicit viscous-characteristic based shock-capturing algorithm, and linearized implicit chemical kinetic algorithm. The resultant model provides spatial marching solutions of the parabolized Navier-Stokes (PNS) equations for supersonic combustion and viscous nozzle flowfields. Specialized procedures are incorporated to deal with the near wall sublayer and with small central pockets of subsonic flow. A parabolic option is available which can be utilized in the direct or inverse mode for design applications. The ability of the model to treat shock waves and wave/mixing layer interactions is assessed via comparisons of predictions with those of well established explicit shock-capturing Euler (SCIPPY) and PNS (SCIPVIS) models. Applications to a variety of supersonic combustion flowfield problems involving tangential or moderately inclined fuel injection, and, to a viscous nozzle flow problem, are presented. This paper serves to exhibit overall capabilities of the model developed and no comparisons with supersonic combustion data are presented. Such comparisons serve mainly to validate the modeling of the turbulence and turbulence/chemical interactions and will be the subject of a future paper.

Sinha, N.↗

A new approximate LU factorization scheme for the Reynolds-averaged Navier-Stokes equations

A new approximate LU factorization scheme is developed to solve the steady state Reynolds-averaged Navier-Stokes (NS) equations. Central differencing is used for both implicit and explicit operator and special care is taken to obtain well-conditioned factors on the implicit side. The scheme is then analyzed and optimized according to a simple linear analysis. It is unconditionally stable for the model hyperbolic equation in both two- and three-dimensions. However, the requirement for well-conditioned factors has essentially limited the effective time step the scheme can achieve. Supersonic and transonic three-dimensinal flows past a hemisphere cylinder are computed to demonstrate the convergence characteristics of the scheme. A good convergence rate is achieved for the inviscid case. Finally, an explicit eigenvector annihilation procedure is adopted successfully to remove the stiffness caused by the fine grid spacing for viscous flows.

Pan, D.↗

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

Calculations of two and three-dimensional transonic cascade flow fields using the Navier-Stokes equations

The multidimensional, ensemble-averaged, compressible, time-dependent Navier-Stokes equations have been used to study the turbulent flow field in two and three-dimensional turbine cascades. The viscous regions of the flow were resolved and non-slip boundary conditions were utilized on solid surfaces. The calculations were performed in a constructive 'O'-type grid which allows representation of the blade rounded trailing edge. Converged solutions were obtained in relatively few time steps (about 80-150) and comparisons for both surface pressure and heat transfer showed good agreement with data. The three-dimensional turbine cascade calculation showed many of the expected flow-field features.

Weinberg, B. C.↗

An implicit flux-difference splitting scheme for three-dimensional, incompressible Navier-Stokes solutions to leading edge vortex flows

A new, implicit finite-difference scheme designed to solve the conservative, flux-difference split Navier-Stokes equations is used to compute incompressible vortex flows around delta wings. The completely vectorizable hybrid algorithm is constructed in delta form for steady state solutions independent of the time-step sizes. The scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. The governing equations are solved in curvilinear, body-fitted coordinates for treating complex geometries. The computed flow field results are compared with other theoretical and experimental data.

Hartwich, P.-M.↗

Navier-Stokes computations of lee-side flows over delta wings

Solutions to the Navier-Stokes equations for the flow over delta wings are computed with emphasis on the separated vortical flows developing on the lee side at high angles of attack. A recently developed implicit algorithm is used which employs upwind differencing for the pressure and convection terms and central differencing for the shear stress and heat transfer terms. Solutions to both the three-dimensional equations and the approximate conical flow equations are compared parametrically with an extensive experimental data base at supersonic speeds. The computations indicate that the conical flow approximation provides results in close agreement with the three-dimensional equations, even to angles of attack as high as 20 degrees. Good agreement with experimentally measured pressures and vapor screen photographs is obtained for the conditions investigated. The method predicts the classical pattern of vortical flow over a delta wing and transition to other flow patterns as the leading edge sweep angle and leading edge normal Mach number are varied.

Thomas, J. L.↗

An upwind algorithm for the parabolized Navier-Stokes equations

A new upwind algorithm based on Roe's scheme has been developed to solve the two-dimensional parabolized Navier-Stokes (PNS) equations. This method does not require the addition of user specified smoothing terms for the capture of discontinuities such as shock waves. Thus, the method is easy to use and can be applied without modification to a wide variety of supersonic flowfields. The advantages and disadvantages of this adaptation are discussed in relation to those of the conventional Beam-Warming scheme in terms of accuracy, stability, computer time and storage, and programming effort. The new algorithm has been validated by applying it to three laminar test cases including flat plate boundary-layer flow, hypersonic flow past a 15 deg compression corner, and hypersonic flow into a converging inlet. The computed results compare well with experiment and show a dramatic improvement in the resolution of flowfield details when compared with the results obtained using the conventional Beam-Warming algorithm.

Lawrence, S. L.↗

Study of hypersonic inlet flow fields with a three-dimensional Navier-Stokes code

A three-dimensional Navier-Stokes code is used to study the flow through an inlet with a rectangular capture area, but with circular or essentially circular internal cross-sections. The inlet may have a center body extending partially or fully over the length of the inlet. Calculations are made under a variety of geometrical constraints imposed on the inlet cross-sections to study the flow quality in the inlet. It is observed that for this class of inlets, the boundary layer has a tendency to thicken and separate readily due to the adverse pressure gradient caused by the continuous compression in the inlet. Since no experimental results are yet available, this paper has presented only the numerical results from the present analysis.

Kumar, A.↗