Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fluid 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

A comparison of numerical solutions to the inviscid equations of fluid motion

The second-order MacCormack method and the third-order Rusanov and Kutler-Warming-Lomax methods are applied to the inviscid 1-d Burgers' equation, wedge flow, and the problem of shock reflection from a rigid boundary. The numerical solution in each case is compared to the exact solution and the quantitative estimates of accuracy are obtained. Results of this study show that the third-order Kutler-Warming-Lomax or the Rusanov methods tuned for minimum dissipation or minimum dispersion provide the most accurate solution in each of the examples considered.

Anderson, D. A.↗
DOE OSTI2022

Assessing physics of ion temperature gradient turbulence via hierarchical reduced-model representations

In this work, the saturation physics of ion temperature gradient (ITG) turbulence is probed by studying how amplitudes and scalings with key parameters vary across a hierarchy of reduced models. The models derive from nonlinear fluid equations for toroidal ITG turbulence under approximations to the mode coupling interactions in wavenumber space and the representation of turbulent decorrelation. Mode coupling approximations include local-in-wavenumber treatments like the spectral density of flux in quasilinear theory, a truncation to three nonlinearly interacting waves, and the interactions in a cascade to high radial wavenumber mediated by a single zonal flow. Turbulent decorrelation treatments are based on the triplet correlation time with and without eddy damping. Model fidelity is assessed by the scalings and magnitudes of the squared amplitudes of unstable mode, stable mode, and zonal flow with respect to the flow-damping rate and temperature gradient. It is shown that all models reproduce fundamental scalings, provided they incorporate the coupling of unstable mode, stable mode, and zonal flow. Accurate amplitude prediction requires eddy damping in the triplet correlation time and proper representation of the zonal-flow drive by interactions associated with the radial wavenumber cascade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Ion Dynamics Model for Collisionless Radio Frequency Sheaths

Full scale reactor model based on fluid equations is widely used to analyze high density plasma reactors. It is well known that the submillimeter scale sheath in front of a biased electrode supporting the wafer is difficult to resolve in numerical simulations, and the common practice is to use results for electric field from some form of analytical sheath model as boundary conditions for full scale reactor simulation. There are several sheath models in the literature ranging from Child's law to a recent unified sheath model [P. A. Miller and M. E. Riley, J. Appl. Phys. 82, 3689 (1997)l. In the present work, the cold ion fluid equations in the radio frequency sheath are solved numerically to show that the spatiotemporal variation of ion flux inside the sheath, commonly ignored in analytical models, is important in determining the electric field and ion energy at the electrode. Consequently, a semianalytical model that includes the spatiotemporal variation of ion flux is developed for use as boundary condition in reactor simulations. This semianalytical model is shown to yield results for sheath properties in close agreement with numerical solutions.

Bose, Deepak↗
DOE OSTI2023

Evolution of the age-included nearest pair distribution in disperse multiphase flows

The age of the nearest particle pair is introduced as the difference between the current time and the most recent time when the nearest particle pair was formed. The evolution equation for the age-included nearest pair distribution function is derived. With the assumption of random destruction of the nearest particle pairs, the evolution equation predicts the exponential probability distribution of the ages of the nearest particle pairs. Particle-resolved numerical simulations with moving particles are performed to verify this prediction. The equation is then used to derive the evolution equation for the particle–fluid–particle (PFP) stress, which is known to be related to hyperbolicity of the two-fluid equations. It is found that the relaxation time of the age probability distribution is also the relaxation time for the PFP stress. Guided by the closure terms in the PFP stress evolution equation, we study kinematics of the nearest particle pairs in the particle-resolved simulations for flows caused by sedimentation of the particles with initially isotropic and homogeneous particle distributions. At the steady states, the particle Reynolds numbers are around 20. Anisotropy and inhomogeneity of particle distributions are seen to develop in these flows. The mean distances to the nearest particles and evolution of the distribution of the Voronoi cell volumes are studied. We also found the PFP stress is closely related to the changes in these inter-particle scale quantities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1986

Perturbation Method for Computational Fluid-Dynamical Equations

Perturbation technique yields accurate flow solutions using as few as one-fourth number of grid points required by finite-difference methods. Technique originally developed to solve Euler equations of two-dimensional, steady, inviscid transonic flow about airfoils, applicable to arbitrary equation sets and higher dimensions. New perturbations scheme used in design cycle where potential solutions generated routinely; Euler perturbation method used in second-cut analysis. Method also used to couple other equation sets.

Chow, L. J.↗
NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗
NASA NTRS1977

Discontinuous transitions in a current-carrying plasma

