Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

UBIQUITOUS NON-THERMALS IN ASTROPHYSICAL PLASMAS: RESTATING THE DIFFICULTY OF MAINTAINING MAXWELLIANS

This paper outlines the rather narrow conditions on a radiatively decoupled plasma where a Maxwell-Boltzmann (MB) distribution can be assumed with confidence. The complementary non-thermal distribution with non-perturbative kurtosis is argued to have a much broader purview than has previously been accepted. These conditions are expressed in terms of the electron Knudsen number, K{sub e} , the ratio of the electron mean free path to the scale length of electron pressure. Rather generally, f(v < v{sub 2}(K{sub e} )) will be Gaussian, so that MB atomic or wave particle effects controlled by speeds v < v{sub 2} {identical_to} w(15/8K{sub e} ){sup 1/4} will remain defensible, where w is the most probable speed. The sufficient condition for Spitzer-Braginskii plasma fluid closure at the energy equation requires globally K{sub e} (s) {<=} 0.01; this global condition pertains to the maximum value of K{sub e} along the arc length s of the magnetic field (to its extremities) provided that contiguous plasma remains uncoupled from the radiation field. The non-thermal regime K{sub e} > 0.01 is common in all main-sequence stellar atmospheres above approximately 0.05 stellar radii from the surface. The entire solar corona and wind are included in this regime where non-thermal distributions with kurtosis are shown to be ubiquitous, heat flux is not well modeled by Spitzer-Braginskii closure, and fluid modeling is qualitative at best.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

AN EXACT SOLUTION FOR MAGNETIC ANNIHILATION IN A CURVED CURRENT SHEET

An exact magnetohydrodynamic solution is presented for steady magnetic annihilation (merging) in an incompressible resistive viscous plasma. The merging, driven by an axisymmetric stagnation flow on a cylinder, takes place in a curved current sheet that is perpendicular to the plane in which the plasma flow stagnates. The new solution extends earlier models of flux pileup merging in a flat current sheet, driven by stagnation-point flows. The new solution remains valid in the presence of both the isotropic and anisotropic (parallel) plasma viscosity. The geometry of the solution may make it useful in modeling the photospheric flux cancellation on the Sun.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Measurement of the dispersion of radiation from a steady cosmological source

The 'missing baryons' of the near universe are believed to be principally in a partially ionized state. Although passing electromagnetic waves are dispersed by the plasma, the effect has hitherto not been utilized as a means of detection because it is generally believed that a successful observation requires the background source to be highly variable, i.e., the class of sources that could potentially deliver a verdict is limited. We argue in two stages that this condition is not necessary. First, by modeling the fluctuations on macroscopic scales as interference between wave packets, we show that, in accordance with the ideas advanced by Einstein in 1917, both the behavior of photons as bosons (i.e., the intensity variance has contributions from Poisson and phase noise) and the van-Cittert-Zernike theorem are a consequence of wave-particle duality. Nevertheless, we then point out that, in general, the variance on some macroscopic timescale τ consists of (1) a main contributing term ∝1/τ, plus (2) a small negative term ∝1/τ{sup 2} due to the finite size of the wave packets. If the radiation passes through a dispersive medium, this size will be enlarged well beyond its vacuum minimum value of Δt ≈ 1/Δν, leading to a more negative (2) term (while (1) remains unchanged), and hence a suppression of the variance wrt the vacuum scenario. The phenomenon, which is typically at a few parts in 10{sup 5} level, enables one to measure cosmological dispersion in principle. Signal-to-noise estimates, along with systematic issues and how to overcome them, will be presented.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

X-ray and EUV spectroscopy of various astrophysical and laboratory plasmas: Collisional, photoionization and charge-exchange plasmas

Several laboratory facilities were used to benchmark theoretical spectral models that are extensively used by astronomical communities. However, there are still many differences between astrophysical environments and laboratory miniatures that can be archived. Here we setup a spectral analysis system for astrophysical and laboratory plasmas to make a bridge between them, and we investigate the effects from non-thermal electrons and the contributions from a metastable level population on level populations and charge stage distribution for coronal-like, photoionized, and geocoronal plasmas. Test applications to laboratory measurement (i.e., electron beam ion trap plasma) and astrophysical observation (i.e., Comet, Cygnus X-3) are presented. A time evolution of the charge stage and level population are also explored for collisional and photoionized plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Forward modeling of gyrosynchrotron intensity perturbations by sausage modes

