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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Disruption of Current Filaments and Isotropization of the Magnetic Field in Counterstreaming Plasmas

In this paper, we study the stability of current filaments produced by the Weibel, or current filamentation, instability in weakly magnetized counterstreaming plasmas. It is shown that a resonance exists between the current-carrying ions and a longitudinal drift-kink mode that strongly deforms and eventually breaks the current filaments. Analytical estimates of the wavelength, growth rate, and saturation level of the resonant mode are derived and validated by three-dimensional particle-in-cell simulations. Finally and furthermore, self-consistent simulations of counterstreaming plasmas indicate that this drift-kink mode is dominant in the slow down of the flows and in the isotropization of the magnetic field, playing an important role in the formation of collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

CRADA Final Report: CRADA Number NFE-22-09330 with General Fusion

General Fusion is developing a magnetized target fusion (MTF) approach that involves compressing an initial magnetically confined plasma inside a cavity formed in liquid metal. This approach builds from concepts initially developed under the Linus program at the U.S. Naval Research Laboratory and combines it with advances from compact toroid experiment (CTX) and sustained spheromak plasma experiment (SSPX) in compact toroid plasmas and coaxial Marshall gun systems. Modeling the tokamak during compression is central to designing a successful MTF device. The plasma is formed by coaxial helicity injection in the General Fusion device. Immediately after formation, the plasma has a diverted tokamak configuration with a single null. As the wall moves inwards, the plasma is repelled from the conducting surface and driven inwards by currents induced by its magnetic field in the liquid metal wall. As the liquid metal closes (or bridges) the opening of the coaxial plasma injector, the magnetic field topology alters to remove the null. Due to this, the plasma moves from a diverted to a wall-limited configuration. The liquid metal liner continues to close in and change shape, reducing in radius by a factor of ten at the peak of plasma compression. A model of the MTF plasma must be able to handle this continually varying geometry, and to be predictive, it must faithfully include the real imperfections arising in the process. In this project, we pursued a Monte Carlo approach to closures for MHD by computing kinetic electron trajectories in an MHD plasma background from simulations of GF devices. This requires enhancing the capabilities of the KORC-T code for running large ensembles of kinetic trajectories by porting it to GPU architectures and enabling workflows for large ensembles on OLCF machines. With these capabilities, it is possible to produce a large library of kinetic calculations of electron orbits evolving in plasma configurations spanning the magnetic configurations and plasma density profiles, including non-axisymmetry, arising in the General Fusion’s existing PI3 spherical tokamak device. Using ensembles will capture particles passing a single point in space in a given magnetic configuration, and the entire dataset will cover a range of global magnetic field geometries. By sampling around many starting points, this dataset will capture the spatial dependence of the plasma parameters. From this large dataset, it is possible to produce a reduced model for the kinetic effects not captured in MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Comparison of Saturation Rules Used for Gyrokinetic Quasilinear Transport Modeling

Theory-based transport modeling has been widely successful and is built on the foundations of quasilinear theory. Specifically, the quasilinear expression of the flux can be used in combination with a saturation rule for the toroidal mode amplitude. Most transport models follow this approach. Saturation rules are heuristic and difficult to rigorously derive. We compare three common saturation rules using a fairly accurate quasilinear expression for the fluxes computed using local linear gyrokinetic simulation. We take plasma parameters from experimental H-mode profiles and magnetic equilibrium and include electrons, deuterium, and carbon species. We find that the various saturation rules provide qualitatively similar behavior. This may help to explain why the different theory-based transport models can all predict core tokamak profiles reasonably well. Comparisons with nonlinear local and global gyrokinetic simulations are discussed.

Parker, Scott E. (ORCID:000000030942630X)↗
DOE OSTI2019

Mean-force scattering potential for calculating optical properties of dense plasmas

Here, we assess the relative importance of ionic structure on the opacity of dense plasmas by using the potential of mean force as a scattering potential within the Kubo-Greenwood formalism. We compare results from the potential of mean force, which includes realistic ionic structure, to results using an average atom potential, which includes only a crude treatment of ionic structure. Comparisons with less approximate but more expensive DFT-MD simulations for aluminum plasma reveal that the mean force generally improves agreement for DC conductivity. We also see improvement when applying the mean force to free-free transitions, whereas for bound-bound and bound-free transitions the mean force leads to poorer agreement on transition energies. Further, we assess the impact of accounting for correlations within the plasma at the temperature and density conditions relevant to iron opacity measurements at Sandia’s Z machine facility and find that these correlations do not account for the discrepancy between the measurements and leading opacity calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Modeling K-DEMO Plasma Power Exhaust with UEDGE

