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,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Mixing and unmixedness in plasma jets 1: Near-field analysis

The flow characteristics in the near-field of a plasma jet are simulated with a two-fluid model. This model accounts for both gradient-diffusion mixing and uni-directional sifting motion resulting from pressure-gradient-body-force imbalance. This latter mechanism is believed to be responsible for the umixedness observed in plasma jets. The unmixedness is considered to be essentially a Rayleigh-Taylor kind instability. Transport equations are solved for the individual plasma and ambient gas velocities, temperatures and volume fractions. Empirical relations are employed for the interface transfers of mass, momentum and heat. The empirical coefficients are first established by comparison of predictions with available experimental data for shear flows. The model is then applied to an Argon plasma jet ejecting into stagnant air. The predicted results show the significant build-up of unmixed air within the plasma gas, even relatively far downstream of the torch. By adjusting the inlet condition, the model adequately reproduces the experimental data.

Ilegbusi, Olusegun J.↗
DOE OSTI2025

Impact of nitrogen molecular breakup on divertor conditions in JET L-mode plasmas using SOLPS-ITER

SOLPS-ITER simulations of nitrogen-seeded, low-confinement mode plasmas in the Joint European Torus (JET) predict that the electron temperature in the low-field side (LFS) divertor leg is reduced locally by up to an order of magnitude when nitrogen is assumed to recycle as molecules (N 2 ) instead of atoms using a fixed nitrogen injection rate. The LFS divertor temperature reduction under the assumption of molecular recycling occurs due to a three-step mechanism: (1) the plasma penetration of nitrogen atoms is increased due to the strong triple bond of the N 2 molecule and the kinetic energy release in the dissociation event, both mechanisms contributing equally, (2) the abundance of (particularly multiply-charged) nitrogen ions in the divertor is increased and (3) the electron temperature is reduced due to the increase in radiation (by up to a factor of 4) from nitrogen ions. Setting the volume-integrated nitrogen radiated power to a constant value (0.6 MW) instead of the nitrogen injection rate, SOLPS-ITER predicts under the molecular nitrogen recycling assumption that the peak line-integrated N II, N III and N IV intensities in the LFS divertor are approximately within 15%, 35% and 5%, respectively, of the reference atomic nitrogen recycling case. The predicted peak N II, N III and N IV intensities under either assumption are within 30%, 65% and 5%, respectively, of measurements using the vertically viewing mirror-link divertor spectrometer (Meigs et al., 2010) in nitrogen-seeded JET L-mode plasmas (Lomanowski et al., 2019). ERO2.0 simulations using a constant nitrogen seeding rate on static background plasma solutions from EDGE2D-EIRENE (previously presented in Mäenpää et al., (2022), revised here to include fast reflections) predict that N II to N IV line emission is increased by 20% to 30% when nitrogen is assumed to recycle as molecules, demonstrating the importance of considering the effect of molecular dissociation reactions on the divertor plasma in a self-consistent manner.

Divertor↗
DOE OSTI2025

Combined Influence of Rotation and Scrape-Off Layer Drifts on Recycling Asymmetries in Tokamak Plasmas

Coupled 2D fluid-kinetic simulations of a DIII-D high confinement tokamak plasma show that plasma rotation coupled with drift effects near the plasma edge play a significant role in the creation of the observed poloidal distribution of neutrals. It is observed that including either drift or rotation effects enhances particle flux at the inner target in the case of ion 𝐵×∇𝐵 drift toward the 𝑋-point. However, the particle flux asymmetry is significantly higher with the combination of drifts and rotation than either effect alone. The heightened particle flux asymmetry allows for improved simulation of the strong in-out asymmetry of the Lyman-𝛼 brightness profiles measured in the experiment. Enhancement of radial transport of parallel momentum changes the upstream scrape-off layer flow pattern, increasing the fraction of deuterium flux that reaches the inboard divertor entrance while lowering that which arrives at the outboard. In conclusion, this Letter indicates that by combining drifts, rotation, and viscous coupling, existing boundary plasma models can achieve a satisfactory agreement with experimentally measured neutral asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Study of the characteristics of reconfigurable plasma antenna array

