SEARCH · Search NASA

Results for “Equations of fluid dynamics”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Stellar convection theory. I - The anelastic modal equations

Methods are developed for dealing with the various dynamical problems that arise because of convective zones in stars. A system of equations for stellar convection is derived from the full equations of compressible fluid dynamics with the aid of two major approximations. The first of these is the anelastic approximation, which involves both the filtering out of acoustic waves and a suitable linearization of the fluctuating thermodynamic variables. The second one approximates the horizontal structure of convection by expanding the motion in a set of horizontal cellular platforms and severely truncating the expansion. The resulting system of partial differential equations, referred to as the anelastic modal equations, is outlined along with suggested boundary conditions and techniques for solving the equations. Ways of assessing the overall validity of the present treatment are discussed.

Latour, J.↗

NASA NTRS1990

A three-dimensional upwind parabolized Navier-Stokes code for chemically reacting flows

A new upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the three-dimensional flow of chemically reacting air around hypersonic vehicles. The code is a modification of the perfect gas, three-dimensional UPS code of Lawrence et al. (1986) which has been extended in the present study to permit the calculation of hypersonic, viscous flows in chemical nonequilibrium. The algorithm solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that has been modified to account for real gas effects. The present code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. The fluid medium is assumed to be a chemically reacting mixture of thermally perfect (but calorically imperfect) gases in thermal equilibrium. Results are presented for the hypersonic laminar flow over a cone at 0- and 10-deg angles of attack and for a generic hypersonic vehicle. Calculations are performed assuming either perfect gas, equilibrium air, or finite-rate chemistry.

Buelow, Philip E.↗

NASA NTRS1992

A new method for transonic static aeroelasticity problems

A new method has been developed to calculate the steady flow and structural deformations for fluid/structure interaction problems. The discretized fluid dynamic and structural equations are regarded as a single set of coupled, nonlinear, algebraic equations. The equilibrium solution is directly obtained using Newton's method. The governing equations used for the fluid flow are the two-dimensional Navier-Stokes equations, and a finite-element model is used to represent the structure. This paper describes the analytical method and presents sample calculations demonstrating the technique. The results show rapid convergence and good agreement with experimental data.

Felker, Fort F.↗

NASA NTRS1985

GMRES acceleration of computational fluid dynamics codes

The generalized minimal residual algorithm (GMRES) is a conjugate-gradient like method that applies directly to nonsymmetric linear systems of equations. In this paper, GMRES is modified to handle nonlinear equations characteristic of computational fluid dynamics. Attention is devoted to the concept of preconditioning and the role it plays in assuring rapid convergence. A formulation is developed that allows GMRES to be preconditioned by the solution procedures already built into existing computer codes. Examples are provided that demonstrate the ability of GMRES to greatly improve the robustness and rate of convergence of current state-of-the-art fluid dynamics codes. Theoretical aspects of GMRES are presented that explain why it works. Finally, the advantage GMRES enjoys over related methods such as conjugate gradients are discussed.

Wigton, L. B.↗

NASA NTRS1975

Kinetic description of ionospheric dynamics in the three-fluid approximation

Conservation equations are developed in the three-fluid approximation for general application problems of ionospheric dynamics in the altitude region 90 km to 800 km for all geographic locations. These equations are applied to a detailed study of auroral E region neutral winds and their relationship to ionospheric plasma motions.

Comfort, R. H.↗

NASA NTRS1977

A block iterative LU solver for weakly coupled linear systems

A hybrid technique, called the block iterative LU solver, is proposed for solving the linear equations resulting from a finite element numerical analysis of certain fluid dynamics problems where the equations are weakly coupled between distinct sets of variables. Either the block Jacobi iterative method or the block Gauss-Seidel iterative solver is combined with LU decomposition.

Cooke, C. H.↗

NASA NTRS2015

A Versatile Nonlinear Method for Predictive Modeling

