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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A solution-adaptive grid procedure for the three-dimensional parabolized Navier-Stokes equations

A solution-adaptive grid procedure for the three-dimensional parabolized Navier-Stokes equations is developed on the basis of the two-dimensional line-by-line solution-adaptive technique of Harvey et al. (1990). It is shown that the new adaptive grid algorithm provides improved shock resolving characteristics over the conventional flow algorithm. The adaptation process has the ability of aligning grid lines with the existing flow-field structure, yielding increasing resolution of high gradient regions.

Harvey, Albert D., III↗
NASA NTRS1991

A grid-independent approximate Riemann solver with applications to the Euler and Navier-Stokes equations

A new two-dimensional approximate Riemann solver has been developed that obtains fluxes on grid faces via wave decomposition. By utilizing information propagation in the velocity-difference directions rather than in the grid-normal directions, this flux function more appropriately interprets and hence more sharply resolves shock and shear waves when they lie oblique to the grid. The model uses five waves to describe the difference in states at a grid face. Two acoustic waves, one shear wave, and one entropy wave propagate in the direction defined by the local velocity difference vector, while the fifth wave is a shear wave that propagates at a right angle to the other four. Test cases presented include a shock reflecting off a wall, a pure shear wave, supersonic flow over an airfoil, and viscous separated airfoil flow. Results using the new model give significantly sharper shock and shear contours than a grid-aligned solver. Navier-Stokes computations over an aifoil show reduced pressure distortions in the separated region as a result of the grid-independent upwinding.

Rumsey, Christopher L.↗
NASA NTRS1991

Jet noise predictions from unsteady Navier-Stokes simulations

Numerical solutions of the Navier-Stokes equations are employed to predict the characteristics of round jets at supersonic speeds. The simulations are performed with a finite volume method which is fourth order accurate in space and second order in time. The overall sound pressure level (OASPL) in the near field of a round free jet is over-predicted by roughly 8 dB relative to an experimental correlation. In an impinging jet, shock motion and vortex stretching are identified as noise generation mechanisms in the impingement zone.

Childs, Robert E.↗
NASA NTRS1990

Navier-Stokes simulation of transonic blade-vortex interactions

Transonic strong blade-vortex interaction is numerically analyzed by solving the unsteady two-dimensional Navier-Stokes equations using an iterative implicit second order scheme. The dominant processes during the interaction are the development of large transverse pressure gradients in the upper leading edge region and the development of disturbances at the root of the lower surface shock wave. As a result of this interaction, high pressure pulses are emitted from the leading edge, and acoustic waves are radiated from the lower surface in a region originally occupied by a supersonic pocket. In addition, severe load variations occur when the vortex is within one chord length of the blade.

Liu, N.-S.↗
NASA NTRS1989

On hypersonic shock layer and its extension beyond the Navier-Stokes level

An extension of the continuum model beyond the Navier-Stokes (NS) level and related issues on problem formulation are examined for a hypersonic shock layer on the basis of Grad's thirteen-moment equations for a Maxwell gas. The 13-moment system, simplified consistently with a fully viscous version of the thin shock-layer approximation, permits correlation with the corresponding NS-based solution. With the exception of pressure and density, several flow properties including normal stress, shear stress and normal heat flux along a streamline are unaffected by translational nonequilibrium and are therefore predicted correctly by the NS solution to the leading order, even in a domain far from translational equilibrium where molecular-transport processes rank equally with the convection.

Cheng, H. K.↗
NASA NTRS1989

Prediction of transitional (laminar-turbulent) hypersonic flows using the parabolized Navier-Stokes equations

Hypersonic transitional flow predictions have been made using the parabolized Navier-Stokes equations with an algebraic transition/turbulence model by appropriately modulating the turbulent viscosity with the available intermittency functions for incompressible and compressible flows. A comparison between the predictions with and without a low Reynolds number correction has also been made. The predictions are compared with the available experimental data and with the theory over a range of Mach number. A simple Re(theta)/M(delta) criterion is shown to satisfactorily predict the meridional variation of the onset location of transition on a cone at a small angle of attack, whereas none of the correlations discussed can do that. Various available correlations are discussed vis-a-vis the predictions as to the locations of the onset and the end of transition.

Kaul, Upender K.↗
NASA NTRS1991

