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.

SourceResearch connectionFrom yearThrough year

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

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS1992

CAN-DO, CFD-based Aerodynamic Nozzle Design and Optimization program for supersonic/hypersonic wind tunnels

A design program is developed which incorporates a modern approach to the design of supersonic/hypersonic wind-tunnel nozzles. The approach is obtained by the coupling of computational fluid dynamics (CFD) with design optimization. The program can be used to design a 2D or axisymmetric, supersonic or hypersonic, wind-tunnel nozzles that can be modeled with a calorically perfect gas. The nozzle design is obtained by solving a nonlinear least-squares optimization problem (LSOP). The LSOP is solved using an iterative procedure which requires intermediate flowfield solutions. The nozzle flowfield is simulated by solving the Navier-Stokes equations for the subsonic and transonic flow regions and the parabolized Navier-Stokes equations for the supersonic flow regions. The advantages of this method are that the design is based on the solution of the viscous equations eliminating the need to make separate corrections to a design contour, and the flexibility of applying the procedure to different types of nozzle design problems.

Korte, John J.↗

NASA NTRS1993

Preconditioned conjugate-gradient methods for low-speed flow calculations

An investigation is conducted into the viability of using a generalized Conjugate Gradient-like method as an iterative solver to obtain steady-state solutions of very low-speed fluid flow problems. Low-speed flow at Mach 0.1 over a backward-facing step is chosen as a representative test problem. The unsteady form of the two dimensional, compressible Navier-Stokes equations is integrated in time using discrete time-steps. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. The new iterative solver is used to solve a linear system of equations at each step of the time-integration. Preconditioning techniques are used with the new solver to enhance the stability and convergence rate of the solver and are found to be critical to the overall success of the solver. A study of various preconditioners reveals that a preconditioner based on the Lower-Upper Successive Symmetric Over-Relaxation iterative scheme is more efficient than a preconditioner based on Incomplete L-U factorizations of the iteration matrix. The performance of the new preconditioned solver is compared with a conventional Line Gauss-Seidel Relaxation (LGSR) solver. Overall speed-up factors of 28 (in terms of global time-steps required to converge to a steady-state solution) and 20 (in terms of total CPU time on one processor of a CRAY-YMP) are found in favor of the new preconditioned solver, when compared with the LGSR solver.

Ajmani, Kumud↗

NASA NTRS1993

Preconditioned Conjugate Gradient methods for low speed flow calculations

An investigation is conducted into the viability of using a generalized Conjugate Gradient-like method as an iterative solver to obtain steady-state solutions of very low-speed fluid flow problems. Low-speed flow at Mach 0.1 over a backward-facing step is chosen as a representative test problem. The unsteady form of the two dimensional, compressible Navier-Stokes equations are integrated in time using discrete time-steps. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. The new iterative solver is used to solve a linear system of equations at each step of the time-integration. Preconditioning techniques are used with the new solver to enhance the stability and the convergence rate of the solver and are found to be critical to the overall success of the solver. A study of various preconditioners reveals that a preconditioner based on the lower-upper (L-U)-successive symmetric over-relaxation iterative scheme is more efficient than a preconditioner based on incomplete L-U factorizations of the iteration matrix. The performance of the new preconditioned solver is compared with a conventional line Gauss-Seidel relaxation (LGSR) solver. Overall speed-up factors of 28 (in terms of global time-steps required to converge to a steady-state solution) and 20 (in terms of total CPU time on one processor of a CRAY-YMP) are found in favor of the new preconditioned solver, when compared with the LGSR solver.

Ajmani, Kumud↗

NASA NTRS1991

Stationary turbulent closure via the Hopf functional equation

