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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Two-Fluid Models and Interfacial Area Transport in Microgravity Condition

The objective of the present study is to develop a two-fluid model formulation with interfacial area transport equation applicable for microgravity conditions. The new model is expected to make a leapfrog improvement by furnishing the constitutive relations for the interfacial interaction terms with the interfacial area transport equation, which can dynamically model the changes of the interfacial structures. In the first year of this three-year project supported by the U.S. NASA, Office of Biological and Physics Research, the primary focus is to design and construct a ground-based, microgravity two-phase flow simulation facility, in which two immiscible fluids with close density will be used. In predicting the two-phase flow behaviors in any two-phase flow system, the interfacial transfer terms are among the most essential factors in the modeling. These interfacial transfer terms in a two-fluid model specify the rate of phase change, momentum exchange, and energy transfer at the interface between the two phases. For the two-phase flow under the microgravity condition, the stability of the fluid particle interface and the interfacial structures are quite different from those under normal gravity condition. The flow structure may not reach an equilibrium condition and the two fluids may be loosely coupled such that the inertia terms of each fluid should be considered separately by use of the two-fluid model. Previous studies indicated that, unless phase-interaction terms are accurately modeled in the two-fluid model, the complex modeling does not necessarily warrant an accurate solution.

Ishii, Mamoru↗

DOE OSTI2020

Theoretical Review of the MFIX Fluid and Two-Fluid Models

MFIX (Multiphase Flow with Interphase eXchanges) is a general purpose code that can be used for describing the hydrodynamics, chemical reactions, and heat transfer of dense or dilute multiphase fluid flows. MFIX calculations give detailed information on pressure, temperature, composition, and velocity distributions in the defined system. This report provides an overview of MFIX’s expressions for single phase flow and continuous dispersed multiphase flow (i.e., the two-fluid model). In addition to the two fluid model (MFIX-TFM), MFIX contains a discrete element model (MFIX-DEM) and a multi-phase particle in cell model (MFIX-PIC) that are not reviewed within this document.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1966

Two-fluid model of the solar wind.

Two-fluid model of solar wind, assuming spherical symmetric flow without rotation, magnetic fields or viscous stress

SPACE RADIATION↗

NASA NTRS1973

Alfven waves in a two-fluid model of the solar wind.

A two-fluid model for the solar wind is presented that includes the presence of Alfven waves which originate at the sun. The effective pressure of these Alfven waves is considered, as well as a model representation for proton heating through nonlinear Alfven wave damping. The effects of rotation in the solar equatorial plane are allowed for. The implications of the computational results obtained are discussed.

Hollweg, J. V.↗

NASA NTRS1968

Two-fluid model of the solar wind.

Solar wind studied with two fluid model, using proton and electron fluid type heat equations, equation of motion and continuity equation

Hartle, R. E.↗

DOE OSTI2020

Benchmark of interfacial area concentration approaches for the two-fluid model in gas-dispersed condensing flow

Accurate predictions of interfacial area concentration (IAC) are required for closure of the condensation phase change term. The interfacial area transport equation (IATE) can dynamically predict interfacial area concentration and is increasingly being used in the benchmarking of phase change flows over static based correlations which require a flow regime map. For higher flow regimes, the two-group two-fluid IATE model can separately model the transport of spherical bubbles and cap/slug bubbles. Recently, a modified formulation of the one-group and two-group two-fluid model with interfacial area transport equation was proposed for gas-dispersed condensing flows. Three datasets have been used for validation consisting of condensing two-phase flow data in the unheated section of a vertical annulus with conditions spanning a range of pressures, inlet subcoolings and mass fluxes for benchmarking the various IAC approaches. For bubbly flows, the coupled void transport - IATE simulation results in good agreement in the predictions of the void fraction and the interfacial area concentration under nearly all sixty conditions. Considering the proposed two-group two-fluid IATE model formulation, the group-2 Nusselt number correlation is investigated to account for the appropriate heat transfer length scale for group-2 bubbles. Validation of the proposed two-group two-fluid model against the new dataset shows good agreement with experimental data with significant improvement over the one-group IATE model for conditions with substantial group-2 void fraction throughout the simulation domain.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1987

