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 1,009 records · Page 56

A second-order closure analysis of turbulent diffusion flames

A complete second-order closure computer program for the investigation of compressible, turbulent, reacting shear layers was developed. The equations for the means and the second order correlations were derived from the time-averaged Navier-Stokes equations and contain third order and higher order correlations, which have to be modeled in terms of the lower-order correlations to close the system of equations. In addition to fluid mechanical turbulence models and parameters used in previous studies of a variety of incompressible and compressible shear flows, a number of additional scalar correlations were modeled for chemically reacting flows, and a typical eddy model developed for the joint probability density function for all the scalars. The program which is capable of handling multi-species, multistep chemical reactions, was used to calculate nonreacting and reacting flows in a hydrogen-air diffusion flame.

Varma, A. K.↗

Numerical simulation of reentry flow around the Shuttle Orbiter including real gas effects

The flight environment of a reentry vehicle is predicted from the numerical solution of fluid dynamics equations for the region between the blunt nose and the swept wings. The inviscid portion of the shock layer is modeled by the Euler equations, but the laminar viscous flow adjacent to the wall is modeled by the approximate parabolic Navier-Stokes equations. The approximations made to the axial gradients of pressure and diffusive fluxes enable the coupled inviscid and viscous equations to be solved efficiently along the body axis. The equilibrium air aftbody code contains significant improvements over its predecessor which only considers a perfect gas model and noncircular configurations. The inclusion of numerical damping either explicitly or implicitly has extended its capabilities for predicting flow field around a winged configuration at higher Machs and greater angles of attack. The results are obtained on the cylindrical coordinates and satisfactory for the Shuttle Orbiter at a free-stream Mach number of 22 and an angle of attack of 40 deg. Also discussed are the inviscid formulation and its application for wind-tunnel conditions.

Li, C. P.↗

Simulation of transonic viscous wing and wing-fuselage flows using zonal methods

The thin-layer Navier-Stokes equations are coupled with a zonal scheme (or domain-decomposition method) to develop the Transonic Navier-Stokes (TNS) wing-alone code. The TNS has a total of 4 zones and is extended to a total of 16 zones for the wing-fuselage version of the code. Results are compared on the Cray X-MP-48 and compared with experimental data.

Flores, Jolen↗

Simulation of transonic viscous wing and wing-fuselage flows using zonal methods

The thin-layer Navier-Stokes equations are coupled with a zonal scheme (or domain-decomposition method) to develop the Transonic Navier-Stokes (TNS) wing-alone code. TNS has a total of four zones and is extended to a total of 16 zones for the wing-fuselage version of the code. Results are computed on the Cray X-MP-48 and compared with experimental data.

Flores, Jolen↗

Aerostructural Vortical Flow Interactions with Applications to F/A-18 and F-117 Tail Buffet

The buffet response of the flexible twin-tail configuration of the F/A-18 and a generic F-111 aircraft are computationally simulated and experimentally validated. The problem is a multidisciplinary one which requires the sequential solution of three sets of equations on a multi-block grid structure. The first set is the unsteady, compressible, full Navier-Stokes equations. The second set is the aeroelastic equations for bending and torsional twin-tail responses. The third set is the grid-displacement equations which are used to update the grid coordinates due to the tail deflections. The computational models consist of a 76 deg. swept back, sharp edged delta wing of aspect ratio of one and a swept-back F/A-18 or F-117 twin-tail. The configuration is pitched at 30 deg. angle-of-attack. The problem is solved for the initial flow conditions with the twin tails kept rigid. Next, the aeroelastic equations of the tails are turned on along with the grid-displacement equations to solve for the bending and torsional tails responses due to the unsteady loads produced by the vortex breakdown flow of the leading-edge vortex cores of the delta wing. Several spanwise locations of the twin tails are investigated. The computational results are validated using several existing experimental data.

Kandil, Osama A.↗

Twin Tail/Delta Wing Configuration Buffet Due to Unsteady Vortex Breakdown Flow

The buffet response of the twin-tail configuration of the F/A-18 aircraft; a multidisciplinary problem, is investigated using three sets of equations on a multi-block grid structure. The first set is the unsteady, compressible, full Navier-Stokes equations. The second set is the coupled aeroelastic equations for bending and torsional twin-tail responses. The third set is the grid-displacement equations which are used to update the grid coordinates due to the tail deflections. The computational model consists of a 76 deg-swept back, sharp edged delta wing of aspect ratio of one and a swept-back F/A-18 twin-tails. The configuration is pitched at 32 deg angle of attack and the freestream Mach number and Reynolds number are 0.2 and 0.75 x 10(exp 6) respectively. The problem is solved for the initial flow conditions with the twin tail kept rigid. Next, the aeroelastic equations of the tails are turned on along with the grid-displacement equations to solve for the uncoupled bending and torsional tails response due to the unsteady loads produced by the vortex breakdown flow of the vortex cores of the delta wing. Two lateral locations of the twin tail are investigated. These locations are called the midspan and inboard locations.

