Search NASA⌕ Search

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.

At least 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Symplectic neural network and its application to charged particle dynamics in electromagnetic fields

Recently, machine learning models have shown many successes in various applications in science and technology. In this work, we focus on the charged particle dynamics, with the development of a class of symplectic neural networks, including a linear version, SympMat, and a nonlinear version, HénonNet. Both are designed to preserve the structure of Hamiltonian systems. We show that they can be used to model relevant Hamiltonian systems of interest in plasma physics and astrophysics, for linear and nonlinear charged particle dynamics, with the potential to bridge multi-scale simulations. These symplectic neural networks are adapted to the applications in plasma simulations and particle-wave interaction with parametric dependence and periodicity, where we have investigated their performance and accuracy. In particular, SympMat is shown to outperform the traditional Boris particle pusher down to the sub-gyroperiod scale in the case of charged particles in uniform magnetic fields. HénonNet successfully predicts the hot electron distribution, which is validated against theoretical results. These results highlight the potential of symplectic neural networks as a trajectory integrator for particle-in-cell simulations or a fast surrogate to replace conventional numerical schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Feedforward equilibrium trajectory optimization with GSPulse

One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.

equilibrium↗
DOE OSTI2023

Toroidal modeling of plasma flow damping and density pump-out by RMP during ELM mitigation in HL-2A

Reduction of both the plasma density and toroidal flow speed, due to application of the predominantly n = 1 (n is the toroidal mode number) resonant magnetic perturbation (RMP) for controlling the edge localized mode in the HL-2A tokamak, is numerically investigated utilizing the quasi-linear initial-value code MARS-Q (Liu et al 2013 Phys. Plasmas 20 042503). Simulation results reveal that the neoclassical toroidal viscosity (NTV) due to three dimensional fields plays the key role in modifying the plasma momentum and particle transport in the HL-2A discharge. By comparing the modeling results with the measured density pump-out in the experiment, the electron NTV particle flux model, in combination with the free-boundary condition for the axisymmetric change of the density at the plasma edge, is found to yield the best agreement in terms of both the pump-out level and the overall time scale. Further sensitivity studies show that the simulated density pump-out level is reasonably robust against variations in the model assumptions, including the particle diffusion model and the non-ambipolar versus ambipolar NTV particle flux. The latter however affects the time scale for reaching the steady state solution. Finally, it is found that the plasma edge-peeling response, the NTV torque, as well as the plasma momentum and particle transport, all are sensitive to the toroidal phase difference between the upper and lower rows of the RMP coil currents in HL-2A, with the 30 ∘ coil phasing producing the minimal side effects on the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Effects of the plasma profiles on photon and pair production in ultrahigh intensity laser solid interaction

The model of photon and pair production in strong field quantum electrodynamics is implemented into our 1D3V particle-in-cell code with Monte Carlo algorithm. Using this code, the evolution of the particles in ultrahigh intensity laser (∼10{sup 23} W/cm{sup 2}) interaction with aluminum foil target is observed. Four different initial plasma profiles are considered in the simulations. The effects of initial plasma profiles on photon and pair production, energy spectra, and energy evolution are analyzed. The results imply that one can set an optimal initial plasma profile to obtain the desired photon distributions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Milestone in predicting core plasma turbulence: successful multi-channel validation of the gyrokinetic code GENE

On the basis of several recent breakthroughs in fusion research, many activities have been launched around the world to develop fusion power plants on the fastest possible time scale. In this context, high-fidelity simulations of the plasma behavior on large supercomputers provide one of the main pathways to accelerating progress by guiding crucial design decisions. When it comes to determining the energy confinement time of a magnetic confinement fusion device, which is a key quantity of interest, gyrokinetic turbulence simulations are considered the approach of choice – but the question, whether they are really able to reliably predict the plasma behavior is still open. The present study addresses this important issue by means of careful comparisons between state-of-the-art gyrokinetic turbulence simulations with the GENE code and experimental observations in the ASDEX Upgrade tokamak for an unprecedented number of simultaneous plasma observables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Changes in peak oxygen uptake and plasma volume in fit and unfit subjects following exposure to a simulation of microgravity

