Engineering Papers⌕ Search

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.

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

At least 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Supersonic flow over three-dimensional ablated nosetips using an unsteady implicit numerical procedure

The three-dimensional supersonic flow over passive, that is, nonablating, indented nosetips of reentry vehicles is determined using an unsteady implicit numerical algorithm which solves either the inviscid Euler equations or the 'thin-layer' Navier-Stokes equations. A nonorthogonal independent variable transformation is used to map the distorted physical domain, containing multiple zones of embedded subsonic flow and separated flow regions into a rectangular computational volume at whose boundaries the required permeable or impermeable boundary conditions are simulated. Use of the implicit algorithm results in faster convergence to the steady state because of a larger allowable time step over conventional explicit schemes. The numerical results obtained compare favorably with existing numerical solutions and experimental data for simple spheres which validates the program. Results are also presented for analytically defined indented bodies for both laminar and turbulent flow conditions that demonstrate the program's capability for computing such flows.

Kutler, P.↗
NASA NTRS1986

Grid adaption for bluff bodies

Methods of grid adaptation are reviewed and a method is developed with the capability of adaptation to several flow variables. This method is based on a variational approach and is an algebraic method which does not require the solution of partial differential equations. Also the method was formulated in such a way that there is no need for any matrix inversion. The method is used in conjunction with the calculation of hypersonic flow over a blunt nose. The equations of motion are the compressible Navier-Stokes equations where all viscous terms are retained. They are solved by the MacCormack time-splitting method and a movie was produced which shows simulataneously the transient behavior of the solution and the grid adaptation. The results are compared with the experimental and other numerical results.

Abolhassani, Jamshid S.↗
NASA NTRS1987

Grid adaption for hypersonic flow

The methods of grid adaption are reviewed and a method is developed with the capability of adaption to several flow variables. This method is based on a variational approach and is an algebraic method which does not require the solution of partial differential equations. Also the method has been formulated in such a way that there is no need for any matrix inversion. The method is used in conjunction with the calculation of hypersonic flow over a blunt nose body. The equations of motion are the compressible Navier-Stokes equations where all viscous terms are retained. They are solved by the MacCormack time-splitting method. A movie has been produced which shows simultaneously the transient behavior of the solution and the grid adaption.

Abolhassani, Jamshid S.↗
NASA NTRS1987

A spectral multi-domain technique for analysis of chemically reacting flows

A technique is presented which models the chemical kinetics and flow kinematics of a nonionized air mixture passing through a fully resolved shock wave, thus alleviating the need for artificial viscosity. The quasi-one-dimensional Navier-Stokes equations and the species conversation equation are the governing equations. The validity of the multidomain Navier-Stokes algorithm is demonstrated via comparison with experiment.

Macaraeg, Michele G.↗
NASA NTRS1988

Application of unsteady aerodynamic methods for transonic aeroelastic analysis

Aerodynamic methods for aeroelastic analysis are applied to various flow problems. These methods include those that solve the three-dimensional transonic small disturbance (TSD) potential equation, the two-dimensional full potential (FP) equation, and the thin-layer Navier-Stokes equations. Flutter analysis performed using TSD aerodynamics shows that such methods can be used to analyze some aeroelastic phenomena. For thicker bodies and larger amplitude motions a nonisentropic full potential method is presented.

Whitlow, Woodrow, Jr.↗
NASA NTRS1988

Instability and sound emission from a flow over a curved surface

The growth and decay of a wavepacket convecting in a boundary layer over a concave-convex surface is studied numerically using direct computations of the Navier-Stokes equations. The resulting sound radiation is computed using the linearized Euler equations with the pressure from the Navier-Stokes solution as a time-dependent boundary condition. It is shown that on the concave portion the amplitude of the wavepacket increases and its bandwidth broadens while on the convex portion some of the components in the packet are stabilized. The pressure field decays exponentially away from the surface and then algebraically exhibits a decay characteristic of acoustic waves in two dimensions. The far-field acoustic pressure exhibits a peak at a frequency corresponding to the inflow instability frequency.

Maestrello, L.↗
NASA NTRS1989

Control of sound radiation from a wavepacket over a curved surface

Active control of acoustic pressure in the far field resulting from the growth and decay of a wavepacket convecting in a boundary layer over a concave-convex surface is investigated numerically using direct computations of the Navier-Stokes equations. The resulting sound radiation is computed using linearized Euler equations with the pressure from the Navier-Stokes solution as a time-dependent boundary condition. The acoustic far field exhibits directivity type of behavior that points upstream to the flow direction. A fixed control algorithm is used where the attenuation signal is synthesized by a filter which actively adapt it to the amplitude-time response of the outgoing acoustic wave.

Maestrello, Lucio↗
NASA NTRS1989

Numerical modeling of a vortex stabilized arcjet thruster

A numerical method to solve the equations governing a vortex stabilized arcjet thruster is being developed. The model will allow the effects of swirling flow and geometry on arcjet thruster performance to be determined. The propellant flow equations, i.e., the axisymmetric, thin layer, Navier-Stokes equations, are solved using a Gauss-Siedel line-relaxation procedure. An implicit FTCS method is used to solve the electromagnetic field equations. A grid generation scheme was developed for an arbitrary arcjet geometry. The model will allow the radial and axial components of velocity and current distributions to be determined from a region upstream of the cathode, through the constrictor, to the exit plane of the nozzle.