The goal of this project is to develop simulation capability for the edge plasma region of K-DEMO in order to assess the heat and particle loads on plasma facing components. The roles of Lawrence Livermore National Laboratory (LLNL) are to provide the UEDGE simulation code for this purpose, to obtain initial edge-plasma solutions for K-DEMO, and to train K-DEMO staff members in both the use of the 2D UEDGE and the code structure such that future modification can be made. A schematic of K-DEMO highlighting the power exhaust is shown in Fig. 1. The major radius of the device is 6.8 m, and the total power to be exhausted is 280 MW in the edge plasma region. Evaluating impact of different divertor geometries and the use of seeded impurity radiation to manage the peak heat flux to divertor-plate surfaces are key aspects that have been addressed during this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Experiment and numerical simulation of RF heating in the Tandem Mirror plasma propulsion device

Recent results of experimental and theoretical studies of ICRF heating of the plasma in the Tendem Mirror rocket are presented. The radial and axial profile of the magnetic field of the wave have been measured, and the data agree with the results from numerical simulation. One very important new finding is that the wave damped when it approached the resonance plane. This is a strong indication of beach heating effect and that RF power is absorbed by the ions in the plasma as expected. This power absorption phenomenon was also further confirmed by the refined analytical study of the wave propagation in a slab model. The electron density in the central cell has been measured by a microwave interferometer. The existence of two types of discharge similar to those of H-alpha emission have also been observed from the density measurements. A preliminary design of a flight system for a 100-day Mars transit has been undertaken. The specific weight for a 10 MW rocket is only 0.04 kg/kW.

Yang, T. F.↗
DOE OSTI2017

Kinetic Simulations of the Interruption of Large-Amplitude Shear-Alfvén Waves in a High- β Plasma

Using two-dimensional hybrid-kinetic simulations, we explore the nonlinear “interruption” of standing and traveling shear-Alfvén waves in collisionless plasmas. Interruption involves a self-generated pressure anisotropy removing the restoring force of a linearly polarized Alfvénic perturbation, and occurs for wave amplitudes δB ⊥ /B 0 ≳β –1/2 (where β is the ratio of thermal to magnetic pressure). We use highly elongated domains to obtain maximal scale separation between the wave and the ion gyroscale. For standing waves above the amplitude limit, we find that the large-scale magnetic field of the wave decays rapidly. The dynamics are strongly affected by the excitation of oblique firehose modes, which transition into long-lived parallel fluctuations at the ion gyroscale and cause significant particle scattering. Traveling waves are damped more slowly, but are also influenced by small-scale parallel fluctuations created by the decay of firehose modes. Our results demonstrate that collisionless plasmas cannot support linearly polarized Alfvén waves above δB ⊥ /B 0 ~β –1/2 . Here, they also provide a vivid illustration of two key aspects of low-collisionality plasma dynamics: (i) the importance of velocity-space instabilities in regulating plasma dynamics at high β, and (ii) how nonlinear collisionless processes can transfer mechanical energy directly from the largest scales into thermal energy and microscale fluctuations, without the need for a scale-by-scale turbulent cascade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Passive deconfinement of runaway electrons using an in-vessel helical coil

A helical coil designed to passively generate non-axisymmetric fields during a plasma disruption is shown (via electromagnetic analysis, linear MHD modeling, and relativistic drift orbit tracing) to be effective at deconfining runaway electrons (REs) on a time scale significantly faster than the plasma current quench. Magnetic equilibria from DIII-D RE-producing scenarios are used to calculate the toroidal electric field generated during the current quench phase of a disruption, which in turn drives current in the proposed n = 1 in-vessel helical coil, without the need for any external power supplies or disruption detection or prediction techniques. Simulations of the plasma evolution using the TokSys GS Evolve code predict the inductive coupling of coil currents up to 12% of the pre-disruption plasma current into the helical coil. The coil geometry is parametrically varied to maximize both the non-resonant and resonant components of the 3D magnetic perturbation, resulting in δB/B ≈ 10 –2 and a vacuumisland overlapwidth of up to 0.7ψ N . The REORBIT module of the MARS-F code is used to model the full non-axisymmetric magnetic field and trace RE drift orbits to determine the effect on RE deconfinement, with up to 70% of the RE orbits lost after 0.2 ms. A two-stage evolution of the RE orbit loss fraction is observed to be caused by resonant trapping between multiple magnetic island chains. Finally, electromagnetic and thermal stresses on the coil are calculated to be within operational limits for installation in DIII-D, and scale favorably to a reactor-size device. Furthermore, these findings motivate future experimental study of the helical coil concept in DIII-D or other tokamaks.

3D coil↗
DOE OSTI2020

Interpenetration and kinetic effects in converging, high-energy plasma jets

