SEARCH · Engineering Papers

Results for “Two-fluid flows”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Inertial Currents in Isotropic Plasma

The magnetospheric convection electric field contributes to Birkeland currents. The effects of the field are to polarize the plasma by displacing the bounce paths of the ions from those of electrons, to redistribute the pressure so that it is not constant along magnetic field lines, and to enhance the pressure gradient by the gradient of the bulk speed. Changes in the polarization charge during the convection of the plasma are neutralized by electrons in the form of field-aligned currents that close through the ionosphere. The pressure drives field-aligned currents through its gradient in the same manner as in quasi-static plasma, but with modifications that are important if the bulk speed is of the order of the ion thermal speed; the variations in the pressure along field lines are maintained by a weak parallel potential drop. These effects are described in terms of the field-aligned currents in steady state, isotropic, MED plasma. Solutions are developed by taking the MHD limit of two-fluid solutions and illustrated in the special case of Maxwellian plasma for which the temperature is constant along magnetic field lines. The expression for the Birkeland current density is a generalization of Vasyliunas' expression for the field-aligned current density in quasi-static plasma and provides a unifying expression when both pressure gradients and ion inertia operate simultaneously as sources of field-aligned currents. It contains a full account of different aspects of the ion flow (parallel and perpendicular velocity and vorticity) that contribute to the currents. Contributions of ion inertia to field-aligned currents will occur in regions of strong velocity shear, electric field reversal, or large gradients in the parallel velocity or number density, and may be important in the low-latitude boundary layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet.

Heinemann, M.↗

NASA NTRS1985

Coronal temperatures, heating, and energy flow in a polar region of the sun at solar maximum

The profiles of resonantly scattered Lyman-alpha coronal radiation have been used to determine the hydrogen kinetic temperature from 1.5 to 4 solar radius from the center of the polar region of the corona observed in 1980 at solar maximum. Hydrogen temperatures derived from the line profiles were found to decrease with height from 1.2 million K at r = 1.5 solar radii to 600,000 K at r = 4 solar radius. Comparison of the measured kinetic temperatures with predictions from a semiempirical two-fluid model showed evidence of a small amount of heating or a nonthermal contribution to the motions of coronal protons between 1.5 and 4 solar radius. The widths of the profiles confirmed an upper limit of 110 + or - 15 km/s on the rms magnitude of the line-of-sight component of velocities between 1.5 and 4 solar radius. Density measurements obtained in situ in the solar wind in the ecliptic were used to locate the sources of low speed and high-speed winds in the polar region. An eclipse photograph of the corona at solar maximum is provided.

Withbroe, G. L.↗

NASA NTRS1986

Characterization of simulated small-droplet fuel sprays

A two-fluid pneumatic atomizer operating at relatively high liquid and gas pressures produced water sprays that simulated small-droplet clouds of liquid fuel for use in studying vaporization and fuel-air mixing effects on combustor performance and emissions. To characterize the sprays, a scattered-light scanning instrument was developed and measurements of volume median or volume mean diameter, D sub V.5, were correlated with D sub O, W sub w, and W sub n, i.e., orifice diameter, water, and nitrogen gas flow rates, respectively, to give the general expression: D sub v.5 approx. (D sub o sup 0.2) (W sub w sup m) (W sub n sup n), which yields D sub v.5 = 45 (D sub o sup 0.2) (W sub w sup 0.2) (W sub w sup - 1.2). Values of D sub o, W sub w, and W sub n are in centimeters and grams/second, respectively. Farther downstream at an axial distance of 6.7 cm, exponent m increased from 0.2 to 0.4 and exponent n decreased from -1.2 to -1.0 and at a distance of 25 cm downstream of the atomizer, n decreased to -0.8. The increase in exponent m and decrease in exponent n was attributed to a loss of very small droplets from the spray due primarily to vaporization and diffusion effects on clouds of small droplets traveling a distance of 25 cm.

Ingebo, R. D.↗

NASA NTRS1986

Characterization of simulated small-droplet fuel sprays

A two-fluid pneumatic atomizer operating at relatively high liquid and gas pressures produced water sprays that simulated small-droplet clouds of liquid fuel for use in studying vaporization and fuel-air mixing effects on combustor performance and emissions. To characterize the sprays, a scattered-light scanning instrument was developed and measurements of volume median or volume mean diameter, D sub V.5, were correlated with D sub O, W sub w, and W sub n, i.e., orifice diameter, water, and nitrogen gas flow rates, respectively, to give the general expression: D sub v.5 approx. (D sub O sup 0.2) (W sub w sup m) (W sub n sup n), which yields D sub v.5 = 45 (D sub o sup 0.2) (W sub w sup 0.2) (W sub w sup -1.2). Values of D sub 0, W sub w, and W sub n are in centimeters and grams/second, respectively. Farther downstream at an axial distance of 6.7 cm, exponent m increased from 0.2 to 0.4 and exponent n decreased from -1.2 to -1.0 and at a distance of 25 cm downstream of the atomizer, n decreased to -0.8. The increase in exponent m and decrease in exponent n was attributed to a loss of very small droplets from the spray due primarily to vaporization and diffusion effects on clouds of small droplets traveling a distance of 25 cm.