As computational fluid dynamics techniques and tools become widely accepted for realworld practice today, it is intriguing to ask: what areas can it be utilized to its potential in the future. Some promising areas include design optimization and exploration of fluid dynamics phenomena (the concept of numerical wind tunnel), in which both have the common feature where some parameters are varied repeatedly and the computation can be costly. We are especially interested in the need for an accurate and efficient approach for handling these applications: (1) capturing complex nonlinear dynamics inherent in a system under consideration and (2) versatility (robustness) to encompass a range of parametric variations. In our previous paper, we proposed to use first-order Taylor expansion collected at numerous sampling points along a trajectory and assembled together via nonlinear weighting functions. The validity and performance of this approach was demonstrated for a number of problems with a vastly different input functions. In this study, we are especially interested in enhancing the method's accuracy; we extend it to include the second-orer Taylor expansion, which however requires a complicated evaluation of Hessian matrices for a system of equations, like in fluid dynamics. We propose a method to avoid these Hessian matrices, while maintaining the accuracy. Results based on the method are presented to confirm its validity.

nonlinear modeling↗

NASA NTRS1987

Adaptive finite difference methods in fluid dynamics

An adaptive method to solve partial differential equations in fluid mechanics is presented. The approach requires internal boundary conditions that must be conservative, data structures for keeping track of several layers of fine grid patches, error estimation, and heuristics for automatic grid generation. In practical calculations gains in computer efficiency up to 10 over nonadaptive methods are observed. The whole procedure takes 3000 lines of FORTRAN code.

Berger, Marsha J.↗

NASA NTRS1996

[?]Nonlinear Issues in the Aerothermochemistry of Gases and Materials and the Associated Physics and Dynamics of Interfaces

Our research and technology are focused on nonlinear issues in the aerothermochemistry of gases and materials and the associated physics and dynamics of interfaces. Our program is now organized to aggressively support the NASA Aeronautics Enterprise so as to: (a) develop a new generation of environmentally compatible, economic subsonic aircraft; (b) develop the technology base for an economically viable and environmentally compatible high-speed civil transport; (c) develop the technology options for new capabilities in high-performance aircraft; (d) develop hypersonic technologies for air-breathing flight; and (e) develop advanced concepts, understanding of physical phenomena, and theoretical, experimental, and computational tools for advanced aerospace systems. The implications from our research for aeronautical and aerospace technology have been both broad and deep. For example, using advanced computational techniques, we have determined exact solutions for the Schrodinger equation in electron-molecule scattering allowing us to evaluate atmospheric models important to reentry physics. We have also found a new class of exact solutions for the Navier Stokes equations. In experimental fluid dynamics, we have found explicit evidence of turbulence modification of droplet sizes in shock tube flow with condensation. We have developed a new diagnostic tool for the direct estimation of flow velocities at MHz sampling rates in quasi-one dimensional turbulent flow. This procedure suggests an unexpected confirmation of the possibility of 'natural' closure in Reynolds stresses with deep implications for the development of turbulent models. A transient increase is observed in both the spectral energy decay rate and the degree of chaotic complexity at the interface of a shock wave and a turbulent ionized gas. Even though the gas is apparently brought to rest by the shock wave, no evidence is found either of the expected relaminarization. A unique diamond-shaped nozzle has been designed for a detailed investigation of the effect of significant streamwise vorticity on the acoustic and IR characteristics of supersonic jets. Our results provide convincing evidence of the significant effect of vorticity on the far-field noise for the diamond jet as compared to the conventional round jets. We have found that the countercurrent shear layer mixes much more efficiently than conventional coflowing shear layers. We also developed the fluid thrust vectoring procedures which use counter flow to vector a jet. Our materials research has shown that the steep stress gradients at the fiber-matrix interface could be the primary cause of interface cracks after the processing of metallic and intermetallic matrix composites. New techniques have been evolved for: the microcharacterization of materials including microplastic strain and, point by point, the misorientation and plasticity for matrix composites; thermally induced stress measurements and load relaxation; the growth and characterization of metallic matrix composite interfaces; and for the growth of ferrite materials by pulsed laser deposition. The FAMU commitment to the HBCU Research Center also continues to be broad and deep.

Johnson, Joseph A., III↗

NASA NTRS1996

Unsteady, Transonic Flow Around Delta Wings Undergoing Coupled and Natural Modes Response: A Multidisciplinary Problem