Kandil, Osama A.↗

Computation and Validation of Fluid/Structure Twin Tail Buffet Response

The buffet response of the flexible twin-tail/delta wing configuration, a multidisciplinary problem, is solved using three sets of equations on a multi-block grid structure. The first set is the unsteady, compressible, full Navier-Stokes equations which are used for obtaining the flow-field vector and the aerodynamic loads on the twin tails. The second set is the coupled aeroelastic equations which are used for obtaining the bending and torsional deflections of the twin tails. The third set is the grid-displacement equations which are used for updating the grid coordinates due to the tail deflections. The computational model is similar to the one used by Washburn et. al. which consists of a delta wing of aspect ratio one and twin tails with taper ratios of 0.23. The vortex of the twin tails are located at the wing trailing edge. The configuration is pitched at 30 deg angle of attack, and the freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. With the twin tails fixed as rigid surfaces, 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 produced by the vortex breakdown flow of the leading-edge vortex cores. The configuration is investigated for three spanwise positions of the twin tails; inboard, midspan and outboard locations. The computational results are validated and are in very good agreement with the experimental data of Washburn, et. al.

Massey, Steven J.↗

Effect of Dynamic Rolling Oscillations on Twin Tail Buffet Response

The effect of dynamic rolling oscillations of delta-wing/twin-tail configuration on twin-tail buffet response is investigated. The computational model consists of a sharp-edged delta wing of aspect ratio one and swept-back flexible twin tail with taper ratio of 0.23. The configuration model is statically pitched at 30 deg. angle of attack and then forced to oscillate in roll around the symmetry axis at a constant amplitude of 4 deg. and reduced frequency of pi and 2(pi). The freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. This multidisciplinary problem is solved using three sets of equations on a dynamic multi-block grid structure. The first set is the unsteady, full Navier-Stokes equations, the second set is the aeroelastic equations for coupled bending and torsion vibrations of the tails, and the third set is the grid-displacement equations. The configuration is investigated for inboard position of the twin tails which corresponds to a separation distance between the twin tails of 33% wing span. The computed results are compared with the results of stationary configuration, which previously have been validated using experimental data. The results conclusively showed that the rolling oscillations of the configuration have led to higher loads, higher deflections, and higher excitation peaks than those of the stationary configuration. Moreover, increasing the reduced frequency has led to higher loads and excitation peaks and lower bending and torsion deflections and acceleration.

Sheta, Essam F.↗

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗

Oscillatory flow with heat transfer in a square cavity

A computational study is presented for the flow inside an oscillatory cavity. The numerical scheme employs a semi-implicit, time-splitting method to integrate the two-dimensional full Navier-Stokes equations satisfying continuity to machine accuracy. The efficient use of direct solvers for the uncoupled momentum and pressure equations is demonstrated. The oscillatory cavity flow is studied considering the effects of heat transfer, Reynolds number, and oscillatory Stokes number.

Biringen, S.↗

Oscillatory flow with heat transfer in a square cavity

A computational study is presented for the flow inside an oscillatory cavity. The numerical scheme employs a semiimplicit, time-splitting method to integrate the two-dimensional full Navier-Stokes equations satisfying continuity to machine accuracy. The efficient use of direct solvers for the uncoupled momentum and pressure equations is demonstrated. The oscillatory cavity flow is studied considering the effects of heat transfer, Reynolds number and oscillatory Stokes number.

Danabasoglu, G.↗

Experiments with explicit filtering for LES using a finite-difference method

The equations for large-eddy simulation (LES) are derived formally by applying a spatial filter to the Navier-Stokes equations. The filter width as well as the details of the filter shape are free parameters in LES, and these can be used both to control the effective resolution of the simulation and to establish the relative importance of different portions of the resolved spectrum. An analogous, but less well justified, approach to filtering is more or less universally used in conjunction with LES using finite-difference methods. In this approach, the finite support provided by the computational mesh as well as the wavenumber-dependent truncation errors associated with the finite-difference operators are assumed to define the filter operation. This approach has the advantage that it is also 'automatic' in the sense that no explicit filtering: operations need to be performed. While it is certainly convenient to avoid the explicit filtering operation, there are some practical considerations associated with finite-difference methods that favor the use of an explicit filter. Foremost among these considerations is the issue of truncation error. All finite-difference approximations have an associated truncation error that increases with increasing wavenumber. These errors can be quite severe for the smallest resolved scales, and these errors will interfere with the dynamics of the small eddies if no corrective action is taken. Years of experience at CTR with a second-order finite-difference scheme for high Reynolds number LES has repeatedly indicated that truncation errors must be minimized in order to obtain acceptable simulation results. While the potential advantages of explicit filtering are rather clear, there is a significant cost associated with its implementation. In particular, explicit filtering reduces the effective resolution of the simulation compared with that afforded by the mesh. The resolution requirements for LES are usually set by the need to capture most of the energy-containing eddies, and if explicit filtering is used, the mesh must be enlarged so that these motions are passed by the filter. Given the high cost of explicit filtering, the following interesting question arises. Since the mesh must be expanded in order to perform the explicit filter, might it be better to take advantage of the increased resolution and simply perform an unfiltered simulation on the larger mesh? The cost of the two approaches is roughly the same, but the philosophy is rather different. In the filtered simulation, resolution is sacrificed in order to minimize the various forms of numerical error. In the unfiltered simulation, the errors are left intact, but they are concentrated at very small scales that could be dynamically unimportant from a LES perspective. Very little is known about this tradeoff and the objective of this work is to study this relationship in high Reynolds number channel flow simulations using a second-order finite-difference method.

