Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plasma kinetic simulation”

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.

At least 829 records · Page 46

Characterization of laser-accelerated proton beams from a 0.5 kJ sub-picosecond laser for radiography applications

Laser-accelerated ion beams show promise for many applications, including high-resolution flash imaging of static or dynamic objects in next-generation radiography to probe materials and plasmas in extreme environments and inertial confinement fusion. To scale up ion beam production for radiography applications, we conducted experiments using sub-picosecond lasers up to 0.5 kJ at the OMEGA-EP facility to characterize proton beams from solid targets, primarily CH/CD sub-micron thin films from which ion beams were also used for static and dynamic radiography for the first time. For standalone sub-micron thin CH films, the highest detected proton energy is in the range of 72–97 MeV. Proton beams with highest energy near or above 60 MeV at full laser energy and similar beam profiles are also measured from low-density CD foams and flat CH foil target of micrometer-scale thickness. The ~ 700–800 nm CH/CD foils achieve the highest ion yield among the targets tested. For sub-micron thin films, the laser prepulse can expand the target and lead to complex interactions, which is simulated using coupled hydrodynamic and two-step kinetic models. Simulations suggest the presence of a micrometer-scale preplasma plateau with near-critical density and further indicate that target normal sheath acceleration, electron heating from Relativistic transparency in the preplasma plateau, and background proton reflection from carbon ion front at the rear side contribute to the resulting proton spectrum from these sub-micron thin targets at various stages. These proton beams show strong potential for radiography and for production of secondary sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Novel SOLPS-ITER simulations of X-point target and snowflake divertors

Abstract The design and understanding of alternative divertor configurations may be crucial for achieving acceptable steady-state heat and particle material loads for magnetic confinement fusion reactors. Multiple X-point alternative divertor geometries such as snowflakes and X-point targets have great potential in reducing power loads, but have not yet been simulated widely in codes with kinetic neutrals. This paper discusses recent changes made to the SOLPS-ITER code to allow for the simulation of X-point target and low-field side snowflake divertor geometries. Snowflake simulations using this method are presented, in addition to the first SOLPS-ITER simulation of the X-point target. Analysis of these results show reasonable consistency with the simple modelling and theoretical predictions, supporting the validity of the methodology implemented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ion-temperature- and density-gradient-driven instabilities and turbulence in Wendelstein 7-X close to the stability threshold

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the δƒ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic simulations verifying reconnection rates measured in the laboratory, spanning the ion-coupled to near electron-only regimes

The rate of reconnection characterizes how quickly flux and mass can move into and out of the reconnection region. In the Terrestrial Reconnection EXperiment (TREX), the rate at which the antiparallel asymmetric reconnection occurs is modulated by the presence of a shock and a region of flux pileup in the high-density inflow. Simulations utilizing a generalized Harris-sheet geometry have tentatively shown agreement with TREX's measured reconnection rate scaling relative to system size, which is indicative of the transition from ion-coupled toward electron-only reconnection. Here, in this study, we present simulations tailored to reproduce the specific TREX geometry, which confirm both the reconnection rate scale and the shock jump conditions previously characterized experimentally in TREX. The simulations also establish an interplay between the reconnection layer and the Alfvénic expansions of the background plasma associated with the energization of the TREX drive coils; this interplay has not yet been experimentally observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Identifying the Growth Phase of Magnetic Reconnection Using Pressure‐Strain Interaction

Abstract Magnetic reconnection often initiates abruptly and then rapidly progresses to a nonlinear quasi‐steady state. While satellites frequently detect reconnection events, ascertaining whether the system has achieved steady‐state or is still evolving in time remains challenging. Here, we propose that the relatively rapid opening of the reconnection separatrices within the electron diffusion region serves as an indicator of the growth phase of reconnection. The opening of the separatrices is produced by electron flows diverging away from the neutral line downstream of the X‐line and flowing around a dipolarization front. This flow pattern leads to characteristic spatial structures in the electron pressure‐strain interaction that could be a useful indicator for the growth phase of a reconnection event. We employ two‐dimensional particle‐in‐cell numerical simulations of anti‐parallel magnetic reconnection to validate this prediction. We find that the signature discussed here, alongside traditional reconnection indicators, can serve as a marker of the growth phase. This signature is potentially accessible using multi‐spacecraft single‐point measurements, such as with NASA's Magnetospheric Multiscale satellites in Earth's magnetotail. Applications to other settings where reconnection occurs are also discussed.