The unsteady, three-dimensional Navier-Stokes equations coupled with the Euler equations of rigid-body dynamics are sequentially solved to simulate and analyze the aerodynamic response of a high angle of attack delta wing undergoing oscillatory motion. The governing equations of fluid flow and dynamics of the multidisciplinary problem are solved using a time-accurate solution of the laminar, unsteady, compressible, full Navier- Stokes equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme and a four-stage Runge-Kutta scheme, respectively. The primary model under consideration consists of a 65 deg swept, sharp-edged, cropped delta wing of zero thickness at 20 deg angle of attack. In a freestream of Mach 0.85 and Reynolds number of 3.23 x 10(exp 6), the flow over the upper surface of the wing develops a complex shock system which interacts with the leading-edge primary vortices producing vortex breakdown. The effect of the oscillatory motion of the wing on the vortex breakdown and overall aerodynamic response is detailed to provide insight to the complicated physics associated with unsteady flows and the phenomenon of wing rock. Forced sinusoidal single and coupled mode rolling and pitching motion is presented for the wing in a transonic freestream. The Reynolds number, frequency of oscillation, and the phase angle are varied. Comparison between the single and coupled mode forced rolling and pitching oscillation cases illustrate the effects of coupling the motion. This investigation shows that even when coupled, forced rolling oscillation at a reduced frequency of 2(pi) eliminates the vortex breakdown which results in an increase in lift. The coupling effect for in phase forced oscillations show that the lift coefficient of the pitching-alone case and the rolling-moment coefficient of the rolling-alone case dominate the resulting response. However, with a phase lead in the pitching motion, the coupled motion results in a non-periodic response of the rolling moment. The second class of problems involve releasing the wing in roll to respond to the flowfield. Two models of sharp-edged delta wings, the previous 65 deg swept model and an 80 deg swept, sharp-edged delta wing, are used to observe the aerodynamic response of a wing free to roll in a transonic and subsonic freestream, respectively. These cases demonstrate damped oscillations, self-sustained limit cycle oscillations, and divergent rolling oscillations. Ultimately, an active control model using a mass injection system was applied on the surface of the wing to suppress the self-sustained limit cycle oscillation known as wing rock. Comparisons with experimental investigations complete this study, validating the analysis and illustrating the complex details afforded by computational investigations.

Menzies, Margaret Anne↗

NASA NTRS1980

Numerical analysis of complex fluid-flow systems

Very flexible computer-assisted numerical analysis is used to solve dynamic fluid-flow equations characterizing computer-controlled heat dissipation system developed for Space lab. Losses caused by bends, ties, fittings, valves, and like are easily included, and analysis can solve both steady-state and transient cases. It can also interact with parallel thermal analysis.

Holland, R. L.↗

NASA NTRS1987

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.↗

NASA NTRS1990

Global convergence of inexact Newton methods for transonic flow

In computational fluid dynamics, nonlinear differential equations are essential to represent important effects such as shock waves in transonic flow. Discretized versions of these nonlinear equations are solved using iterative methods. In this paper an inexact Newton method using the GMRES algorithm of Saad and Schultz is examined in the context of the full potential equation of aerodynamics. In this setting, reliable and efficient convergence of Newton methods is difficult to achieve. A poor initial solution guess often leads to divergence or very slow convergence. This paper examines several possible solutions to these problems, including a standard local damping strategy for Newton's method and two continuation methods, one of which utilizes interpolation from a coarse grid solution to obtain the initial guess on a finer grid. It is shown that the continuation methods can be used to augment the local damping strategy to achieve convergence for difficult transonic flow problems. These include simple wings with shock waves as well as problems involving engine power effects. These latter cases are modeled using the assumption that each exhaust plume is isentropic but has a different total pressure and/or temperature than the freestream.

Young, David P.↗

NASA NTRS1991

Efficient use of direct solvers for the calculation of compressible flows

While the direct solution of systems of linear equations resulting from fluid dynamic problems has generally not been practical in the past, it is presently demonstrated that the direct method is often more efficient than the most popular iterative schemes when constructed in such a way as to take advantage of presently available vector processing capabilities and large memory. The vertical line Gauss-Seidel algorithm was chosen as the iterative method to be compared with the direct method. It is fond that the direct method becomes efficient only when large residual reductions are desired.

Riggins, David W.↗

NASA NTRS1996

A Parallel Multilevel Spectral Element Scheme

A parallel multilevel strategy is developed using spectral (p) finite elements. Hierarchic bases are particularly well suited since the element matrices and vectors are nested and the multilevel projections easiliy performed. Since the basis degree is used to specify the multigrid level, an EBE strategy is natural br the multilevel technique. Results are presented for two candidate nonlinear elliptic transport problems: the augmented drift-diffusion equations of semiconductor device modeling and the stream function-vorticity equations of incompressible fluid dynamics.

Davis, M. B.↗

NASA NTRS1999