Explicit solutions of the stationary Hopf equation are discussed and their computational possibilities are explored. The motivation is to circumvent the infinite hierarchy of coupled equations for the velocity moments and obtain an exact closure of the steady-state 3D Navier-Stokes equations, without modeling assumptions or truncation. The Hopf formulation of the Navier-Stokes equation is reviewed. A stationary homogeneous solution for 2D flow is displayed and discussed. It is shown how depletion of nonlinearity may arise for 3D forced homogeneous flow. The general 3D forced case is considered and a method for closing the 3D unforced equations with arbitrary boundary conditions is derived.

Shen, Hubert H.↗

NASA NTRS1985

A mixed pseudospectral/finite difference method for a thermally driven fluid in a nonuniform gravitational field

A numerical study of the steady, axisymmetric flow in a heated, rotating spherical shell is conducted to model the Atmospheric General Circulation Experiment (AGCE) proposed to run aboard a later shuttle mission. The AGCE will consist of concentric rotating spheres confining a dielectric fluid. By imposing a dielectric field across the fluid a radial body force will be created. The numerical solution technique is based on the incompressible Navier-Stokes equations. In the method a pseudospectral technique is based on the incompressible Navier-Stokes equations. In the method a pseudospectral technique is used in the latitudinal direction, and a second-order accurate finite difference scheme discretizes time and radial derivatives. This paper discusses the development and performance of this numerical scheme for the AGCE which has been modelled in the past only by pure FD formulations. In addition, previous models have not investigated the effect of using a dielectric force to simulate terrestrial gravity. The effect of this dielectric force on the flow field is investigated as well as a parameter study of varying rotation rates and boundary temperatures. Among the effects noted are the production of larger velocities and enhanced reversals of radial temperature gradients for a body force generated by the electric field.

Macaraeg, M. G.↗

NASA NTRS1986

A verification of unsteady Navier-Stokes solutions around oscillating airfoils

A finite difference solution code for the two dimensional Navier-Stokes equations was combined with a moving-grid system. The thin layer Navier-Stokes equations with a turbulence model are solved in a time-accurate manner in order to study the unsteady aerodynamics around airfoils undergoing small amplitude pitching or heaving motions in the transonic regime. The accuracy of the solutions obtained by the use of the present moving-grid technqiue is investigated. The effects of the minimum grid size and the integrating time-step size on the solutions are also checked. Some of the solutions obtained by the present method are compared with experimental results. It is demonstrated that the unsteady aerodynamics around oscillating airfoils can be predicted fairly well by the present code for cases in which the dynamic angle of attack or displacement is small.

Nakamichi, J.↗

NASA NTRS1990

Combination Of Techniques For Computing Incompressible Flow

Pseudocompressibility, upwind differencing, and other techniques used to solve Navier-Stokes equations. Scheme for finite-difference numerical solution of two-dimensional Navier-Stokes equations of incompressible flow combines several recently developed methods, each developed to increase speed and/or accuracy of computations of this kind.

Kwak, D.↗

NASA NTRS1994

Development of advanced Navier-Stokes solver

The objective of research was to develop and validate new computational algorithms for solving the steady and unsteady Euler and Navier-Stokes equations. The end-products are new three-dimensional Euler and Navier-Stokes codes that are faster, more reliable, more accurate, and easier to use. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible/incompressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. Convergence rates and the robustness of the codes are enhanced by the use of an implicit full approximation storage multigrid method.

Yoon, Seokkwan↗

NASA NTRS1991

A Fourier-Chebyshev pseudospectral method for solving steady 3-D Navier-Stokes and heat equations in cylindrical cavities

A Fourier-Chebyshev pseudospectral method for solving steady 3D Navier-Stokes equations in cylindrical cavities is presented and discussed. The general method is pseudo-unsteady and uses a semi-implicit finite difference scheme for the time discretization. The generalized ADI (Alternating Direction Implicit) procedure is then applied to reduce the problem to a successive solution of one-dimensional problems. The spatial approximation uses a Fourier-Galerkin approximation in the periodic direction and a Chebyshev-collocation approximation in the other directions. Difficulties related to the pressure are surmounted by using the artifical compressibility method. A suitable variable change was chosen to avoid the problem of singularity at the axis generated by cylindrical coordinates. The method is first tested on an advection-diffusion equation and then on the Navier-Stokes equations. Finally, the method is illustrated by a convection problem of a differentially heated fluid.