Barbhuiya, M. Hasan [Department of Physics and Ast↗

Ion and Electron Acoustic Bursts during Anti-Parallel Reconnection Driven by Lasers

Magnetic reconnection converts magnetic energy into thermal and kinetic energy in plasma. Among the numerous candidate mechanisms, ion acoustic instabilities driven by the relative drift between ions and electrons (or equivalently, electric current) have been suggested to play a critical role in dissipating magnetic energy in collisionless plasmas. However, their existence and effectiveness during reconnection have not been well understood due to ion Landau damping and difficulties in resolving the Debye length scale in the laboratory. We report a sudden onset of ion acoustic bursts measured by collective Thomson scattering in the exhaust of anti-parallel magnetically driven reconnection using high-power lasers. The ion acoustic bursts are followed by electron acoustic bursts with electron heating and bulk acceleration. We reproduce these observations with one- and two-dimensional particle-in-cell simulations in which an electron outflow jet drives ion acoustic instabilities, forming double layers. These layers induce electron two-stream instabilities that generate electron acoustic bursts and energize electrons. Our results demonstrate the importance of ion and electron acoustic dynamics during reconnection when ion Landau damping is ineffective, a condition applicable to a range of astrophysical plasmas including near-Earth space, stellar flares and black hole accretion engines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Inferring thermal ion temperature and residual kinetic energy from nuclear measurements in inertial confinement fusion implosions

In inertial con finement fusion implosion experiments, the presence of residual anisotropic fluid motion within the stagnating hot spot leads to signifi cant variations in ion-temperature measurements using neutron time- of-flight detectors along different lines of sight. The minimum ion-temperature measurement is typically used as representative of the thermal temperature. In the presence of isotropic flows, however, even the minimum DT neutron-inferred ion temperature can be well above the plasma thermal temperature. Using both DD and DT neutron-inferred ion-temperature measurements, we show that it is possible to determine the contribution of isotropic flows and infer the DT burn-averaged thermal ion temperature. The contribution of large isotropic flows on driving the ratio of DD to DT neutron-inferred ion temperatures well below unity and approaching the lower bound of 0.8 is demonstrated in multimode simulations. The minimum DD neutron- inferred ion temperature is determined from the velocity variance analysis, accounting for the presence of isotropic flows. Being close to the DT burn-averaged thermal ion temperature, the inferred DD minimum ion temperatures demonstrate a strong correlation with the experimental yields in the OMEGA implosion database. An analytical expression is also derived to explain the effect of mode ℓ = 1 ion-temperature measurement asymmetry on yield degradations caused by the anisotropic flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-Consistent Magnetosphere-Ionosphere Coupling

A self-consistent ring current (RC) model has been developed that couples electron and ion magnetospheric dynamics with the calculation of the electric field. Two new features were taken into account in order to close the self-consistent magnetosphere-ionosphere coupling loop. First, in addition to the RC ions, we have solved an electron kinetic equation in our model. Second, using the relation of Galand and Richmond, we have calculated the height integrated ionospheric conductances as a function of the precipitated high energy magnetospheric electrons and ions that are produced by our model. To validate the results of our model we simulate the magnetic storm of May 2, 1986, a storm that has been comprehensively studied by Fok et al., and have compared our results with different theoretical approaches. The self-consistent inclusion of the hot electrons and their effect on the conductance results in deeper penetration of the magnetospheric electric field. In addition, a slight westward rotation of the potential pattern (compared to previous self-consistent results) is evident in the inner magnetosphere. These effects change the hot plasma distribution, especially by allowing increased access of plasma sheet ions and electrons to low L shells.

Six, N. Frank↗

Physics of the implosion up until the time of ignition in a revolver (triple-shell) capsule

Here, a revolver is designed to have two dynamical fuel implosion stages to reach ignition: a shock phase in which the fuel is pre-heated up until the shock collapses at the capsule center, followed by an adiabatic compression phase to ignition. This picture assumes ideal hydrodynamics. Employing the state-of-the-art, hybrid (kinetic-ion/fluid electron), multi-ion Vlasov-Fokker-Planck code, iFP, along with semi-analytic predictions from ideal hydrodynamics in spherical geometry, we confirm this two phase picture—resolving some outstanding questions in Revolver design along the way. Although we find that shock kinetic effects and non-ideality are present in the course of the implosion, these effects do not change the overall dynamics (which is well described by ideal hydrodynamics theory). Additionally, we put the assumption of adiabaticity, for the post-collapse phase, to the test. Finally, we verify the ideal hydro-conjecture that the shock collapse (fuel) convergence is ~2.12 given a single shock, and we confirm that the convergence is about 2.5 in the case of two shocks—in agreement with previous radiation-hydrodynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Perturbed Ion Temperature and Toroidal Flow Profile Measurements in Rotating Neoclassical Tearing Mode Magnetic Islands

The perturbed ion temperature and toroidal flow were measured in rotating Neoclassical Tearing Modes (NTM) in a tokamak for the first time. These toroidally and radially resolved profiles were obtained by impurity ion spectroscopy in a 2,1 NTM in DIII-D. In agreement with drift-kinetic simulations, the electron temperature profile is flat, while the ion temperature gradient is restored across the magnetic island O-point in the presence of fast ions; the perturbed flow has minima in the O-points and maxima at the X-points. Furthermore, these measurements provide the first confirmation of the theoretically expected ion temperature and flow response to a magnetic island needed to predict the NTM on set threshold scaling for ITER and other future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of mid- Z gas fill on dynamics and performance of shock-driven implosions at the OMEGA laser

Shock-driven implosions with 100% deuterium (D 2 ) gas fill compared to implosions with 50:50 nitrogen-deuterium (N 2 ⁢D 2 ) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-Ζ fill gas on implosion performance. Ion temperature (Τ ion ) as inferred from the width of measured DD-neutron spectra is seen to be 34%±6% higher for the N 2⁢ D 2 implosions than for the D 2 -only case, while the DD-neutron yield from the D 2 -only implosion is 7.2±0.5 times higher than from the N 2⁢ D 2 gas fill. The T ion enhancement for N 2 ⁢D 2 is observed in spite of the higher Z, which might be expected to lead to higher radiative loss, and higher shock strength for the D 2 -only versus N 2 ⁢D 2 implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-Z material. Further, this picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the N 2 ⁢D 2 implosions relative to the pure D 2 -fill implosion is determined to result from the reduced amount of D 2 in the fill (fourfold effect on yield) combined with a lower fraction of the D 2 fuel being hot enough to burn in the N 2 ⁢D 2 case. The experimental yield and T ion ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot D 2 in the N 2 ⁢D 2 relative to the D 2 -only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D XRAGE simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased T ion and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the N 2 ⁢D 2 versus D 2 -only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

47 OTHER INSTRUMENTATION↗

Dusty-Plasma Particle Accelerator

A dusty-plasma apparatus is being investigated as means of accelerating nanometer- and micrometer-sized particles. Applications for the dusty-plasma particle accelerators fall into two classes: Simulation of a variety of rapidly moving dust particles and micrometeoroids in outer-space environments that include micrometeoroid streams, comet tails, planetary rings, and nebulae and Deposition or implantation of nanoparticles on substrates for diverse industrial purposes that could include hardening, increasing thermal insulation, altering optical properties, and/or increasing permittivities of substrate materials. Relative to prior apparatuses used for similar applications, dusty-plasma particle accelerators offer such potential advantages as smaller size, lower cost, less complexity, and increased particle flux densities. A dusty-plasma particle accelerator exploits the fact that an isolated particle immersed in plasma acquires a net electric charge that depends on the relative mobilities of electrons and ions. Typically, a particle that is immersed in a low-temperature, partially ionized gas, wherein the average kinetic energy of electrons exceeds that of ions, causes the particle to become negatively charged. The particle can then be accelerated by applying an appropriate electric field. A dusty-plasma particle accelerator (see figure) includes a plasma source such as a radio-frequency induction discharge apparatus containing (1) a shallow cup with a biasable electrode to hold the particles to be accelerated and (2) a holder for the substrate on which the particles are to impinge. Depending on the specific design, a pair of electrostatic-acceleration grids between the substrate and discharge plasma can be used to both collimate and further accelerate particles exiting the particle holder. Once exposed to the discharge plasma, the particles in the cup quickly acquire a negative charge. Application of a negative voltage pulse to the biasable electrode results in the initiation of a low-current, high-voltage cathode spot. Plasma pressure associated with the cathode spot as well as the large voltage drop at the cathode spot accelerates the charged particles toward the substrate. The ultimate kinetic energy attained by particles exiting the particle holder depends in part on the magnitude of the cathode spot sheath potential difference, which is proportional to the magnitude of the voltage pulse, and the on the electric charge on the dust. The magnitude of the voltage pulse can be controlled directly, whereas the particle s electric charge can be controlled indirectly by controlling the operating parameters of the plasma apparatus.

Foster, John E.↗

Causal explicit algorithm for heat conduction in a plasma

Hyperbolic heat conduction extends standard Spitzer-Harm heat conduction by including a term proportional to the time derivative of the heat flux. The new term arises from a kinetic derivation of the heat flux that includes higher order corrections. Here we present a causal explicit numerical algorithm for solving the nonlinear hyperbolic heat conduction equation in an unmagnetized plasma. The maximum stable timestep for the causal explicit algorithm scales linearly with the cell size, owing to the hyperbolic nature of the problem. This is in contrast to the quadratic scaling of the maximum stable timestep with the cell size for the parabolic forward time centered space algorithm. The favorable scaling of the timestep with the cell size enables a practical explicit implementation of heat conduction in high-performance massively parallel plasma codes. In particular, we have implemented the causal explicit algorithm in the laser plasma interaction code pF3D. We verify the CE algorithm and analyze its convergence rate by simulating a harmonic mode, which has an analytic solution within the context of the HHC model. We also compare simulations using the CE algorithm to those using the forward time centered space algorithm on a pair of test problems: evolution in time of a Gaussian temperature perturbation in a uniform plasma and heat transport in the presence of inverse bremsstrahlung heating by a Gaussian laser speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Near-Simultaneous POLAR and DMSP Measurements of Topside Ionospheric Up and Down Flows at High Latitudes

We will present what we believe to be the first reported observations of up- and downflows of topside ionospheric thermal plasmas from multiple near-simultaneous tracks through the high-latitude topside ionosphere. From several Southern polar passes, it has been possible to construct plots of field-aligned flows of 0+ observed by the Thermal Ion Dynamics Experiment(TIDE) on the POLAR spacecraft near 5000 km altitude together with vertical ion flow observations from one or more DMSP spacecraft near 800 km altitude. These observations provide a glimpse of the wide-spread upward and downward ionospheric ion flows over the broad polar region along multiple distinct satellite tracks. For the instances where DMSP and POLAR cross nearly the same field lines at the 800 and 5000 kin altitudes, we will show ion parameter measurements at these to altitudes in comparison to Dynamic Fluid-Kinetic(DyFK) transport simulations of anticipated altitude profiles of these parameters.

Zeng, W.↗

Inertial and anisotropic pressure effects on cross-field electron transport in low-temperature magnetized plasmas

In this paper, a one-dimensional (1D) particle-in-cell Monte Carlo collision (PIC-MCC) model is developed to investigate the effects of anisotropic pressure and inertial terms due to non-Maxwellian velocity distribution functions on cross-field electron transport. The conservation of momentum is evaluated by taking the moments of the first-principles gas-kinetic equation. A steady-state discharge is obtained without any low-frequency ionization oscillations by considering an anomalous electron scattering profile. Further, the results obtained from the 1D PIC-MCC model are compared with fluid models, including the quasi-neutral drift-diffusion (DD), non-neutral DD, and full fluid moment models. The discharge current obtained from the PIC-MCC model is in good agreement with the fluid models. The cross-field electron transport due to the inertial terms, i.e. the gradient of axial and azimuthal drift, is evaluated. Moreover, PIC-MCC simulation results show non-zero, anisotropic, off-diagonal pressure tensor terms due to asymmetric non-Maxwellian electron velocity distribution function, potentially contributing to cross-field electron transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Coalescence of two current loops with a kink instability simulated by a three-dimensional electromagnetic particle code

We have studied the dynamics of a coalescence of current loops using three-dimensional electromagnetic (EM) particle simulation code. Our focus is the investigation of such kinetic processes as energy trasnfer, heating particles, and electromagnetic emissions associated with a current loop coalescence which cannot be studied by MHD simulations. First, the two loops undergo a pinching oscillation due to a pressure imbalance between the inside and outside of the current loop. During the pinching oscillation, a kinetic kink instability is excited and electrons in the loops are heated perpendicularly to an ambient magnetic field. Next, the two current loops collide and coalesce, while at the same time a helical structure grows further. Subsequently, the perturbed current, which is due to these helically bunched electrons, can drive a whistler instability. It should be noted in this case that the whistler wave is excited by the kinetic kink instability and not a beam instability. After the coalescence of two helical loops, tilting motions can be observed in the direction of left-hand rotation, and the helical structure will relax resulting in strong plasma heating mostly in the direction perpendicular to the ambient magnetic field. It is also shown that high-frequency electromagnetic waves can be emitted from the region where the two loops coalesce and propagate strongly in the direction of the electron drift velocity. These processes may be important in understanding heating mechansims for coronal loops as well as radio wave emission mechanisms from active regions of solar plasmas.

Nishikawa, K.-I.↗

Comparing XRISM Cluster Velocity Dispersions with Predictions from Cosmological Simulations: Are Feedback Models Too Ejective?

The dynamics of the intracluster medium (ICM), the hot plasma that fills galaxy clusters, are shaped by gravity-driven cluster mergers and feedback from supermassive black holes (SMBHs) in the cluster cores. XRISM measurements of ICM velocities in several clusters offer insights into these processes. We compare XRISM measurements for nine galaxy clusters (Virgo, Perseus, Centaurus, Hydra-A, PKS 0745–19, A2029, Coma, A2319, and Ophiuchus) with predictions from three state-of-the-art cosmological simulation suites, TNG-Cluster, the Three Hundred Project GADGET-X, and GIZMO-SIMBA, that employ different models of feedback. In cool cores, XRISM reveals systematically lower velocity dispersions than the simulations predict, with all 10 measurements below the median simulated values by a factor of 1.5–1.7 on average and all falling within the bottom 10% of the predicted distributions. The observed kinetic-to-total pressure ratio is also lower, with a median value of 2.2%, compared to the predicted 5.0%–6.5% for the three simulations. Outside the cool cores and in non-cool-core (NCC) clusters, simulations show better agreement with XRISM measurements, except for the outskirts of the relaxed, cool-core cluster A2029, which exhibits an exceptionally low kinetic pressure support (<1%), with none of the simulated systems in either of the three suites reaching such low levels. The NCC Coma and A2319 exhibit dispersions at the lower end but within the simulated spread. Our comparison suggests that the three numerical models may overestimate the kinetic effects of SMBH feedback in cluster cores. Additional XRISM observations of NCC clusters will clarify if there is a systematic tension in the gravity-dominated regime as well.

galaxy clusters↗