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 667 records · Page 37

Turbulence modeling for hypersonic flows

Turbulence modeling for high speed compressible flows is described and discussed. Starting with the compressible Navier-Stokes equations, methods of statistical averaging are described by means of which the Reynolds-averaged Navier-Stokes equations are developed. Unknown averages in these equations are approximated using various closure concepts. Zero-, one-, and two-equation eddy viscosity models, algebraic stress models and Reynolds stress transport models are discussed. Computations of supersonic and hypersonic flows obtained using several of the models are discussed and compared with experimental results. Specific examples include attached boundary layer flows, shock wave boundary layer interactions and compressible shear layers. From these examples, conclusions regarding the status of modeling and recommendations for future studies are discussed.

Marvin, J. G.↗

Turbulence modeling for hypersonic flows

Turbulence modeling for high speed compressible flows is described and discussed. Starting with the compressible Navier-Stokes equations, methods of statistical averaging are described by means of which the Reynolds-averaged Navier-Stokes equations are developed. Unknown averages in these equations are approximated using various closure concepts. Zero-, one-, and two-equation eddy viscosity models, algebraic stress models and Reynolds stress transport models are discussed. Computations of supersonic and hypersonic flows obtained using several of the models are discussed and compared with experimental results. Specific examples include attached boundary layer flows, shock wave boundary layer interactions and compressible shear layers. From these examples, conclusions regarding the status of modeling and recommendations for future studies are discussed.

Marvin, J. G.↗

Turbulence modeling for hypersonic flows

Turbulence modeling for high-speed compressible flows is described and discussed. Starting with the compressible Navier-Stokes equations, methods of statistical averaging are described by means of which the Reynolds-averaged Navier-Stokes equations are developed. Unknown averages in these equations are approximated using various closure concepts. Zero-, one-, and two-equation eddy viscosity models, algebraic stress models, and Reynolds stress transport models are discussed. Computations of supersonic and hypersonic flows obtained using several of the models are discussed and compared with experimental results. Specific examples include attached boundary-layer flows, shock-wave boundary-layer interactions, and compressible shear layers. From these examples, conclusions regarding the status of modeling and recommendations for future studies are discussed.

Marvin, J. G.↗

Computational Investigation and Validation of Twin-Tail Buffet Response Including Dynamics and Control

Multidisciplinary tools for prediction of single rectangular-tail buffet are extended to single swept-back-tail buffet in transonic-speed flow, and multidisciplinary tools for prediction and control of twin-tail buffet are developed and presented. The configuration model consists of a sharp-edged delta wing with single or twin tails that are oriented normal to the wing surface. The tails are treated as cantilevered beams fixed at the root and allowed to oscillate in both bending and torsion. This complex multidisciplinary problem is solved sequentially using three sets of equations on a dynamic single or multi-block grid structure. The first set is the unsteady, compressible, Reynolds-averaged Navier-Stokes equations which are used for obtaining the flow field vector and the aerodynamic loads on the tails. The Navier-Stokes equations are solved accurately in time using the implicit, upwind, flux-difference splitting, finite volume scheme. The second set is the coupled bending and torsion aeroelastic equations of cantilevered beams which are used for obtaining the bending and torsion deflections of the tails. The aeroelastic equations'are solved accurately in time using, a fifth-order-accurate Runge-Kutta scheme. The third set is the grid-displacement equations and the rigid-body dynamics equations, which are used for updating the grid coordinates due to the tail deflections and rigid-body motions. The tail-buffet phenomenon is predicted for highly-swept, single vertical tail placed at the plane of geometric symmetry, and for highly-swept, vertical twin tails placed at three different spanwise separation distances. The investigation demonstrates the effects of structural inertial coupling and uncoupling of the bending and torsion modes of vibration, spanwise positions of the twin-tail, angle of attack, and pitching and rolling dynamic motions of the configuration model on the tail buffet loading and response. The fundamental issue of twin-tail buffet alleviation is addressed using two active flow-control methods. These methods are the tangential leading-edge blowing and the flow suction from the leading-edge vortex cores along their paths. Qualitative and quantitative comparisons with the available experimental data are presented. The comparisons indicate that the present multidisciplinary aeroelastic analysis tools are robust, accurate and efficient.

Kandil, Osama A.↗

Development of an integrated BEM (Boundary Element Mesh) for hot fluid-structure interaction

A boundary integral representation for a coupled approach to fluid flow and solid deformation problems associated with the design of hot-section components such as those in the Space Shuttle Main Engine is discussed. The formulation is based on the fundamental analytical solution of the Navier-Stokes equation for fluid velocity in an infinite domain. This fundamental solution was obtained by decomposing a Navier-Stokes equation into vorticity and dilation transport equations. A boundary integral involving convolutions in time was then constructed in which the convective terms appear in the volume integral.

