Comparison of Reactive Flow Equilibrium Closure Assumptions in CTH .
Abstract not provided.
Ruggirello, Kevin↗
SEARCH · Search NASA
Results for “EQUILIBRIUM FLOW”
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.
Screening of resonant magnetic perturbation fields assuming various plasma flow models
A recently updated version of the MARS-F code [Y. Q. Liu et al., Phys. Plasmas 7, 3681 (2000); L. Li et al., Phys. Plasmas 25, 082512 (2018); and G. L. Xia et al., Nucl. Fusion 59, 126035 (2019)] is utilized to numerically investigate the plasma screening effect on the applied resonant magnetic perturbation (RMP) field, assuming various equilibrium flow models, including the toroidal flow, the parallel flow and their combinations, and poloidal and toroidal projections of the parallel flow. A parallel equilibrium flow with a uniform radial profile is found to have no effect on plasma screening of the RMP field. A sheared parallel flow, however, does change plasma screening. The poloidal projection of the parallel flow weakens plasma screening in the resistive-inertial regime. The effect on the favorable average curvature regime is found, however, to be non-monotonic. With the increasing flow speed, the poloidal projection first weakens Glasser-Green-Johnson (GGJ)-screening. Further increase in the flow speed results in enhanced GGJ-screening again. This non-monotonic behavior is related to the perturbed parallel shielding current, which appears also off the mode rational surface at fast flow due to additional resonances between the RMP perturbation and the sound wave continuum. These results indicate that flow induced plasma screening to the RMP field can have complicated characteristics, which, in turn, can have implications on the RMP field penetration into the plasma in experiments for controlling the edge localized modes.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Formulation of a General Collisional-Radiative Model for NO to Study Non-Equilibrium, Hypersonic Flows
Explore the source record for details and available documents.
Levin, Deborah A↗
Laser absorption spectroscopy diagnostics for non-equilibrium hypersonic flows
Explore the source record for details and available documents.
Goldenstein, Christopher [Purdue University]↗
Invariant regimes of Spencer scaling law for magnetic compression of rotating FRC plasma
Abstract The scaling laws for the magnetic compression of a toroidally rotating field reversed configuration (FRC) have been investigated in this work. The magnetohydrodynamics (MHD) simulations of the magnetic compression on rotating FRCs employing the NIMROD code (Sovinec et al 2004 J. Comput. Phys. 195 355), are compared with the Spencer’s one-dimensional (1D) theory (Spencer et al 1983 Phys. Fluids 26 1564) for a wide range of initial flow speeds and profiles. The toroidal flow can influence the scalings directly through the alteration of the compressional work as also evidenced in the 1D adiabatic model, and indirectly by reshaping the initial equilibrium. However, in comparison to the static initial FRC equilibrium cases, the pressure and the radius scalings remain invariant for the magnetic compression ratio B w 2 / B w 1 up to 6 in presence of the initial equilibrium flow, suggesting a broader applicable regime of the Spencer scaling law for FRC magnetic compression. The invariant scaling has been proven a natural consequence of the conservation of angular momentum of both fluid and magnetic field during the dynamic compression process.
Ma, Yiming↗
Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Final Report)
This research has focused on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aimed to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
A quasi-linear model of electromagnetic turbulent transport and its application to flux-driven transport predictions for STEP
A quasi-linear reduced transport model is developed from a database of high-β electromagnetic nonlinear gyrokinetic simulations performed with spherical tokamak for energy production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid kinetic ballooning mode turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat top of the conceptual design, but do not demonstrate how this state can be accessed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Effects of NBI-induced energetic ions on internal kink stability in HL-2M
A systematic investigation is carried out, studying the effect of the neutral beam injection induced energetic particles (EPs) on the n = 1 (n is the toroidal mode number) internal kink (IK) instability in the HL-2M tokamak, utilizing the MARS-F/K code [Liu et al., Phys. Plasmas 7, 3681 (2000) and 15, 112503 (2008)]. A high-beta sawteething HL-2M scenario, simulated by the TRANSP code [Breslau et al. Computer Software (2018)], is chosen for this study. Compared to the fluid model, non-perturbative magnetohydrodynamic (MHD)-kinetic hybrid computations with MARS-K show a generally stabilization effect on the IK, due to drift kinetic resonances associated with EPs. Further, the bounce resonance of trapped EPs has minor influence on the mode stability. In the absence of the plasma equilibrium flow and with the assumed particle pitch distribution, the transit resonance of co-current (countercurrent) passing EPs destabilizes (stabilizes) the IK. With plasma flow, both co- and countercurrent passing EPs tend to stabilize the mode, but the effect is stronger with the countercurrent particles. These modeling results provide useful guidance for interpreting MHD instabilities in the future high-performance experiments in HL-2M.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Spatially resolved measurements of plasma ion velocity distributions in a dipole magnetic field
The equilibrium flows of a plasma discharge in a dipole magnetic field are a topic of interest in low temperature plasma physics. Experimentalists typically rely on probe-based and line-integrated diagnostic techniques in these environments to describe plasma behavior. Presented here are measurements of argon ion dynamics with laser induced fluorescence techniques to provide insight into plasma dynamics in dipole magnetic fields with nonperturbative, spatially localized measurements. Simulation results from a Lagrangian approach to track particle orbits are compared to measured density profiles and provide evidence to support the mechanism distinguishing experimental configurations is the initial approach of particles. Applying a negative DC bias to the magnet induces strong E×B flows around the magnet, even exceeding the ion acoustic speed as measured far from the magnet. A strong enough bias also produces two distinct ion populations and provides a method for controlling the density gradient on the equator.
McLaughlin, Jacob W. (ORCID:0000000152661888)↗
Magnetorotational instability in a swirling partially ionized gas
ABSTRACT The magnetorotational instability (MRI) has been proposed as the method of angular momentum transport that enables accretion in astrophysical discs. However, for weakly ionized discs, such as protoplanetary discs, it remains unclear whether the combined non-ideal magnetohydrodynamic (MHD) effects of Ohmic resistivity, ambipolar diffusion, and the Hall effect make these discs MRI stable. While much effort has been made to simulate non-ideal MHD MRI, these simulations make simplifying assumptions and are not always in agreement with each other. Furthermore, it is difficult to directly observe the MRI astrophysically because it occurs on small scales. Here, we propose the concept of a swirling gas experiment of weakly ionized argon gas between two concentric cylinders threaded with an axial magnetic field that can be used to study non-ideal MHD MRI. For our proposed experiment, we derive the hydrodynamic equilibrium flow and a dispersion relation for MRI that includes the three non-ideal effects. We solve this dispersion relation numerically for the parameters of our proposed experiment. We find it should be possible to produce a non-ideal MRI in such an experiment because of the Hall effect, which increases the MRI growth rate when the vertical magnetic field is anti-aligned with the rotation axis. As a proof of concept, we also present experimental results for a hydrodynamic flow in an unmagnetized prototype. We find that our prototype has a small, but non-negligible, α-parameter that could serve as a baseline for comparison to our proposed magnetized experiment, which could be subject to additional turbulence from the MRI.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Three-beam hybrid fs/ps coherent anti-Stokes Raman scattering of rotation–vibration non-equilibrium
In this work, we demonstrate time-resolved rotational and vibrational temperature measurements to probe plasma non-equilibrium through a simple three-beam hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (fs/ps CARS) system. A single pump/Stokes pair is employed to generate both the pure-rotational and ro-vibrational Raman coherence. A novel phase-matching scheme, to our knowledge, is employed to spatially overlap the two CARS signals, simplifying the optical design and allowing signal selection by tuning a single spectrometer grating. Measurements were performed near the electrodes in a N 2 DC glow discharge. Both the rotational temperature and the vibrational populations up to $v$ = 8 were calculated from the separately measured single-shot pure-rotational and ro-vibrational CARS spectra. Strong rotation–vibration non-equilibrium was observed at both electrodes, with the cathode showing higher vibrational and rotational temperatures. In addition, non-Boltzmann behaviors were observed at both electrodes. This simplified approach enables measurements of rotational and vibrational temperatures with high spatial resolution in non-equilibrium flows.
42 ENGINEERING↗
Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport
A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Rapid depletion analysis of flowing-pebble reactor systems at equilibrium using SCALE
Several high-temperature gas-cooled reactor concepts (and more recently, salt-cooled designs such as the fluoride salt-cooled high-temperature reactor) feature core designs employing continuously circulating fuel pebbles. These reactor designs permit both continuous online refueling of fuel elements as well as higher overall achievable discharge burnups. However, rapid calculation of time-dependent fuel isotopic inventories proves challenging for this class of dynamic systems with current analysis tools. While iterative approaches employing coupled neutron transport have been developed to solve this issue, rapid depletion analysis techniques are needed to calculate time-dependent inventories for individual pebbles and batches (and thus the construction of full- core inventory at equilibrium). We propose a depletion analysis strategy for this type of system for cores at equilibrium. Drawing upon previous neutronic analysis of the PBMR-400 equilibrium core, we demonstrate the viability of developing collapsed one-group cross section libraries suitable for performing rapid depletion analyses with SCALE. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Bead–Spring Simulation of Ionomer Melts—Studying the Effects of Chain-Length and Associating Group Fraction on Equilibrium Structure and Extensional Flow Behavior
Ionomers are associative polymers with diverse applications ranging from selective membranes and high-performance adhesives to abrasion- and chemical-resistant coatings, insulation layers, vacuum packaging, and foamed sheets. Within equilibrium melt, the ionic or associating groups are known to form thermally reversible, associative clusters whose presence can significantly affect the system’s mechanical, viscoelastic, and transport properties. It is, thus, of great interest to understand how to control such clusters’ size distribution, shape, and stability through the designed choice of polymer architecture and the ionic groups’ fraction, arrangement, and interaction strength. In this work, we represent linear associating polymers using a Kremer–Grest type bead–spring model and perform large-scale MD simulations to explore the effect of polymer chain-length (l) and fraction (f s ) of randomly placed associating groups on the size distribution and stability of formed clusters. We consider different chain-lengths (below and above entanglement), varying fractions of associating groups (represented by ‘sticky’ beads) between 5 and 20%, and a fixed sticky–sticky nonbond interaction strength of four times that between regular non-associating beads. For all melts containing associating groups the equilibrium structure factor S(q) displays a signature ionomer peak at low wave vector q whose intensity increases with increasing f s and l. The average cluster size Nc increases with f s . However, the effect of chain-length on Nc appears to be pronounced only at higher values of f s . Under extensional flows, the computed stress (and viscosity) is higher at higher f s and l regardless of strain rate. Beyond a critical strain rate, we observe fragmentation of the associative clusters, which has interesting effects on the stress/viscous response.
36 MATERIALS SCIENCE↗
Stable Non-equilibrium Structures in Chiral Nematics under Microfluidic Flow
Cholesteric liquid crystals (CLCs) are compelling responsive materials with applications in next-generation sensing, imaging, and display technologies. While electric fields and surface treatments have been used to manipulate the molecular organization and, subsequently, the optical properties of CLCs, their response to controlled fluid flow has remained largely unexplored. Here, in this study, we investigate the influence of microfluidic flow on the structure of thermotropic CLCs that can exhibit structural coloration. We demonstrate that the shear forces that arise from microfluidic flow align the helical axis of CLCs; alignment is a prerequisite for harnessing the promising photonic properties of CLCs. Moreover, we show that microfluidic flow can generate non-equilibrium structures exhibiting photonic band gaps that are inaccessible in the stationary cholesteric phase. Our findings have implications for the use of CLCs in applications involving flow processing such as additive manufacturing.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Nonequilibrium design strategies for functional colloidal assemblies
Here, we use a nonequilibrium variational principle to optimize the steady-state, shear-induced interconversion of self-assembled nanoclusters of DNA-coated colloids. Employing this principle within a stochastic optimization algorithm allows us to identify design strategies for functional materials. We find that far-from-equilibrium shear flow can significantly enhance the flux between specific colloidal states by decoupling trade-offs between stability and reactivity required by systems in equilibrium. For isolated nanoclusters, we find nonequilibrium strategies for amplifying transition rates by coupling a given reaction coordinate to the background shear flow. We also find that shear flow can be made to selectively break detailed balance and maximize probability currents by coupling orientational degrees of freedom to conformational transitions. For a microphase consisting of many nanoclusters, we study the flux of colloids hopping between clusters. We find that a shear flow can amplify the flux without a proportional compromise on the microphase structure. This approach provides a general means of uncovering design principles for nanoscale, autonomous, functional materials driven far from equilibrium.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Nonequilibrium kinetic freeze-out properties in relativistic heavy ion collisions from energies employed at the RHIC beam energy scan to those available at the LHC
In this paper, we investigate the kinetic freeze-out properties in relativistic heavy ion collisions at different collision energies. We present a study of standard Boltzmann-Gibbs blast-wave (BGBW) fits and Tsallis blast-wave (TBW) fits performed on the transverse momentum spectra of identified hadrons produced in Au + Au collisions at collision energies of $\sqrt{s_{\rm{NN}}}=$7.7 –200 GeV at the Relativistic Heavy Ion Collider (RHIC), and in Pb + Pb collisions at collision energies of $\sqrt{s_{\rm{NN}}}=$ 2.76 and 5.02 TeV at the Large Hadron Collider (LHC). The behavior of strange and multistrange particles is also investigated. We found that the TBW model describes data better than the BGBW one overall, and the contrast is more prominent as the collision energy increases as the degree of nonequilibrium of the produced system is found to increase. From TBW fits, the kinetic freeze-out temperature at the same centrality shows a weak dependence of collision energy between 7.7 and 39 GeV, while it decreases as collision energy continues to increase up to 5.02 TeV. The radial flow is found to be consistent with zero in peripheral collisions at RHIC energies but sizable at LHC energies and central collisions at all RHIC energies. We also observed that the strange hadrons, with higher temperature and similar radial flow, approach equilibrium more quickly from peripheral to central collisions than light hadrons. The dependence of temperature and flow velocity on nonequilibrium parameter (q - 1) is characterized by two second-order polynomials. Both a and dξ from the polynomials fit, related to the influence of the system bulk viscosity, increase toward lower RHIC energies.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry
A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.
Mechanics↗