Ingebo, R. D.↗

NASA NTRS1992

Effect of vaporization on cryogenic spray dropsize measurement

The fluid mechanics of multi-phase flow breakup of liquid nitrogen, LN2, jets injected into sonic velocity nitrogen gasflow, was experimentally investigated. A scattered-light scanning instrument was used to measure the characteristic dropsize, D(sub v.5), of LN2 sprays and to determine the effect of droplet vaporization on experimental dropsize measurements. Under sonic gas-velocity conditions, liquid-jet breakup occurred in the regime of aerodynamic stripping. As a result, the following correlation of volume-median drop diameter, D(sub v.5), with atomizing gas flowrate, W(sub g), was derived for two-fluid atomizers; with atomizing gas flowrate, W(sub g), was derived for two-fluid atomizers; (D(sub v.5))(exp -1) = k(sub c)(W(sub g))(sup n), where proportionally constant k(sub c) and exponent n are functions of droplet vaporization rate. Partially vaporized sprays were investigated and it was found that n = 1.11, which is considerably less than the value of 1.33 that is predicted by atomization theory. This was attributed to the evaporative loss of very small droplets. As a result, the following expression was obtained experimentally: (D(sub v.5e))(exp -1) = 301(W(sub g))(sup 1.11). Values of D(sub v.5), that existed prior to partial vaporization of the LN2 sprays, were calculated and the following expression was derived for originally unvaporized LN2 sprays: (D(sub v.5))(exp -1) = 285(W(sub g))(sup 1.33). This expression agrees well with atomization theory that predicts n = 1.33, for liquid jet breakup in high-velocity gasflow.

Ingebo, Robert D.↗

NASA NTRS1993

Effect of vaporization on cryogenic spray dropsize measurement

The fluid mechanics of multi-phase flow breakup of liquid nitrogen, LN2, jets injected into sonic velocity nitrogen gasflow, was experimentally investigated. A scattered-light scanning instrument was used to measure the characteristic dropsize, D(sub v.5), of LN2 sprays and to determine the effect of droplet vaporization on experimental dropsize measurements. Under sonic gas-velocity conditions, liquid-jet breakup occurred in the regime of aerodynamic stripping. As a result, the following correlation of volume-median drop diameter, D(sub v.5), with atomizing gas flowrate, W(sub g), was derived for two-fluid atomizers: (D(sub v.5)) (exp -1) = k(sub c) (W(sub g)) (sup n), where proportionally constant k(sub c) and exponent n are functions of droplet vaporization rate. Partially vaporized sprays were investigated and it was found that n = 1.11, which is considerably less than the value of 1.33 that is predicted by atomization theory. This was attributed to the evaporative loss of very small droplets. As a result, the following expression was obtained experimentally: (D(sub v.5e)) (exp -1) = 301 (W(sub g)) (sub 1.11). Values of D(sub v.5), that existed prior to partial vaporization of the LN2 sprays, were calculated and the following expression was derived for originally unvaporized LN2 sprays: (D(sub v.5)) (exp -1) = 285 (W(sub g)) (sub 1.33). This expression agrees well with atomization theory that predicts n = 1.33, for liquid jet breakup in high-velocity gasflow.

Ingebo, Robert D.↗

NASA NTRS1993

Particle injection and the structure of energetic-particle-modified shocks

A macroscopic 'self-consistent' nonlinear two-fluid model is developed for energetic-particle-modeled shocks. The model incorporates particle exchange between a thermal gas and an energetic population via a 'thermal leakage' mechanism. In the model, injection is regarded as a process whereby particles in the cooler thermal core of a total particle distribution are energized sufficiently to cross a 'momentum boundary' into the wings of the distribution. Particles in the wings are identified as 'energetic particles' or 'cosmic rays' and are assumed to propagate according to the diffusive transport equation of cosmic ray theory. Thermal particles are energized either as a result of adiabatic compression in a decelerating flow or directly via anomalous heating within a subshock. It is shown that the nature and model of particle injection directly affects the nature and dynamics of the shock, and hence the efficiency of particle acceleration, and that the shock itself regulates particle injection dynamically.

Zank, G. P.↗

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 NTRS2009

Ultimate Temperature of Pulse Tube Cryocoolers