Dargush, G. F.↗

Navier-Stokes analysis of three-dimensional unsteady flows inside turbine stages

This study presents a numerical method for solving the 3D Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully 3D Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms, while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flowfield of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flowfields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗

Turbulent flow simulation - A large eddy simulator's viewpoint

The state of the art of the data base on turbulent structures and numerical turbulence modelling is assessed. The failure of statistical methods to produce accurate results revealed that turbulence is more than simply random noise. Both randomness and coherence may be characteristics of turbulence, which may also feature noise that is a remnant of previous coherence. Five methods are presently used for numerically predicting turbulent flow behavior: correlations, integral methods, the k-epsilon model, and phenomenological and exact solutions to the Navier-Stokes equations. The higher level models include large eddy simulation and full simulation with numerical solutions to the Navier-Stokes equations; the two methods are applicable to experimental, rather than design work.

Ferziger, J. H.↗

Numerical study of vortex breakdown

The incompressible axisymmetric steady Navier-Stokes equations and the Euler equations are solved numerically to model the breakdown of a vortex. Although the solutions obtained for the Euler equations show a 'vortex breakdown-like' structure, their behavior is very different from that of the Navier-Stokes equations. The details of the numerical algorithmms used are presented, and the results obtained are compared to those in the literature.

Hafez, M.↗