In this work, we report on numerical simulations of laser-driven convergent plasma fusion targets. These “inverted corona” fusion targets are useful for the study of counter-streaming and converging rarefied plasma flows, and previous experiments have demonstrated their potential as neutron sources. The scheme consists of a fuel layer lined along the interior surface of a hollow plastic shell that is laser-ablated and expands inward towards the target center. The plasma streams generated in these targets are initially nearly collisionless as they converge, leading to wide interaction length scales and long interaction time scales as the jets interpenetrate. Such kinetic effects impact mixing of constituent ions - a phenomenon not properly captured by single-fluid hydrodynamic simulations. Here we conduct numerical simulations using two different methods: (1) single-fluid simulations in HYDRA, and (2) kinetic-ion, fluid-electron hybrid particle-in-cell (PIC) simulations in the code Chicago. It is shown that the initially nearly collisionless plasma fronts interpenetrate deeply and lead to broader interaction regions in space and time resulting in significant beam-beam fusion. The two approaches make different, testable predictions for the effect of fuel-layer thickness on neutron yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Recent progress in plasma modelling at INFN-LNS

At Istituto Nazionale di Fisica Nucleare - Laboratori Nazionali del Sud (INFN-LNS), the development of intense ion and proton sources has been supported by a great deal of work on the modelling of microwave generated plasmas for many years. First, a stationary version of the particle-in-cell code was developed for plasma modelling starting from an iterative strategy adopted for the space charge dominated beam transport simulations. Electromagnetic properties of the plasma and full-waves simulations are now affordable for non-homogenous and non-isotropic magnetized plasma via “cold” approximation. The effects of Coulomb collisions on plasma particles dynamics was implemented with the Langevin formalism, instead of simply applying the Spitzer 90° collisions through a Monte Carlo technique. A wide database of different cross sections related to reactions occurring in a hydrogen plasma was implemented. The next step consists of merging such a variety of approaches for retrieving an “as-a-whole” picture of plasma dynamics in ion sources. The preliminary results will be summarized in the paper for a microwave discharge ion source designed for intense and high quality proton beams production, proton source for European Spallation Source project. Even if the realization of a predictive software including the complete processes involved in plasma formation is still rather far, a better comprehension of the source behavior is possible and so the simulations may support the optimization phase.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Development of a high charge 10 GeV laser electron accelerator

Low-density meter-scale plasma waveguides produced in meter-scale supersonic gas jets have paved the way for recent demonstrations of all-optical multi-gigaelectronvolt laser wakefield acceleration (LWFA). This paper reviews recent advances by the University of Maryland, which have enabled these results, focusing on the development of elongated supersonic gas jets up to ∼1 m in length, experimental and simulation studies of plasma waveguide formation, and a new three-stage model for relativistic pulse propagation dynamics in these waveguides. We also present results from recent LWFA experiments conducted at the Laboratory for Advanced Lasers and Extreme Photonics at Colorado State University demonstrating high charge, low divergence electron bunches to ∼10 GeV, with laser-to-electron beam efficiency of at least ∼30%.

Rockafellow, E. [Univ. of Maryland, College Park, ↗
DOE OSTI2025

MHD Analysis of Dual-Coolant Lead-Lithium Blanket for Spherical Tokamak Advanced Reactor

The tritium breeding blanket is vital for future fusion power plants, with the Spherical Tokamak Advanced Reactor (STAR) project highlighting the dual-coolant lead-lithium (DCLL) design. The DCLL blanket performs shielding, energy exhaust, and tritium breeding using a lead-lithium alloy , with lithium as the breeder and lead as the neutron multiplier. It also serves as the primary coolant, with helium providing supplemental cooling. Reduced-activation ferritic/martensitic steel is used for the blanket structure. Magnetohydrodynamic (MHD) phenomena influence the liquid metal flow in a magnetic field, affecting heat transfer in the breeder affected by energetic neutrons. Understanding key flow parameters in such conditions is critical for efficient DCLL design. This study uses three-dimensional thermofluid MHD analysis with ANSYS CFX software, modified at Princeton Plasma Physics Laboratory, to simulate high Hartmann flows. The neutronics code MCNP, coupled with plasma equilibrium, provides heat source distribution. In conclusion, we examine electromagnetic interactions in adjacent fluid domains and analyze the magnetic field’s impact on flow distribution in the inboard and outboard blanket layout, using detailed mesh generation for accurate results.

DCLL↗
NASA NTRS1987

Optical efficiencies of lightning in planetary atmospheres

The first simulations of lightning in planetary atmospheres by laser-induced plasmas are reported. These simulations show that the fraction of the energy in lightning discharge channels that is radiated in the visible spectrum is similar for earth, Venus, and Titan, but quite different for Jupiter. One implication of these results is that the amount of trace gases produced by lightning in the Jovian atmosphere must be larger than previously estimated.

Borucki, W. J.↗
NASA NTRSDate not supplied

Nanoflare Heating Frequency of an X-Ray Bright Point Observed By MaGIXS

Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multimillion kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission from an X-ray bright point observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with concurrent observations from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO) and X-Ray Telescope (XRT) onboard Hinode observatory. The length and magnetic field strength of the coronal loops are derived from the potential field extrapolation of the observed photospheric magnetogram by Helioseismic and Magnetic Imager (HMI) onboard SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the loop length and magnetic field strength. The simulation results are then compared and matched against the measured intensity from AIA, XRT, and MaGIXS. Our model results indicate the loop morphology and emissions from the XBP under study could be well matched by a distribution of nanoflares with average delay times 400 s to 800 s, which strongly suggest that the heating is dominated by high-frequency events. Further, we demonstrate the high sensitivity of MaGIXS and XRT to diagnose the heating frequency using this method, while AIA passbands are found to be the least sensitive.

