Search NASA⌕ Search

SEARCH · Search NASA

Results for “two fluid model”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Chaotic behaviour of high Mach number flows

The stability of the super-Alfvenic flow of a two-fluid plasma model with respect to the Mach number and the angle between the flow direction and the magnetic field is investigated. It is found that, in general, a large scale chaotic region develops around the initial equilibrium of the laminar flow when the Mach number exceeds a certain threshold value. After reaching a maximum the size of this region begins shrinking and goes to zero as the Mach number tends to infinity. As a result high Mach number flows in time independent astrophysical plasmas may lead to the formation of 'quasi-shocks' in the presence of little or no dissipation.

Varvoglis, H.↗
NASA NTRS1987

Hydrodynamical constraints on cosmic-ray acceleration in relativistic shocks

A two-fluid hydrodynamical model governing the transport of cosmic rays in a relativistically moving background plasma is developed. The equations are used to discuss the time asymptotic structure of a relativistic, plane-parallel shock wave significantly modified by first-order Fermi acceleration of cosmic rays. The model allows for an anisotropic cosmic-ray pressure tensor with pressure components parallel and perpendicular to the shock normal. It is surmised that the well-known energy threshold problem for diffusive shock acceleration that downstream particles have a large enough velocity to overtake the shock and subsequently scatter in the upstream medium before returning to the downstream medium may be expressed hydrodynamically by the condition V less than V(cr), where V is the fluid velocity relative to the shock and V(cr) is the relativistic cosmic-ray sound speed. Astrophysical implications of the results are briefly discussed.

Webb, G. M.↗
NASA NTRS1988

Modified non-linear Burgers' equations and cosmic ray shocks

A reductive perturbation scheme is used to derive a generalized non-linear Burgers' equation, which includes the effects of dispersion, in the long wavelength regime for the two-fluid hydrodynamical model used to describe cosmic ray acceleration by the first-order Fermi process in astrophysical shocks. The generalized Burger's equation is derived for both relativistic and non-relativistic cosmic ray shocks, and describes the time evolution of weak shocks in the theory of diffusive shock acceleration. The inclusion of dispersive effects modifies the phase velocity of the shock obtained from the lower order non-linear Burger's equation through the introduction of higher order terms from the long wavelength dispersion equation. The travelling wave solution of the generalized Burgers' equation for a single shock shows that larger cosmic ray pressures result in broader shock transitions. The results for relativistic shocks show a steepening of the shock as the shock speed approaches the relativistic cosmic ray sound speed. The dependence of the shock speed on the cosmic ray pressure is also discussed.

Zank, G. P.↗
NASA NTRS1989

On wave stability in relativistic cosmic-ray hydrodynamics

Wave stability of a two-fluid hydrodynamical model describing the acceleration of cosmic rays by the first-order Fermi mechanism in relativistic, cosmic-ray-modified shocks is investigated. For a uniform background state, the short- and long-wavelength wave speeds are shown to interlace, thus assuring wave stability in this case. A JWKB analysis is performed to investigate the stability of short-wavelength thermal gas sound waves in the smooth, decelerating supersonic flow upstream of a relativistic, cosmic-ray-modified shock. The stability of the waves is assessed both in terms of the fluid velocity and density perturbations, as well as in terms of the wave action. The stability and interaction of the short-wavelength cosmic-ray coherent mode with the background flow is also studied.

Webb, G. M.↗
NASA NTRS1989

Collisional and inertial dynamics of the ionospheric interchange instability

The dynamical properties and spectral characteristics of the ionospheric interchange instability in the collisional and inertial regimes are investigated numerically and analytically, using a two-dimensional fluid model. This model, which describes the collisional and inertial dynamics of the F region ionosphere, is similar to the models of Ott (1978) and Kintner and Seyler (1985). By numerically integrating the model equations, the nonlinear evolution of three interchange unstable configurations in the collisional and inertial regimes are followed, and their dynamical and spectral properties are studied. The solutions indicate that the dynamical properties of a collision-dominated regime are distinctly different from those of an inertia-dominated regime. The qualitative characteristics of a plasma bubble in both regime-types are described.