Time-accurate Navier-Stokes calculations with multigrid acceleration

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is implemented by modifying the multigrid-multiblock version of the steady Navier-Stokes equations solver, TLNS3D. The scheme is fully implicit in time and uses TLNS3D to iteratively invert the equations at each physical time step. The design objective of the scheme is unconditional stability (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability which is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. An analysis of the iterative procedure and the implementation of this procedure in TLNS3D are discussed. Numerical results are presented to show both the capabilities of the scheme and its speed up relative to the use of global minimum time stepping. Reductions in computational times of an order of magnitude are demonstrated.

Melson, N. Duane↗

On the breakup of viscous liquid threads

A one-dimensional model evolution equation is used to describe the nonlinear dynamics that can lead to the breakup of a cylindrical thread of Newtonian fluid when capillary forces drive the motion. The model is derived from the Stokes equations by use of rational asymptotic expansions and under a slender jet approximation. The equations are solved numerically and the jet radius is found to vanish after a finite time yielding breakup. The slender jet approximation is valid throughout the evolution leading to pinching. The model admits self-similar pinching solutions which yield symmetric shapes at breakup. These solutions are shown to be the ones selected by the initial boundary value problem, for general initial conditions. Further more, the terminal state of the model equation is shown to be identical to that predicted by a theory which looks for singular pinching solutions directly from the Stokes equations without invoking the slender jet approximation throughout the evolution. It is shown quantitatively, therefore, that the one-dimensional model gives a consistent terminal state with the jet shape being locally symmetric at breakup. The asymptotic expansion scheme is also extended to include unsteady and inerticial forces in the momentum equations to derive an evolution system modelling the breakup of Navier-Stokes jets. The model is employed in extensive simulations to compute breakup times for different initial conditions; satellite drop formation is also supported by the model and the dependence of satellite drop volumes on initial conditions is studied.

Papageorgiou, Demetrios T.↗

Viscous aerodynamic design using the adjoint variable approach

The use of classical optimal control methods, in particular variational methods, to solve the airfoil optimization problem, by deriving a set of adjoint (costate) equations and boundary conditions has already been done for inviscid (potential and Euler flows) and two dimensional, steady state, incompressible flow governed by the Navier-Stokes equations. The interior and boundary terms of the volume integral have been derived (in this work) for the steady Navier-Stokes equations in three dimensions for a viscous, compressible heat conducting fluid. This can be used to derive the adjoint equations and numerical boundary conditions for general classes of problems and hence paves the way for a solution to the aerodynamic optimization problem for compressible viscous flows. The next steps to the realization of that goal are projected as below. The usual square integral pressure functional as an objective function is being replaced by a more realistic drag functional subject to a lift constraint. The feasibility of attempting the more difficult time dependent problem is being investigated. It remains to get the full system of adjoint equations and boundary conditions with the new functional. The state and adjoint equations must be discretized and coded. An appropriate optimization program must be used (steepest descents seems inadequate) and various known airfoil shapes should be recovered in test cases of the computer program.

DeRise, George↗

Numerical investigations in the backflow region of a vacuum plume

The objective of this research is to numerically simulate the vacuum plume flow field in the backflow region of a low thrust nozzle exit. In space applications, the low thrust nozzles are used as a propulsion device to control the vehicle attitude, or to maneuver the vehicle flight trajectory. When the spacecraft is deployed in the orbit or cruising in a planetary mission, the vacuum plume is created behind the nozzle exit (so called backflow region), by the exhausting gas of the propulsion system or by venting internal gas to the extremely low density ambient. The low density vacuum plume flow regions cover the continuum, transitional and free molecular flow regimes, which were characterized by the Knudsen number K(sub n), K(sub n) = lambda(sub m)/L where lambda(sub m) is the mean free path of the gas molecules and L is the characteristic length of the flow field. The transitional regime is defined by 0.01 is less than or equal to K(sub n) is less than or equal to 10. The conventional Navier-Stokes equations are valid only in the flow region close to the nozzle exit since the validity of the Navier-Stokes equations fails asymptotically as the Knudsen number increases. The vacuum plume characteristics prediction is primarily a problem of transitional aerodynamics.

Liaw, Goang-Shin↗

Effective treatments of the singular line boundary problem for three dimensional grids

The axis singularity problem for 3D blunt body flows is examined. Two techniques of eliminating the axis-singularity difficulty are presented for utilization with finite-difference codes: (1) the 3D Navier-Stokes equations are reformulated utilizing a redefined Jacobian that is nonzero at the singular line, and (2) the 3D Navier-Stokes equations are solved using the Roe flux-difference splitting technique with a finite-volume based method for evaluation of the grid Jacobian and the metrics, and appropriate boundary conditions. Real-gas and viscous computations are also presented, demonstrating the generality of the two methods.

Palmer, Grant↗

Boundary conditions for direct simulations of compressible viscous flows

The present consideration of procedures for the definition of boundary conditions for the Navier-Stokes equations emphasizes the derivation of boundary conditions that are compatible with nondissipative algorithms applicable to direct simulations of turbulent flows. A novel formulation for the Euler equations is derived on the basis of characteristic wave relations through boundaries; this formulation is generalized to the Navier-Stokes equations. The method, which applies to both sub- and supersonic flows, is used in reflecting and nonreflecting boundary-condition treatments. Attention is given to practical implementations involving inlet and outlet boundaries and slip and nonslip walls, as well as the test cases of a ducted shear layer, vortices propagating through boundaries, and Poiseuille flow.

Poinsot, T. J.↗

Unsteady aerodynamic flow phenomena in a transonic compressor stage

A three-dimensional unsteady, viscous aerodynamic analysis has been developed for the flow inside a transonic, high-through-flow, single stage compressor. The compressor stage is comprised of a low-aspect-ratio rotor and a closely coupled stator. The analysis is based on a numerical method for solving the three-dimensional Navier-Stokes equation for unsteady viscous flow through multiple turbomachinery blade rows. The method solves the fully three-dimensional Navier-Stokes equation with an implicit scheme. A two-equation turbulence model with a low-Reynolds-number modification is applied for the turbulence closure. A third-order accurate upwinding scheme is used to approximate convection terms while a second-order accurate central difference scheme is used for the discretization of the viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow field inside a transonic, high-through-flow, axial compressor stage. The numerical results are compared with available experimental data.

Hah, C.↗

Three-dimensional flow analysis inside turbomachinery stages with steady and unsteady Navier-Stokes method

This study presents a numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed in the present study. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow field of a subsonic turbine stage and the unsteady flow field inside a transonic, high-through-flow, axial compressor stage. The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid.

Copenhaver, W. W.↗

Calculation of transonic flow over supercritical airfoil sections

Two numerical methods for calculating the transonic flow, including viscous effects over lifting airfoil sections and experimental data are compared for turbulent flow over a supercritical airfoil. In addition, results for a NACA 64A010 airfoil at nonzero angle of attack are compared to demonstrate the applicability of the numerical methods to classical, lifting airfoils. One numerical method is a solution to the time-averaged Navier-Stokes equations throughout the entire flow field. The other is a hybrid method that combines inviscid, boundary-layer, and Navier-Stokes equations in appropriate regions of the flow field. Both methods adequately predict the surface pressures and flow field about the 64A010 airfoil at M = 0.8 and alpha = 2 deg when an appropriate turbulence model is used. The methods are not as successful for the supercritical airfoil.

Rose, W. C.↗

Efficiency trade-offs of steady-state methods using FEM and FDM

The efficiency characteristics of finite element and finite difference approximations for the steady-state solution of the Navier-Stokes equations are examined. The finite element method discussed is a standard Galerkin formulation of the incompressible, steady-state Navier-Stokes equations. The finite difference formulation uses simple centered differences that are O(delta x-squared). Operation counts indicate that a rapidly converging Newton-Raphson-Kantorovitch iteration scheme is generally preferable over a Picard method. A split NOS Picard iterative algorithm for the finite difference method was most efficient.

Gartling, D. K.↗