Keith, Theo G., Jr.↗
NASA NTRS1987

Application of sound and temperature to control boundary-layer transition

The growth and decay of a wave packet convecting in a boundary layer over a concave-convex surface and its active control by localized surface heating are studied numerically using direct computations of the Navier-Stokes equations. The resulting sound radiations are computed using linearized Euler equations with the pressure from the Navier-Stokes solution as a time-dependent boundary condition. It is shown that on the concave portion the amplitude of the wave packet increases and its bandwidth broadens while on the convex portion some of the components in the packet are stabilized. The pressure field decays exponentially away from the surface and then algebraically, exhibiting a decay characteristic of acoustic waves in two dimensions. The far-field acoustic behavior exhibits a super-directivity type of behavior with a beaming downstream. Active control by surface heating is shown to reduce the growth of the wave packet but have little effect on acoustic far field behavior for the cases considered. Active control by sound emanating from the surface of an airfoil in the vicinity of the leading edge is experimentally investigated. The purpose is to control the separated region at high angles of attack. The results show that injection of sound at shedding frequency of the flow is effective in an increase of lift and reduction of drag.

Maestrello, Lucio↗
NASA NTRS1991

Distributed Minimal Residual (DMR) method for acceleration of iterative algorithms

A new method for enhancing the convergence rate of iterative algorithms for the numerical integration of systems of partial differential equations was developed. It is termed the Distributed Minimal Residual (DMR) method and it is based on general Krylov subspace methods. The DMR method differs from the Krylov subspace methods by the fact that the iterative acceleration factors are different from equation to equation in the system. At the same time, the DMR method can be viewed as an incomplete Newton iteration method. The DMR method was applied to Euler equations of gas dynamics and incompressible Navier-Stokes equations. All numerical test cases were obtained using either explicit four stage Runge-Kutta or Euler implicit time integration. The formulation for the DMR method is general in nature and can be applied to explicit and implicit iterative algorithms for arbitrary systems of partial differential equations.

Lee, Seungsoo↗
NASA NTRS1993

Analytical description of the breakup of liquid jets in air

A viscous or inviscid cylindrical jet with surface tension in a vacuum tends to pinch due to the mechanism of capillary instability. Similarity solutions are constructed which describe this phenomenon as a critical time is encountered, for two physically distinct cases: inviscid jets governed by the Euler equations and highly viscous jets governed by the Stokes equations. In both cases the only assumption imposed is that at the time of pinching the jet shape has a radial length scale which is smaller than the axial length scale. For the inviscid case, we show that our solution corresponds exactly to one member of the one-parameter family of solutions obtained from slender jet theories and the shape of the jet is locally concave at breakup. For highly viscous jets our theory predicts local shapes which are monotonic increasing or decreasing indicating the formation of a mother drop connected to the jet by a thin fluid tube. This qualitative behavior is in complete agreement with both direct numerical simulations and experimental observations.

Papageorgiou, Demetrios T.↗
NASA NTRS1995

A Zonal Approach for Prediction of Jet Noise

A zonal approach for direct computation of sound generation and propagation from a supersonic jet is investigated. The present work splits the computational domain into a nonlinear, acoustic-source regime and a linear acoustic wave propagation regime. In the nonlinear regime, the unsteady flow is governed by the large-scale equations, which are the filtered compressible Navier-Stokes equations. In the linear acoustic regime, the sound wave propagation is described by the linearized Euler equations. Computational results are presented for a supersonic jet at M = 2. 1. It is demonstrated that no spurious modes are generated in the matching region and the computational expense is reduced substantially as opposed to fully large-scale simulation.

Shih, S. H.↗
NASA NTRS1961

Analysis of Multipoint-Multitime Correlations and Diffusion in Decaying Homogeneous Turbulence

Two-point, two-time correlation equations are obtained by considering the Navier-Stokes equations for two points in a fluid at two time. By neglecting the triple correlations in the equations, a solution is obtained for the final period of decay. The analysis is extended to earlier times by considering three points at three different times. The set of equations is made determinate by neglecting the quadruple correlations in comparison with the triple correlations. The diffusion of particles from a source in a decaying turbulent field is calculated approximately by assuming that the velocity fluctuations are small.

Deissler, Robert G.↗
NASA NTRS1981

Calculation of Unsteady Transonic Flow over an Airfoil

An implicit finite-difference solver for either the Euler equations or the "thin-layer" Navier-Stokes equations was used to calculate a transonic flow over the NACA 64A010 airfoil pitching about its one-quarter chord. An unsteady automatic grid-generation procedure that will improve significantly the computational efficiency of various unsteady flow problems is described. The calculated results for both inviscid and viscous flows at Much number 0.8 over the airfoil oscillating with reduced frequency referenced to one-half chord, 0.2, are compared with experimental data measured in the Ames 11 x 11-ft Transonic Wind Tunnel. Nonlinear, unsteady effects of the flow on the surface pressure variations, shock-wave excursions, and overall airloads are examined. Good agreements between the results of computations and experiments were obtained. In the shock-wave region, however, the results of the viscous-flow computations showed closer agreement with the experimental data.

