Unstructured mesh technologies for massively parallel simulation and data analysis of magnetically confined plasmas
Explore the source record for details and available documents.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
SEARCH · Search NASA
Results for “cell fusion”
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.
Radial drift of plasma blobs in a toroidal magnetic field with fully kinetic and reduced fluid models
Abstract In curved magnetic geometries, field-aligned regions of enhanced plasma pressure and density, termed ‘blobs,’ move as coherent filaments across the magnetic field lines. Coherent blobs account for a significant fraction of transport at the edges of magnetic fusion experiments and arise in naturally-occurring space plasmas. This work examines the dynamics of blobs with a fully kinetic electromagnetic particle-in-cell code and with a drift-reduced fluid code. In low-beta regimes with moderate blob speeds, good agreement is found in the maximum blob velocity between the two simulation schemes and simple analytical estimates. The fully kinetic code demonstrates that blob speeds saturate near the initial sound speed, which is a regime outside the validity of the reduced fluid model.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Radial drift of plasma blobs in a toroidal magnetic field with fully kinetic and reduced fluid models
In curved magnetic geometries, field-aligned regions of enhanced plasma pressure and density, termed ‘blobs,’ move as coherent filaments across the magnetic field lines. Coherent blobs account for a significant fraction of transport at the edges of magnetic fusion experiments and arise in naturally-occurring space plasmas. This work examines the dynamics of blobs with a fully kinetic electromagnetic particle-in-cell code and with a drift-reduced fluid code. In low-beta regimes with moderate blob speeds, good agreement is found in the maximum blob velocity between the two simulation schemes and simple analytical estimates. The fully kinetic code demonstrates that blob speeds saturate near the initial sound speed, which is a regime outside the validity of the reduced fluid model.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Integrative Quantitative-Phase and Airy Light-Sheet Imaging
Light-sheet microscopy enables considerable speed and phototoxicity gains, while quantitative-phase imaging confers label-free organelle recognition and metabolic information that are inaccessible by conventional methods. We report the fusion of these two modalities onto a standard inverted microscope that retains compatibility with microfluidics. We describe the utilization of an accelerating Airy-beam light-sheet yielding identical imaging areas with interferometry, and an application in unmasking the effects of cellular noise on metabolic compartmentalization.
Biological sciences, Biological techniques, Micros↗
A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering
Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Modeling of convective cells, turbulence, and transport induced by a radio-frequency antenna in the tokamak boundary plasma
The edge turbulence model Hermes (Dudson et al 2017 Plasma Phys. Control. Fusion 59 05401) is set up for plasma boundary simulations with an radiofrequency (RF) antenna, using parameters characteristic of a tokamak edge. Cartesian slab geometry is used with thin plate limiters representing the ion cyclotron range of frequency (ICRF) antenna side-wall limiters. Ad-hoc DC electric biasing of the limiters, motivated by calculations with VSim (Nieter et al 2004 J. Comput. Phys. 196 448), represents an induced RF sheath rectified potential in the plasma turbulence model. Flux-driven turbulence simulations demonstrate a realistic distribution of plasma profiles and fluctuations. There is a clear effect of the antenna sheath voltage leading to formation of convective cells; bias-induced convective transport flattens the scrape-off layer density profile and fluctuations penetrate into the shadow region of the limiters as the bias voltage increases. Turbulent transport for impurity ions is inferred by following ion trajectories in the simulated plasma turbulence fields, showing Bohm-like effective diffusion rates. All in all, the model elucidates the key physical phenomena governing the effects of ICRF-induced antenna biasing on the tokamak boundary plasma.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Hybrid-VPIC: An open-source kinetic/fluid hybrid particle-in-cell code
Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview of applications of the code to space and laboratory plasma physics problems. Particular choices in how the hybrid solvers were implemented are documented for reference by users. A few solutions for handling numerical complications particular to hybrid codes are also described. Finally, special emphasis is given to the computationally taxing problem of modeling mix in collisional high-energy-density regimes, for which more accurate electron fluid transport coefficients have been implemented for the first time in a hybrid PIC code.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Observation of weakly damped modes using high resolution measurement of turbulence in a dipole confined plasma
When plasma is confined by a levitated superconducting current ring, the turbulent radial particle flux changes the direction depending upon the gradient of the particle number within tubes of equal magnetic flux [Garnier et al., Phys. Plasmas 24, 012506 (2017)]. When the density profile is highly peaked, the particle flux is outward. When the density profile is less peaked, turbulence drives an inward “particle pinch.” This change in the direction of particle flux coincides with a reversal of the mean toroidal propagation direction of turbulent fluctuations. Previously, only the mean wavenumber of the fluctuations was reported. The mean wavenumber was computed from the ensemble correlation between two probes. Here, the full frequency-wavenumber power spectrum of the saturated turbulence is calculated by applying Capon’s “maximum likelihood method” to floating potential measurements from an array of probes. The power spectrum shows that the highest intensity modes are due to unstable entropy modes, and the less intense modes coincide with weakly damped convective cells, which rotate toroidally in the direction opposite to the unstable modes. Because the weakly damped modes make up a significant fraction of the fluctuating power, these measurements imply that saturated turbulence in a dipole-confined plasma involves the excitation of weakly damped rotating convective cells.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Energy transfer in reconnection and turbulence
Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady energy transfer rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. Finally, these results provide explicit evidence that reconnection dynamics fundamentally involves turbulence-like energy transfer.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS
Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Accurate simulation of direct laser acceleration in a laser wakefield accelerator
In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The traditional particle-in-cell (PIC) algorithm, although often the tool of choice to study DLA, contains inherent errors due to numerical dispersion and the time staggering of the electric and magnetic fields. Furthermore, conventional PIC implementations cannot reliably disentangle the fields of the plasma wave and laser pulse, which obscures interpretation of the dominant acceleration mechanism. Here, a customized field solver that reduces errors from both numerical dispersion and time staggering is used in conjunction with a field decomposition into azimuthal modes to perform PIC simulations of DLA in an LWFA. Comparisons with traditional PIC methods, model equations, and experimental data show improved accuracy with the customized solver and convergence with an order-of-magnitude fewer cells. Furthermore, the azimuthal-mode decomposition reveals that the most energetic electrons receive comparable energy from DLA and LWFA.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Nonlocal, diamagnetic electromagnetic effects in magnetically insulated transmission lines
We identify the time-dependent physics responsible for the critical reduction of current losses in magnetically insulated transmission lines (MITLs) due to uninsulated space charge-limited currents of electrons emitted by field stress. A drive current of sufficiently short pulse length introduces a strong enough time dependence that steady-state results alone become inadequate for the complete understanding of current losses. The time-dependent physics can be described as a nonlocal, diamagnetic electromagnetic response of space charge limited currents. As the pulse length is increased or equivalently, the MITL length reduced, these time-dependent effects diminish and current losses converge to those predicted by the well-known Child–Langmuir law in the external (vacuum) fields. We present a simple one-dimensional (1D) model that encapsulates the essence of this physics. We find excellent agreement with 2D particle-in-cell simulations for two MITL geometries, Cartesian parallel plate and azimuthally symmetric straight coaxial. Based on the 1D model, we explore various scaling dependencies of MITL losses with relevant parameters, e.g., peak current, pulse length, geometrical dimensions, etc. We propose an improved physics model of magnetic insulation in the form of a Hull curve, which could also help improve predictions of current losses by common circuit element codes, such as BERTHA. Finally, we describe how to calculate the temperature rise due to electron impact within the 1D model.
Computer simulation↗
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↗
Particle-in-cell simulations of expanding high energy density plasmas with laser ray tracing
The design and analysis of high energy density (HED) laser experiments typically rely on radiation hydrodynamics simulations. However, some laser–plasma interaction regimes are not collisional and cannot be adequately modeled with hydrodynamics. For example, strongly driven magnetic reconnection and magnetized collisionless shock experiments possess extended hydrodynamic or even kinetic properties, necessitating first-principles kinetic simulations. In this paper, we present the benchmarking and first results obtained with a laser-ray-tracing and inverse bremsstrahlung absorption module implemented in the particle-in-cell code PSC. The simulation results are compared to radiation hydrodynamic simulations using the FLASH code as well as analytical estimates. We successfully benchmark the energy deposition model and overall hydrodynamic evolution of the systems. We also consider possible kinetic effects that may be expected from laser-target ablation in the HED regime, including non-local transport and two-temperature effects.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Effects of ion trapping and fluctuations of electron temperature and plasma flow on cross-beam energy transfer
The influences of ion trapping and fluctuations of electron temperature and plasma flow on cross-beam energy transfer (CBET) are examined using two- and three-dimensional particle-in-cell simulations in parameter regimes relevant to recent CBET experiments at the OMEGA laser facility. In mid-Z plasma irradiated by an intense pump beam and weaker probe beam, ion trapping, collisional de-trapping, and plasma flow induced by thermal effects are shown to affect the CBET gain. Ion trapping can enhance or detune the CBET resonance [Nguyen et al., Phys. Plasmas 28, 082705 (2021)]. Collisional de-trapping can affect the CBET gain at low seed beam intensity near the onset threshold for ion trapping. Thermal-effects-induced flow can also detune the CBET resonance at a level comparable to that from trapping at low seed beam intensity. As a consequence, the CBET gain is sensitive to collisions and dimensionality at low seed beam intensity where ion trapping is weak but is insensitive to collisions and dimensionality at high seed beam intensity where ion trapping is strong.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Particle-in-cell simulation of the neutrino fast flavor instability
Not Available
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Liner implosion experiments driven by a dynamic screw pinch
This paper expands upon recent experimental results [Campbell et al., Phys. Rev. Lett. 125, 035001 (2020)], where thin-foil liner implosions were driven by a dynamic screw pinch (DSP) and found to have magneto-Rayleigh–Taylor instability (MRTI) amplitudes up to three times smaller than in implosions driven by a standard z-pinch (SZP). Here, the expanded discussion presented herein includes: (1) a detailed comparison of the MRTI growth measured in the experiment with that calculated from theory; (2) measurements of axial magnetic field injection into the liner interior prior to the implosion, as well as the subsequent compression of this field during the implosion; (3) an in-depth description of how the helical geometry of the DSP can result in earlier implosion and stagnation times relative to the SZP; and (4) particle-in-cell simulations showing different electron drift behavior in the anode–cathode gap of the DSP relative to the SZP, and how this difference may be related to the different current waveforms recorded during the experiments.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Time-dependent saturation and physics-based nonlinear model of cross-beam energy transfer
The nonlinear physics of cross-beam energy transfer (CBET) for multi-speckled laser beams is examined using large-scale particle-in-cell simulations for a range of laser and plasma conditions relevant to indirect-drive inertial confinement fusion (ICF) experiments. The time-dependent growth and saturation of CBET involve complex, nonlinear ion and electron dynamics, including ion trapping-induced enhancement and detuning, ion acoustic wave (IAW) nonlinearity, oblique forward stimulated Raman scattering (FSRS), and backward stimulated Brillouin scattering (BSBS) in a CBET-amplified seed beam. Ion-trapping-induced detuning of CBET is captured in the kinetic linear response by a new δf-Gaussian-mixture algorithm, enabling an accurate characterization of trapping-induced non-Maxwellian distributions. Ion trapping induces nonlinear processes, such as changes to the IAW dispersion and nonlinearities (e.g., bowing and self-focusing), which, together with pump depletion, FSRS, and BSBS, determine the time-dependent nature and level of CBET gain as the system approaches a steady state. Using VPIC simulations at intensities at and above the onset threshold for ion trapping and the insight from the time-dependent saturation analyses, we construct a nonlinear CBET model from local laser and plasma conditions that predicts the CBET gain and the energy deposition into the plasma. This model is intended to provide a more accurate, physics-based description of CBET saturation over a wide range of conditions encountered in ICF hohlraums compared with linear CBET gain models with ad hoc saturation clamps often used in laser ray-based methods in multi-physics codes.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