Efficacy of the radial pair potential approximation for molecular dynamics simulations in warm dense plasmas.
Abstract not provided.
Stanek, Lucas James↗
SEARCH · Search NASA
Results for “PLASMA DYNAMICS”
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.
3D modeling of the Staged Z-pinch with the FLASH code
The Staged Z-pinch (SZP) fusion concept is a magneto-inertial compression scheme developed by Magneto-Inertial Fusion Technologies, Inc. (MIFTI), in which small amounts of fusion fuel are brought to fusion-relevant conditions by passing multi-million amperes strong current through a cylindrical shell of high atomic number material. One- and two-dimensional modeling performed by MIFTI with the MACH2 code, suggests that net fusion energy gain can be achieved when currents in the 10 million amperes range compress a 50%-50% mixture of deuterium and tritium gas. In this project we enhance and use the FLASH code to execute high-fidelity, simulations of various SZP configurations, a collaboration between the Flash Center for Computational Science at the University of Rochester and the Laboratory for Laser Energetics, and MIFTI. FLASH is a high-performance computing, multi-physics, radiation magnetohydrodynamic (MHD) code with extended physics capabilities, which is developed at the Flash Center. The goal of this project is to assess the shell/fuel stability of the pinch to two- and three-dimensional MHD instabilities and to utilize FLASH’s extended physics capabilities to understand how extended-MHD effects impact implosion dynamics and plasma conditions at stagnation. The project is a natural extension of an ongoing collaborative effort between MIFTI and the Flash Center through the U.S. DOE Advanced Research Projects Agency-Energy (ARPA-E) BETHE program and can provide MIFTI with simulation capabilities that are currently beyond their reach with MACH2. Ultimately, MIFTI wants FLASH to become one of their simulation workhorses for reliable, high-fidelity SZP platform design. During the project, the Flash Center team developed a suite of in one- and two-dimensional FLASH simulations to model different variants of the SZP platform, performed code-to-code comparisons with MACH2, and used the FLASH code to design and model SZP experiments like the Double Eagle experiments that MIFTI recently performed. Also, the Flash Center team developed the capabilities of FLASH to be able to do three-dimensional simulations of pulsed-power experiments with the code for the first time. Two publications from this effort are currently under review and two more are in preparation. MIFTI has also committed to the use of FLASH for future SZP simulation efforts.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Non-equilibrium Dynamics of the Quark Gluon Plasma (Final Report)
During this funding period I continued my work on far-from-equilibrium formulations of hydrodynamics and bottomonium suppression in the quark-gluon plasma. I published 23 peer-reviewed articles and 9 proceedings contributions. The publications produced during this period have thus far received 705 citations. I delivered 22 talks, seminars, and colloquia about my work, most of which were invited. In addition to the work on hydrodynamics and heavy quarkonia, I also worked on high-order perturbative calculations of N=4 supersymmetric Yang-Mills thermodynamics and contributed to the efforts of the MUSES collaboration whose focus is developing a modular equation-of-state framework for use in modeling neutron stars and neutron star mergers.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
LLNL FESP Theory Highlights: October 2024
I. Novikau, I. Y. Dodin, E. A. Startsev, I. Joseph, Quantum algorithms for simulating dissipative linear and nonlinear dynamics of plasmas. Invited talk at the 66th Annual Meeting of the APS Division of Plasma Physics, Atlanta, Georgia. Novikau I., Dodin I.Y., Startsev E.A., Encoding of linear kinetic plasma problems in quantum circuits via data compression, Journal of Plasma Physics. 2024;90(4):805900401, doi:10.1017/S0022377824000795. We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation to be solved by using the quantum signal processing algorithm. The latter requires encoding of a matrix in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
MMMnet: A Neural Network Surrogate for Real-Time Transport Prediction Based on the Updated Multi-Mode Model
The Multi-Mode Model (MMM) is a physics-based anomalous transport model integrated into TRANSP for predicting electron and ion thermal transport, electron and impurity particle transport, and toroidal and poloidal momentum transport. While MMM provides valuable predictive capabilities, its computational cost, although manageable for standard simulations, is too high for real-time control applications. MMMnet, a neural network-based surrogate model, is developed to address this challenge by significantly reducing computation time while maintaining high accuracy. Trained on TRANSP simulations of DIII-D discharges, MMMnet incorporates an updated version of MMM (9.0.10) with enhanced physics, including isotopic effects, plasma shaping via effective magnetic shear, unified correlation lengths for ion-scale modes, and a new physics-based model for the electromagnetic electron temperature gradient mode. A key advancement is MMMnet’s ability to predict all six transport coefficients, providing a comprehensive representation of plasma transport dynamics. MMMnet achieves a two-order-of-magnitude speed improvement while maintaining strong correlation with MMM diffusivities, making it well-suited for real-time tokamak control and scenario optimization.
DIII-D↗
L-H Transition Dynamics in ITER-Similar Plasmas with Applied n=3 Magnetic Perturbations*
In ITER-similar shape (ISS), low rotation plasmas [1] in DIII-D [ =1.5-5x1019m-3, Bt=1.9-2T, Ip=1.5MA, q95~3.6, balanced neutral beam injection (NBI)], the L-H threshold power PLH with applied n=3 Resonant Magnetic Perturbations (RMP) [2,3] is found to increase strongly with decreasing collisionality. This is a concern for H-mode access in primarily ECH-heated ITER plasmas since RMP may be applied before the L-H transition in ITER to safely suppress the first ELM. The L-H power threshold increase with RMP is determined by initially accessing H-mode via applying balanced NBI and increasing steps in ECH power at high RMP field (applied via the DIII-D I-coil system consisting of six internal coils above and six coils below the outboard tokamak midplane [2]). NBI and ECH power are then removed to let the plasma return to L-mode. A second L-H transition is induced later in the same shot by applying NBI and stepping up ECH power again to obtain the threshold power with no applied RMP as a reference. Using this method to bracket PLH, a scan of PLH vs. RMP perturbation strength is then obtained via stepping down the I-coil current (RMP perturbation strength) at constant beam and ECH power, as shown in Fig. 1.
Schmitz, L.↗
Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition
We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Linear Theory of Electron-Plasma Waves at Arbitrary Collisionality
The dynamics of electron-plasma waves are described at arbitrary collisionality by considering the full Coulomb collision operator. The description is based on a Hermite- Laguerre decomposition of the velocity dependence of the electron distribution function. The damping rate, frequency, and eigenmode spectrum of electron-plasma waves are found as functions of the collision frequency and wavelength. A comparison is made between the collisionless Landau damping limit, the Lenard-Bernstein and Dougherty collision operators, and the electron-ion collision operator, finding large deviations in the damping rates and eigenmode spectra. A purely damped entropy mode, characteristic of a plasma where pitch-angle scattering effects are dominant with respect to collisionless effects, is shown to emerge numerically, and its dispersion relation is analytically derived. It is shown that such a mode is absent when simplified collision operators are used, and that like-particle collisions strongly influence the damping rate of the entropy mode.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Enhancing Fundamental Understanding of Magnetic Confinement in the Presence of Flows, Non-Conservative Forces, Dynamics, Collisions, Relaxation and Kinetic Effects
Pulsed high-power technology is used to create MHD-governed, dynamically-evolving, plasma-filled magnetic flux tubes. Advanced diagnostics measure the dynamical evolution and the measurements are then related to theoretical models. The dynamics are relevant to solar coronal plasma loops, astrophysical jets, and other situations where magnetic field lines intercept a boundary and current flows from the boundary into the plasma along the field lines. The morphology is time-dependent and determined by complex physics involving magnetic forces, hydrodynamic forces, and boundary conditions. Unlike magnetic fusion plasma experiments which are in approximate MHD equilibrium, fast accelerated plasma flows are a dominant feature. Typical parameters are: electric current 50-150 kA, applied voltage 3-6 kV, magnetic field 0.3 to 3 kG, experiment duration 10 μs, jet velocity 10’s of km/s, flux tube dimensions a few cm in radius and 10’s of cm in length. A scheme for measuring the wavelength of a quasi-neutral magnetic plasma wave using a single point of measurement was demonstrated both by an actual spacecraft (MMS) and in the lab. The experimental work is supported by a strong theoretical component. Theoretical activities included proposing a model for driving astrophysical jets (model motivated by lab experiments), and providing new insights into fast collisionless magnetic reconnection by demonstration of a new, more fundamental parameter than previously used.
47 OTHER INSTRUMENTATION↗
Micropinch formation dynamics in X pinches
High temporal resolution x-ray streak camera studies of micropinch formation in Cu hybrid x pinches reveal key plasma conditions. Analysis of Ne-like Cu lines indicate an average electron temperature of about 200 eV and 4.5×10 28 m -3 electron density. Here, the spectra suggest that the electron temperature jumps to about 1 keV, inferred from the continuum and the postcontinuum line emission that includes Li-like Cu lines. There is no sign of a rapid temperature change or a substantial surge in radiation emission during the 200 ps precontinuum x-ray burst, suggesting that the radiative collapse process does not play a major role in micropinch formation. Two-dimensional extended Magnetohydrodynamic (MHD) simulations, coupled to a collisional-radiative spectral analysis code, suggest the significance of the rapid radial implosion of high-temperature, low-density plasma, the axial outflow, and the dynamic plasma pressure in micropinch formation.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Subharmonic Shapiro steps in depinning dynamics of a two-dimensional solid dusty plasma modulated by one-dimensional nonlinear deformed periodic substrates
Langevin dynamical simulations are performed to investigate the depinning dynamics of a two-dimensional (2D) solid dusty plasma, which is modulated by one-dimensional (1D) nonlinear deformed periodic substrates, and also driven by the combination of the DC and AC forces. As the DC driving force increases gradually, pronounced subharmonic and harmonic Shapiro steps are discovered under the combined modulation of the deformed substrate and the external AC drive. Here, these observed subharmonic and harmonic Shapiro steps are attributed to the dynamic mode locking. The data analysis results indicate that the nonlinear deformed substrate strongly influences these subharmonic Shapiro steps, which can be accurately diagnosed using the fraction of sixfold coordinated particles. Furthermore, the diagnostic of the kinetic temperature clearly indicates the difference between harmonic and subharmonic Shapiro steps, i.e., the particle motion is slightly less synchronized at subharmonic Shapiro steps, caused by the unstable locations of the deformed substrate.
36 MATERIALS SCIENCE↗
Correlation and transport properties for mixtures at constant pressure and temperature
Transport properties of mixtures of elements in the dense plasma regime play an important role in natural astrophysical and experimental systems, e.g., inertial confinement fusion. In this paper, we present a series of orbital-free molecular dynamics simulations on dense plasma mixtures with comparison to a global pseudo ion in jellium model. Hydrogen is mixed with elements of increasingly high atomic number (lithium, carbon, aluminum, copper, and silver) at a fixed temperature of 100 eV and constant pressure set by pure hydrogen at 2g/cm 3 , namely, 370 Mbars. We compute ionic transport coefficients, such as self-diffusion, mutual diffusion, and viscosity for various concentrations. Small concentrations of the heavy atoms significantly change the density of the plasma and decrease the transport coefficients. The structure of the mixture evidences a strong Coulomb coupling between heavy ions and the appearance of a broad correlation peak at short distances between hydrogen atoms. Finally, the concept of an effective one component plasma is used to quantify the overcorrelation of the light element induced by the admixture of a heavy element.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Helium flux effects on bubble growth and surface morphology in plasma-facing tungsten from large-scale molecular dynamics simulations
We investigate helium flux effects on helium transport and surface evolution in plasma-facing tungsten using molecular dynamics. The simulations span two orders of magnitude, from ITER-relevant levels to those more typical of simulations published to date. Simulation times of up to 2.5 µs (corresponding to actual fluences of 1.4 x 10 21 m -2 ) are achieved, revealing concerted bubble-bursting events that are responsible for significant and very sudden changes in surface morphology. The depth distribution of helium depends very strongly on helium flux: helium self-trapping becomes more probable near the surface at high flux, and a layer of near-surface bubbles forms. Helium retention prior to the onset of bubble bursting is also substantially lower at low flux than it is at high flux. Surface features at low fluence are correlated with the positions of bubbles, but at high fluence, bubbles tend to coalesce, venting to the surface at one or more locations and leaving large interconnected cavities below the surface. Ruptured bubbles may serve as pathways deeper into the material, allowing helium to bypass the layer of near-surface bubbles and fill deeper, potentially much larger, bubbles that can produce more substantial surface features. Deeper bubbles also emit prismatic dislocation loops that can fill in cavities closer to the surface. Our results suggest that nearly all molecular dynamics simulations published to date are hampered by finite-size effects, and that helium flux is a very important parameter in determining the behavior of helium in plasma-facing components.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Nonlinear dynamics of cold magnetized non-relativistic plasma in the presence of electron-ion collisions
A numerical study is presented of the nonlinear dynamics of a magnetized, cold, non-relativistic plasma, in the presence of electron-ion collisions. The ions are considered to be immobile while the electrons move with non-relativistic velocities. The primary interest is to study the effects of the collision parameter, external magnetic field strength, and the initial electromagnetic polarization on the evolution of the plasma system.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Dynamical control in a prethermalized molecular ultracold plasma: Local dissipation drives global relaxation
Prethermalization occurs as an important phase in the dynamics of isolated many-body systems when coupled degrees of freedom relax in a subspace separated from a state of complete thermodynamic equilibrium by a gap in energy or other conserved quantity. Slow equilibration from a prethermal phase can localize the dynamics of natural and model systems despite high dimensionality and limited disorder. Here, we report the signature of an enduring prethermal regime of arrested relaxation in the molecular ultracold plasma that forms following the avalanche of a state-selected Rydberg gas of nitric oxide. For a wide range of initial conditions, this system enters a critical phase in which a density of NO + and electrons balances a population of Rydberg molecules. Electron collisions mix orbital angular momentum, scattering Rydberg electrons to states of very high ℓ . The rapid predissociation of molecules in states of low- ℓ purifies this high- ℓ ensemble. The angular momentum barrier that separates NO + ions and electrons creates an extraordinary gap between the plasma states of n ≈ ℓ , with measured n > 200 and penetrating states of ℓ = 0 , 1 , and 2. Evolution to a statistically equilibrated state of N and O atoms cannot occur without Rydberg electron penetration, and this gap blocks relaxation for a millisecond or more. A weak radiofrequency (RF) field drives ℓ -mixing electron collisions that erase this gap, causing wholesale dissipation. Remarkably, a local quantum-state transition promoting an exceedingly small fraction of the molecules in the prethermalized ensemble to a predissociative state also acts with a global effect, driving the entire system to a dissipative equilibrium. Using the Lindblad master equation, we illustrate qualitatively similar dynamics for a toy model of an open quantum system that consists of a localized set of spins on which dissipation acts at a single site as a gateway to equilibrium. Published by the American Physical Society 2025
Wang, Ruoxi↗
Phase space analysis for dynamics of three vortices of pure electron plasma trapped with Penning trap
A trilinear phase space analysis is applied for dynamics of three electron clumps confined with a Penning-Malmberg trap. We show that the Aref’s concept of phase space describe the observed features of the dynamics of three point vortices qualitatively. In vacuum, phase point P moves to physical region boundary in phase space, i.e. triangular configuration cannot be kept. With the addition of a low level background vorticity distribution (BGVD), the excursion of the clumps is reduced and the distance between P and stable point does not extend in the phase space.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Statistical properties of three-dimensional two-fluid plasma model
The nonlinear dynamics of incompressible non-dissipative two-fluid plasma model is investigated through classical Gibbs ensemble methods. Liouville's theorem of phase space for each wave number is proved, and the absolute equilibrium spectra for Galerkin truncated two-fluid model are calculated. In two-fluid theory, the equilibrium is built on the conservation of three quadratic invariants: the total energy and the self-helicities for ions and electrons fluid, respectively. The implications of statistic equilibrium spectra with arbitrary ratios of conserved invariants are discussed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Dynamic formation of stable current-driven plasma jets
Abstract Instabilities play a prominent role in determining the inherent structure and properties of magnetized plasma jets spanning both laboratory and astrophysical settings. The manner in which prominent unstable modes dynamically evolve remains key to understanding plasma behavior and control. In astrophysical phenomena, self-similar jets are observed to propagate over vast distances while avoiding breakup caused by unstable mode growth. However, the production of stable dense plasma jets in the laboratory has been limited by the onset of unstable modes that restrict jet lifetime, collimation, and scalability. In this work, we visualize the formation of stable laboratory-generated, dense, super-magnetosonic plasma jets in real time, and we identify an underlying mechanism that contributes to this behavior. The current-driven plasma jets generated in our experiments form a flowing Z-pinch, which is generally unstable to the m = 1 kink instability. Our results indicate that a stable dense plasma jet can be maintained for timescales over which a steady pinch current can be sustained, even at levels which would otherwise lead to rapid unstable mode growth and resultant pinch disassembly.
97 MATHEMATICS AND COMPUTING↗