A two-fluid model for solar wind fluid with higher order moments

A spherically symmetric two-fluid model for the solar wind with higher-order moments is presented. In this model, continuity, momentum, temperature and heat flux equations for two components (electrons and protrons) in steady solar wind states are simultaneously solved by using a time-dependent method. This work is used to compare solutions of the steady-state solar wind with and without higher order moments, and to study the effects of thermal conduction. The coupling between electrons and protrons is also given special attention. The numerical solutions of the steady-state solar wind in both subsonic and supersonic regions between the Sun and 1 AU are obtained and graphically illustrated.

Wang, Jingfang↗

NASA NTRS1994

Flow properties of the solar wind obtained from white light data and a two-fluid model

The flow properties of the solar wind from 1 R(sub s) to 1 AU were obtained using a two fluid model constrained by density and scale height temperatures derived from white light observations, as well as knowledge of the electron temperature in coronal holes. The observations were obtained with the white light coronographs on SPARTAN 201-1 and at Mauna Loa (Hawaii), in a north polar coronal hole from 1.16 to 5.5 R(sub s) on 11 Apr. 1993. By specifying the density, temperature, Alfven wave velocity amplitude and heating function at the coronal base, it was found that the model parameters fit well the constraints of the empirical density profiles and temperatures. The optimal range of the input parameters was found to yield a higher proton temperature than electron temperature in the inner corona. The results indicate that no preferential heating of the protons at larger distances is needed to produce higher proton than electron temperatures at 1 AU, as observed in the high speed solar wind.

Habbal, Shadia Rifai↗

NASA NTRS1995

Flow properties of the solar wind obtained from white light data, Ulysses observations and a two-fluid model

Using the empirical constraints provided by observations in the inner corona and in interplanetary space. we derive the flow properties of the solar wind using a two fluid model. Density and scale height temperatures are derived from White Light coronagraph observations on SPARTAN 201-1 and at Mauna Loa, from 1.16 to 5.5 R, in the two polar coronal holes on 11-12 Apr. 1993. Interplanetary measurements of the flow speed and proton mass flux are taken from the Ulysses south polar passage. By comparing the results of the model computations that fit the empirical constraints in the two coronal hole regions, we show how the effects of the line of sight influence the empirical inferences and subsequently the corresponding numerical results.

Habbal, Shadia Rifai↗

NASA NTRS1995