Chyu, W. J.↗
NASA NTRS1981

Calculation of Unsteady Transonic Flow Over an Airfoil

An implicit finite-difference solver for either the Euler equations or the "thin-layer" Navier-Stokes equations was used to calculate a transonic flow over the NACA 64A010 airfoil pitching about its one-quarter chord. An unsteady automatic grid-generation procedure that will improve significantly the computational efficiency of various unsteady flow problems is described. The calculated results for both inviscid and viscous flows at Mach number 0.8 over the airfoil oscillating with reduced frequency referenced to one-half chord, 0.2, are compared with experimental data measured in the Ames 11 x 11 ft Transonic Wind Tunnel. Nonlinear, unsteady effects of the flow on the surface pressure variations, shock-wave excursions, and overall airloads are examined. Good agreements between the results of computations and experiments were obtained. In the shock-wave region, however, the results of the viscous-flow computations showed closer agreement with the experimental data.

Chyu, W. J.↗
NASA NTRS1998

Effect of Configuration Pitching Motion on Twin Tail Buffet Response

The effect of dynamic pitch-up motion of delta wing on twin-tail buffet response is investigated. The computational model consists of a delta wing-twin tail configuration. The computations are carried out on a dynamic multi-block grid structure. This multidisciplinary problem is solved using three sets of equations which consists of the unsteady Navier-Stokes equations, the aeroelastic equations, and the grid displacement equations. The configuration is pitched-up from zero up to 60 deg. angle of attack, and the freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. With the twin tail fixed as rigid surfaces and with no-forced pitch-up motion, the problem is solved for the initial flow conditions. Next, the problem is solved for the twin-tail response for uncoupled bending and torsional vibrations due to the unsteady loads on the twin tail and due to the forced pitch-up motion. The dynamic pitch-up problem is also solved for the flow response with the twin tail kept rigid. The configuration is investigated for inboard position of the twin tail which corresponds to a separation distance between the twin tail of 33% wing chord. The computed results are compared with the available experimental data.

Sheta, Essam F.↗
NASA NTRS1996

Aerodynamic Influence Coefficient Computations Using Euler/Navier-Stokes Equations on Parallel Computers

Modem design requirements for an aircraft push current technologies used in the design process to their limit or sometimes require more advanced technologies to meet the requirement. New design requirements always demand to improve the operational performance. Accurate prediction of aerodynamic coefficients is essential to improve the performance. For example, in the design of an advanced subsonic civil transport, since the fluid flow at transonic regime shows strong nonlinearities, high fidelity equations, such as the Euler or Navier-Stokes equations predict flow characteristics more accurately than the linear aerodynamics, which are widely used in the current design process However, high fidelity flow equations are computationally expensive and require an order of magnitude longer time to obtain aerodynamic coefficients required in the design. Parallel computing is one possibility to cut down the computational turn-around time in using high fidelity equations so that high fidelity equations would be incorporated into the design process. By doing so, high fidelity equations would be used in the routine design process. This work will demonstrate the feasibility of using high fidelity flow equations in a design process by computing aerodynamic influence coefficients of a wing-body-empennage configuration on a multiple-instruction, multiple-data parallel computer.

Byun, Chansup↗
NASA NTRS2004

A Variable Turbulent Schmidt Number Formulation for Scramjet Application

In high speed engines, thorough turbulent mixing of fuel and air is required to obtain high performance and high efficiency. Thus, the ability to predict turbulent mixing is crucial in obtaining accurate numerical simulation of an engine and its performance. Current state of the art in CFD simulation is to assume both turbulent Prandtl number and Schmidt numbers to be constants. However, since the mixing of fuel and air is inversely proportional to the Schmidt number, a value of 0.45 for the Schmidt number will produce twice as much diffusion as that with a value of 0.9. Because of this, current CFD tools and models have not been able to provide the needed guidance required for the efficient design of a scramjet engine. The goal of this investigation is to develop the framework needed to calculate turbulent Prandtl and Schmidt numbers as part of the solution. This requires four additional equations: two for the temperature variance and its dissipation rate and two for the concentration variance and its dissipation rate. In the current investigation emphasis will be placed on studying mixing without reactions. For such flows, variable Prandtl number does not play a major role in determining the flow. This, however, will have to be addressed when combustion is present. The approach to be used is similar to that used to develop the k-zeta model. In this approach, relevant equations are derived from the exact Navier-Stokes equations and each individual correlation is modeled. This ensures that relevant physics is incorporated into the model equations. This task has been accomplished. The final set of equations have no wall or damping functions. Moreover, they are tensorially consistent and Galilean invariant. The derivation of the model equations is rather lengthy and thus will not be incorporated into this abstract, but will be included in the final paper. As a preliminary to formulating the proposed model, the original k-zeta model with constant turbulent Prandtl and Schmidt numbers is used to model the supersonic coaxial jet mixing experiments involving He, O2 and air.

Xiao, X.↗