coronal heating↗
NASA NTRSDate not supplied

Determining The Nanoflare Heating Frequency of an X-Ray Bright Point Observed by MaGIXS

Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multi-million kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission from an X-ray bright point (XBP) observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with nearest-available observations from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO) and X-Ray Telescope (XRT) onboard Hinode observatory. The length and magnetic field strength of the coronal loops are derived from the linear-force-free extrapolation of the observed photospheric magnetogram by Helioseismic and Magnetic Imager (HMI) onboard SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the loop length and magnetic field strength. The simulation results are then compared and matched against the measured intensity from AIA, XRT, and MaGIXS. Our model results indicate the observed emissions from the XBP under study could be well matched by a distribution of nanoflares with average delay times 1500 s to 3000 s, which suggest that the heating is dominated by high-frequency events. Further, we demonstrate the high sensitivity of MaGIXS and XRT to diagnose the heating frequency using this method, while AIA passbands are found to be the least sensitive.

coronal heating↗
DOE OSTI2017

Numerical study of core formation of asymmetrically driven cone-guided targets

Compression of a directly driven fast ignition cone-sphere target with a finite number of laser beams is numerically studied using a three-dimensional hydrodynamics code IMPACT-3D. The formation of a dense plasma core is simulated for 12-, 9-, 6-, and 4-beam configurations of the GEKKO XII laser. The complex 3D shapes of the cores are analyzed by elucidating synthetic 2D x-ray radiographic images in two orthogonal directions. Finally, the simulated x-ray images show significant differences in the core shape between the two viewing directions and rotation of the stagnating core axis in the top view for the axisymmetric 9- and 6-beam configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2009

Efficacy of Compression Garments to Simulate Fluid Shifts during Lunar Bedrest

The effectiveness of JOBST(Registered TradeMark) compression stockings for fluid redistribution was examined in the lunar bed rest analog based on Digital Astronaut computer model predictions. NASA's future goals include missions to the moon which will require extended exposure to the lunar gravity environment. To prepare for these missions, physiological adaptations in various systems must be resolved. A bed rest at a 10 degree head-up tilt was used to simulate lunar gravity. Plasma volume (PV) loss may cause some of the cardiovascular adaptations which occur during space flight so PV is measured in the bed rest analog to evaluate fluid loss. PV varies from individual to individual so PV index (PVI) is used to determine the magnitude and time course of fluid shift and cardiovascular adaptation to 1/6 g. The Digital Astronaut, a computer simulation tool, predicts a 6% PV loss during an extended simulated lunar mission for a male with a body surface area of 1.95 m(exp 2). Simple geometry calculations suggest that 10 degree head-up tilt is most useful for simulated measurement of deconditioning in bone and muscle, however, 2 degree head-up tilt may best imitate cardiovascular fluid shifts. In order to reconcile these different models, compression stocking must be used in the 10 degree paradigm to better approximate expected cardiovascular changes.

Brinley, Alaina↗
DOE OSTI2019

Real-time capable modeling of neutral beam injection on NSTX-U using neural networks

A new model of heating, current drive, torque and other effects of neutral beam injection on NSTX-U that uses neural networks has been developed. The model has been trained and tested on the results of the Monte Carlo code NUBEAM for the database of experimental discharges taken during the first operational campaign of NSTX-U. By projecting flux surface quantities onto empirically derived basis functions, the model is able to efficiently and accurately reproduce the behavior of both scalars, like the total neutron rate and shine through, and profiles, like beam current drive and heating. With this being said, the model has been tested on the NSTX-U real-time computer, demonstrating a rapid execution time orders of magnitude faster than the Monte Carlo code that is well suited for the iterative calculations needed to interpret experimental results, optimization during scenario development activities, and real-time plasma control applications. Simulation results of a proposed design for a nonlinear observer that embeds the neural network calculations to estimate the poloidal flux profile evolution, as well as Z eff and fast ion diffusivity, are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