Pulicani, J. P.↗

NASA NTRS1985

Simulation of Multistage Turbine Flows

The numerical simulation of turbine flows serves to enhance the understanding of the flow phenomena within multistage turbomachinery components. The direct benefit of this activity is improved modeling capability, which can be used to improve component efficiency and durability. A hierarchy of equations was formulated to assess the difficulty in analyzing the flow field within multistage turbomachinery components. The Navier-Stokes equations provides the most complete description. The simplest description is given by a set of equations that govern the quasi-one-dimensional flow. The number of unknowns to be solved for increases monotonically above the number of equations. The development of the additional set of equations needed to mathematically close the system of equations forms the closure problem associated with that level of description. For the Navier-Stokes equation there is no closure problem. For the quasi-one-dimensional equation set random flow fluctuations, unsteady fluctuations, nonaxisymmetric flow variations, and hub-to-shroud variations on the quasi-one-dimensional flow must be accounted for.

Celestina, M. L.↗

NASA NTRS1991

Computation of compressible quasi-axisymmetric slender vortex flow and breakdown

The unsteady, compressible Navier-Stokes equations are used to compute and analyze compressible quasi-axisymmetric isolated vortices. The Navier-Stokes equations are solved using an implicit, upwind, flux difference splitting finite volume scheme. The developed three dimensional solver was verified by comparing its solution profiles with those of a slender, quasi-axisymmetric vortex solver for a subsonic, quasi-axisymmetric vortex in an unbounded domain. The Navier-Stokes solver is then used to solve for a supersonic, quasi-axisymmetric vortex flow in a configured circular duct. Steady and unsteady vortex-shock interactions and breakdown were captured. The problem was also calculated using the Euler solver of the same code; the results were compared with those of the Navier-Stokes solver. The effect of the initial swirl was investigated.

Kandil, Osama A.↗

NASA NTRS1991

Computation of steady and unsteady compressible quasi-axisymmetric vortex flow and breakdown

The unsteady, compressible Navier-Stokes equations are used to compute and analyze compressible quasi-axisymmetric isolated vortices. The Navier-Stokes equations are solved using an implicit, upwind, flux-difference splitting finite-volume scheme. The developed three-dimensional solver has been verified by comparing its solution profiles with those of a slender, quasi-axisymmetric vortex solver for a subsonic, isolated quasi-axisymmetric vortex in an unbounded domain. The Navier-Stokes solver is then used to solve for a supersonic quasi-axisymmetric vortex flow in a configured circular duct. Steady and unsteady vortex-shock interactions and breakdown have been captured. The problem has also been calculated using the Euler solver of the same code and the results are compared with those of the Navier-Stokes solver. The effect of the initial swirl has been tentatively studied.

Kandil, Osama A.↗

NASA NTRS1990

An implicit finite element method for high speed flows

A fast algorithm is presented for constructing continuous lines, consisting of element sides, on general unstructured two-dimensional triangular meshes. The lines must pass through each node of the mesh once and only once. The discussion focuses on the use of these lines in a relaxation method for the solution of the equation system arising from an implicit algorithm for the solution of two-dimensional Euler and Navier-Stokes equations on general unstructured grids. It is also shown that the method can be used for solving three-dimensional Navier-Stokes equations on a grid composed of both structured and unstructured regions.

Hassan, O.↗

NASA NTRS1992

Bleed-boundary conditions for numerically simulated mixed-compression supersonic inlet flow