The properties of 'discontinuous' transitions in a current-carrying plasma are investigated by adopting an approach based on the two-fluid equations for a plasma. Braginskii's (1965) transport equations for a two-fluid plasma carrying a current are presented, and allowance is made for collisions with a background neutral gas as well as for effects of reflected particles. The speed of sound in a current-carrying plasma is derived so that the Mach number can be defined correctly as the ratio of the flow speed to the wave phase velocity in such a plasma. The various types of discontinuities permitted by the relevant jump relations in the absence of reflected particles are classified, the effect of reflected particles on such discontinuities is examined, and it is shown that the degree of indeterminacy of a discontinuity can be changed by the presence of reflected particles. Consideration of shocks containing a population of reflected particles and adiabatic transitions indicates that the pressure of reflected particles alone can give rise to both weak and strong transitions. The construction of discontinuous periodic flow solutions is also discussed.

Lee, H.-J.↗
NASA NTRS1987

The onset of Alfvenic turbulence

An investigation is conducted on how low-frequency MHD oscillations in a warm plasma may undergo a transition from a coherent state to one of turbulence. A driven/dissipative derivative nonlinear Schroedinger equation is derived from the fluid equations. The time evolution of an arbitrary spectrum of waves is analyzed in the case where one k-mode is unstable, with the rest damped. It is found that the transition from order to chaos in the driven/dissipative system is correlated with the existence or absence of 'breathing' solitons in the associated conservative system.

Ghosh, S.↗
NASA NTRS1982

Magnetohydrodynamic equilibrium. III - Helically symmetric fields. IV - Nonequilibrium of nonsymmetric hydrodynamic topologies

It is pointed out that plasma confinement in stable equilibrium states constitutes a fundamental and still unresolved question in plasma astrophysics and thermonuclear fusion research. The problem has two parts related to the equilibrium states themselves and their mechanical stability. The question of the existence of general solutions of the field and fluid equations for the steady dynamical interaction of inviscid compressible fluids of high electrical conductivity with magnetic and gravity fields is considered. In the absence of fluid motions, the presented equations become the familiar equations of magnetostatics supplemented by an equation of state. Starting from this simplest case of magnetostatic equilibrium, the investigation proceeds to the more complex case of magnetohydrodynamic equilibrium. Examples of helically symmetric fields are presented to illustrate the use of the formulation for treating the dynamics of helically symmetric hydromagnetic flows.

Tsinganos, K. C.↗
NASA NTRS1975

Computational aspects of the prediction of multidimensional transonic flows in turbomachinery

The analytical prediction and description of transonic flow in turbomachinery is complicated by three fundamental effects: (1) the fluid equations describing the transonic regime are inherently nonlinear, (2) shock waves may be present in the flow, and (3) turbomachine blading is geometrically complex, possessing large amounts of curvature, stagger, and twist. A three-dimensional computation procedure for the study of transonic turbomachine fluid mechanics is described. The fluid differential equations and corresponding difference operators are presented, the boundary conditions for complex blade shapes are described, and the computational implementation and mapping procedures are developed. Illustrative results of a typical unthrottled transonic rotor are also presented.

Oliver, D. A.↗
DOE OSTI2026

HygroThermFEM v1.0

HygroThermFEM is a Finite Element Method-based numerical calculation engine for solving 2-D heat and moisture transfer problems. This numerical engine is used in the THERM software tool, and its primary purpose is for the analysis of building envelopes (e.g., windows, walls, roofs, foundations, etc.). However, the engine can also be used for any heat and moisture transfer problems that require solving fundamental 2-D energy and mass transfer equations. Fluid flow solutions (Navier-Stokes momentum equations) are not included, but the correlations for various convection heat transfer situations are provided, including the translation of complex cavity geometries into those for which correlations are applicable. The calculation engine is written in C++ and includes an API for connecting to third-party tools.