To determine the observable radio signatures of the fast sausage standing wave, we examine gyrosynchrotron (GS) emission modulation using a linear three-dimensional magnetohydrodynamic model of a plasma cylinder. Effects of the line-of-sight angle and instrumental resolution on perturbations of the GS intensity are analyzed for two models: a base model with strong Razin suppression and a low-density model in which the Razin effect was unimportant. Our finding contradicts previous predictions made with simpler models: an in-phase variation of intensity between low (f < f {sub peak}) and high (f > f {sub peak}) frequencies is found for the low-density model and an anti-phase variation for the base model in the case of a viewing angle of 45°. The spatially inhomogeneous character of the oscillating emission source and the spatial resolution of the model are found to have a significant effect on the resulting intensity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Onset of electron acceleration in a flare loop

We carried out a detailed analysis of X-ray and radio observations of a simple flare loop that occurred on 2002 August 12, with the impulsive hard X-ray (HXR) light curves dominated by a single pulse. The emission spectra of the early impulsive phase are consistent with an isothermal model in the coronal loop with a temperature reaching several keV. A power-law high-energy spectral tail is evident near the HXR peak time, in accordance with the appearance of footpoints at high energies, and is well correlated with the radio emission. The energy content of the thermal component keeps increasing gradually after the disappearance of this nonthermal component. These results suggest that electron acceleration only covers the central period of a longer and more gradual energy dissipation process and that the electron transport within the loop plays a crucial role in the formation of the inferred power-law electron distribution. The spectral index of power-law photons shows a very gradual evolution, indicating that the electron accelerator is in a quasi-steady state, which is confirmed by radio observations. These results are consistent with the theory of stochastic electron acceleration from a thermal background. Advanced modeling with coupled electron acceleration and spatial transport processes is needed to explain these observations more quantitatively, which may reveal the dependence of the electron acceleration on the spatial structure of the acceleration region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Global structure of isothermal diffuse X-ray emission along the Fermi bubbles

In our previous works, we found absorbed thermal X-ray plasma with kT ≃ 0.3 keV observed ubiquitously near the edges of the Fermi bubbles and interpreted this emission as weakly shock-heated Galactic halo gas. Here we present a systematic and uniform analysis of archival Suzaku (29 pointings; 6 newly presented) and Swift (68 pointings; 49 newly presented) data within Galactic longitudes |l| < 20° and latitude 5°≲ |b| < 60°, covering the whole extent of the Fermi bubbles. We show that the plasma temperature is constant at kT ≃ 0.30 ± 0.07 keV, while the emission measure (EM) varies by an order of magnitude, increasing toward the Galactic center (i.e., low |b|) with enhancements at the North Polar Spur (NPS), SE-claw, and NW-clump features. Moreover, the EM distribution of kT ≃ 0.30 keV plasma is highly asymmetric in the northern and southern bubbles. Although the association of the X-ray emission with the bubbles is not conclusive, we compare the observed EM properties with simple models assuming (i) a filled halo without bubbles, whose gas density follows a hydrostatic isothermal model (King profile), and (ii) a bubble-in-halo in which two identical bubbles expand into the halo, forming thick shells of swept halo gas. We argue that the EM profile in the north (b > 0°) favors (ii), whereas that of the south (b < 0°) is rather close to (i), but a weak excess signature is clearly detected also in the south like NPS (South Polar Spur). Such an asymmetry, if due to the bubbles, cannot be fully understood only by the inclination of bubbles’ axis against the Galactic disk normal, thus suggesting asymmetric outflow due to different environmental/initial conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A SIMPLE METHOD FOR MODELING COLLISION PROCESSES IN PLASMAS WITH A KAPPA ENERGY DISTRIBUTION

We demonstrate that a nonthermal distribution of particles described by a kappa distribution can be accurately approximated by a weighted sum of Maxwell–Boltzmann distributions. We apply this method to modeling collision processes in kappa-distribution plasmas, with a particular focus on atomic processes important for solar physics. The relevant collision process rate coefficients are generated by summing appropriately weighted Maxwellian rate coefficients. This method reproduces the rate coefficients for a kappa distribution to an estimated accuracy of better than 3%. This is equal to or better than the accuracy of rate coefficients generated using “reverse-engineering” methods, which attempt to extract the needed cross sections from the published Maxwellian rate coefficient data and then reconvolve the extracted cross sections with the desired kappa distribution. Our approach of summing Maxwellian rate coefficients is easy to implement using existing spectral analysis software. Moreover, the weights in the sum of the Maxwell–Boltzmann distribution rate coefficients can be found for any value of the parameter κ, thereby enabling one to model plasmas with a time-varying κ. Tabulated Maxwellian fitting parameters are given for specific values of κ from 1.7 to 100. We also provide polynomial fits to these parameters over this entire range. Several applications of our technique are presented, including the plasma equilibrium charge state distribution (CSD), predicting line ratios, modeling the influence of electron impact multiple ionization on the equilibrium CSD of kappa-distribution plasmas, and calculating the time-varying CSD of plasmas during a solar flare.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Physics and technology in the ion-cyclotron range of frequency on Tore Supra and TITAN test facility: implication for ITER