This paper presents a design and simulation of a reconfigurable array of plasma antenna. The plasma column is used as radiating elements instead of metal to create an antenna. The advantages of the plasma antenna over the conventional antenna are its possible to change the operating parameters, such as the working pressure, input power, radius of the discharge tube, resonant frequency, and length of the plasma column. In addition, plasma antenna can be reconfigurable with respect to shape, frequency and radiation parameters in a very short time. The plasma discharge tube was designed with a length of 200 mm, the radius of the plasma column was 2.5 mm and the coupling sleeve was connected to the SMA as the ground. This simulation was performed by using the simulation software Computer Simulation Technology (CST). The frequency is set in the range of 1 GHz to 10 GHz. The performance of the designed antenna was analyzed in term of return loss, gain and radiation pattern. For reconfigurable plasma antenna array, it shows that the gain is increase when the number of antenna element is increase. The combination of the discharge tube and metal rod as an antenna array has been done, and the result shows that an array with the plasma element can achieve higher gain.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Multiple scattering theory for dense plasmas

Dense plasmas occur in stars, giant planets, and in inertial fusion experiments. Accurate modeling of the electronic structure of these plasmas allows for prediction of material properties that can in turn be used to simulate these astrophysical objects and terrestrial experiments. But modeling them remains a challenge. Here we explore the Korringa-Kohn-Rostoker Green's function (KKR-GF) method for this purpose. We find that it is able to predict equation of state in good agreement with other state-of-the-art methods, where they are accurate and viable. In addition, it is shown that the computational cost does not significantly change with temperature, in contrast with other approaches. Moreover, the method does not use pseudopotentials—core states are calculated self consistently. We conclude that KKR-GF is a very promising method for dense plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Impact of the plasma geometry on divertor power exhaust: experimental evidence from $\mathrm{TCV}$ and simulations with SolEdge2D and $\mathrm{TOKAM3X}$

A deep understanding of plasma transport at the edge of magnetically confined fusion plasmas is needed for the handling and control of heat loads on the machine first wall. Experimental observations collected on a number of tokamaks over the last three decades taught us that heat flux profiles at the divertor targets of X-point configurations can be parametrized by using two scale lengths for the scrape-off layer (SOL) transport, separately characterizing the main SOL (λ q ) and the divertor SOL (S q ). Here in this work we challenge the current interpretation of these two scale lengths as well as their dependence on plasma parameters by studying the effect of divertor geometry modifications on heat exhaust in the Tokamak à Configuration Variable. In particular, a significant broadening of the heat flux profiles at the outer divertor target is diagnosed while increasing the length of the outer divertor leg in lower single null, Ohmic, L-mode discharges. Efforts to reproduce this experimental finding with both diffusive (SolEdge2D-EIRENE) and turbulent (TOKAM3X) modelling tools confirm the validity of a diffusive approach for simulating heat flux profiles in more traditional, short leg, configurations while highlighting the need of a turbulent description for modified, long leg, ones in which strongly asymmetric divertor perpendicular transport develops.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Non-perturbative measurement of cross-field thermal diffusivity reduction at the O-point of 2/1 neoclassical tearing mode islands in the DIII-D tokamak

Neoclassical tearing modes (NTMs) often lead to the decrease of plasma performance and can lead to disruptions, which makes them a major impediment in the development of operating scenarios in present toroidal fusion devices. Recent gyrokinetic simulations predict a decrease of plasma turbulence and cross-field transport at the O-point of the islands, which in turn affects the NTM dynamics. In this paper, a heat transport model of magnetic islands employing spatially non-uniform cross-field thermal diffusivity (χ{sub ⊥}) is presented. This model is used to derive χ{sub ⊥} at the O-point from electron temperature data measured across 2/1 NTM islands in DIII-D. It was found that χ{sub ⊥} at the O-point is 1 to 2 orders of magnitude smaller than the background plasma transport, in qualitative agreement with gyrokinetic predictions. As the anomalously large values of χ{sub ⊥} are often attributed to turbulence driven transport, the reduction of the O-point χ{sub ⊥} is consistent with turbulence reduction found in recent experiments. Finally, the implication of reduced χ{sub ⊥} at the O-point on NTM dynamics was investigated using the modified Rutherford equation that predicts a significant effect of reduced χ{sub ⊥} at the O-point on NTM saturation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Impact of electron-scale turbulence and multi-scale interactions in the JET tokamak