An ideal pulse tube cryocooler using an ideal gas can operate at any temperature. This is not true for real gases. The enthalpy flow resulting from the real gas effects of He-3, He-4, and their mixtures in ideal pulse tube cryocoolers puts limits on the operating temperature of pulse tube cryocoolers. The discussion of these effects follows a previous description of the real gas effects in ideal pulse tube cryocoolers and makes use of models of the thermophysical properties of He-3 and He-4. Published data is used to extend the analysis to mixtures of He-3 and He-4. The analysis was done for pressures below 2 MPa and temperatures below 2.5 K. Both gases and their mixtures show low temperature limits for pulse tube cryocoolers. These limits are in the 0.5-2.2 K range and depend on pressure and mixture. In some circumstances, even lower temperatures may be possible. Pulse tube cryocoolers using the two-fluid properties of dilute 3He in superfluid He-4 appear to have no limit.

Kittel, Peter↗

NASA NTRS1977

A reexamination of two-fluid solar wind models

The two-fluid solar-wind equations have been solved by a method which is approximately 50 times faster than any previously developed, through the use of asymptotic expansions which are self-consistently iterated upon to find a solution that passes through the critical point. The energy assumptions in two-fluid solar-wind models are reexamined, and the conclusions are as follows: (1) proton thermal conduction may not be neglected, (2) the Coulomb logarithm must be calculated as a function of radius, and (3) the electron and proton temperatures at the base need not be equal, even when the time scale for energy exchange between the species is an order of magnitude smaller than the expansion time at the base. It is possible to reproduce reasonable quiet-time solar-wind parameters at 1 AU, but only if the proton temperature is approximately twice the electron temperature at 1 solar radius. This may indicate that extended proton heating is important in the outer solar corona. Winds with velocities at 1 AU of 450 km/s are generated without nonthermal energy deposition but require high proton temperatures as well as very low densities at the base. Higher-velocity solutions are not possible in a spherically symmetric geometry for reasonable particle fluxes at 1 AU, and it is suggested that these higher-velocity states probably require additional heating, acceleration mechanisms, or nonradial flow.

Nerney, S.↗

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 NTRS1973

Fundamental equations of a mixture of gas and small spherical solid particles from simple kinetic theory.

The fundamental equations of a mixture of a gas and pseudofluid of small spherical solid particles are derived from the Boltzmann equation of two-fluid theory. The distribution function of the gas molecules is defined in the same manner as in the ordinary kinetic theory of gases, but the distribution function for the solid particles is different from that of the gas molecules, because it is necessary to take into account the different size and physical properties of solid particles. In the proposed simple kinetic theory, two additional parameters are introduced: one is the radius of the spheres and the other is the instantaneous temperature of the solid particles in the distribution of the solid particles. The Boltzmann equation for each species of the mixture is formally written, and the transfer equations of these Boltzmann equations are derived and compared to the well-known fundamental equations of the mixture of a gas and small solid particles from continuum theory. The equations obtained reveal some insight into various terms in the fundamental equations. For instance, the partial pressure of the pseudofluid of solid particles is not negligible if the volume fraction of solid particles is not negligible as in the case of lunar ash flow.

Pai, S. I.↗

NASA NTRS1982

Statistical mechanics of dilute liquid mixtures of He3 in He4

The Woo-Tan-Massey quantum-mechanical theory of dilute mixtures of He-3 in liquid He-4 is extended to incorporate superfluid flow. Then statistical-mechanical principles are employed to find explicit formulas for thermodyanmic and hydrodynamic properties of mixtures, such as internal energy, free energy, entropy, pressure, chemical potential, and normal-fluid and superfluid densities. This paper is the first in a planned series that applies methods originally developed for pure liquid He-4 to dilute mixtures as a means of establishing a comprehensive theory that includes a theoretical bridge between a microscopic description of the mixtures and a hydrodynamic two-fluid model. Ultimately effects of applied electric and magnetic fields on dilute solutions are also planned to be treated.

Jackson, H. W.↗

NASA NTRS1989

Numerical simulation of forced convection in a two fluid layered system in a floating zone configuration

Forced convection due to rotational effects has been numerically studied in a stable two-fluid layered system in which the lighter fluid occupies the upper portion of the computational and physical domain. Results are presented for three different fluid height combinations for conditions prevailing in a microgravity environment. It is found that the heavier fluid tends to dominate the flow field, and that equal corotation and counterrotation cases produce nearly identical behavior in the flow and thermal fields.

Ramachandran, N.↗

NASA NTRS1994

Analysis of phase distribution phenomena in microgravity environments

