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 487 records · Page 27

Spatial decay of electromagnetic waves from hypervelocity impact plasmas

Due to the high kinetic energy of impactors in space, hypervelocity impacts on spacecraft can ablate spacecraft material and produce plasmas. These plasmas have been observed in ground-based experiments to produce electromagnetic pulses (EMPs) and are suspected to be a cause of spacecraft anomalies. The initial EMP that is produced after plasma formation is replicated using a 2D electromagnetic plasma simulation, consisting of a discontinuous Galerkin, particle-in-cell (DG-PIC) framework. The spontaneous formation of EMPs from an initial quasineutral plasma is consistent with previous works and show behavior similar to that of waves emitted by antennas. The EMP emits from the plasma with a frequency less than the plasma frequency associated with the point of maximum density, showing the importance of the density gradient, and the far-field region of the “plasma antenna” is estimated. From the proportionality between plasma frequency and EMP frequency, there is a direct correlation between plasma density and radiation. Finally, as the plasma density decreases, the EMP frequency decreases, which suggests that several observable phenomena, such as optical flashes, radio frequency emission, and microwave emission, all result from the same plasma-related mechanism within a single hypervelocity impact.

42 ENGINEERING↗

Simulating ‘Two Ribbon’ Type Solar Flares [Slides]

The underlying mechanism in solar flares in a process known as magnetic reconnection. This astrophysical process is when magnetic field lines, running anti-parallel, break and reconnect. The change in configuration of field lines results in an explosive release of energy as the leftover magnetic energy is converted to kinetic and thermal energies. In this project, we are using an astrophysical magnetohydrodynamic (MHD) simulation code known as Athena++ with a reconnection problem generator file created by Li et al. 2018. Our current research is to successfully implement a radiative cooling term into the MHD equations as it could have important physical effects on the plasma. We use Klimchuk et al. 2008 and SPEX_DM as our chosen cooling functions. The results from the SPEX_DM function show there is a condensation forming that we had predicted but had not seen before with our simulation code. We will continue to analyze the SPEX_DM function and potentially implement thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗

Analytic quasi-steady evolution of marginally unstable kinetic plasma instabilities in the presence of drag and scattering

The 1D bump-on-tail problem is studied in order to determine the influence of drag on quasi-steady solutions near marginal stability (1-γ d /γ L << 1) when effective collisions are much larger than the instability growth rate (ν>>γ). In this common tokamak regime, it is rigorously shown that the paradigmatic Berk–Breizman cubic equation for the nonlinear mode evolution reduces to a much simpler differential equation, dubbed the time-local cubic equation, which can be solved directly. It is found that in addition to increasing the saturation amplitude, drag introduces a shift in the apparent oscillation frequency by modulating the saturated wave envelope. Furthermore, excellent agreement is found between the analytic solution for the mode evolution and both the numerically integrated Berk–Breizman cubic equation and fully nonlinear 1D Vlasov simulations. Experimentally isolating the contribution of drag to the saturated mode amplitude for verification purposes is explored but complicated by the reality that the amount of drag cannot be varied independently of other key parameters in realistic scenarios. While the influence of drag is modest when the ratio of effective drag to effective scattering α/ν is very small, it can become substantial when α/ν ≳ 0.5, suggesting that drag should be accounted for in quantitative models of fast-ion-driven instabilities in fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic Sorbent for the Removal of Selenium(IV) from Simulated Industrial Wastewaters: Determination of Column Kinetic Parameters

A novel meso- and microporous tire-derived-carbon support with magnetic iron oxide nanoparticle adsorbents that selectively adsorbs Se(IV) ions from simulated contaminated water has been developed. In this work, the physicochemical characteristics of the composite adsorbent are characterized with respect to porosity and surface area, chemical composition, and microstructure morphology. The kinetics of this composite adsorbent in a fixed-bed setting has been determined. Several column runs were conducted and analyzed by inductively coupled plasma-optical emission spectroscopy (ICP-OES) to determine the concentration gradient vs time. These results were then fit to a pseudo-second order rate law to obtain equilibrium values. Combining calculated equilibrium values with effluent concentration data, enabled the application of the Adams–Bohart model to determine reaction constants and column coefficients. Column parameters obtained from different flow rates and fittings of the Adams–Bohart model were remarkably consistent. These findings enable the application of this sorbent to fixed-bed column systems and opens up further research into mixed pollutants tests with real wastewater and scaling of selenium pollutant removal.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Viscosity of the magnetized strongly coupled one-component plasma

