ANALYSIS OF A SKIMMER FOR A HIGH-INTENSITY MOLECULAR BEAM USING A THREE-FLUID MODEL
Analysis of skimmer for high intensity molecular beam using three-fluid model
SEARCH · Search NASA
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.
Analysis of skimmer for high intensity molecular beam using three-fluid model
It is commonly assumed in models of the solar wind-interstellar neutral hydrogen interaction that the ionized interstellar particles are quickly assimilated into the solar wind proton population and 'become indistinguishable' from the original solar wind. This assumption leads to the prediction that the solar wind proton temperature should increase with radius in the outer heliosphere. This temperature increase has not been observed. It is pointed out that assimilation of the interstellar particles to the point of indistinguishability takes place on the very long Coulomb collision time scale, and is not expected to occur within the heliosphere. Results are presented of a three-fluid model of the solar wind which consists of comoving thermal populations of protons of solar origin, protons produced by ionization of interstellar hydrogen, and electrons. The steady-state results yield a solar wind with a 'core' proton distribution which cools adiabatically, and a 'halo' of interstellar pickup protons which is maintained near 10 to the 7th K by the energy input of continued ionization and pickup. Such a distribution will not be observed to manifest the temperature increase at large heliocentric distances which is predicted from a one-fluid analysis. Further time-dependent calculations show a strong correlation between the densities of the solar wind and the interstellar pickup protons. It is suggested that the interstellar pickup population may be observable by the Voyager plasma instruments in low resolution mode during periods of high solar wind density and low solar wind temperature.
High energy, high intensity molecular beam flow is analyzed, using a skimmer which is a small inlet where the beam is separated from the core flow
Solar wind model with electrons, protons and alpha particles coupled by electric field and expanding due to pressure gradient and solar gravitation
Observations show that, when alpha particles and other minor ions in the solar wind plasma encounter fast shocks, they are heated more than protons and their bulk motion is decelerated less than protons. These effects have been studied using a three-fluid model, and the model predictions have been compared with observations. The comparison indicates that, for supercritical fast shocks, the three-fluid model can explain cross-shock minor ion heating which is significantly greater than that of protons. When the ratio of specific heats for minor ions, gamma (alpha), equals 2, both the lesser cross-shock deceleration and the greater heating of minor ions than of protons can be predicted by the model; thus, the minor ion heating through the shock transition region is consistent with the involvement of two degrees of freedom. Because the analysis is formulated in the de Heffmann-Teller frame of reference, the method is not valid for perpendicular shocks or when the angle is large. These results agree with the few extant observations and might be confirmed by further observations at the earth's bow shock.
Explore the source record for details and available documents.
Wave equations for disturbances in partly ionized gas due to current sources determined, using three-fluid model
This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.
A detailed model for a photospheric dynamo has been presented by Henoux and Somov (1987), who used the three-fluid model to treat the properties of the weakly ionized plasma. Only the equations for the two ionized components were solved. The equation for the neutral component is considered, and it is argued that the model is unacceptable becaused of an implied impossibly large unbalanced stress on the neutral gas. It is argued more generally that all existing photospheric dynamo models are untenable.
The solution was obtained without approximating the dispersion relation. This was made possible by the use of a computer method for carrying out the extremely large number of algebraic manipulations involved. The three-fluid model employed accounts for the interaction of the electron, ion, and neutral species. The solution for the e-folding distance describes the damping characteristics of each wave mode. Each of the e-folding solutions corresponds to one of the phase velocity solutions.
A self-consistent, three-fluid model of plasma transport and implanted ion acceleration in the unshocked solar wind is presented. The solar wind plasma is depleted by charge exchange with the expanding cometary exosphere, while implanted protons and heavy ions are produced by photoionization and charge transfer and lost by charge exchange. A generalized transport equation describing convection, adiabatic and diffusive velocity change, and the appropriate production terms is used to describe the evolution of the two cometary ion components, while the moments of the Boltzmann equation are used to calculate the solar wind density and pressure. The flow velocity is obtained self-consistently by combining the conservation equations of the three ion species. The results imply that second-order Fermi acceleration can explain the implanted spectra observed in the unshocked solar wind. Comparison of measured and calculated distribution indicates that spatial diffusion of implanted ions probably plays an important role in forming the energetic particle environment in the shock vicinity.
The coronal slow shock was predicted to exist embedded in large coronal holes at 4 to 10 solar radii. A three-fluid model was used to study the jumps in minor ions propertes across the coronal slow shock. The jump conditions were formulated in the de Hoffmann-Teller frame of reference. The Rankine-Hugoniot solution determines the MHD flow and the magnetic field across the shocks. For each minor ion species, the fluid equations for the conservation of mass, momentum, and energy can be solved to determine the velocity and the temperature of the ions across the shock. A simularity solution was also obtained for heavy ions. The results show that on the downstream side of the coronal slow shock the ion temperatures are nearly proportional to the ion masses for He, O, Si, and Fe in agreement with observed ion temperatures in the inner solar wind. This indicates that the possibly existing coronal slow shock can be responsible for the observed heating of minor ions in the solar wind.
The coronal slow shock has been predicted to exist embedded in large coronal holes at 4-10 solar radii. In this paper, a three-fluid model is used to study the jumps in minor ion properties across a slow shock such as the coronal slow shock. The jump conditions are formulated in the de Hoffmann-Teller frame of reference. The Rankine-Hugoniot solution determines the MHD flow and the magnetic field across the shocks. For each minor ion species, the fluid equations for the conservation of mass, momentum, and energy can be solved to determine the velocity and the temperature of the ions across the shock. Also obtained is a similarity solution for heavy ions. The results show that, on the downstream side of the slow shock, the ion temperatures are nearly proportional to the ion masses for He, O, Si, and Fe, in agreement with observed ion temperatures in the inner solar wind. This indicates that the possibly existing coronal slow shock can be responsible for the observed heating of minor ions in the solar wind.
Mathematical models of wave propagation in three-fluid plasma and numerical solutions to dispersion and propagation properties
In heavy-ion collision experiments, the global collectivity of final-state particles can be quantified by anisotropic flow coefficients (𝑣 𝑛 ) . The first-order flow coefficient, also referred to as the directed flow (𝑣 1 ), describes the collective sideward motion of produced particles and nuclear fragments in heavy-ion collisions. It carries information on the very early stage of the collision, especially at large pseudorapidity (𝜂), where it is believed to be generated during the nuclear passage time. Directed flow therefore probes the onset of bulk collective dynamics during thermalization, providing valuable experimental guidance to models of the pre-equilibrium stage. In 2018, the Event Plane Detector (EPD) was installed in STAR and used for the Beam Energy Scan phase-II (BES-II) data taking. The combination of EPD (2.1 < |𝜂| < 5.1) and high-statistics BES-II data enables us to extend the 𝑣 1 measurement to the forward and backward 𝜂 regions. In this paper, we present the measurement of 𝑣 1 over a wide 𝜂 range in Au + Au collisions at $\sqrt{s_{NN}}$ = 19.6 and 27 GeV using the STAR EPD. The results of the analysis at $\sqrt{s_{NN}}$ = 19.6 GeV exhibit excellent consistency with the previous PHOBOS measurement, while elevating the precision of the overall measurement. The increased precision of the measurement also revealed finer structures in heavy-ion collisions, including a potential observation of the first-order event-plane decorrelation. Multiple physics models were compared to the experimental results. Only a transport model and a three-fluid hybrid model can reproduce a sizable 𝑣 1 at large 𝜂 as was observed experimentally. The model comparison also indicates 𝑣 1 at large 𝜂 might be sensitive to the QGP phase transition.
Continuum source molecular beam facility, examining effect of source-skimmer distance variation on flux and spatial distribution of beam in collimation chamber downstream
As the nuclear industry explores new fuel designs to accommodate increased burnup, studying fuel behavior during loss-of-coolant accidents is essential for ensuring the safe operation of light-water reactors. Here, this study presents a computational model to simulate the intricate three-phase flow dynamics of fuel dispersal following cladding breach, where solid fuel particles and carrier fission gases interact with the surrounding liquid or gas medium. The model utilizes a Eulerian three-fluid framework to depict the phase interaction within interpenetrating continua, treating fragmented fuel as granular material. Closure relationships for bulk-flow properties of the solid phase and detailed interfacial exchange terms in high-temperature, high-speed gas-solid-liquid flow scenarios are proposed. The model was validated using adiabatic experimental results in the literature, focusing on characterizing fuel dispersal behavior during loss-of-coolant accidents conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched probability density functions reported in experimental studies. Key highlights of this study include the theoretical description of a coupled three-phase flow with phase change, interpretation of the high-pressure boundary conditions and insights into the transient behavior of gas-solid-liquid phase dynamics during a simulated high-speed dispersal event.
This work presents an investigation of the properties of the one-dimensional, two-fluid turbulence-driven solar wind model introduced by Hollweg and Johnson (1988). It is found that the model has serious difficulties in reproducing the observed high-speed wind at 1 AU. In particular, the model proton temperatures are lower than those observed by a factor of 2 or more, and the highest temperature models yield excessive wave intensities at 1 AU. It appears that the problem stems from the specific spatial distribution of heat deposition in the model. Thus this study does not rule out a turbulence-driven fast solar wind, since other forms of the turbulent evolution could probably achieve better results. A three-fluid version of the model is also presented to show that the addition of alpha particles does not significantly reduce the extreme proton temperatures displayed near the sun by the Hollweg and Johnson work. Finally, it is suggested that the additional heating needs to be located well beyond the critical point, implying that the heating mechanism for the fast solar wind is likely not the same as that heating the solar corona.