The purpose of the research presented in this paper is to demonstrate the ability of multidimensional two-fluid models for bubbly two-phase flow to accurately predict lateral phase distribution phenomena in microgravity environments. If successful, this research should provide NASA with mechanistically-based analytical methods which can be used for multiphase space system design and evaluation, and should be the basis for future shuttle experiments for model verification.

Lahey, Richard T., Jr.↗

NASA NTRS1994

Local shear instabilities in weakly ionized, weakly magnetized disks

We extend the analysis of axisymmetric magnetic shear instabilities from ideal magnetohydrodynamic (MHD) flows to weakly ionized plasmas with coupling between ions and neutrals caused by collisions, ionization, and recombination. As part of the analysis, we derive the single-fluid MHD dispersion relation without invoking the Boussinesq approximation. This work expands the range of applications of these instabilities from fully ionized accretion disks to molecular disks in galaxies and, with somewhat more uncertainty, to protostellar disks. Instability generally requires the angular velocity to decrease outward, the magnetic field strengths to be subthermal, and the ions and neutrals to be sufficiently well coupled. If ionization and recombination processes can be neglected on an orbital timescale, adequate coupling is achieved when the collision frequency of a given neutral with the ions exceeds the local epicyclic freqency. When ionization equilibrium is maintained on an orbital timescale, a new feature is present in the disk dynamics: in contrast to a single-fluid system, subthermal azimuthal fields can affect the axisymmetric stability of weakly ionized two-fluid systems. We discuss the underlying causes for this behavior. Azimuthal fields tend to be stabilizing under these circumstances, and good coupling between the neutrals and ions requires the collision frequency to exceed the epicyclic frequency by a potentially large secant factor related to the magnetic field geometry. When the instability is present, subthermal azimuthal fields may also reduce the growth rate unless the collision frequency is high, but this is important only if the field strengths are very subthermal and/or the azimuthal field is the dominant field component. We briefly discuss our results in the context of the Galactic center circumnuclear disk, and suggest that the shear instability might be present there, and be responsible for the observed turbulent motions.

Blaes, Omer M.↗

NASA NTRS1995

Observations and physical interpretations of the solar wind flow properties as obtained from white light coronagraph aboard SPARTAN 201-01

The solar corona was observed with an externally occulted White Light Coronagraph (WLC) carried on the SPARTAN 201-1 spacecraft on 11-12 Apr. 1993. With observations from WLC and the ground based Mauna Loa White Light Coronagraph, a large number of polar plumes both in the north and south polar holes were traced from 1.16 to 5.5 Rs. Flow properties of the solar wind in coronal holes have been determined (Habbal et al., 1995) by using a two fluid model constrained by density profiles and scale height temperatures from the white light observations, and interplanetary measurements of the flow speed and proton mass flux from Ulysses' south polar passage. Provisions for acceleration by Alfven waves, as well as electron and proton heating, are included in the momentum and the energy equations respectively. The model computations fit remarkably well the empirical constraints of the two different density structures (plumes and coronal holes) for a range of input parameters. In this study we investigate the physical nature of the heating function used in the two-fluid model. Alfven waves have been suggested as the possible source of heating that accelerates the solar wind (Ofman and Davila, 1995). We utilize the density contrast observed in WLC data in the plume and ambient coronal hole region to estimate the Alfven wave frequencies responsible for heating these structures. The source heating function utilized in the two fluid model of the solar wind acceleration will be compared with the resonant Alfven wave heating function.

Guhathakurta, Madhulika↗

NASA NTRS1999

Effects of Gravity on Sheared Turbulence Laden with Bubbles or Droplets

The objective of this numerical/experimental study is to improve the understanding of the effects of gravity on the two-way interaction between dispersed particles (bubbles or liquid droplets) and the carrier turbulent flow. The first phase of the project considers isotropic turbulence. Turbulent homogeneous shear flows laden with droplets/bubbles will be studied in the next phase. The experiments reported here are concerned with the dispersion of liquid droplets by homogeneous turbulence under various gravitational conditions and the effect of these droplets on the evolution of the turbulence of the carrier fluid (air). Direct numerical simulations (DNS) of bubble - laden isotropic decaying turbulence are performed using the two-fluid approach (TF) instead of the Eulerian-Lagrangian approach (EL). The motivation for using the TF formulation is that EL requires considerable computational resources especially for the case of two-way coupling where the instantaneous trajectories of a large number of individual bubbles need to be computed. The TF formulation is developed by spatially averaging the instantaneous equations of the carrier flow and bubble phase over a scale of the order of the Kolmogorov length scale which, in our case, is much larger than the bubble diameter. On that scale, the bubbles are treated as a continuum (without molecular diffusivity) characterized by the bubble phase velocity field and concentration (volume fraction). The bubble concentration, C, is assumed small enough to neglect the bubble-bubble interactions.

Elghobashi, Said↗