The viscosity tensor of the magnetized one-component plasma, consisting of five independent shear viscosity coefficients, a bulk viscosity coefficient, and a cross coefficient, is computed using equilibrium molecular dynamics simulations and the Green-Kubo relations. A broad range of Coulomb coupling and magnetization strength conditions are studied. Magnetization is found to strongly influence the shear viscosity coefficients when the gyrofrequency exceeds the Coulomb collision frequency. Three regimes are identified as the Coulomb coupling strength and magnetization strength are varied. Here, the Green-Kubo relations are used to separate kinetic and potential energy contributions to each viscosity coefficient, showing how each contribution depends upon the magnetization strength. The shear viscosity coefficient associated with the component of the pressure tensor parallel to the magnetic field, and the two coefficients associated with the component perpendicular to the magnetic field, are all found to merge to a common value at strong Coulomb coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation and Control of Self-Organized Nonlinear Kinetic Structures in High Energy Density Plasmas in the Presence of Intense Magnetic Fields and Ultrashort Laser Pulses

In the project we investigated mechanisms to sculpt phase space by controlling plasma instabilities. We found that linear theory, which applies in the initial growth phase of an instability, allows for initial conditions that can suppress instabilities. This provides a possible mechanism to release the energy available to an instability on demand. We developed various computational techniques that can reduce computational cost of simulations while maintaining or improving physical fidelity. We considered the Buneman instability in the presence of a ponderomotively driven plasma wave to examine the interplay between an instability and steady-state plasma response.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Role of 3D neoclassical particle flux in density pump-out during ELM control by RMP in DIII-D

Application of three-dimensional (3-D) resonant magnetic perturbations, for the purpose of controlling edge localized modes (ELMs) in high confinement (H-mode) tokamak fusion operations, routinely leads to significantly reduction of the plasma density during the discharge, a phenomenon termed 'density pump-out'. Understanding the density pump-out phenomenon due to 3-D fields is a longstanding challenging issue. This work discusses a self-consistent toroidal model that allows initial value simulation of the density pump-out effect. The model successfully reproduces pump-out phenomenology in DIII-D ELM suppression experiments. The radial particle flux, associated with the 3-D neoclassical effects (the neoclassical toroidal viscosity, NTV), is discovered to play an important role in producing the observed large density pump- out. The key physics is associated with the drift kinetic Landau resonances between the toroidal precession of trapped thermal ions and electrons from one side, and the magnetic field perturbation from the other side. Such resonances, occurring near the top of the pedestal region, substantially enhances the radial particle flux associated with NTV. As a general consequence of the resonant NTV flux, larger density pump- out occurs in plasmas with lower collisionality and slower toroidal rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Omnigenous stellarators with improved ideal and kinetic ballooning stability

Omnigenity is a property of a magnetic field which ensures confinement of trapped particles. It is a necessary requirement for any high-performance stellarator. After creating an omnigenous equilibrium, one must also ensure reduced transport resulting from kinetic and magnetohydrodynamic (MHD) instabilities. To this end, we leverage the GPU-accelerated DESC optimization suite, which is used to design stable, finite-β omnigenous equilibria with poloidal, toroidal, and helical symmetry, achieving Mercier, ideal ballooning, and as a consequence, improved kinetic ballooning stability. We discover stellarators with second stability, a regime of large pressure gradient where an equilibrium becomes ideal ballooning stable, and demonstrate and explore both using theory and gyrokinetic simulations the connection between ideal and kinetic ballooning stability.

optimization↗

Laboratory Demonstration of Collisionless Blob Formation via Laser-Produced Plasma Self-Focusing

Strongly localized, propagating plasma density structures that are capable of crossing magnetic field lines are known as “blobs.” Here we demonstrate a novel mechanism for the formation and propagation of an ion gyroradius-scale blob-cavity structure at the interface between a super-Alfvénic laser-produced plasma (LPP) and an ambient magnetized plasma. The LPP self-focuses along the edge of the diamagnetic cavity which results in a dense, jetlike structure as compared to ballistic motion. This collimated flow couples momentum to the ambient plasma through a collisionless process known as Larmor coupling. The Larmor electric fields locally displace the ambient ions forming a blob above the LPP flow. In the region between a gyrating blob and collimated LPP flow, a secondary cavity of expelled magnetic field forms. In conclusion, these findings are supported by particle-in-cell simulations that replicate the blob formation mechanism and provide insight to similar processes in space, astrophysical, and laboratory settings characterized by ion kinetic scales.

Fluorescence spectroscopy↗

The effects of laser polarization and wavelength on injection dynamics of a laser wakefield accelerator