Zargham, S.↗
NASA NTRS1993

Three-dimensional solidification and melting using magnetic field control

A new two-fluid mathematical model for fully three dimensional steady solidification under the influence of an arbitrary acceleration vector and with or without an arbitrary externally applied steady magnetic field have been formulated and integrated numerically. The model includes Joule heating and allows for separate temperature dependent physical properties within the melt and the solid. Latent heat of phase change during melting/solidification was incorporated using an enthalpy method. Mushy region was automatically captured by varying viscosity orders of magnitude between liquidus and solidus temperature. Computational results were obtained for silicon melt solidification in a parallelepiped container cooled from above and from a side. The results confirm that the magnetic field has a profound influence on the solidifying melt flow field thus changing convective heat transfer through the boundaries and the amount and shape of the solid accrued. This suggests that development of a quick-response algorithm for active control of three dimensional solidification is feasible since it would require low strength magnetic fields.

Dulikravich, George S.↗
NASA NTRS2002

Fundamentals of Plasma Sails Propulsion Concept

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000], is based on the two-fluid plasma model and requires a 15-km frontal standoff distance (or 20-km cross-sectional diameter) in order for the magnetic bubble to absorb sufficient momentum from the SW to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Winglee et al. [2000] derived the above size requirement based on an extrapolation of their simulated results in which a system much smaller than a M2P2 was used (p. 21,074 of their study). We submit, however, that a fluid model has no validity for such a small scale size-even in the region near the plasma source! It is assumed in the MHD fluid model, normally applied to the magnetosphere, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the SW. In the case of the M2P2, however, the size of the magnetic bubble is actually less than or, at best, comparable to, the scale of these characteristic parameters and, therefore, a kinetic approach, which addresses the smallscale physical mechanisms involved, must be used. A fully three-dimensional version of the hybrid code is used in our M2P2 (Plasma Sails) studies was originally developed by Delamere et al. [1999]. The M2P2 plasma sail is an excellent application for this hybrid code. The primary advantage of this code is the seamless interface between fluid and kinetic descriptions of the ion populations. A kinetic description is not necessary for the dense inner regions of the magnetic bubble and tremendous computational savings can be realized by treating this dense, magnetized ion population with the fluid description. It is essential, however, that the outer bubble regions be treated kinetically as well as the SW protons. Comparison of full size M2P2 simulation based on 3D MHD and kinetic models show that kinetic treatment introduces much more asymmetry to the considering problem and the possibility of kinetic instabilities development.

Khazanov, G. V.↗
NASA NTRS2018

Computational Investigation of the Near-Field Plasma Plume in Ion-Ion Propulsion

A two-fluid computational model of plasma flows was developed to investigate the plume of an ion-ion propulsion system. The densities of positive and negative ions, along with the associated values of net charge, electric field, and electric potential were calculated throughout the domain. The computational domain was chosen to be large enough (25 thruster diameters downstream of the accelerating grids) to examine the neutralization of the plume. The resulting plasma electric potential and charge neutrality at the downstream end of the domain are shown. The results from this simulation are compared to existing literature on ion-ion plasma thrusters.

Sankaran, K.↗
NASA NTRS2018

Computational Investigation of the Near-Field Plasma Plume in Ion-Ion Propulsion

A two-fluid numerical model of plasma flows was developed to investigate the plume of an ion-ion propulsion system. The densities of positive and negative ions, and the associated values of net charge, electric field, and electric potential were calculated as a function of time throughout the domain. The computational domain was chosen to be large enough (25 thruster diameters downstream of the exit plane) to allow for examining the neutralization of the plume. The resulting plasma electric potential and charge neutrality at the downstream end of the domain are shown and they indicate that it is possible to alternatively accelerate oppositely charged ions without the need for an electron-emitting neutralizer and without facing any electric potential hills that could cause stagnation. However, compared to existing literature on ion-ion plasma thrusters, the results from this simulation predict a longer length-scale for voltage decay.