Navier-Stokes simulation of a close-coupled canard-wing-body configuration

The thin-layer Navier-Stokes equations are solved for the flow about a coplanar close-coupled canard-wing-body configuration at a transonic Mach number of 0.90 and at angles of attack ranging from 0 to 12 degrees. The influence of the canard on the wing flowfield, including canard-wing vortex interaction and wing vortex breakdown, is investigated. A study of canard downwash and canard leading-edge vortex effects, which are the primary mechanisms of the canard-wing interaction, is emphasized. Comparisons between the computations and experimental measurements of surface pressure coefficients, lift, drag and pitching moment data are favorable. A grid refinement study for configurations with and without canard shows that accurate results are obtained using a refined grid for angles of attack where vortex burst is present. At an angle of attack of approximately 12 deg, favorable canard-wing interaction which delays wing vortex breakdown is indicated by the computations and is in good agreement with experimental findings.

Tu, Eugene L.↗
NASA NTRS1991

The prediction of viscous hypersonic flows about complex configurations using an upwind parabolized Navier-Stokes code

A three-dimensional parabolized Navier-Stokes computer code which employs an upwind algorithm is used to conduct a numerical study of an advanced maneuvering reentry vehicle configuration. Comparisons between numerical solutions and experimental data are presented for surface pressure, wall heat flux, and overall forces and moments. The effects of angle of attack, angle of yaw, and surface mass injection are investigated. Good agreement is observed between the calculated and measured data. The results of this investigation demonstrate the accuracy and efficiency of an upwind scheme in predicting the hypersonic flow field characteristics about a complex configuration.

Narain, J. P.↗
NASA NTRS1991

Navier-Stokes computation of wing/rotor interaction for a tilt rotor in hover

A method has been developed to analyze the wing/rotor interaction of tilt rotor aircraft in hover. The unsteady, thin-layer compressible Navier-Stokes equations are solved using an implicit, finite difference scheme that employs LU-ADI factorization. The rotor is modeled as an actuator disk which imparts a radial and azimuthal distribution of pressure rise and swirl to the flowfield. The 'chimera' approach of grid point blanking is used to update the rotor boundary conditions. Results are presented for both a rotor alone and for wing/rotor interaction where the thrust coefficient is 0.0164 and wing flap deflection is 67 degrees. Many of the complex flow features are captured including the fountain effect, leading and trailing edge separation, and the unsteady wake beneath the wing. Wing surface pressures compare fairly well with experimental data although the time-averaged download is about twenty percent higher than the measured value. This discrepancy is due to a combination of factors that are discussed.

Fejtek, Ian↗
NASA NTRS1990

Prediction of vortical flows on wings using incompressible Navier-Stokes equations

Grid-density and Reynolds-number effects on integral values and static pressure distributions are presently studied in the steady-state Navier-Stokes solutions obtained for the flow around a round-edged double delta wing using an upwind-relaxation finite-difference algorithm. It is found that while the computed longitudinal aerodynamic coefficients are in good agreement with extant experimental data, the magnitudes of suction-pressure peaks are underpredicted in the vicinity of the trailing edge. It is judged that additional computations employing finer-grid solutions in the vortical flow region are required.

Hsu, C.-H.↗
NASA NTRS1990

An isoparametric finite-element Stokes flow test problem

Two-dimensional Stokes flow of an incompressible Newtonian fluid about a horizontal cylinder situated next to a moving wall is presented and suggested as a nontrivial benchmark problem. Both the primitive variable Galerkin finite-element and the exact solutions are generated on a sequence of refined meshes. Error norms and convergence rates are presented for the C(o) bilinear-velocity/C(-1) piecewise constant-pressure and the C(o) biquadratic-velocity/C(-1) linear-pressure quadrilateral elements.

Maslanik, Miriam K.↗
NASA NTRS1991

Newton's method applied to finite-difference approximations for the steady-state compressible Navier-Stokes equations

Finite-difference approximations for steady-state compressible Navier-Stokes equations, whose two spatial dimensions are written in generalized curvilinear coordinates and strong conservation-law form, are presently solved by means of Newton's method in order to obtain a lifting-airfoil flow field under subsonic and transonnic conditions. In addition to ascertaining the computational requirements of an initial guess ensuring convergence and the degree of computational efficiency obtainable via the approximate Newton method's freezing of the Jacobian matrices, attention is given to the need for auxiliary methods assessing the temporal stability of steady-state solutions. It is demonstrated that nonunique solutions of the finite-difference equations are obtainable by Newton's method in conjunction with a continuation method.