Here, we investigate the effects of laser polarization and wavelength on electron injection dynamics in a laser wakefield accelerator. During the ionization process, electrons gain residual momentum and kinetic energy via above threshold ionization, which has a strong dependence on laser polarization. A circularly polarized laser pulse results in a much higher residual momentum and kinetic energy gain for the ionized electrons compared with the linearly polarized case. This residual momentum results in particle injection because of the sensitivity of particle trapping to the initial conditions and enhanced the total injected beam charge in both experiments and particle-in-cell simulations. Due to the strong correlation of above threshold ionization with laser wavelength, in this work we extended the investigation to long wavelength (up to 20 μm) drive pulses using particle-in-cell simulations. Owing to the gain in kinetic energy, it may be expected that the charge trapped would consistently increase for circular polarization with increasing laser wavelength, but this was not observed. Instead, there are oscillations with wavelength in the relative trapped charge between linear and circular polarization cases, which arise because of ionization and heating effects on the plasma. Our studies highlight the complex interplay between several different physical effects, including injection regimes—above threshold ionization assisted injection, wave-breaking injection by carrier-envelope-phase effects and ionization injection—ionization gradient induced laser pulse evolution, and thermal modifications to the wake structure that need considering when extrapolating laser wakefield acceleration to different wavelength regimes.

Ma, Y. (ORCID:0000000239194856)↗

Interplay between the non-resonant streaming instability and self-generated pressure anisotropies

The non-thermal particles escaping from collision-less shocks into the surrounding medium can trigger a non-resonant streaming instability that converts parts of their drift kinetic energy into large amplitude magnetic field perturbations, and promote the confinement and acceleration of high energy cosmic rays. We present simulations of the instability using an hybrid-Particle-in-Cell approach including Monte Carlo collisions, and demonstrate that the development of the non-resonant mode is associated with important ion pressure anisotropies in the background plasma. Depending on the initial conditions, the anisotropies may act on the instability by lowering its growth and trigger secondary micro-instabilities. Introducing collisions with neutrals yield a strong reduction of the magnetic field amplification as predicted by linear fluid theory. In contrast, Coulomb collisions in fully ionized plasmas are found to mitigate the self-generated pressure anisotropies and promote the growth of the magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗

Phase-space Energization of Ions in Oblique Shocks

Abstract Examining energization of kinetic plasmas in phase space is a growing topic of interest, owing to the wealth of data in phase space compared to traditional bulk energization diagnostics. Via the field-particle correlation (FPC) technique and using multiple means of numerically integrating the plasma kinetic equation, we have studied the energization of ions in phase space within oblique collisionless shocks. The perspective afforded to us with this analysis in phase space allows us to characterize distinct populations of energized ions. In particular, we focus on ions that reflect multiple times off the shock front through shock-drift acceleration, and how to distinguish these different reflected populations in phase space using the FPC technique. We further extend our analysis to simulations of three-dimensional shocks undergoing more complicated dynamics, such as shock ripple, to demonstrate the ability to recover the phase-space signatures of this energization process in a more general system. This work thus extends previous applications of the FPC technique to more realistic collisionless shock environments, providing stronger evidence of the technique’s utility for simulation, laboratory, and spacecraft analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Dispersion relation for gauge-free electromagnetic drift kinetics

Recently, a new approach to gyrokinetics, invariant under electromagnetic gauge transformations, was developed. The gyrocenter equations of motion are now expressed in terms of the perturbed fields instead of the potentials, in a form suitable for numerical simulations and analytic studies. In this paper, we verify that the long-wavelength limit, i.e., the drift-kinetic limit of the new gyrokinetic theory, is in line with existing work, providing a solid foundation for simulations. Here, we compute the dispersion relation of the new drift-kinetic theory in slab geometry and find agreement with a long-wavelength limit of the full Vlasov–Maxwell model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modelling of sawtooth-induced fast ion transport in positive and negative triangularity in TCV

Internal kinks are a common magneto hydro-dynamic (MHD) instability observed in tokamak operation when the q profile in the plasma core is close to unity. This MHD instability impacts both the transport of the bulk plasma (current, particle and energy transport) and minority species, such as fast ions. In TCV (R 0 /a = 0.88 m/0.25 m) the fast ion population is generated in the plasma by neutral beam tangential injection of energies up to 28 keV. TCV features 16 active shaping coils permitting a great flexibility in plasma shape, including negative triangularity (δ) configurations that show surprisingly high confinement. This study focuses on the transport of fast ions induced by sawteeth, by comparing two triangularity cases and simulation results with experimental data. Comparison of two equilibria with opposite δ shows that the fast ion drifts are larger for δ < 0. Furthermore, the sawtooth-induced transport in this case is larger than δ > 0 in similar conditions. Additionally, comparison with experimental data confirms the dominance of the modification of thermal kinetic profiles following the sawtooth crash in explaining drops in the neutron rates and FIDA signals. Additional fast ion diffusion, however, improves the interpretation of the experimental data. For δ < 0, the amplitude of the perturbation better representing the experimental data is larger. Finaly, an exploratory study for 50 keV particles (soon available in TCV) shows that the situation does not worsen for such particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of net direct current on the properties of radio frequency sheaths: simulation and cross-code comparison

