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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The role of focusing geometry in MeV x-ray production from petawatt laser–solid interaction

Relativistic laser–plasma interactions provide a compact and flexible route to generating bright, ultrashort pulses of MeV x rays, with applications in high-energy density science, nuclear physics, and radiography. Despite extensive study of intensity scaling in laser-driven electron acceleration, the role of focusing geometry and focal-volume effects in MeV radiation production remains insufficiently understood. Here, we present experimental and kinetic simulation results from the Texas Petawatt Laser (120 J, 140 fs) in which the focusing geometry (f/3 or f/1.5) is varied, while the laser energy and pulse duration are fixed. Experimentally, the f/3 geometry produces approximately three times more MeV radiation than the f/1.5 geometry, despite its twofold lower nominal vacuum intensity. This result is based on the time-integrated radiation per unit solid angle, as measured along the diagnostic line of sight. Three-dimensional particle-in-cell simulations reproduce this trend when a modest (10s of μm) effective focal plane shift is introduced, demonstrating that relativistic laser–plasma coupling is highly sensitive to focal geometry in the presence of a preplasma. This behavior is consistent with differences in interaction length and effective intensity at the critical surface between the two configurations. The sensitivity of the tightly focused f/1.5 configuration ($z_R ≈$⁠ 5 μ m) reflects the combined influence of thermal lensing and an extended preplasma, while the f/3 geometry (⁠$z_R ≈$ 80 μ m) remains comparatively robust. These results demonstrate a breakdown of conventional intensity scaling and identify focusing geometry as a critical control parameter for MeV electron and x-ray generation at petawatt powers. Implications for laser–plasma-based radiographic facilities are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Role of Parallel Solenoidal Electric Field on Energy Conversion in 2.5D Decaying Turbulence with a Guide Magnetic Field

We perform 2.5D particle-in-cell simulations of decaying turbulence in the presence of a guide (out-of-plane) background magnetic field. The fluctuating magnetic field initially consists of Fourier modes at low wavenumbers (long wavelengths). With time, the electromagnetic energy is converted to plasma kinetic energy (bulk flow+thermal energy) at the rate per unit volume of J· E for current density J and electric field E. Such decaying turbulence is well known to evolve toward a state with strongly intermittent plasma current. Here we decompose the electric field into components that are irrotational, E ir , and solenoidal (divergence-free), E so . E ir is associated with charge separation, and J · E ir is a rate of energy transfer between ions and electrons with little net change in plasma kinetic energy. Therefore, the net rate of conversion of electromagnetic energy to plasma kinetic energy is strongly dominated by J · E so , and for a strong guide magnetic field, this mainly involves the component E so,∥ parallel to the total magnetic field B. We examine various indicators of the spatial distribution of the energy transfer rate J ∥ · E so,∥ , which relates to magnetic reconnection, the best of which are (1) the ratio of the out-of-plane electric field to the in-plane magnetic field, (2) the out-of-plane component of the nonideal electric field, and (3) the magnitude of the estimate of current helicity.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

The interaction of collisionless shocks in astrophysical plasmas

The interaction between collisionless shocks, such as may arise in a number of situations in space plasmas, is examined by means of hybrid numerical simulations. The production of energetic particles by shock collisions is investigated. When the shocks are quasi-perpendicular, ions with energies greater than 10E(0) are produced, where E(0) is the kinetic energy of the unshocked ions. The effect is optimized when the two colliding shocks have roughly equal strengths. As the shocks become closer to quasi-parallel, the number and energy of the accelerated ions increases. Energies in excess of 30E(0) are obtained in the quasi-parallel regime, with a small fall-off when the shocks are exactly parallel. For quasi-parallel collisions, the ion energization can be attributed to the interaction of hot ions with both strong electromagnetic waves and the motional electric field located at the shocks. These results are applied to shock collisions in the distant solar wind, at planetary bow shocks, as well as to the impulsive phase of solar flares.

Cargill, P. J.↗
NASA NTRS2004