Time-dependent simulation of oblique MHD cosmic-ray shocks using the two-fluid model

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibria consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(sub f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(sub c). We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, results in enhanced particle accretion through drift acceleration.

Frank, Adam↗

NASA NTRS1985

Two-fluid models of turbulence

The defects of turbulence models are summarized and the importance of so-called nongradient diffusion in turbulent fluxes is discussed. The mathematical theory of the flow of two interpenetrating continua is reviewed, and the mathematical formulation of the two fluid model is outlined. Results from plane wake, axisymmetric jet, and combustion studies are shown.

Spalding, D. B.↗

NASA NTRS1971

Two fluid model for thermal stratification in the Apollo cryogenic oxygen tanks

It is observed on Apollo flights that for oxygen densities greater than 42 lb/cu ft the heater operation time to achieve a desired pressure rise is half that which would be required if the energy were uniformly dispersed. A two fluid model was developed to describe this effect of stratification on the tank's pressurization behavior. The increase in pressure rise rate is attributed to compression of the bulk fluid in response to expansion of the heated fluid adjacent to the heater. It is assumed that during the heating operation the adjacent fluid retains nearly all the heater energy and has no mass interaction with the bulk fluid. The hot fluid temperatures predicted by the model are in reasonable correspondence with the observed heater temperatures. The model indicates that the fraction of the fluid interacting with the heater increases with decreasing density, from less than 0.1% at full tank to roughly 1% at 50% quantity, the lowest quantity at which stratification effects on pressure are observed.

Saxton, J. A.↗

DOE OSTI2021

Two-fluid model of rf current condensation in magnetic islands

The stabilization of tearing modes with rf waves is subject to a nonlinear effect, termed rf current condensation, that has the potential to greatly enhance and localize current driven within magnetic islands. Here we extend previous investigations of this effect with a two fluid model that captures the balance of diffusive and thermal equilibration processes within the island. We show that the effective power and resulting strength of the condensation effect can be greatly enhanced by avoiding collisional heat loss to the ions. The relative impact of collisions on the overall power balance within the island depends on the ratio of the characteristic diffusion timescale and the electron–ion equilibration time, rather than the latter alone. Although relative heat loss to ions increases with island size, the heating efficiency does as well. In particular, we show that the latter safely dominates for large deposition profiles, as is typically the case for lower hybrid current drive. This supports the possibility of passive stabilization of neoclassical tearing modes without the precise aiming of the rf waves required for electron cyclotron current drive stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

A two-fluid model of the solar wind

A method is presented for the integration of the two-fluid solar-wind equations which is applicable to a wide variety of coronal base densities and temperatures. The method involves proton heat conduction, and may be applied to coronal base conditions for which subsonic-supersonic solar wind solutions exist.

Sandbaek, O.↗

DOE OSTI2021

Nonlinear two-fluid modeling of plasma response to RMPs for the ELM control in the ITER baseline

Numerical modeling, combining the toroidal ideal MHD code GPEC and the nonlinear two-fluid MHD code TM1, was used for comprehensive studies of the plasma response to resonant magnetic perturbations (RMPs) with toroidal mode number n = 1–5 for controlling edge-localized modes (ELMs) in ITER for the standard operation scenario (15 MA Q = 10). Several issues related to RMP ELM control are investigated, including the optimization of the RMP coils configuration, the evaluation of the magnitude of density pump-out and the q95 windows of ELM suppression. GPEC calculates the magnetic response, which consistently includes the very important edge kink/peeling response to static magnetic perturbations. Furthermore, GPEC two-dimensional scans of the relative coil current phasing among the three rows of internal coils, at fixed coil current amplitude, reveal the optimal phasing for the RMP coil configuration with n = 1–5, respectively. The poloidal half wavelength of resonant mode at the edge of plasma calculated by GPEC indicates that the midplane row coils have the best resonant coupling with the plasma for n = 2, while the upper and lower row coils have the best resonant coupling with the plasma for n = 3. Based on the plasma kinetic equilibrium and the GPEC calculations of the magnetic response, TM1 was used to simulate the conditions for RMP field penetration in the ITER pedestal. TM1 shows magnetic island formation at the foot of ITER pedestal with RMP coil current threshold ranging from 4 kAt to 8 kAt with n = 2 to 4. These magnetic islands at the pedestal-foot lead to density pump-out, the magnitude of which scales as ${I}_{\text{RMP}}^{0.5}$ and ranges from 5% to 20% at the pedestal-top when scanning the coil current from 4 to 60 kAt. The density pump-out is found to be weaker for higher n RMP. The nonlinear TM1 simulations also show field penetration at the pedestal-top, where the threshold of RMP coil current depends on the q 95 . The alignment of the magnetic island and the location of the pedestal-top decreases the height and width of the pedestal to suppress ELMs. Simulations by two-dimensional scans of RMP coil current and q 95 reveal the accessible q 95 windows of ELM suppression for both n = 3 and 4 RMPs. The predicted q 95 windows of ELM suppression are very similar to the ones in currently operating tokamaks and the required RMP coil current for ELM suppression is less than 40–50 kAt, which is well within the designed capability for ITER. In addition, the simulations indicate that wide q 95 windows of ELM suppression may be accessible in ITER by operating with dominant n = 4 (or n = 5) RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Development of five-moment two-fluid modeling for Z-pinch physics

The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is γ τ A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of γ τ A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1974

Two-fluid model of the solar corona

A simple model of the lower corona which allows for a possible difference in the electron and proton temperatures is analyzed. With the introduction of a phenomenological heating term, temperature and density profiles are calculated for several different cases. It is found that, under certain circumstances, the electron and proton temperatures may differ significantly.

Knight, J. W.↗

NASA NTRS1974

Two-fluid model of the solar corona

Knight, J. W.↗