To support the design of an ITER ion-cyclotron range of frequency heating (ICRH) system and to mitigate risks of operation in ITER, CEA has initiated an ambitious Research & Development program accompanied by experiments on Tore Supra or test-bed facility together with a significant modelling effort. The paper summarizes the recent results in the following areas: Comprehensive characterization (experiments and modelling) of a new Faraday screen concept tested on the Tore Supra antenna. A new model is developed for calculating the ICRH sheath rectification at the antenna vicinity. The model is applied to calculate the local heat flux on Tore Supra and ITER ICRH antennas. Full-wave modelling of ITER ICRH heating and current drive scenarios with the EVE code. With 20 MW of power, a current of 400 kA could be driven on axis in the DT scenario. Comparison between DT and DT(3He) scenario is given for heating and current drive efficiencies. First operation of CW test-bed facility, TITAN, designed for ITER ICRH components testing and could host up to a quarter of an ITER antenna. R&D of high permittivity materials to improve load of test facilities to better simulate ITER plasma antenna loading conditions.

Litaudon, X↗
DOE OSTI2016

Suppression of type-I ELMs with reduced RMP coil set on DIII-D

Recent experiments on DIII-D have demonstrated that having a toroidally-monochromatic spectral content of edge-resonant magnetic perturbations (RMPs) is not a necessary condition for suppression of edge localized modes (ELMs). Robust ELM suppression has been reproducibly obtained on DIII-D during experiments in which various non-axisymmetric coil loops were turned off pseudo-randomly producing a variety of n=1 , n=2, and n=3 spectral contributions. It was demonstrated that RMP ELM suppression could be achieved with as few as 5 out of 12 internal coil loops (I-coils) on DIII-D at similar coil currents and with good plasma confinement. Linear MHD plasma response (M3DC1, IPEC, MARS) and vacuum (SURFMN, TRIP3D) modelling have been performed in order to understand the effects of the perturbation spectrum on the plasma response and ELM suppression. The results suggest that reduction of the dominant n=3 perturbation field is compensated by increased n=2 field in the plasma that may lead to RMP ELM suppression at lower levels of n=3 perturbative magnetic flux from the I-coils. In conclusion, these results provide additional confidence that ITER may be capable of RMP ELM suppression in the event of multiple internal coil failures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Opacity calculations for solar mixtures

The OPAS opacity model has been used to calculate the radiative opacity of stellar plasmas in local thermodynamic equilibrium. According to the recent chemical composition revision of the solar photosphere, opacities have been computed for various hydrogen and metallic element mass fractions. Calculations have been performed toward solar interior modeling for temperatures between log[T(K)]=6 and log[T(K)]=7.2, and for electron densities between 10{sup 20} and 10{sup 26} cm{sup −3}. We discuss possible sources of uncertainty in the calculations. We also compare Rosseland opacities to OPAL and OP data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Collisionless plasma transport mechanisms in stochastic open magnetic field lines in tokamaks

The collisionless plasma transport in stochastic magnetic fields connecting to wall boundary has been studied to understand the underlying mechanisms of thermal quench physics in tokamak disruption. In this study, we present a comprehensive understanding of the plasma dynamics in three-dimensional stochastic open magnetic field lines, taking into account the consistent coupling of electron and ion dynamics through the ambipolar electric field. The open stochastic field lines act as a 3D magnetic mirror consisting of the magnetic well and hill regions. It was found that the magnetic hill regions play a critical role in enhancing electron thermal transport. When the plasma collapses to the wall, the 3D ambipolar potential arises in the stochastic layer to maintain the quasi-neutrality of the plasma. The E × B vortices induced by the 3D ambipolar potential mix the plasma across stochastic field lines and enhance the radial transport. Particularly, the E × B mixing between the magnetic well and hill regions provides a collisionless detrapping mechanism that plays a major role in the loss of high-energy trapped electrons. As a result, the electron temperature steadily decreases at a rate of about -0.5 keV ms -1 , comparable to the rate observed in thermal quench experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A kinetic line-driven radiation operator and its application to Gyrokinetics