In order to understand, predict and control ion cyclotron range of frequency (ICRF) interactions with tokamak scrape-off layer plasmas, computational tools which can model radio frequency (RF) sheaths are needed. In particular, models for the effective surface impedance and DC rectified sheath potentials may be coupled with full wave RF simulation codes to predict self-consistent wave fields near surfaces and the resulting power dissipation and plasma-material interactions from ion sputtering. In this study, previous work assuming zero net DC current flow through the sheath is generalized to allow the surface to collect net positive or negative current, as is often observed in experiments. The waveforms, DC potential and RF admittance are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) codes. Cross-code comparisons provide detailed model verification and elucidate the roles of ion and electron kinetics. When the sheath draws negative (positive) DC current, the voltage rectification is reduced (increased) compared with the zero-current case, and the magnitude of both the real and imaginary parts of the admittance are increased (reduced). In conclusion, a previous four-input parametrization of the sheath rectification and admittance properties is generalized to include a fifth parameter describing the DC sheath current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic-ballooning-mode turbulence in low-average-magnetic-shear equilibria

In this work, kinetic-ballooning-mode (KBM) turbulence is studied via gyrokinetic flux-tube simulations in three magnetic equilibria that exhibit small average magnetic shear: the Helically Symmetric eXperiment (HSX), the helical-axis Heliotron-J and a circular tokamak geometry. For HSX, the onset of KBM being the dominant instability at low wavenumber occurs at a critical value of normalized plasma pressure $\beta ^{\rm KBM}_{\rm crit}$ that is an order of magnitude smaller than the magnetohydrodynamic (MHD) ballooning limit $\beta ^{\rm MHD}_{\rm crit}$ when a strong ion temperature gradient (ITG) is present. However, $\beta ^{\rm KBM}_{\rm crit}$ increases and approaches the MHD ballooning limit as the ITG tends to zero. For these configurations, $\beta ^{\rm KBM}_{\rm crit}$ also increases as the magnitude of the average magnetic shear increases, regardless of the sign of the normalized magnetic shear. Simulations of Heliotron-J and a circular axisymmetric geometry display behaviour similar to HSX with respect to $\beta ^{\rm KBM}_{\rm crit}$ . Despite large KBM growth rates at long wavelengths in HSX, saturation of KBM turbulence with $\beta > \beta _{\rm crit}^{\rm KBM}$ is achievable in HSX and results in lower heat transport relative to the electrostatic limit by a factor of roughly five. Nonlinear simulations also show that KBM transport dominates the dynamics when KBMs are destabilized linearly, even if KBM growth rates are subdominant to ITG growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global gyrokinetic study of shaping effects on electromagnetic modes at NSTX aspect ratio with ad hoc parallel magnetic perturbation effects

Plasma shaping may have a stronger effect on global turbulence in tight-aspect-ratio tokamaks than in conventional-aspect-ratio tokamaks due to the higher toroidicity and more acute poloidal asymmetry in the magnetic field. In addition, previous local gyrokinetic studies have shown that it is necessary to include parallel magnetic field perturbations in order to accurately compute growth rates of electromagnetic modes in tight-aspect-ratio tokamaks. In this work, the effects of elongation and triangularity on global, ion-scale, linear electromagnetic modes are studied at National Spherical Torus Experiment (NSTX) aspect ratio and high plasma β using the global gyrokinetic particle-in-cell code XGC. The effects of compressional magnetic perturbations are approximated via a well-known modification to the particle drifts that was developed for flux-tube simulations [Joiner et al., Phys. Plasmas 17, 072104 (2010)], without proof of its validity in a global simulation, with the gyrokinetic codes GENE and GEM being used for local verification and global cross-verification. Magnetic equilibria are re-constructed for each distinct plasma profile that is used. Coulomb collision effects are not considered. Within the limitations imposed by the present study, it is found that linear growth rates of electromagnetic modes (collisionless microtearing modes and kinetic ballooning modes) are significantly reduced in a high-elongation and high-triangularity NSTX-like geometry compared to a circular NSTX-like geometry. For example, growth rates of kinetic ballooning modes at high- β are reduced to the level of that of collisionless trapped electron modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