Sankaran, K.↗
NASA NTRS1994

A two-fluid approximation for calculating the cosmic microwave background anisotropies

We present a simplified treatment for calculating the cosmic microwave background anisotropy power spectrum in adiabatic models. It consists of solving for the evolution of a two-fluid model until the epoch of recombination and then integrating over the sources to obtain the cosmic microwave background (CMB) anisotropy power spectrum. The approximation is useful both for a physical understanding of CMB anisotropies as well as for a quantitative analysis of cosmological models. Comparison with exact calculations shows that the accuracy is typically 10%-20% over a large range of angles and cosmological models, including those with curvature and cosmological constant. Using this approximation we investigate the dependence of the CMB anisotropy on the cosmological parameters. We identify six dimensionless parameters that uniquely determine the anisotropy power spectrum within our approximation. CMB experiments on different angular scales could in principle provide information on all these parameters. In particular, mapping of the Doppler peaks would allow an independent determination of baryon mass density, matter mass density, and the Hubble constant.

Seljak, Uros↗
NASA NTRSDate not supplied

Validation of a Two-Phase CFD Model for Predicting Tank Self-Pressurization in the Ground-Based K-Site Experiment

A two-phase CFD model for self-pressurization of a cryogenic storage tank partially filled with liquid hydrogen is presented using the Volume-Of-Fluid approach for modeling two-phase flow, as well as interfacial heat, mass and momentum transfer between the liquid and vapor regions. The CFD model is validated against self-pressurization experiment performed using the K-site flightweight hydrogen storage tank at NASA Glenn Research Center 1 . Laminar and turbulent simulations are performed together with conjugate heat transfer analysis. Effects of turbulence, the value of accommodation coefficient used for predicting phase change rates, as well as tank wall geometry are presented and discussed. Predicted tank pressures, fluid and wall temperatures are compared with the experimental data at 49% fill level and two different heat loads for validating this CFD model.

Computational Fluid Dynamics↗
NASA NTRS1991

Particle flows to shape and voltage surface discontinuities in the electron sheath surrounding a high voltage solar array in LEO

This paper discusses the numerical modeling of electron flows from the sheath surrounding high positively biased objects in LEO (Low Earth Orbit) to regions of voltage or shape discontinuity on the biased surfaces. The sheath equations are derived from the Two-fluid, Warm Plasma Model. An equipotential corner and a plane containing strips of alternating voltage bias are treated in two dimensions. A self-consistent field solution of the sheath equations is outlined and is pursued through one cycle. The electron density field is determined by numerical solution of Poisson's equation for the electrostatic potential in the sheath using the NASCAP-LEO relation between electrostatic potential and charge density. Electron flows are calculated numerically from the electron continuity equation. Magnetic field effects are not treated.

Metz, Roger N.↗
NASA NTRS1992

Ablative subduction - A two-sided alternative to the conventional subduction model

The plausibility of a two-sided fluid-based model of lithospheric subduction that is based upon current views of lithospheric structure is examined. In this model the viscous lower lithosphere flows downward, and the brittle upper lithosphere deforms in passive response. This process is potentially double-sided, since it is found that even a buoyant plate can be dragged downward by a dense descending neighbor. Thus an apparent overriding plate may be worn away by a process of viscous ablation, with the rate of ablation a function of plate buoyancy. This process, called 'ablative subduction,' makes it possible to simply interpret observations concerning slab profiles, interplate seismicity, back arc tectonics, and complex processes such as double subduction and subduction polarity reversal. When experiments modeling the evolution of simple fluid 'slabs' are performed, slab profile is found to be strongly influenced by ablation in the overriding plate. When ablation is weak, as when a buoyant continent borders the trench, deformable slabs adopt shallow Andean-style profiles.

Tao, Winston C.↗
NASA NTRS1996

Simulation of Two-Fluid Flows by the Least-Squares Finite Element Method Using a Continuum Surface Tension Model