A velocity dependent, kinetic model for line radiation is developed for continuum kinetic codes. It has been implemented in the full-f gyrokinetic code Gkeyll. The total radiation for a charge state is modeled as an advection in velocity space with a form of $\nabla_v \cdot(v\nu(v)f(v))$, guaranteeing particle conservation. The velocity dependence (in the form of an effective frequency $\nu(v)$) is found through fitting the energy loss of the operator, i.e. the second velocity moment, to the radiation data in the OpenADAS database. Therefore, each individual transition does not need to be evaluated every time step, significantly reducing the computational cost of including line radiation in a kinetic model. The dependence on velocity instead of the usual, temperature, allows the radiation to be computed from non-Maxwellian electron distribution functions: We benchmark the model against a collisional radiative model using isotropic non-Maxwellian distribution functions. A velocity dependent model of radiation can more accurately describe the radiation in the more kinetic regimes expected in reactor-scale devices. The velocity dependence qualitatively captures the quantum mechanical need for a minimum velocity before any radiation occurs.

kinetic↗
DOE OSTI2026

Autoregressive long-horizon prediction of plasma edge dynamics *

Accurate modeling of scrape-off layer (SOL) and divertor-edge dynamics is vital for designing plasma-facing components in fusion devices. High-fidelity edge fluid/neutral codes such as SOLPS-ITER capture SOL physics with high accuracy, but their computational cost limits broad parameter scans and long transient studies. We present transformer-based, autoregressive surrogates for efficient prediction of 2D, time-dependent plasma edge state fields. Trained on SOLPS-ITER spatiotemporal data for the KSTAR tokamak, the surrogates forecast electron temperature, electron density, and radiated power over extended horizons. We evaluate model variants trained with increasing autoregressive horizons (1–100 steps) on short- and long-horizon prediction tasks. Longer-horizon training systematically improves rollout stability and mitigates error accumulation, enabling stable predictions over hundreds to thousands of steps and reproducing key dynamical features such as the motion of high-radiation regions. Measured end-to-end wall-clock times show the surrogate is orders of magnitude faster than SOLPS-ITER, enabling rapid parameter exploration. Prediction accuracy degrades when the surrogate enters physical regimes not represented in the training dataset, motivating future work on data enrichment and physics-informed constraints. Overall, this approach provides a fast, accurate surrogate for computationally intensive plasma edge simulations, supporting rapid scenario exploration, control-oriented studies, and progress toward real-time applications in fusion devices.

autoregressive deep learning↗
DOE OSTI2013

COSMIC RAY MODULATION BEYOND THE HELIOPAUSE: A HYBRID MODELING APPROACH

Results from a newly developed hybrid cosmic ray (CR) modulation model are presented. In this approach, the transport of CRs is computed by incorporating the plasma flow from a magnetohydrodynamic model for the heliospheric environment, resulting in representative CR transport. The model is applied to the modulation of CRs beyond the heliopause (HP) and we show that (1) CR modulation persists beyond the HP, so it is unlikely that the Voyager spacecraft will measure the pristine local interstellar spectra of galactic CRs when crossing the HP. (2) CR modulation in the outer heliosheath could maintain solar-cycle-related changes. (3) The modulation of CRs in the outer heliosheath is primarily determined by the ratio of perpendicular to parallel diffusion, so that the value of the individual diffusion coefficients cannot be determined uniquely using this approach. (4) CRs can efficiently diffuse between the nose and tail regions of the heliosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Implementing a microphysics model in hydrodynamic simulations to study the initial plasma formation in dielectric ablator materials for direct-drive implosions

A microphysics model to describe the photoionization and impact ionization processes in dielectric ablator materials like plastic has been implemented into the 1-D hydrodynamic code LILAC for planar and spherical targets. Currently, the initial plasma formation during the early stages of a laser drive are modeled in an ad hoc manner, until the formation of a critical surface. Implementation of the physics-based models predict higher values of electron density, electron temperature, and pressure than the ad hoc model. Moreover, the numerical predictions are consistent with previous experimental observations of the shinethrough mechanism in plastic ablators. For planar targets, a decompression of the rear end of the target was observed, that is similar to recent experiments. An application of this model is to understand the laser-imprint mechanism that is caused by nonuniform laser irradiation due to single beam speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