To test the hypothesis that the magnitude of reduction in plasma volume and work capacity following exposure to simulated microgravity is dependent on the initial level of aerobic fitness, peak oxygen uptake (VO2peak) was measured in a group of physically fit subjects and compared with VO2peak in a group of relatively unfit subjects before and after 10 days of continuous 6 degrees head-down tilt (HDT). Ten fit subjects (40 +/- 2 year) with mean +/- SE VO2peak = 48.9 +/- 1.7 mL kg-1 min-1 were matched for age, height, and lean body weight with 10 unfit subjects (VO2peak = 37.7 +/- 1.6 mL kg-1 min-1). Before and after HDT, plasma, blood, and red cell volumes and body composition were measured and all subjects underwent a graded supine cycle ergometer test to determine VO2peak period needed. Reduced VO2peak in fit subjects (-16.2%) was greater than that of unfit subjects (-6.1%). Similarly, reductions in plasma (-18.3%) and blood volumes (-16.0%) in fit subjects were larger than those of unfit subjects (blood volume = -5.6%; plasma volume = -6.6%). Reduced plasma volume was associated with greater negative body fluid balance during the initial 24 h of HDT in the fit group (912 +/- 154 mL) compared with unfit subjects (453 +/- 200 mL). The percentage change for VO2peak correlated with percentage change in plasma volume (r = +0.79). Following exposure to simulated microgravity, fit subjects demonstrated larger reductions in VO2peak than unfit subjects which was associated with larger reductions in plasma and blood volume. These data suggest that the magnitude of physical deconditioning induced by exposure to microgravity without intervention of countermeasures was influenced by the initial fitness of the subjects.

Non-NASA Center↗
DOE OSTI2023

The Origin of Power-law Spectra in Relativistic Magnetic Reconnection