In this paper a numerical procedure for simulating two-fluid flows is presented. This procedure is based on the Volume of Fluid (VOF) method proposed by Hirt and Nichols and the continuum surface force (CSF) model developed by Brackbill, et al. In the VOF method fluids of different properties are identified through the use of a continuous field variable (color function). The color function assigns a unique constant (color) to each fluid. The interfaces between different fluids are distinct due to sharp gradients of the color function. The evolution of the interfaces is captured by solving the convective equation of the color function. The CSF model is used as a means to treat surface tension effect at the interfaces. Here a modified version of the CSF model, proposed by Jacqmin, is used to calculate the tension force. In the modified version, the force term is obtained by calculating the divergence of a stress tensor defined by the gradient of the color function. In its analytical form, this stress formulation is equivalent to the original CSF model. Numerically, however, the use of the stress formulation has some advantages over the original CSF model, as it bypasses the difficulty in approximating the curvatures of the interfaces. The least-squares finite element method (LSFEM) is used to discretize the governing equation systems. The LSFEM has proven to be effective in solving incompressible Navier-Stokes equations and pure convection equations, making it an ideal candidate for the present applications. The LSFEM handles all the equations in a unified manner without any additional special treatment such as upwinding or artificial dissipation. Various bench mark tests have been carried out for both two dimensional planar and axisymmetric flows, including a dam breaking, oscillating and stationary bubbles and a conical liquid sheet in a pressure swirl atomizer.

Wu, Jie↗
NASA NTRS1988

Numerical models of thermally and mechanically coupled two-layer convection of highly viscous fluids

Thermal convection in the earth mantle is investigated by means of numerical simulations. The mantle models comprise two horizontal layers of viscous incompressible fluid (with identical or differing properties) separated by a fixed horizontal interface; an isothermal, fixed-heat-flux, or insulating lower boundary; and heat supplied either from below or internally (in the lower layer only, equally in both layers, or primarily in the upper layer). The mathematical formulation of the models is explained, and results for linear stability and finite-amplitude convection are presented in extensive tables and graphs and discussed in detail. Particular attention is given to the presence of thermal coupling without interface distortion in many cases, the predominance of long-wavelength cells in the finite-amplitude models, and the large (50 percent) temperature difference across the interface in all cases.

Ellsworth, Kirk↗
NASA NTRS2004

Aeolian Sediment Transport Pathways and Aerodynamics at Troughs on Mars

Interaction between wind regimes and topography can give rise to complex suites of aeolian landforms. This paper considers aeolian sediment associated wit11 troughs on Mars and identifies a wider range of deposit types than has previously been documented. These include wind streaks, falling dunes, "lateral" dunes, barchan dunes, linear dunes, transverse ridges, sand ramps, climbing dunes, sand streamers, and sand patches. The sediment incorporated into these deposits is supplied by wind streaks and ambient Planitia sources as well as originating within the trough itself, notably from the trough walls and floor. There is also transmission of sediment between d~~neTsh. e flow dynamics which account for the distribution of aeolian sediment have been modeled using two-dimensional computational fluid dynamics. The model predicts flow separation on the upwind side of the trough followed by reattachment and acceleration at the downwind margin. The inferred patterns of sediment transport compare well with the distribution of aeolian forms. Model data indicate an increase of wind velocity by approx. 30 % at the downwind trough margin. This suggests that the threshold wind speed necessary for sand mobilization on Mars will be more freq~~entmlye t in these inclined locations.

Bourke, Mary C.↗
NASA NTRS2017

Two-Pendulum Model of Propellant Slosh in Europa Clipper PMD Tank

Model propellant slosh for Europa Clipper using two pendulums such that controls engineers can predict slosh behavior during the mission. Importance of predicting propellant slosh; (1) Sloshing changes CM (center of mass) of spacecraft and exerts forces and torques on spacecraft. (2) Avoid natural frequencies of structures. (3) Size ACS (Attitude Control Systems) thrusters to counteract forces and torques. Can model sloshing fluid as two pendulums with specific parameters (mass, length, damping),

Two-Pendulum Model of Propellant Slosh in Europa C↗