Coupling Kinetic PIC to Multi-Fluid Plasma Models using Compatible Discretizations
Explore the source record for details and available documents.
Cyr, Eric Christopher↗
SEARCH · Search NASA
Results for “Multi-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.
Gauss' Law Preserving Methods for the Multi-Fluid Plasma Model
Explore the source record for details and available documents.
Logemann, Caleb↗
Closure models for the feedback of energetic particles on plasma turbulence
Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.
79 ASTRONOMY AND ASTROPHYSICS↗
Speed of Sound Measurements of Binary Mixtures of trans -1,2-Difluoroethylene (R-1132(E)) with Difluoromethane (R-32) or 2,3,3,3-Tetrafluoropropene (R-1234yf)
The speed of sound of two binary mixtures containing 0.336 mole fraction trans-1,2- difluoroethylene (R-1132(E)) with difluoromethane (R-32) and 0.435 mole fraction R-1132(E) with 2,3,3,3-tetrafluoropropene (R-1234yf) were measured with a dual-path pulse-echo instrument. The speed of sound was measured along several pseudo-isochores for each blend at temperatures ranging from 230 K to 325 K for blends with R-32 and to 342 K for blends with R- 1234yf. Measurements started at pressures just above each mixtures bubble point pressure and were limited to 8 MPa to avoid potential disproportionation reactions of the R-1132(E). The data were compared to multi-fluid models incorporating Helmholtz-energy-explicit equations of state (EOS) for each pure fluid.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Collisionless ion-electron energy exchange in magnetized shocks
Energy partition between ions and electrons in collisionless shocks has been a long-standing unsolved fundamental physical question. Here, we show that kinetic simulations of moderate Alfv´enic Mach number, magnetized, collisionless shocks reveal rapid, faster-than-Coulomb, energy exchange between ions and electrons when the plasma is sufficiently magnetized. Using kinetic and multi-fluid models with counter-streaming ions, we identify resonances between electron whistler and ion magnetohydrodynamic waves that account for this rapid energy exchange.
High-energy-density plasmas↗
Ten-moment fluid model for low-temperature magnetized plasmas
In this paper, a one-dimensional 10-moment multi-fluid plasma model is developed and applied to low-temperature magnetized plasmas. The 10-moment model solves for six anisotropic pressure terms, in addition to density and three components of fluid momentum, which allows the model to capture finite kinetic effects. The results are benchmarked against a 5-moment model, which assumes that the gas constituents follow a Maxwellian velocity distribution function (VDF), and a particle-in-cell Monte Carlo collision model that allows for arbitrary non-Maxwellian VDFs. The models are compared in a one-dimensional, low-temperature, partially magnetized plasma test case. The 10-moment results accurately reproduce the anisotropic temperature profile in low-temperature magnetized plasmas, where shear gradients exist due to the E×B drift. We discuss the mechanisms by which the anisotropic pressure can be generated in low-temperature magnetized plasmas. In addition, the importance of a self-consistent heat flux closure to the 10-moment model is studied, showing consistency with other models only when the assumptions of the underlying model are met. The 10-moment model allows for study of electron inertia effects and non-Maxwellian VDFs without the need for kinetic methods that are more computationally expensive.
Kuldinow, Derek Amur (ORCID:0000000319730196)↗
A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics
The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.
54 ENVIRONMENTAL SCIENCES↗
Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency
Here, we consider a kinetic description of multi-species gas mixture modeled with Bhatnagar-Gross-Krook (BGK) collision operators, in which the collision frequency varies not only in time and space but also with the microscopic velocity. In this model, the Maxwellians typically used in standard BGK operators are replaced by a generalization of such target functions, which are defined by a variational procedure. In this paper we present a numerical method for simulating this model, which uses an Implicit-Explicit (IMEX) scheme to minimize a certain potential function, mimicking the Lagrange functional that appears in the theoretical derivation. We show that theoretical properties such as conservation of mass, total momentum and total energy as well as positivity of the distribution functions are preserved by the numerical method, and illustrate its usefulness and effectiveness with numerical examples.
97 MATHEMATICS AND COMPUTING↗
Boiling and cavitation caused by transient heat transfer in superfluid helium-4
Superfluid helium-4 (He II) has been widely utilized as a coolant in various scientific and engineering applications due to its superior heat transfer capability. An important parameter required in the design of many He II based cooling systems is the peak heat flux $q^*$, which refers to the threshold heat flux above which boiling spontaneously occurs in He II. Past experimental and numerical studies showed that $q^*$ increases when the heating time $t_h$ is reduced, which leads to an intuitive expectation that very high $q^*$ may be achievable at sufficiently small $t_h$. Knowledge on how $q^*$ actually behaves at small $t_h$ is important for applications such as laser ablation in He II. Furthermore we present a numerical study on the evolution of the thermodynamic state of the He II in front of a planar heater by solving the He II two-fluid equations of motion. For an applied heat flux, we determine the heating time beyond which the He II near the heater transits to the vapor phase. As such, a curve correlating $q^*$ and $t_h$ can be obtained, which nicely reproduces some relevant experimental data. Surprisingly, we find that there exists a critical peak heat flux $q^*_c$, above which boiling occurs nearly instantaneously regardless of $t_h$. We reveal that the boiling in this regime is essentially cavitation caused by the combined effects of the first-sound and the second-sound waves in He II. Based on this physical picture, an analytical model for $q^*_c$ is developed, which reproduces the simulated $q^*_c$ values at various He II bath temperatures and hydrostatic head pressures. This work represents a major progress in our understanding of transient heat transfer in He II.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Boiling peak heat flux for steady inhomogeneous heat transfer in superfluid He 4
Superfluid helium-4 (He II) is a widely adopted coolant in scientific and engineering applications owing to its exceptional heat transfer capabilities. However, boiling can spontaneously occur on a heating surface in He II when the heat flux exceeds a threshold value $q^*$, referred to as the peak heat flux. While the parameter $q^*$ holds paramount importance in the design of He II based cooling systems, extensive research has primarily focused on its behavior in steady homogeneous heat transfer from a flat heating surface. For inhomogeneous heat transfer from curved surfaces, $q^*$ exhibits intricate dependance on parameters such as the He II bath temperature $T_b$, the immersion depth $h$, and the curvature radius $R_0$ of the heating surface. A comprehensive understanding on how $q^*$ depends on these parameters remains elusive. In this paper, we report our systematic study on $q^*$ for steady heat transfer from cylindrical and spherical heaters in He II. We compute $q^*$ for a wide range of parameter combinations $(T_b, h, R_0)$ by solving the He II two-fluid equations of motion. The generated data have allowed us to develop a robust correlation that accurately reproduces $q^*$ for all the parameter combinations we explored. Furthermore, our findings, particularly the establishment of the correlation, carry valuable implications for emergent applications that involve steady inhomogeneous heat transfer in He II systems.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Physical regimes of electrostatic wave-wave nonlinear interactions generated by an electron beam propagating in a background plasma
Not Available
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Three-dimensional modeling of a plasma in a strong azimuthal magnetic field
Three-dimensional magnetohydrodynamic (MHD) simulations are able to model the generation of disc-shaped plasma, driven by laser ablation from a current carrying rod in a pulsed-power machine producing azimuthal magnetic fields of 2–3 MG. The plasma at such extreme conditions is unique in that the parameter space for the plasma β and Hall parameter χ transition from below unity to greater than unity at different stages of the plasma generation. In simulations, the formation of the plasma disc in the azimuthal direction, is driven by heat flux from the laser spot and depends on the set of transport coefficients used in simulations. The most recent set of transport coefficients leads to the formation of plasma ejecta at the back end of the rod, that qualitatively matches experiments. Specifically, the cross-gradient Nernst effect, which twists magnetic field, is shown to have a large effect on the shape of the back end ejecta. In the direction along the axis of the rod there is propagation of perturbations from the disc as observed in experiments. In simulations the period of temperature perturbations is in good agreement with experimental results. Furthermore, an instability due to coupling of heat flux and the magnetic field advection provides a possible explanation for perturbation growth along the axis of the rod, and the instability growth rate is consistent with experimental results.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Nonequilibrium effects in high-gain inertial confinement fusion
Recent experimental demonstrations of ignition and target gain in inertial confinement fusion (ICF) have stimulated interest in exploring the fundamental physics of violent deuterium-tritium (DT) burn in high-gain ICF targets. A significant DT-burn fraction is a necessary condition for high energy gain and large neutron yields (>100MJ). Using classical molecular-dynamics (MD) simulations and a hybrid fluid-kinetic model, we examine how a large fraction of low-energy 𝛼 particles can kick D and T ions out of equilibrium in high-gain ICF targets. The MD results suggest that (1) temperatures of 𝑇 𝐷 and 𝑇 𝑇 can differ by as much as ∼20% of their mean temperature and (2) the deviation of the DT energy distribution from the Maxwell-Boltzmann function can exceed ∼30%. Some of these MD observations, such as the preferential heating of D ions by low-energy 𝛼 particles and the temperature separation, can be explained by a proposed hybrid fluid-kinetic model. Furthermore, the implication of such nonequilibrium effects on the DT reactivity is also discussed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Higher Harmonics in Multipactor Induced Plasma Ionization Breakdown near a Dielectric Surface
In this Letter, the novel physics of higher harmonic (HH) generation of the normal electric field near a dielectric surface is reported in multipactor induced plasma ionization breakdown, as determined by kinetic particle-in-cell simulations. The observed HH frequency is around ten times the fundamental rf driving frequency, but lower than the electron plasma frequency. A theory is constructed which indicates that stream plasma interaction-induced instability is the mechanism of HH generation in the collisional regime. The HH frequency and its corresponding growth rate of the HH oscillation amplitude from the theory are in good agreement with kinetic particle-in-cell simulations.
60 APPLIED LIFE SCIENCES↗
Data-driven discovery of reduced plasma physics models from fully kinetic simulations
Not Available
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Data-driven modeling of Landau damping by physics-informed neural networks
Not Available
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Combined Influence of Rotation and Scrape-Off Layer Drifts on Recycling Asymmetries in Tokamak Plasmas
Coupled 2D fluid-kinetic simulations of a DIII-D high confinement tokamak plasma show that plasma rotation coupled with drift effects near the plasma edge play a significant role in the creation of the observed poloidal distribution of neutrals. It is observed that including either drift or rotation effects enhances particle flux at the inner target in the case of ion 𝐵×∇𝐵 drift toward the 𝑋-point. However, the particle flux asymmetry is significantly higher with the combination of drifts and rotation than either effect alone. The heightened particle flux asymmetry allows for improved simulation of the strong in-out asymmetry of the Lyman-𝛼 brightness profiles measured in the experiment. Enhancement of radial transport of parallel momentum changes the upstream scrape-off layer flow pattern, increasing the fraction of deuterium flux that reaches the inboard divertor entrance while lowering that which arrives at the outboard. In conclusion, this Letter indicates that by combining drifts, rotation, and viscous coupling, existing boundary plasma models can achieve a satisfactory agreement with experimentally measured neutral asymmetries.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Ponderomotive electron physics captured in a single-fluid extended magnetohydrodynamics model
The ponderomotive force, arising from the interaction between electromagnetic waves and plasma, plays a critical role in laser fusion, astrophysical plasmas, and laser diagnostics. Traditionally, modeling this force requires multi-fluid or particle-in-cell simulations due to its strong coupling to electron-scale dynamics. In this work, we demonstrate that a one-fluid, two-temperature extended magnetohydrodynamics (XMHD) model—augmented with a generalized Ohm's law (GOL) including electron inertia—can accurately reproduce key ponderomotive effects. We derive the ponderomotive force within this framework using a phasor-based approach and then validate its nonlinear manifestations through direct numerical simulations in the PERSEUS code, where steepening and density modulation phenomena typically associated with kinetic-scale models are reproduced. These results and prior work establish XMHD as a robust and efficient alternative for modeling nonlinear laser–plasma dynamics, bridging the gap between ideal MHD and fully kinetic approaches.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