Magnetic-Nozzle Studies for Fusion Propulsion Applications: Gigawatt Plasma Source Operation and Magnetic Nozzle Analysis

This project has been a multiyear effort to assess the feasibility of a key process inherent to virtually all fusion propulsion concepts: the expansion of a fusion-grade plasma through a diverging magnetic field. Current fusion energy research touches on this process only indirectly through studies of plasma divertors designed to remove the fusion products from a reactor. This project was aimed at directly addressing propulsion system issues, without the expense of constructing a fusion reactor. Instead, the program designed, constructed, and operated a facility suitable for simulating fusion reactor grade edge plasmas, and to examine their expansion in an expanding magnetic nozzle. The approach was to create and accelerate a dense (up to l0(exp 20)/m) plasma, stagnate it in a converging magnetic field to convert kinetic energy to thermal energy, and examine the subsequent expansion of the hot (100's eV) plasma in a subsequent magnetic nozzle. Throughout the project, there has been a parallel effort between theoretical and numerical design and modelling of the experiment and the experiment itself. In particular, the MACH2 code was used to design and predict the performance of the magnetoplasmadynamic (MPD) plasma accelerator, and to design and predict the design and expected behavior for the magnetic field coils that could be added later. Progress to date includes the theoretical accelerator design and construction, development of the power and vacuum systems to accommodate the powers and mass flow rates of interest to out research, operation of the accelerator and comparison to theoretical predictions, and computational analysis of future magnetic field coils and the expected performance of an integrated source-nozzle experiment.

Gilland, James H.↗
DOE OSTI2025

Three-dimensional reconstruction of inertial confinement fusion hot-spot plasma from x-ray and nuclear diagnostics on OMEGA

Multidimensional effects degrade the neutron yield and the compressed areal density of laser-direct-drive inertial confinement fusion implosions of layered deuterium–tritium cryogenic targets on the OMEGA Laser System with respect to 1D radiation-hydrodynamic simulation predictions. A comprehensive physics-informed 3D reconstruction effort is under way to infer hot-spot and shell conditions at stagnation from four x-ray and seven neutron detectors distributed around the OMEGA target chamber. Neutron diagnostics, providing measurements of the neutron yield, hot-spot flow velocity, and apparent ion-temperature distribution, are used to infer the mode-1 perturbation at stagnation. The x-ray imagers record the shape of the hot-spot plasma to diagnose mode-1 and mode-2 perturbations. A deep-learning convolutional neural network trained on an extensive set of 3D radiation-hydrodynamic simulations is used to interpret the x-ray and nuclear measurements to infer the 3D profiles of the hot-spot plasma conditions and the amount of laser energy coupled to the hot-spot plasma. A 3D simulation database shows that larger mode-1 asymmetries are correlated with higher hot-spot flow velocities and reduced laser-energy coupling and neutron yield. Three-dimensional hot-spot reconstructions from x-ray measurements indicate that higher amounts of residual kinetic energy are correlated with higher measured hot-spot flow velocities, consistent with 3D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Hybrid simulations of sub-cyclotron compressional and global Alfvén eigenmode stability in spherical tokamaks

A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron-frequency compressional Alfvén eigenmodes (CAEs) and global Alfvén eigenmodes (GAEs) in low-aspect-ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam powers in NSTX experiments, prompting a further examination of the conditions necessary for their excitation. Linear simulations have been performed with the hybrid MHD–kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Additionally, three distinct types of modes were found in the simulations—co-CAEs, cntr-GAEs, and co-GAEs—with differing spectral and stability properties. The simulations revealed that unstable GAEs are more ubiquitous than unstable CAEs, which is consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from the simulations and analytically estimated for relevant sources absent from the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Laser-driven collisionless shock acceleration of ions from near-critical plasmas

Here, this paper overviews experimental and numerical results on acceleration of narrow energy spread ion beams by an electrostatic collisionless shockwave driven by 1 μm (Omega EP) and 10 μm (UCLA Neptune Laboratory) lasers in near critical density CH and He plasmas, respectively. Shock waves in CH targets produced high-energy ~50 MeV protons (ΔE/E of ≤30%) and 314 MeV C 6+ ions (ΔE/E of ≤10%). Observation of acceleration of both protons and carbon ions to similar velocities is consistent with reflection of particles off the moving potential of a shock front. For shocks driven by CO 2 laser in a gas jet, ~30 MeV peak in He ion spectrum was detected. Particle-in-cell simulations indicate that regardless of the target further control over its density profile is needed for optimization of accelerated ion beams in part of energy spread, yield and maximum kinetic energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

An improved ten-moment closure for reconnection and instabilities

The integration of kinetic effects in fluid models is important for global simulations of Earth's magnetosphere. The use of the two-fluid ten-moment model, which includes the pressure tensor for both species, has had some success in simulating Ganymede and Mercury with a simple closure model. Here, we discuss a heat flux closure which accounts for some limitations of the earlier work while remaining computationally tractable. Comparisons with kinetic simulations for magnetic reconnection and lower-hybrid drift instabilities show good agreement with kinetic results and improvements on previous closure models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

SOLPS-ITER validation with TCV L-mode discharges

This work presents a quantitative test of SOLPS-ITER simulations against tokamak a configuration variable (TCV) L-mode experiments. These simulations account for drifts, currents, kinetic neutrals, and carbon impurities providing the most complete edge transport simulations for TCV to date. The comparison is performed on nominally identical discharges carried out to assess the effectiveness of TCV’s divertor baffles in the framework of the European Plasma Exhaust program and employs numerous edge diagnostics providing a detailed code-experiment benchmark for TCV. The simulations show a qualitative consistency, but the quantitative differences remain, which are assessed herein. It is found that, for a given separatrix density, the simulations most notably yield a colder, and denser, divertor state with a higher divertor neutral pressure than measured.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Plasma heating at collisionless shocks due to the kinetic cross-field streaming instability

Heating at collisionless shocks due to the kinetic cross-field streaming instability, which is the finite beta (ratio of plasma to magnetic pressure) extension of the modified two stream instability, is studied. Heating rates are derived from quasi-linear theory and compared with results from particle simulations to show that electron heating relative to ion heating and heating parallel to the magnetic field relative to perpendicular heating for both the electrons and ions increase with beta. The simulations suggest that electron dynamics determine the saturation level of the instability, which is manifested by the formation of a flattop electron distribution parallel to the magnetic field. As a result, both the saturation levels of the fluctuations and the heating rates decrease sharply with beta. Applications of these results to plasma heating in simulations of shocks and the earth's bow shock are described.

Winske, D.↗
DOE OSTI2023

Improved melt model for power flow

Accelerators that drive z-pinch experiments transport current densities in excess of 1 MA/cm 2 in order to melt or ionize the target and implode it on axis. These high current densities stress the transmission lines upstream from the target, where rapid electrode heating causes plasma formation, melt, and possibly vaporization. These plasmas negatively impact accelerator efficiency by diverting some portion of the current away from the target, referred to as “current loss”. Simulations that are able to reproduce this behavior may be applied to improving the efficiency of existing accelerators and to designing systems operating at ever higher current densities. The relativistic particle-in-cell code CHICAGO ® is the primary code for modeling power flow on Sandia National Laboratories’ Z accelerator. We report here on new algorithms that incorporate vaporization and melt into the standard power-flow simulation framework. Taking a hybrid approach, the CHICAGO® kinetic/multi-fluid treatment has been expanded to include vaporization while the quasi-neutral equation-of-motion has been updated for melt at high current-densities. For vaporization, a new one-dimensional substrate model provides a more accurate calculation of electrode thermal, mass, and magnetic field diffusion as well as a means of emitting absorbed contaminants and vaporized metal ions. A quasi-fluid model has been implemented expressly to mimic the motion of imploding liners for accurate inductance histories. For melt, a multi-ion Hall-MHD option has been implemented and benchmarked against Alegra MHD. This new model is described with sufficient detail to reproduce these algorithms in any hybrid kinetic code. Physics results from the new code are also presented. A CHICAGO ® Hall-MHD simulation of a radial transmission line demonstrates that Hall physics, not included in Alegra, has no significant impact on the diffusion of electrode material. When surface contaminant desorption is mocked in as a hydrogen surface plasma, both the surface and bulk-material plasmas largely compress under the influence of the j × B force. Similar results are seen in Alegra, which also shows magnetic and material diffusion scaling with peak current. Test vaporization simulations using MagLIF and a power-flow experimental geometry show Fe + ions diffuse only a few hundred µm from the electrodes, so present models of Z power flow remain valid.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

ARC physics basis–magnetohydrodynamics

ARC is designed to produce 400 ⁢MW of net electricity and prove the commercial feasibility of a fusion power plant. In order to achieve this goal ARC has to operate with optimal core performance in a stationary scenario that minimises wear on the first wall and divertor. This requires avoiding or mitigating magnetohydrodynamic (MHD) instabilities which have the potential to not only degrade the plasma core but also lead to deleterious transient heat loads on plasma facing components. Therefore, this work aims at characterising the MHD stability of the high performance ARC scenario and inform the design of error field correction coils. Firstly, simulations of vertical displacement events show that an in-vessel coil is not needed and instead the poloidal shaping coils can be used to control vertical stability. These simulations also inform the demands on the corresponding coil power supplies. Stability analysis of the ideal kink mode with or without a conducting wall and kinetic effects suggests that the ARC baseline scenario operates deeply in the stable region. Using RDCON, tearing modes at the 𝑚/𝑛 =2/1 and 3/2 surfaces (with poloidal mode number 𝑚, and toroidal mode number 𝑛) are shown to be linearly stable, and including thermal transport effects in the rational surfaces lead to further stabilisation. However, other transient plasma instabilities can seed neoclassical tearing modes (NTMs). The marginally stable width of NTMs in ARC strongly depends on the internal inductance and can fall below 0.1% of the normalised poloidal flux. Furthermore, an empirical cross-machine model of the 𝑛 =1 error field leading to a disruption predicts a critical error field larger than SPARC but smaller than ITER. Three-dimensional coils can be designed with the Generalised Purturbed Equilbium Code based on a simple model that calculates the maximum correctable error field that is limited by the neoclassical toroidal viscosity torque. Broad scans of different coil geometries identify a set of 2 rows of off-midplane coils to be a suitable solution. It is also determined that such a set of three-dimensional coils is capable of correcting 𝑛 =2 error fields to some degree and creating strong enough 𝑛 =2 or 𝑛 =3 edge resonant perturbation fields for the suppression of edge-localised modes at reasonable coil currents. The final design of the first ARC will be further informed by results from SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Integrated control of individual plasma scalars with simultaneous neoclassical tearing-mode suppression

A novel integrated-control architecture has been tested in nonlinear, one-dimensional simulations using the control-oriented transport simulator (COTSIM©) and in DIII-D experiments. Integrated architectures that can perform continuous-mission control while also handling off-normal events will be vital in future reactor-grade tokamaks. Continuous-mission controllers for individual magnetic and kinetic scalars (thermal stored-energy (W), volume-average toroidal rotation (Ω Φ ), and safety factor profile (q) at different spatial locations) have been integrated in this work with event-triggered neoclassical tearing-mode (NTM) suppression controllers by combining them into an architecture augmented by a supervisory and exception handling (S&EH) system and an actuator management (AM) system. Here, the AM system, which enables the integration of competing controllers, solves in real time a nonlinear optimization problem that takes into account the high-level control priorities dictated by the S&EH system. The resulting architecture offers a high level of integration and some of the functionalities that will be required to fulfill the advanced-control requirements anticipated for ITER. Initial simulations using COTSIM suggest that the plasma performance and its MHD stability may be improved under integrated feedback control. In addition, the integrated-control architecture has been implemented in the DIII-D plasma control system and tested experimentally for the first time ever in DIII-D in a high-q min scenario, which is a candidate for steady-state operation in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Solitary Magnetic Structures at Quasi–Parallel Collisionless Shocks: Formation

Solitary magnetic structures known as SLAMS (short large-amplitude magnetic structures) have been considered as essential elements of collisionless shocks with quasi-parallel geometries. Yet the physics underlying their formation remains an open question. In this paper, we use measurements from the magnetospheric multiscale mission combined with fully kinetic simulations to study the formation of SLAMS. We find that gyro-resonance between solar wind ions and right-hand circularly polarized electromagnetic waves results in magnetic field amplification. Gyro-trapping by the growing magnetic field builds up the plasma density that further enhances the current and field. The solitary nature of SLAMS stems from a beat-like magnetic field envelope where the maximum sets the initial location for nonlinear growth. Our results present a conceptual advance on SLAMS, and may shed new light on the open question of magnetic field amplification at astrophysical shocks.

58 GEOSCIENCES↗
DOE OSTI2022

Solitary Magnetic Structures Developed From Gyro-Resonance With Solar Wind Ions at Mars and Earth

We investigate solitary magnetic structures in the foreshock of Mars and Earth. The structures exhibit pulse-like magnetic field and density enhancements along with plasma heating and local solar wind (SW) ion reflection. The structures at Mars resemble the foreshock structures developed from ultra-low-frequency electromagnetic waves gyro-resonant with SW ions at Earth, and propagate toward the magnetosphere. We perform fully kinetic simulations to reproduce the solitary structures and the nonlinear evolution of ion distribution functions, illustrating their resonance with SW ions. The structures present self-induced foreshock turbulence that can have space weather effects. Our results advance the fundamental understanding of how SW interacts with planetary magnetospheres, and have potential impact on the current picture of planet-origin ion escape.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Development of magnetohydrodynamic turbulence in coronal loops

It is proposed that a stationary spectrum of magnetohydrodynamic turbulence is generated in solar coronal loops by the interweaving of magnetic field lines driven by the turbulent velocity field in the convective region. Starting from the set of incompressible magnetohydrodynamic equations, two ordinary differential equations describing the stationary turbulent energy spectrum and the frequency of direct energy cascade for a given kinetic energy spectrum of the subphotospheric convective motions are derived. An estimate of the typical dissipation length scale, at which magnetic energy efficiently heats the plasma by Joule dissipation is obtained. The impact of recent results of two-dimensional magnetohydrodynamic numerical simulations on the coronal heating problem is also discussed.

Gomez, Daniel O.↗
NASA NTRS2002

Phase Transition in Dusty Plasmas: A Microphysical Description

Dust grains immersed in plasma discharges acquire a large negative charge and settle into a dust cloud at the edge of the sheath. In this region, the plasma ions stream toward the electrode at a velocity u approx. cs=(T(sub e)/m(sub i))(exp 1/2). Experimentally at sufficiently high gas pressure P, the random kinetic energy of the grains is damped by gas friction, and the grains are strongly coupled and self-organize into a crystalline configuration. For lower pressures despite the dissipation of grain kinetic energy to gas friction, the dust grains reach a steady-state kinetic temperature T(sub d) which is much larger than the temperature of any other component in the plasma. T(sub d) is so large that the dust acts like a fluid. We have used the dynamically shielded dust (DSD) model to simulate these physical processes. We find that the known experimental features are nicely reproduced in the simulations, and that additional features are revealed. In the figure we plot the variation of T(sub d) as P is continuously varied in a DSD code run. A marked difference is evident between the critical pressure P(sub m) for the melting transition as P is decreased, and the critical pressure P(sub c) for the condensation transition as P is increased. For P(sub m) is less than P is less than P(sub c), mixed phase states are seen. This hysteresis occurs because the instability which triggers melting is different from the instability that heats the dust in the fluid phase and inhibits freezing. At low pressure, the dust is subject to a two-stream instability with the ions. This instability is responsible for the high temperature of the dust at low pressure. The basic physics underlying the melting transition has been elucidated in a series of papers. We are developing a first-principles analytic approach to the melting transition, which embodies the same physics that is present in the DSD code.

Joyce, Glenn↗