Abstract Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model—supported and benchmarked with large-scale three-dimensional kinetic particle-in-cell simulations in electron–positron plasmas—that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factorγ≳ 3σ(here,σis the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time is t acc ∼ γ η rec − 1 ω c − 1 ≃ 20 γ ω c − 1 , whereη rec ≃ 0.06 is the inflow speed in units of the speed of light andω c =eB 0 /mcis the gyrofrequency in the upstream magnetic field. They leave the region of active energization aftert esc , when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measuret esc in our simulations and find thatt esc ∼t acc forσ≳ few. This leads to a universal (i.e.,σ-independent) power-law spectrum dN free / d γ ∝ γ − 1 for the particles undergoing active acceleration, and dN / d γ ∝ γ − 2 for the overall particle population. Our results help to shed light on the ubiquitous presence of power-law particle and photon spectra in astrophysical nonthermal sources.

Astronomy & Astrophysics↗
NASA NTRS1987

A possible explanation of the electron temperature enhancement in the wake of a satellite

A numerical simulation was performed in an attempt to explain the enhancement of electron temperature in the wake of an ionospheric satellite in terms of counterstreaming plasma expansions. The simulation was conducted using full kinetic treatments (Vlasov-Poisson code) for both the plasma electrons and ions. The calculations indicated that at a distance of a few Debye lengths from the surface of the body the wake/ambient electron-temperature ratio reaches a value of about 2. This result is in good agreement with some in situ and laboratory experimental results.

Singh, Nagendra↗
DOE OSTI2023

Core-Pedestal Plasma Configurations in Advanced Tokamaks

Here, several configurations for the core and pedestal plasma are examined for a predefined tokamak design by implementing multiple heating/current drive (H/CD) sources to achieve an optimum configuration of high fusion power in a noninductive operation while maintaining an ideally magnetohydrodynamic (MHD) stable core plasma using the IPS-FASTRAN framework. IPS-FASTRAN is a component-based lightweight coupled simulation framework that is used to simulate magnetically confined plasma by integrating a set of high-fidelity codes to construct the plasma equilibrium (EFIT, TOQ, and CHEASE), calculate the turbulent heat and particle transport fluxes (TGLF), model various H/CD systems (TORIC, TORAY, GENRAY, and NUBEAM), model the pedestal pressure and width (EPED), and estimate the ideal MHD stability (DCON). The TGLF core transport model and EPED pedestal model are used to self-consistently predict plasma profiles consistent with ideal MHD stability and H/CD (and bootstrap) current sources. In order to evaluate the achievable and sustainable plasma beta, varying configurations are produced ranging from the no-wall stability to with-wall stability regimes, simultaneously subject to the self-consistent TGLF, EPED, and H/CD source profile predictions that optimize configuration performance. The pedestal density, plasma current, and total injected power are scanned to explore their impact on the target plasma configuration, fusion power, and confinement quality. A set of fully noninductive scenarios are achieved by employing ion-cyclotron, neutral beam injection, helicon, and lower-hybrid H/CDs to provide a broad profile for the total current drive in the core region for a predefined tokamak design. These noninductive scenarios are characterized by high fusion gain (Q ~ 4) and power (P fus ~ 600 MW), optimum confinement quality (H 98 ~ 1.1), and high bootstrap current fraction (f BS ~ 0.7) for Greenwald fraction below unity. The broad current profile configurations identified are stable to low-n kink modes either because the normalized pressure β N is below the no-wall limit or a wall is present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Plasma physical parameters along CME-driven shocks. II. Observation–simulation comparison

In this work, we compare the spatial distribution of the plasma parameters along the 1999 June 11 coronal mass ejection (CME)-driven shock front with the results obtained from a CME-like event simulated with the FLIPMHD3D code, based on the FLIP-MHD particle-in-cell method. The observational data are retrieved from the combination of white-light coronagraphic data (for the upstream values) and the application of the Rankine–Hugoniot equations (for the downstream values). The comparison shows a higher compression ratio X and Alfvénic Mach number M{sub A} at the shock nose, and a stronger magnetic field deflection d toward the flanks, in agreement with observations. Then, we compare the spatial distribution of M{sub A} with the profiles obtained from the solutions of the shock adiabatic equation relating M{sub A}, X, and θ{sub Bn} (the angle between the upstream magnetic field and the shock front normal) for the special cases of parallel and perpendicular shock, and with a semi-empirical expression for a generically oblique shock. The semi-empirical curve approximates the actual values of M{sub A} very well, if the effects of a non-negligible shock thickness δ{sub sh} and plasma-to magnetic pressure ratio β{sub u} are taken into account throughout the computation. Moreover, the simulated shock turns out to be supercritical at the nose and sub-critical at the flanks. Finally, we develop a new one-dimensional Lagrangian ideal MHD method based on the GrAALE code, to simulate the ion-electron temperature decoupling due to the shock transit. Two models are used, a simple solar wind model and a variable-γ model. Both produce results in agreement with observations, the second one being capable of introducing the physics responsible for the additional electron heating due to secondary effects (collisions, Alfvén waves, etc.).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

OBSERVATIONAL SIGNATURES OF CORONAL LOOP HEATING AND COOLING DRIVEN BY FOOTPOINT SHUFFLING

The evolution of a coronal loop is studied by means of numerical simulations of the fully compressible three-dimensional magnetohydrodynamic equations using the HYPERION code. The footpoints of the loop magnetic field are advected by random motions. As a consequence, the magnetic field in the loop is energized and develops turbulent nonlinear dynamics characterized by the continuous formation and dissipation of field-aligned current sheets: energy is deposited at small scales where heating occurs. Dissipation is nonuniformly distributed so that only a fraction of the coronal mass and volume gets heated at any time. Temperature and density are highly structured at scales that, in the solar corona, remain observationally unresolved: the plasma of our simulated loop is multithermal, where highly dynamical hotter and cooler plasma strands are scattered throughout the loop at sub-observational scales. Numerical simulations of coronal loops of 50,000 km length and axial magnetic field intensities ranging from 0.01 to 0.04 T are presented. To connect these simulations to observations, we use the computed number densities and temperatures to synthesize the intensities expected in emission lines typically observed with the Extreme Ultraviolet Imaging Spectrometer on Hinode. These intensities are used to compute differential emission measure distributions using the Monte Carlo Markov Chain code, which are very similar to those derived from observations of solar active regions. We conclude that coronal heating is found to be strongly intermittent in space and time, with only small portions of the coronal loop being heated: in fact, at any given time, most of the corona is cooling down.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Plasma Science and Innovation Center (PSI-Center) at Washington, Wisconsin, and Utah State, ARRA Supplement

The objective of the Plasma Science and Innovation Center (PSI-Center) is to develop and deploy computational models that simulate conditions in smaller, concept-exploration plasma experiments. The PSIC group at the University of Wisconsin-Madison, led by Prof. Carl Sovinec, uses and enhances the Non-Ideal Magnetohydrodynamics with Rotation, Open Discussion (NIMROD) code, to simulate macroscopic plasma dynamics in a number of magnetic confinement configurations. These numerical simulations provide information on how magnetic fields and plasma flows evolve over all three spatial dimensions, which supplements the limited access of diagnostics in plasma experiments. The information gained from simulation helps explain how plasma evolves. It is also used to engineer more effective plasma confinement systems, reducing the need for building many experiments to cover the physical parameter space. The ultimate benefit is a more cost-effective approach to the development of fusion energy for peaceful power production. The supplemental funds provided by the American Recovery and Reinvestment Act of 2009 were used to purchase computer components that were assembled into a 48-core system with 256 Gb of shared memory. The system was engineered and constructed by the group's system administrator at the time, Anthony Hammond. It was successfully used by then graduate student, Dr. John O'Bryan, for computing magnetic relaxation dynamics that occur during experimental tests of non-inductive startup in the Pegasus Toroidal Experiment (pegasus.ep.wisc.edu). Dr. O'Bryan's simulations provided the first detailed explanation of how the driven helical filament of electrical current evolves into a toroidal tokamak-like plasma configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Turbulence in strongly coupled dusty plasmas using generalized hydrodynamic description

The properties of decaying turbulence have been studied with the help of a Generalized Hydrodynamic (GHD) fluid model in the context of strongly coupled dusty plasma medium in two dimensions. The GHD model treats the strongly coupled dusty plasma system as a visco-elastic medium. The incompressible limit of the GHD model is considered here. The studies carried out here are, however, applicable to a wider class of visco-elastic systems, and are not merely confined to the dusty plasma medium. Our simulations studies show that an initial spectrum that is confined in a limited domain of wave numbers becomes broad, even when the Reynold's number is much less than the critical value required for the onset of turbulence in Newtonian fluids. This is a signature of elastic turbulence, where Weissenberg's number also plays an important role on the onset of turbulence. This feature has been observed in several experiments. It is also shown that the existence of memory relaxation time parameter and the transverse shear wave inhibit the normal process (for 2-D systems) of inverse spectral cascade in this case. A detailed simulation study has been carried out for the understanding of this inhibition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Molecular dynamics simulations of soliton-like structures in a dusty plasma medium

The existence and evolution of soliton-like structures in a dusty plasma medium are investigated in a first principles approach using molecular dynamic (MD) simulations of particles interacting via a Yukawa potential. These localized structures are found to exist in both weakly and strongly coupled regimes with their structures becoming sharper as the correlation effects between the dust particles get stronger. A surprising result, compared to fluid simulations, is the existence of rarefactive soliton-like structures in our non-dissipative system, a feature that arises from the charge conjugation symmetry property of the Yukawa fluid. Our simulation findings closely resemble many diverse experimental results reported in the past.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Location of the first plasma response to resonant magnetic perturbations in DIII-D H-mode plasmas

The resonant location of the first plasma response to periodic toroidal phase flips of a Resonant Magnetic Perturbation (RMP) field is experimentally identified in the DIII-D tokamak using phase minima of the modulated plasma density and toroidal rotation relative to the RMP field. Furthermore, the plasma response coincides with the q=3 rational surface and electron fluid velocity null, which is consistent with simulations of the plasma response to the RMP field from a resistive Magnetohydrodynamics modeling. We also observe an asymmetric propagation of the particle and the momentum from the resonant location of the plasma response to the RMP into to core and into the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Final Report for Center for Tokamak Transient Simulations at USU

Providing plasma fluid codes like NIMROD with continuum drift kinetic (CDK) physics that is quantitatively valid and computationally feasible throughout the spatial domain is difficult. Work at Utah State University (USU), in collaboration with the Center for Tokamak Transient Simulations (CTTS), focused on applying CDK closures in disruption-related calculations. Three examples where kinetic physics is paramount are (1) the electron stress tensor closure in Ohms law for accurately describing neoclassical tearing mode (NTM) evolution, (2) runaway electron (RE) density (nRE) and current (jRE) moments in NIMROD’s extended MHD model for self-consistent evolution of RE populations during disruptions and, (3) energetic ion effects on a myriad of MHD instabilities. While NTM simulations and continuum and PIC approaches to energetic ions in NIMROD have been a major goals of USU’s closure work for several years, the development of self-consistent CDK RE capability in NIMROD was started and extended considerably during the CTTS effort. Some goals of CDK RE in NIMROD are to explore the effects of the 2D relativistic phase space in 4D simulations and compare with NIMROD’s fluid RE model. Four publications and two PhD theses came out of the USU CTTS effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