Control of Meridional Flow by a Non-Uniform Rotational Magnetic Field

The diffusive mass transfer of species during crystal growth in vertical ampoules is significantly affected by fluid flow in the liquid mother phase (melt). For electrically conductive melts, an elegant way of remotely inducing and controlling this flow is by utilizing a uniform rotational magnetic field (RMF) in the transverse direction. It induces an azimuthal flow which tends to homogenize the thermal and solutal fields. The rotating field also reduces the diffusion boundary layer, stabilizes temperature fluctuations, and promotes better overall crystal growth. For moderate strengths of the applied magnetic field (2-20 m Tesla) with frequencies of up to 400 Hz, the induced secondary meridional flow becomes significant. It typically consists of one roll at the bottom of the liquid column and a second roll (vortex) at the top. The flow along the centerline (ampoule axis) is directed from the growing solid (interface) towards the liquid (melt). In case of convex interfaces (e.g. in floating zone crystal growth) such flow behavior is beneficial since it suppresses diffusion at the center. However, for concave interfaces (e.g. vertical Bridgman crystal growth) such a flow tends to exacerbate the situation in making the interface shape more concave. It would be beneficial to have some control of this meridional flow- for example, a single recirculating cell with controllable direction and flow magnitude will make this technique even more attractive for crystal growth. Such flow control is a possibility if a non-uniform PNE field is utilized for this purpose. Although this idea has been proposed earlier, it has not been conclusively demonstrated so far. In this work, we derive the governing equations for the fluid dynamics for such a system and obtain solutions for a few important cases. Results from parallel experimental measurements of fluid flow in a mercury column subjected to non-uniform RMF will also be presented.

Mazuruk, Konstantin↗

NASA NTRS2001

Investigation of Conjugate Heat Transfer in Turbine Blades and Vanes

We report on work carried out to develop a 3-D coupled Finite Volume/BEM-based temperature forward/flux back (TFFB) coupling algorithm to solve the conjugate heat transfer (CHT) which arises naturally in analysis of systems exposed to a convective environment. Here, heat conduction within a structure is coupled to heat transfer to the external fluid which is convecting heat into or out of the solid structure. There are two basic approaches to solving coupled fluid structural systems. The first is a direct coupling where the solution of the different fields is solved simultaneously in one large set of equations. The second approach is a loose coupling strategy where each set of field equations is solved to provide boundary conditions for the other. The equations are solved in turn until an iterated convergence criterion is met at the fluid-solid interface. The loose coupling strategy is particularly attractive when coupling auxiliary field equations to computational fluid dynamics codes. We adopt the latter method in which the BEM is used to solve heat conduction inside a structure which is exposed to a convective field which in turn is resolved by solving the NASA Glenn compressible Navier-Stokes finite volume code Glenn-HT. The BEM code features constant and bi-linear discontinuous elements and an ILU-preconditioned GMRES iterative solver for the resulting non-symmetric algebraic set arising in the conduction solution. Interface of flux and temperature is enforced at the solid/fluid interface, and a radial-basis function scheme is used to interpolated information between the CFD and BEM surface grids. Additionally, relaxation is implemented in passing the fluxes from the conduction solution to the fluid solution. Results from a simple test example are reported.

Kassab, A. J.↗

NASA NTRS2005

Predicting Flows of Rarefied Gases

DSMC Analysis Code (DAC) is a flexible, highly automated, easy-to-use computer program for predicting flows of rarefied gases -- especially flows of upper-atmospheric, propulsion, and vented gases impinging on spacecraft surfaces. DAC implements the direct simulation Monte Carlo (DSMC) method, which is widely recognized as standard for simulating flows at densities so low that the continuum-based equations of computational fluid dynamics are invalid. DAC enables users to model complex surface shapes and boundary conditions quickly and easily. The discretization of a flow field into computational grids is automated, thereby relieving the user of a traditionally time-consuming task while ensuring (1) appropriate refinement of grids throughout the computational domain, (2) determination of optimal settings for temporal discretization and other simulation parameters, and (3) satisfaction of the fundamental constraints of the method. In so doing, DAC ensures an accurate and efficient simulation. In addition, DAC can utilize parallel processing to reduce computation time. The domain decomposition needed for parallel processing is completely automated, and the software employs a dynamic load-balancing mechanism to ensure optimal parallel efficiency throughout the simulation.

LeBeau, Gerald J.↗