A novel machine learning based lumping approach for the reduction of large kinetic mechanisms for plasma-assisted combustion applications
Not Available
Rekkas-Ventiris, Georgios (ORCID:0000000226586586)↗
SEARCH · Search NASA
Results for “PLASMA 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.
Scaling of Magnetic Reconnection in Relativistic Collisionless Pair Plasmas
Using fully kinetic simulations, we study the scaling of the inflow speed of collisionless magnetic reconnection in electron-positron plasmas from the non-relativistic to ultra-relativistic limit. In the anti-parallel configuration, the inflow speed increases with the upstream magnetization parameter sigma and approaches the speed of light when sigma is greater than O(100), leading to an enhanced reconnection rate. In all regimes, the divergence of the pressure tensor is the dominant term responsible for breaking the frozen-in condition at the x-line. The observed scaling agrees well with a simple model that accounts for the Lorentz contraction of the plasma passing through the diffusion region. The results demonstrate that the aspect ratio of the diffusion region, modified by the compression factor of proper density, remains approximately 0.1 in both the non-relativistic and relativistic limits.
plasma↗
Coupling Microstructural Evolution Simulations to Material Property Degradation Predictions for Plasma-Facing Materials
Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500 C to 1500 C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 D/m^2-s. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.
36 - MATERIALS SCIENCE↗
Effect of Inductive Coil Geometry and Current Sheet Trajectory of a Conical Theta Pinch Pulsed Inductive Plasma Accelerator
Results are presented demonstrating the e ect of inductive coil geometry and current sheet trajectory on the exhaust velocity of propellant in conical theta pinch pulsed induc- tive plasma accelerators. The electromagnetic coupling between the inductive coil of the accelerator and a plasma current sheet is simulated, substituting a conical copper frustum for the plasma. The variation of system inductance as a function of plasma position is obtained by displacing the simulated current sheet from the coil while measuring the total inductance of the coil. Four coils of differing geometries were employed, and the total inductance of each coil was measured as a function of the axial displacement of two sep- arate copper frusta both having the same cone angle and length as the coil but with one compressed to a smaller size relative to the coil. The measured relationship between total coil inductance and current sheet position closes a dynamical circuit model that is used to calculate the resulting current sheet velocity for various coil and current sheet con gura- tions. The results of this model, which neglects the pinching contribution to thrust, radial propellant con nement, and plume divergence, indicate that in a conical theta pinch ge- ometry current sheet pinching is detrimental to thruster performance, reducing the kinetic energy of the exhausting propellant by up to 50% (at the upper bound for the parameter range of the study). The decrease in exhaust velocity was larger for coils and simulated current sheets of smaller half cone angles. An upper bound for the pinching contribution to thrust is estimated for typical operating parameters. Measurements of coil inductance for three di erent current sheet pinching conditions are used to estimate the magnetic pressure as a function of current sheet radial compression. The gas-dynamic contribution to axial acceleration is also estimated and shown to not compensate for the decrease in axial electromagnetic acceleration that accompanies the radial compression of the plasma in conical theta pinches.
Hallock, Ashley K.↗
Machine learning aided line intensity ratio method for helium–hydrogen mixed recombining plasmas
The helium line intensity ratio (LIR) with the help of a collisional radiative (CR) model has long been used to measure the electron density, n e , and temperature, T e , and its potential and limitations for fusion applications have been discussed. However, it has been reported that the CR model approach leads to deviations in helium–hydrogen mixed plasmas and/or recombining plasmas. In this study, a machine learning (ML) aided LIR method is used to measure n e and T e from spectroscopic data of helium–hydrogen mixed recombining plasmas in the divertor simulator Magnum-PSI. To analyze mixed plasmas, which have more complex spectral shapes, the spectroscopy data were used directly for training instead of separating the intensities of each line. Finally, it is shown that the ML approach can provide a robust and simpler analysis method to deduce n e and T e from the visible emissions in helium–hydrogen mixed plasmas.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Temperature effects in plasma-based positron acceleration schemes using electron filaments
Preserving the quality of a positron beam in a plasma-based accelerator, where a wakefield suitable for positron transport and acceleration is generated by means of an electron filament, is challenging. This is due to the nature of the wakefields, characterized by focusing fields that vary nonlinearly in the transverse direction, and by accelerating fields that are non-uniform. These fields also change slice-by-slice along the beam. Maintaining a high beam quality is pivotal for application of positron beams in a plasma-based collider. In this paper, we show that an initial background plasma temperature can help mitigate the positron beam quality degradation in plasma-based accelerators that rely on electron filaments. We show that temperature effects broaden the electron filament and smooth radially both the non-linear transverse and the non-uniform longitudinal wakefields. Using warm plasmas opens up new possibilities to improve beam quality in several plasma-based positron acceleration concepts.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Axial motion of collector plasma in a relativistic backward wave oscillator
In this paper, it is proposed that plasma formed at the collector may drift back to the cathode and cause pulse shortening of the relativistic backward wave oscillator. Theoretical analysis shows that the axial drift velocity of plasma ions can be up to 5 mm/ns due to the presence of space charge potential provided by an intense relativistic electron beam. Particle-in-cell simulations indicate that the plasma electrons are initially trapped around the collector surface. With the accumulation of the plasma ions, a large electrostatic field forms and drives the plasma electrons to overcome the space charge potential and enter the beam-wave interaction region along the magnetic field lines. As a result, the beam current modulation is disturbed and the output microwave power falls rapidly. The plasma ions move in the beam-wave interaction region with an average axial velocity of 5–8 mm/ns. After the plasma ions reach the diode region, the emitted current at the cathode rises due to the charge neutralizations by the ions. The impedance collapse leads to further decrease of the microwave power. In experiments, when the diode voltage and beam current were 850 kV and 9.2 kA, and the collector radius was 2.15 cm, the output microwave power was 2.4 GW with a pulse width of less than 20 ns. The ion drift velocity was estimated to be about 5 mm/ns. After an improved collector with 3.35 cm radius was adopted, the pulse width was prolonged to more than 30 ns.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Advances in 3D transient plasma dynamics and control through MHD and hybrid fluid-kinetic simulations with JOREK
Transient phenomena and their control are of high relevance in magnetic confinement fusion plasmas to guarantee a stable and safe plasma operation. Interpretative simulations can maximize the insights gained from experiments on present machines and predictive simulations can help in the preparation of design, mitigation techniques and operational scenarios for future devices. In this article, we provide an overview of recent advances and novel scientific results obtained with the 3D non-linear hybrid fluid-kinetic code JOREK, covering physics of plasma transients from the core to the scrape-off layer (SOL) both for tokamak and stellarator devices. Substantial progress was made in the physics understanding, model validation with experiments and experiment interpretation, thus, giving confidence for predictions to devices like DTT, ITER and DEMO. The topics addressed comprise a wide range: the edge physics of new operation scenarios and edge localized mode suppression; major disruptions with a focus on runaway electrons and vertical displacement events as well as disruption mitigation by shattered pellet injection; the physics mechanisms and operational limits of the flux pumping regime for sawtooth control; MHD limits of stellarators and work towards incorporating advanced edge/SOL/exhaust dynamics; continuing improvements of the code for more efficient hybrid simulations on conventional and accelerated high performance computing architectures.
disruptions↗
CHEQUP v0.1
CHEQUP (Castro-based Hofi Expansion with QUasineutral Plasma) is a simulation code for modeling the formation of hydrodynamic optical-field-ionized (HOFI) plasma channels, which are used as waveguides in laser-plasma acceleration experiments. This includes experiments performed at LBNL's BELLA facility as well as other laser facilities across the world. CHEQUP extends the open-source Castro hydrodynamics framework with physics modules tailored for modeling HOFI plasma channels -- including multi-species ionization and three-body recombination for mixtures of hydrogen, nitrogen, helium, and argon ; a two-temperature model tracking electron and heavy-species temperatures separately ; and coupling with other codes of the BLAST ecosystem (https://blast.lbl.gov/) such as WarpX, via the openPMD standard. CHEQUP inherits from Castro the ability to run on modern GPU architectures (NVIDIA CUDA, AMD HIP) and supports adaptive mesh refinement (AMR) for efficient multi-scale resolution. Compared to existing tools, CHEQUP would be, to our knowledge, the first open-source code implementing the full HOFI channel formation physics, and the first implementation capable of running on GPUs. This enables significantly faster, large-scale parameter scans critical for the design of next-generation LPA-based accelerators and light sources.
Lehe, Remi [Lawrence Berkeley National Laboratory ↗
IMAGE EUV Observations and Modeling of the Plasmaspheric Density Trough Associated with the 24 May 2000 Geomagnetic Storm
The IMAGE EUV imager observed a plasmaspheric density trough in association with a geomagnetically active period on 24 May 2000. At EUV wavelengths, this density trough appeared as an Archimedes spiral extending from Earth's shadow to approximately 1800 MLT. We present an analysis of this density trough using simulated EUV images. Observational EUV images are subjected to edge analysis to establish the plasmapause L-shell and the location of the density trough in terms of L-shell, MLT extent, and radial width. The plasmaspheric density distribution is modeled using both static and dynamic models for the plasmasphere. The background plasmasphere is then numerically simulated using the 4-parameter plasmaspheric density model contained within the Global Core Plasma Model (GCPM) [Gallagher et al., 20001 and the Dynamic Global Core Plasma Model (DGCPM). Simulated EUV images of the model plasmasphere are produced once an artificial density depletion, matching the observed MLT extent and width, has been removed. Once the azimuthal extent and width of the trough have been simulated, the depth of the artificial density depletion is iteratively adjusted to produce simulated EUV images that approximate observation. The results of this analysis and discussion of possible origins for this density trough will be presented.
Adrian, M.L.↗
ETHOS: An automated framework to generate multi-fidelity constitutive data tables and propagate uncertainties to hydrodynamic simulations
Accurate constitutive data, such as equations of state and plasma transport coefficients, are necessary for reliable hydrodynamic simulations of plasma systems such as fusion targets, planets, and stars. Here, we develop a framework for automatically generating transport-coefficient tables using a parameterized model that incorporates data from both high-fidelity sources (e.g., density functional theory calculations and reference experiments) and lower-fidelity sources (e.g., average-atom and analytic models). The framework incorporates uncertainties from these multi-fidelity sources, generating ensembles of optimally diverse tables that are suitable for uncertainty quantification of hydrodynamic simulations. We illustrate the utility of the framework with magnetohydrodynamic simulations of magnetically launched flyer plates, which are used to measure material properties in pulsed-power experiments. We explore how changes in the uncertainties assigned to the multi-fidelity data sources propagate to changes in simulation outputs and find that our simulations are most sensitive to uncertainties near the melting transition. The presented framework enables computationally efficient uncertainty quantification that readily incorporates new high-fidelity measurements or calculations and identifies plasma regimes where additional data will have high impact.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Data on erosion and hydrogen fuel retention in Beryllium plasma-facing materials
ITER will use beryllium as a plasma-facing material in the main chamber, covering a total surface area of about 620 m 2 . Given the importance of beryllium erosion and co-deposition for tritium retention in ITER, significant efforts have been made to understand the behaviour of beryllium under fusion-relevant conditions with high particle and heat loads. This paper provides a comprehensive report on the state of knowledge of beryllium behaviour under fusion-relevant conditions: the erosion mechanisms and their consequences, beryllium migration in JET, fuel retention and dust generation. The paper reviews basic laboratory studies, advanced computer simulations and experience from laboratory plasma experiments in linear simulators of plasma–wall interactions and in controlled fusion devices using beryllium plasma-facing components. A critical assessment of analytical methods and simulation codes used in beryllium studies is given. The overall objective is to review the existing set of data with a broad literature survey and to identify gaps and research needs to broaden the database for ITER.
Erosion–deposition↗
Interpretation of electromagnetic modes in the sub-TAE frequency range in JET plasmas with elevated monotonic q-profiles
Recent JET deuterium experiments with an advanced tokamak scenario using an internal transport barrier (ITB) have been performed to clearly observe destabilized toroidicity-induced Alfvén eigenmodes (TAEs) by fast ions; interestingly, these also exhibit unstable electromagnetic (EM) perturbations in the sub-TAE frequency range. We identify such EM perturbations to be beta-induced ion temperature gradient (BTG) eigenmodes and not beta-induced Alfvén eigenmodes nor beta-induced Alfvén acoustic eigenmodes, which are often unstable in such high-beta plasmas with high-power neutral beam injection (NBI). The BTG modes are the most unstable modes due to the high thermal ion temperature gradient related to the ITB, high thermal ion temperature compared to thermal electron temperature (high T i /T e ), and a high ion beta regime. BTG mode experimental characteristics match the analytical theory, i.e., location in the vicinity of a rational magnetic surface with a low magnetic shear, mode frequency scaling with the ion diamagnetic frequency (ω$^*_i$), and a coupling among Alfvén and drift waves. We also perform linear gyrokinetic simulations with validated plasma profiles and equilibrium, and find a mode kinetically driven by thermal ions with similar characteristics as the experimental BTG modes.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
System-size Convergence of Nonthermal Particle Acceleration in Relativistic Plasma Turbulence
We apply collisionless particle-in-cell simulations of relativistic pair plasmas to explore whether driven turbulence is a viable high-energy astrophysical particle accelerator. We characterize nonthermal particle distributions for varying system sizes up to L/2πρ e0 = 163, where L/2π is the driving scale and ρ e0 is the initial characteristic Larmor radius. We show that turbulent particle acceleration produces power-law energy distributions that, when compared at a fixed number of large-scale dynamical times, slowly steepen with increasing system size. We demonstrate, however, that convergence is obtained by comparing the distributions at different times that increase with system size (approximately logarithmically). We suggest that the system-size dependence arises from the time required for particles to reach the highest accessible energies via Fermi acceleration. The converged power-law index of the energy distribution, α ≈ 3.0 for magnetization σ = 3/8, makes turbulence a possible explanation for nonthermal spectra observed in systems such as the Crab Nebula.
Zhdankin, Vladimir↗
Temperature-dependent Saturation of Weibel-type Instabilities in Counter-streaming Plasmas
We present the first 2X2V continuum Vlasov–Maxwell simulations of interpenetrating, unmagnetized plasmas to study the competition between two-stream, Oblique, and filamentation modes in the weakly relativistic regime. We find that after nonlinear saturation of the fastest-growing two-stream and Oblique modes, the effective temperature anisotropy, which drives current filament formation via the secular Weibel instability, has a strong dependence on the internal temperature of the counter-streaming plasmas. The effective temperature anisotropy is significantly more reduced in colder than in hotter plasmas, leading to orders of magnitude lower magnetization for colder plasmas. Here, a strong dependence of the energy conversion efficiency of Weibel-type instabilities on internal beam temperature has implications for determining their contribution to the observed magnetization of many astrophysical and laboratory plasmas.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Fluid simulation of relativistic electron beam driven wakefield in a cold plasma
Excitation of wakefield in a cold homogeneous plasma, driven by an ultra-relativistic electron beam is studied in one dimension using fluid simulation techniques. For a homogeneous rigid beam having density (n{sub b}) less than or equal to half the plasma density (n{sub 0}), simulation results are found to be in good agreement with the analytical work of Rosenzweig [Phys. Rev. Lett. 58, 555 (1987)]. Here, Rosenzweig's work has been analytically extended to regimes where the ratio of beam density to plasma density is greater than half and results have been verified using simulation. Further in contrast to Rosenzweig's work, if the beam is allowed to evolve in a self-consistent manner, several interesting features are observed in simulation viz. splitting of the beam into beam-lets (for l{sub b} > λ{sub p}) and compression of the beam (for l{sub b} < λ{sub p}), l{sub b} and λ{sub p}, respectively, being the initial beam length and plasma wavelength.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
An MHD simulation study of the poloidal mode field line resonance in the Earth's dipole magnetosphere
The poloidal mode field line resonance in the Earth's dipole magnetic field is investigated using cold plasma ideal MHD simulations in dipole geometry. In order to excite the poloidal mode resonance, we use either an initial or a continuous velocity perturbation to drive the system. The perturbation is localized at magnetic shell L = 7 with plasma flow in the radial direction (electric field component in the azimuthal direction). It is found that with the initial perturbation alone, no polodial mode resonance can be obtained and the initially localized perturbation spreads out across all magnetic L shells. With the continuous perturbation, oscillating near the poloidal resonance frequency, a global-scale poloidal cavity mode can be obtained. For the first time, a localized guided poloidal mode resonance is obtained when a radial component of electric field is added to the initial perturbation such that the curl of the electric field is everywhere perpendicular to the background dipole magnetic field. During the localized poloidal resonance, plasma vortices parallel/antiparallel to the background dipole magnetic field B(sub 0). This circular flow, elongated radially, results in twisting of magnetic field flux tubes, which, in turn, leads to the slowdown of the circular plasma flow and reversal of the plasma vortices. The energy associated with the localized poloidal resonance is conserved as it shifts back and forth between the oscillating plasma vortices and the alternately twisted magnetic flux tubes. In the simulations the eigenfunctions associated with the localized poloidal resonance are grid-scale singular functions. This result indicates that ideal MHD is inadequate to describe the underlying problem and nonideal MHD effects are needed for mode broadening.
Ding, D. Q.↗
Time-dependent modeling of dust outburst into tokamak divertor plasma
The first self-consistent simulation studies of the impact of transient injection of tungsten dust on ITER-like edge plasma are presented. The simulations are performed using DUSTT and UEDGE codes coupled in a fully time-dependent manner. The evolution of the divertor and upstream plasma conditions, following the injection of tungsten dust of various quantities and sizes at the outer strike point, is investigated. The estimates of the core plasma impurity fraction and the core impurity screening factor associated with the dust injection are also obtained. Here, the simulations show that the bursts of even modest amounts of tungsten dust injected into the divertor plasmas can result in large perturbations of the divertor power load and unacceptable levels of the core plasma contamination with tungsten impurities.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