Bailey, Harry E.↗
NASA NTRS1991

Role of artificial viscosity in Euler and Navier-Stokes solvers

A method is proposed to determine directly the amount of artificial viscosity needed for stability using an eigenvalue analysis for a finite difference representation of the Navier-Stokes equations. The stability and growth of small perturbations about a steady flow over airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small time-dependent perturbation about a steady inviscid flow over an NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 deg. The method has been applied to inviscid flows here, but as discussed is also applicable to viscous flows. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundaries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time-marching scheme are also presented. The eigenvalue procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of artificial viscosity.

Mahajan, Aparajit J.↗
NASA NTRS1991

Steady and unsteady solutions of the incompressible Navier-Stokes equations

An algorithm for the solution of the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates is presented. The algorithm can be used to compute both steady-state and time-dependent flow problems. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The equations are solved with a line-relaxation scheme that allows the use of very large pseudotime steps leading to fast convergence for steady-state problems as well as for the subiterations of time-dependent problems. The steady-state solution of flow through a square duct with a 90-deg bend is computed, and the results are compared with experimental data. Good agreement is observed. Computations of unsteady flow over a circular cylinder are presented and compared to other experimental and computational results. Finally, the flow through an artificial heart configuration with moving boundaries is calculated and presented.

Rogers, Stuart E.↗
NASA NTRS1991

A fractional step solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

The time-dependent, three-dimensional incompressible Navier-Stokes equations are presently solved in generalized coordinate systems by means of a fractional-step method whose primitive variable formulation uses as dependent variables, in place of the Cartesian components of the velocity: (1) pressure (defined at the center of the computational cell), and (2) volume fluxes across the faces of the cells. The momentum equations are solved by means of an approximate factorization method. A novel 'ZEBRA' scheme incorporating four-color ordering efficiently solves the Poisson equation. Illustrative two- and three-dimensional laminar flow test cases are computed and evaluated relative to extant numerical and experimental results, and good agreement is obtained.

Rosenfeld, Moshe↗
NASA NTRS1991

Airfoil design method using the Navier-Stokes equations

An airfoil design procedure is described that was incorporated into an existing 2-D Navier-Stokes airfoil analysis method. The resulting design method, an iterative procedure based on a residual-correction algorithm, permits the automated design of airfoil sections with prescribed surface pressure distributions. The inverse design method and the technique used to specify target pressure distributions are described. It presents several example problems to demonstrate application of the design procedure. It shows that this inverse design method develops useful airfoil configurations with a reasonable expenditure of computer resources.

Malone, J. B.↗
NASA NTRS1991

Effect of a multi-dimensional flux function on the monotonicity of Euler and Navier-Stokes computations

A limiting method has been devised for a grid-independent flux function for use with the two-dimensional Euler and Navier-Stokes equations. This limiting is derived from a monotonicity analysis of the model and allows for solutions with reduced oscillatory behavior while still maintaining sharper resolution than a grid-aligned method. In addition to capturing oblique waves sharply, the grid-independent flux function also reduces the entropy generated over an airfoil in an Euler computation and reduces pressure distortions in the separated boundary layer of a viscous-flow airfoil computation. The model has also been extended to three dimensions, although no angle-limiting procedure for improving monotonicity characteristics has been incorporated.

Rumsey, Christopher L.↗
NASA NTRS1991

An implicit three-dimensional Navier-Stokes solver for compressible flow

A three-dimensional numerical method based on the lower-upper symmetric-Gauss-Seidel implicit scheme in conjunction with the flux-limited dissipation model is developed for solving the compressible Navier-Stokes equations. A new computer code which is based on this method requires only 9 microsec per grid-point per iteration on a single processor of a Cray YMP computer and executes at the sustained rate of 170 MFLOPS. A reduction of 4 orders of magnitude in the residual for a high Reynolds number flow using 230 K grid points is obtained in 24 minutes. The computational results compare well with available experimental data.

Yoon, Seokkwan↗