Experimental observations in JET tokamak plasmas and gyrokinetic simulations point to an important role, for electron heat transport, of electron-scale instabilities and of their interaction with ion-scale instabilities. Since these effects are maximized for strong electron heating and ion-scale modes close to marginal stability, these findings are of high relevance for ITER plasmas, featuring both conditions. Gyrokinetic and quasi-linear transport models accounting for multi-scale effects are assessed against JET experimental results.

Physics↗
NASA NTRS1987

Solar transition region and coronal response to heating rate perturbations

Observations of Doppler shifts in UV emission lines formed in the solar transition region show continual plasma downflows and impulsive plasma upflows. Using numerical simulations, the authors examine the conjecture that areas of downflowing plasma are the base regions of coronal loops in which the heating is gradually decreasing and that areas of upflowing plasma are the base regions of coronal loops in which the heating rate is gradually increasing. Beginning with a coronal loop in equilibrium, the heating rate is reduced on time scales of 100, 1000, and 2000 s to 10 percent and 1 percent of the initial value, and the loop is allowed to evolve to a new equilibrium. The heating rate for the cooled models is then increased back to the initial value on the same time scales. While significant mass motions do develop in the simulations, both the emission measure and the velocity at 100,000 K do not show the characteristics present in UV observations.

Mariska, John T.↗
DOE OSTI2015

Testing the high turbulence level breakdown of low-frequency gyrokinetics against high-frequency cyclokinetic simulations

This paper presents numerical simulations of the nonlinear cyclokinetic equations in the cyclotron harmonic representation [R. E. Waltz and Zhao Deng, Phys. Plasmas 20, 012507 (2013)]. Simulations are done with a local flux-tube geometry and with the parallel motion and variation suppressed using a newly developed rCYCLO code. Cyclokinetic simulations dynamically follow the high-frequency ion gyro-phase motion which is nonlinearly coupled into the low-frequency drift-waves possibly interrupting and suppressing gyro-averaging and increasing the transport over gyrokinetic levels. By comparing the more fundamental cyclokinetic simulations with the corresponding gyrokinetic simulations, the breakdown of gyrokinetics at high turbulence levels is quantitatively tested over a range of relative ion cyclotron frequency 10 < Ω*{sup }< 100 where Ω*{sup }= 1/ρ*, and ρ* is the relative ion gyroradius. The gyrokinetic linear mode rates closely match the cyclokinetic low-frequency rates for Ω*{sup }> 5. Gyrokinetic transport recovers cyclokinetic transport at high relative ion cyclotron frequency (Ω*{sup }≥ 50) and low turbulence level as required. Cyclokinetic transport is found to be lower than gyrokinetic transport at high turbulence levels and low-Ω* values with stable ion cyclotron (IC) modes. The gyrokinetic approximation is found to break down when the density perturbations exceed 20%. For cyclokinetic simulations with sufficiently unstable IC modes and sufficiently low Ω*{sup }∼ 10, the high-frequency component of cyclokinetic transport level can exceed the gyrokinetic transport level. However, the low-frequency component of the cyclokinetic transport and turbulence level does not exceed that of gyrokinetics. At higher and more physically relevant Ω*{sup }≥ 50 values and physically realistic IC driving rates, the low-frequency component of the cyclokinetic transport and turbulence level is still smaller than that of gyrokinetics. Thus, the cyclokinetic simulations do not account for the so-called “L-mode near edge short fall” seen in some low-frequency gyrokinetic transport and turbulence simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Numerical simulation of positive-potential conductors in the presence of a plasma and a secondary-emitting dielectric