A numerical study was performed to investigate boundary conditions for bleed openings that are needed in the numerical simulation of critical flow through turbojet inlets with bleed. Nine different boundary conditions (BCs) were devised and examined. Of the BCs examined, three were able to stabilize the terminal shock near the inlet's throat. This numerical study was based on the ensemble-averaged, unsteady, three-dimensional, compressible Navier-Stokes equations closed by the Baldwin-Lomax algebraic turbulence model. Solutions to the Navier-Stokes equations were obtained by using a hybrid implicit-explicit method (Beam and Warming/MacCormack). The multiblock grid system used in the solution procedure was generated by solving a set of elliptic partial differential equations. Solutions obtained were compared with available experimental data.

Chyu, W. J.↗

NASA NTRS1994

Unsteady aerodynamic analyses for turbomachinery aeroelastic predictions

Applications for unsteady aerodynamics analysis in this report are: (1) aeroelastic: blade flutter and forced vibration; (2) aeroacoustic: noise generation; (3) vibration and noise control; and (4) effects of unsteadiness on performance. This requires that the numerical simulations and analytical modeling be accurate and efficient and contain realistic operating conditions and arbitrary modes of unsteady excitation. The assumptions of this application contend that: (1) turbulence and transition can be modeled with the Reynolds averaged and using Navier-Stokes equations; (2) 'attached' flow with high Reynolds number will require thin-layer Navier-Stokes equations, or inviscid/viscid interaction analyses; (3) small-amplitude unsteady excitations will need nonlinear steady and linearized unsteady analyses; and (4) Re to infinity will concern inviscid flow. Several computer programs (LINFLO, CLT, UNSVIS, AND SFLOW-IVI) are utilized for these analyses. Results and computerized grid examples are shown. This report was given during NASA LeRC Workshop on Forced Response in Turbomachinery in August of 1993.

Verdon, Joseph M.↗

NASA NTRS1988

Comparison of shock structure solutions using independent continuum and kinetic theory approaches

A vehicle traversing the atmosphere will experience flight regimes at high altitudes in which the thickness of a hypersonic shock wave is not small compared to the shock standoff distance from the hard body. When this occurs, it is essential to compute accurate flow field solutions within the shock structure. In this paper, one-dimensional shock structure is investigated for various monatomic gases from Mach 1.4 to Mach 35. Kinetic theory solutions are computed using the Direct Simulation Monte Carlo method. Steady-state solutions of the Navier-Stokes equations and of a slightly truncated form of the Burnett equations are determined by relaxation to a steady state of the time-dependent continuum equations. Monte Carlo results are in excellent agreement with published experimental data and are used as bases of comparison for continuum solutions. For a Maxwellian gas, the truncated Burnett equations are shown to produce far more accurate solutions of shock structure than the Navier-Stokes equations.

Fiscko, Kurt A.↗

NASA NTRS1978

The computation of incompressible, three-dimensional flows in containers with extensible walls

The theoretical study of flow in biological vascular systems is made very difficult in connection with local and temporal changes of the cross section. Experimental investigations with the aid of model tests are, however, not enough for a solution of the problems, and numerical solutions are more and more employed for a description of such flows. A description is presented of a difference procedure for the solution of the complete Navier-Stokes equations with curvilinear coordinates for three-dimensional flows in containers. The integration of the Navier-Stokes equations for flows in containers with rigid and moving walls is considered.

Liu, N. S.↗

NASA NTRS1988

Time-accurate simulations of a shear layer forced at a single frequency

Calculations are presented for the forced shear layer studied experimentally by Oster and Wygnanski, and Weisbrot. Two different computational approaches are examined: Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES). The DNS approach solves the full three dimensional Navier-Stokes equations for a temporally evolving mixing layer, while the LES approach solves the two dimensional Navier-Stokes equations with a subgrid scale turbulence model. While the comparison between these calculations and experimental data was hampered by a lack of information on the inflow boundary conditions, the calculations are shown to qualitatively agree with several aspects of the experiment. The sensitivity of these calculations to factors such as mesh refinement and Reynolds number is illustrated.

Claus, R. W.↗