Vidanovic, Dragan [Lawrence Berkeley National Labo↗
NASA NTRS1984

Dynamics of two-phase face seals

An analytic study is presented of the effects of phase change on load support for parallel and tapered face seals. Consideration is given to an adiabatic model for low Reynolds number flow. Numerical integration is carried out of the descriptive fluid equations, giving the opening force due to fluid film pressure. The loci of steady-state solutions are then plotted for water to provide curves of load support as a function of film thickness. For axial excursions of the seal rings, a quasi-steady transient analysis is made. It is found that the load support generated by fluid pressure can be multivalued for a given film thickness. Another finding is that axial disturbances of the seal rings may lead to sudden drops in load support generated by fluid pressure with three possible results. The first is that sufficient damping may permit the seal to return to the previous equilibrium operating position. The second is that the seal may collapse to an equilibrium position of smaller film thickness where face contact is more likely and a significantly higher operating temperature is assured. The third is that a limit cycle of self-sustained oscillation in the axial direction may occur if damping is sufficiently low.

Beeler, R. M.↗
DOE OSTI2022

Fluid and hybrid simulations of the ionization instabilities in Hall thruster

Low-frequency axial oscillations in the range of 5–50 kHz stand out as a pervasive feature observed in many types of Hall thrusters. While it is widely recognized that the ionization effects play the central role in this mode, as manifested via the large-scale oscillations of neutral and plasma density, the exact mechanism(s) of the instabilities remain unclear. To gain further insight into the physics of the breathing mode and evaluate the role of kinetic effects, a one-dimensional time-dependent full nonlinear low-frequency model describing neutral atoms, ions, and electrons is developed in full fluid formulation and compared to the hybrid model in which the ions and neutrals are kinetic. Both models are quasi-neutral and share the same electron fluid equations that include the electron diffusion, mobility across the magnetic field, and the electron energy evolution. The ionization models are also similar in both approaches. Further, the predictions of fluid and hybrid simulations are compared for different test cases. Two main regimes are identified in both models: one with pure low-frequency behavior and the other one, where the low-frequency oscillations coexist with high-frequency oscillations in the range of 100–200 kHz, with the characteristic time scale of the ion channel fly-by time, 100–200 kHz. The other test case demonstrates the effect of a finite temperature of injected neutral atoms, which has a substantial suppression effect on the oscillation amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Anti-symmetric representation of the extended magnetohydrodynamic equations

We introduce the anti-symmetric representation of the extended magnetohydrodynamic (MHD) equations. In this representation, the use of the anti-symmetric flux operator (∇·ν+ν·∇) results in conservation theorems with discrete analogs. Inherently robust numerical applications are achieved with little effort, and conservation to machine precision is possible with simple numerical schemes. Starting from the two-fluid equations, we construct a single-fluid MHD model based on generalized center-of-mass variables for the mass (√ρ), momentum (√ρν), and pressure (√ p ). This model is shown to possess identical conservation properties to the two-fluid system, with the only restriction being the use of a single temperature. Common approximations to the Braginskii heat fluxes and to the gyroviscous stress tensor are cast into our representation for convenience. The discrete conservation properties are verified using the classic Orszag–Tang vortex problem. In addition to the favorable mass, momentum, and energy conservation properties, the time reversibility of the simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Calculation of stiffness and damping coefficients for elastically supported gas foil bearings

The stiffness and damping coefficients of an elastically supported gas foil bearing are calculated. A perfect gas is used as the lubricant, and its behavior is described by the Reynolds equation. The structural model consists only of an elastic foundation. The fluid equations and the structural equations are coupled. A perturbation method is used to obtain the linearized dynamic coefficient equations. A finite difference formulation has been developed to solve for the four stiffness and the four damping coefficients. The effect of the bearing compliance on the dynamic coefficients is discussed in this paper.

Peng, J.-P.↗
NASA NTRS1994

Time-dependent evolution of cosmic-ray-modified shock structure: Transition to steady state

Steady state solutions to the two-fluid equations of cosmic-ray-modified shock structure were investigated first by Drury and Volk (1981). Their analysis revealed, among other properties, that there exist regions of upstream parameter space where the equations possess three different downstream solutions for a given upstream state. In this paper we investigate whether or not all these solutions can occur as time-asymptotic states in a physically realistic evolution. To do this, we investigate the time-dependent evolution of the two-fluid cosmic-ray equations in going from a specified initial condition to a steady state. Our results indicate that the time-asymptotic solution is strictly single-valued, and it undergoes a transition from weakly to strongly cosmic-ray-modified at a critical value of the upstream cosmic ray energy density. The expansion of supernova remnant shocks is considered as an example, and it is shown that the strong to weak transition is in fact more likely. The third intermediate solution is shown to influence the time-dependent evolution of the shock, but it is not found to be a stable time-asymptotic state. Timescales for convergence to these states and their implications for the efficiency of shock acceleration are considered. We also investigate the effects of a recently introduced model for the injection of seed particles into the shock accelerated cosmic-ray population. The injection is found to result in a more strongly cosmic-ray-dominated shock, which supports our conclusion that for most classes of intermediate and strong cosmic-ray-modified shocks, the downstream cosmic-ray pressure component is at least as large as the thermal gas pressure, independent of the upstream state. As a result, cosmic rays almost always play a significant role in determining the shock structure and dissipation and they cannot be regarded as test particles.

Donohue, D. J.↗