Lund, T. S.↗

Multistage Schemes with Multigrid for Euler and Navier-Strokes Equations: Components and Analysis

A class of explicit multistage time-stepping schemes with centered spatial differencing and multigrids are considered for the compressible Euler and Navier-Stokes equations. These schemes are the basis for a family of computer programs (flow codes with multigrid (FLOMG) series) currently used to solve a wide range of fluid dynamics problems, including internal and external flows. In this paper, the components of these multistage time-stepping schemes are defined, discussed, and in many cases analyzed to provide additional insight into their behavior. Special emphasis is given to numerical dissipation, stability of Runge-Kutta schemes, and the convergence acceleration techniques of multigrid and implicit residual smoothing. Both the Baldwin and Lomax algebraic equilibrium model and the Johnson and King one-half equation nonequilibrium model are used to establish turbulence closure. Implementation of these models is described.

Swanson, R. C.↗

Numerical study of the effect of an embedded surface-heat source on the separation bubble of supersonic flow

A numerical study of the conjugated problem of a separated supersonic flow field and a conductive solid wall with an embedded heat source is presented. Implicit finite-difference schemes were used to solve the two-dimensional time-dependent compressible Navier-Stokes equations and the time-dependent heat-conduction equation for the solid in both general coordinate systems. A detailed comparison between the thin-layer and Navier-Stokes models was made for steady and unsteady supersonic flow and showed insignificant differences. Steady-state and transient cases were computed and the results show that a temperature pulse at the solid-fluid interface can be used to detect the flow direction near the wall in the vicinity of separation without significant distortion of the flow field.

Degani, D.↗

Large-scale computations in fluid mechanics; Proceedings of the Fifteenth Summer Seminar on Applied Mathematics, University of California, La Jolla, CA, June 27-July 8, 1983. Parts 1 & 2

Papers are presented on such topics as the use of semi-Lagrangian advective schemes in meteorological modeling; computation with high-resolution upwind schemes for hyperbolic equations; dynamics of flame propagation in a turbulent field; a modified finite element method for solving the incompressible Navier-Stokes equations; computational fusion magnetohydrodynamics; and a nonoscillatory shock capturing scheme using flux-limited dissipation. Consideration is also given to the use of spectral techniques in numerical weather prediction; numerical methods for the incorporation of mountains in atmospheric models; techniques for the numerical simulation of large-scale eddies in geophysical fluid dynamics; high-resolution TVD schemes using flux limiters; upwind-difference methods for aerodynamic problems governed by the Euler equations; and an MHD model of the earth's magnetosphere.

Engquist, B. E.↗

Applications of classical and zero-total-pressure-loss sets of Euler equations to delta wings

Classical and zero-total pressure-loss sets of Euler equations were applied to sharp- and round-edge delta wings. The origin of the total pressure was explained in the classical set. For sharp-edged delta wings, all sets of Euler equations produce the same separated flow solutions. For round-edged delta wings and for coarse grids, the solution depends on the level of dissipation, the accuracy of the surface boundary condition, and the type of Euler equations set. For round-edged delta wings and for fine grids, attached flow solutions are obtained. Also presented were three dimensional flow solutions and asymmetric flow solutions including unsteady flow for sharp-edged delta wings. Euler equations should be restricted to sharp-edged wings for real flow solutions. For roung-edged wings, Navier-Stokes equations must be used.

Kandil, Osama A.↗

Boundary-layer equations in generalized curvilinear coordinates

A set of higher-order boundary-layer equations is derived valid for three-dimensional compressible flows. The equations are written in a generalized curvilinear coordinate system, in which the surface coordinates are nonorthogonal; the third axis is restricted to be normal to the surface. Also, higher-order viscous terms which are retained depend on the surface curvature of the body. Thus, the equations are suitable for the calculation of the boundary layer about arbitrary vehicles. As a starting point, the Navier-Stokes equations are derived in a tensorian notation. Then by means of an order-of-magnitude analysis, the boundary-layer equations are developed. To provide an interface between the analytical partial differentiation notation and the compact tensor notation, a brief review of the most essential theorems of the tensor analysis related to the equations of the fluid dynamics is given. Many useful quantities, such as the contravariant and the covariant metrics and the physical velocity components, are written in both notations.

Panaras, Argyris G.↗