A 2 1/2-dimensional, cylindrically symmetric particle in cell program is used to simulate the interaction of a plasma with a positively biased disk covered by a dielectric material (Kapton). Several runs using different combinations of bias voltages and secondary electron yields are discussed. Special attention is paid to the role of the secondaries as the dielectric surface changes from a negative potential to become increasingly positive. This is termed the 'snapover' process. The contribution due to secondaries is examined by plotting the distribution (number versus position) of the particles emitted by the impacting on the dielectric. These plots show that secondaries are drawn back into the dielectric as well as the conductor in the case of snapover. Mechanisms and experimental implications for these particular distributions are discussed.

Kessel, R. L.↗
DOE OSTI2021

Modeling of carbon pellets disruption mitigation in an NSTX-U plasma

Single carbon pellet disruption mitigation simulations using M3D-C1 were conducted for an NSTX-U-like plasma to support the electromagnetic pellet injection concept (EPI). A carbon ablation model has been implemented in M3D-C1 and tested with available data. 2D simulations were conducted in order to estimate the amount of carbon needed to quench the plasma, finding that the content in a 1 mm radius vitreous carbon pellet (~3.2x10E20 atoms) would be enough if it is entirely ablated. 3D simulations were performed, scanning over pellet velocity and parallel thermal conductivity, as well as different injection directions and pellet concepts (solid pellets and shell pellets). The sensitivity of the thermal quench and other related quantities to these parameters has been evaluated. A 1 mm radius solid pellet only partially ablates at velocities of 300 m/s or higher, thus being unable to fully quench the plasma. To further enhance the ablation, approximations to an array of pellets and the shell pellet concept were also explored. 3D field line stochastization plays an important role in both quenching the center of the plasma and in heat flux losses, thus lowering the amount of carbon needed to mitigate the plasma when compared to the 2D case. This study constitutes an important step forward in `predict-first' simulations for disruption mitigation in NSTX-U and other devices, such as ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The Influence of Environment on Post-Detonation Chemistry and Debris Formation (Abbreviated Final Report: 20-SI-006)

Predicting, responding to, or interpreting the chemical record preserved in debris derived from nuclear events can be challenging due to chemical fractionation. Chemical fractionation is where different species of the evolving radionuclide inventory segregate and/or are lost from the system over the timescales of debris formation. Both historic data and recent research suggest that the interaction and character of the local environment may exert controls on chemical fractionation by influencing the cooling and evolution of the associated fireball as well as the composition of the vapor term and resultant speciation. Prior to this work, an integrated platform permitting dynamic and concurrent consideration of physical and chemical evolution of early time post-detonation event environments did not exist. Our work merged historic data and experimental approaches to support development of a computational framework able to simulate fundamental processes (e.g., entrainment of local environment, oxidation chemistry, and cooling time scales) that may perturb the radionuclide inventory captured in post-detonation debris. Work with historic debris confirmed that entrained environmental material affect debris composition, structure, and radionuclide incorporation. Complementary work utilizing a readily controllable and tunable benchtop setup (a plasma flow reactor) simulated the late cooling of a nuclear fireball (e.g., T < 6000 K) and bounded the sensitivity of actinide speciation and particle size distribution to variations in oxygen concentration and cooling rates. Concurrent laser ablation and laser heating experiments were used to investigate the chemistry and physics of processes occurring in vaporized and/or rapidly heated actinides and other elements in the presence of oxygen. A more computationally efficient microphysical model was developed for predicting and evolving size distributions of particles forming from mixed vapor terms and simulating particle formation processes under a variety of extreme conditions. Continued study of historic nuclear event film confirmed that shockwave data and physics codes agree to within the uncertainty of the data. Good agreement was achieved for thermal emission from an airburst, however the paucity of low-temperature molecular opacity data for mixtures of air, bomb debris, entrained dirt, and water vapor complicate agreement for more elaborate scenarios. A multiphysics code (ALE3D) was modified to bring the necessary physics and chemistry, including these new data and insights, onto a single platform. Code development included improved initialization of large physical systems, modernization of chemistry capabilities, and modifications to enable inclusion of particle transport.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

NONTHERMALLY DOMINATED ELECTRON ACCELERATION DURING MAGNETIC RECONNECTION IN A LOW-β PLASMA

By means of fully kinetic simulations, we investigate electron acceleration during magnetic reconnection in a nonrelativistic proton–electron plasma with conditions similar to solar corona and flares. We demonstrate that reconnection leads to a nonthermally dominated electron acceleration with a power-law energy distribution in the nonrelativistic low-β regime but not in the high-β regime, where β is the ratio of the plasma thermal pressure and the magnetic pressure. The accelerated electrons contain most of the dissipated magnetic energy in the low-β regime. A guiding-center current description is used to reveal the role of electron drift motions during the bulk nonthermal energization. We find that the main acceleration mechanism is a Fermi-type acceleration accomplished by the particle curvature drift motion along the electric field induced by the reconnection outflows. Although the acceleration mechanism is similar for different plasma β, low-β reconnection drives fast acceleration on Alfvénic timescales and develops power laws out of thermal distribution. The nonthermally dominated acceleration resulting from magnetic reconnection in low-β plasma may have strong implications for the highly efficient electron acceleration in solar flares and other astrophysical systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Simulation of whistler waves excited in the presence of a cold plasma cloud - Implications for the CRRES mission

A one-dimensional electromagnetic particle simulation model is constructed to study the excitation of whistler waves in the presence of a cold plasma cloud for conditions representative of those after the release of lithium in the inner plasma sheet during the Combined Release and Radiation Effect Satellite mission. The results indicate that a standing-wave pattern with discrete wave frequencies is formed within the cloud. The magnetic wave amplitude inside the cloud, which is limited by quasi-linear diffusion, is of the order of several nanoteslas. Assuming a magnetospheric loss cone of 5 deg, the observed pitch angle diffusion produced by the whistler waves is sufficient to put the electrons on strong diffusion.

Pritchett, P. L.↗
NASA NTRS2008

A Hot Microflare Observed with RHESSI and HINODE

RHESSI and Hinode observations of a GOES B-class flare are combined to investigate the origin of 15 MK plasma. The absence of any detectable hard X-ray emission coupled with weak blueshifted emission lines (indicating upward velocities averaging only 14 km/s) suggests that this was a result of direct heating in the corona, as opposed to nonthermal electron precipitation causing chromospheric evaporation. These findings are in agreement with a recent hydrodynamical simulation of microflare plasmas that found that higher temperatures can be attained when less energy is used to accelerate electrons out of the thermal distribution. In addition, unusual redshifts in the 2 MK Fe xv line (indicating downward velocities of approx.14 km/s) were observed cospatial with one of the flare ribbons during the event. Downflows of such high-temperature plasma are not predicted by any common flare model.

Milligan, Ryan O.↗
DOE OSTI2015

First-principles equation of state of polystyrene and its effect on inertial confinement fusion implosions

Obtaining an accurate equation of state (EOS) of polystyrene (CH) is crucial to reliably design inertial confinement fusion (ICF) capsules using CH/CH-based ablators. Thus, with first-principles calculations, we have investigated the extended EOS of CH over a wide range of plasma conditions (ρ = 0.1 to 100 g/cm 3 and T = 1,000 to 4,000,000 K). When compared with the widely used SESAME -EOS table, the first-principles equation of state (FPEOS) of CH has shown significant differences in the low-temperature regime, in which strong coupling and electron degeneracy play an essential role in determining plasma properties. Hydrodynamic simulations of cryogenic target implosions on OMEGA using the FPEOS table of CH have predicted ~5% reduction in implosion velocity and ~30% decrease in neutron yield in comparison with the usual SESAME simulations. This is attributed to the ~10% lower mass ablation rate of CH predicted by FPEOS. Simulations using CH-FPEOS show better agreement with measurements of Hugoniot temperature and scattered lights from ICF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