Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma boundary physics”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

INGRID: An interactive grid generator for 2D edge plasma modeling

A fusion boundary-plasma domain is defined by axisymmetric magnetic surfaces where the geometry is often complicated by the presence of one or more X-points; and modeling boundary plasmas usually relies on computational grids that account for the magnetic field geometry. The new grid generator INGRID (Interactive Grid Generator) presented here is a Python-based code for calculating grids for fusion boundary plasma modeling, for a variety of configurations with one or two X-points in the domain. INGRID first performs partitioning over the domain consisting of a small number of patches conforming to the magnetic field and wall geometry; then it generates a subgrid on each of the patches and joins them into a global grid. This domain partitioning strategy makes possible a uniform treatment of various configurations with one or two X-points in the domain. This includes single-null, double-null, and other configurations with two X-points in the domain. The INGRID design allows generating grids either interactively, via a parameterfile driven GUI, or using a non-interactive script-controlled workflow. Results of testing demonstrate that INGRID is a flexible, robust, and user-friendly grid-generation tool for fusion boundary-plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1979

Numerical methods for large-scale, time-dependent partial differential equations

A survey of numerical methods for time dependent partial differential equations is presented. The emphasis is on practical applications to large scale problems. A discussion of new developments in high order methods and moving grids is given. The importance of boundary conditions is stressed for both internal and external flows. A description of implicit methods is presented including generalizations to multidimensions. Shocks, aerodynamics, meteorology, plasma physics and combustion applications are also briefly described.

Turkel, E.↗
DOE OSTI2025

Q2 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q2: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $\circ$ Run CESOL with BOUT++/Hermes-3 and immersed boundary condition to directly map to wall: • Run BOUT++/Hermes-3 through the IPS workflow to find radial particle and energy diffusivities to match either the Eich or the physics-based scaling of the SOL heat flux width, and • Expand source of first wall heat flux to include charged particles, neutrals, and radiation from the core+edge. 2. Generate medium fidelity parametrized CAD: $\circ$ Develop the TRACER tool to read an existing CAD, regenerate the geometry based on vertex location and connectivity information, define vertex translation and parameters needed for scaling the CAD, and $\circ$ Utilize the FreeGS code to determine CAT PF coil placement, including minimizing the number of coils, coil current, and electromechanical stresses. 3. Utilize plasma loading for engineering analysis: $\circ$ Couple the plasma loading to input for OpenFOAM and demonstrate initial test of thermal analysis of CAT first wall loading with typical DCLL blanket component cooling boundary conditions. 4. Demonstrate nuclear analysis: $\circ$ Apply initial analysis of tritium transport in DCLL blanket by evaluating spatially resolved tritium generation rates, tritium diffusion and convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Development of Langmuir and Spatially Resolved OES Diagnostics for Arc Jet Shock and Boundary Layer Measurements Relevant to Atmospheric Re-Entry Environments

The second-generation miniature Arc Jet Research Chamber (mARC) ground testing facility at NASA Ames requires physical characterization of its supersonic plasma flow to verify numerical models and material sample testing. Two types of plasma diagnostics will be utilized: Langmuir probes and optical emission spectroscopy (OES). Several single and triple Langmuir probes have been constructed and will be inserted parallel to the flow via a sweep arm to measure post-shock electron temperature and number density along radial profiles at different axial distances. Due to the high heat flux (200-3000 W/cm2) sweep durations are necessarily short (< 0.1 s). Measuring these properties near the probe surface provide characterization of the post-shock environment. This is an important metric for validating coupled CFD-material codes, and will provide insight into longstanding thermocouple anomalies attributed to charging of heatshield surfaces. An OES system has also been assembled to obtain spatially resolved density and temperature for various flow species. The system consists of a convex lens that forms an image of the plasma onto a 16-channel linear fiber array, which will simultaneously capture the spectra from sixteen different locations. The light emission from the plasma is then focused by a spherical mirror into a spectrometer and recorded with a camera. These measurements non-intrusively provide a new spatially resolved flow characterization capability. This will initially be used to resolve issues with shock radiation interference in pyrometer measurements of material surface temperature. The design of each diagnostic and any initial results will be discussed.

Ethan Main-Loam Leong↗
DOE OSTI2020

Mechanisms of energetic-particle transport in magnetically confined plasmas

Super-thermal ions and electrons occur in both space and fusion plasmas. Because these energetic particles (EP) have large velocities, EP orbits necessarily deviate substantially from magnetic surfaces. Orbits are described by conserved constants of motion that define topological boundaries for different orbit types. Electric and magnetic field perturbations produced by instabilities can disrupt particle orbits, causing the constants of motion to change. The statistics of the “kicks” associated with these perturbations determines the resulting cross field transport. A unifying theme of this tutorial is the importance of the perturbation’s phase at the particle’s position Θ = k • r – ωt, where k and ω are the wavevector and frequency of the perturbation, r is the EP position, and t is the time. A distinction is made between field perturbations that resonate with an aspect of the orbital motion and those that do not. Resonance occurs when the wave phase returns to its initial value in an integer multiple of an orbital period. Convective transport occurs when resonant particles experience an unvarying wave phase. Alternatively, multiple wave-particle resonances usually decorrelate the phase, resulting in diffusive transport. Large orbits increase the number of important resonances and can cause chaotic orbits even for relatively small amplitude waves. In contrast, in the case of non-resonant perturbations, orbital phase averaging reduces transport. Here, large field perturbations introduce additional effects, including nonlinear resonances at fractional values of the orbital motion. In summary, large orbits are a blessing and a curse: For non-resonant modes, orbit-averaging reduces transport but, for resonant transport, large orbits facilitate jumps across topological boundaries and enhance the number of important resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Plasma-droplet interaction study to assess transport limitations and the role of • OH, O • ,H • ,O 2 (a 1 Δ g ),O 3 , He(2 3 S) and Ar(1s 5 ) in formate decomposition

Plasmas interacting with liquid microdroplets are gaining momentum due to their ability to significantly enhance the reactivity transfer from the gas phase plasma to the liquid. This is, for example, critically important for efficiently decomposing organic pollutants in water. In this contribution, the role of •OH as well as non-•OH-driven chemistry initiated by the activation of small water microdroplets in a controlled environment by diffuse RF glow discharge in He with different gas admixtures (Ar, O 2 and humidified He) at atmospheric pressure is quantified. The effect of short-lived radicals such as O• and H• atoms, singlet delta oxygen (O 2 (a 1 Δ g )), O 3 and metastable atoms of He and Ar, besides •OH radicals, on the decomposition of formate dissolved in droplets was analyzed using detailed plasma diagnostics, droplet characterization and ex situ chemical analysis of the treated droplets. The formate decomposition increased with increasing droplet residence time in the plasma, with ~70% decomposition occurring within ~15 ms of the plasma treatment time. The formate oxidation in the droplets is shown to be limited by the gas phase •OH flux at lower H 2 O concentrations with a significant enhancement in the formate decomposition at the lowest water concentration, attributed to e – /ion-induced reactions. However, the oxidation is diffusion limited in the liquid phase at higher gaseous •OH concentrations. Here, the formate decomposition in He/O 2 plasma was similar, although with an order of magnitude higher O• radical density than the •OH density in the corresponding He/H 2 O plasma. Using a one-dimensional reaction–diffusion model, we showed that O 2 (a 1 Δ g ) and O 3 did not play a significant role and the decomposition was due to O•, and possibly •OH generated in the vapor containing droplet-plasma boundary layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Validation of the imaging neutral particle analyzer in nearly MHD quiescent plasmas using injected beam ions on DIII-D

The imaging neutral particle analyzer (INPA) is a scintillator-based diagnostic that provides energy and radially resolved measurements of confined fast ions. To verify operation, understanding, and modeling of the diagnostic performance, distributions of deuterium beam ions in low density, nearly MHD-quiescent plasmas are measured on the DIII-D tokamak using the INPA. Here, the distribution is altered by changing the pitch-angle scattering rate with electron cyclotron heating (ECH). Experimental images are mostly reproduced with simulation of core fast ions charge exchanging with injected beam and edge cold neutrals in classical cases. Image comparisons between the experiment and simulation show agreement with less than 25% discrepancies near the beam injection energy of 50 keV. Other strong passive signals, observed in a localized region on the image, are suggested to come from charge exchange between cold neutrals and fast ions around the plasma boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1979

Interplanetary gas. XXIV - Are cometary plasma tail disconnections caused by sector boundary crossings or by encounters with high-speed streams

The paper discusses and compares the original sector boundary model (Niedner and Brandt) and the alternative high-speed stream model (Ip and Mendis) suggested for explaining cometary plasma tail disconnection events (DE) within a largely observational framework not dependent on detailed plasma physics. Among the findings are: (1) the strong solar cycle phase dependence of the amplitude of the Rosenberg-Coleman effect (Svalgaard and Wilcox) yields inferred maximum latitudes of all of the DEs in the original survey to a one-time Rosenberg-Coleman effect measurement made by Pioneer 11, which indicated a disappearance of sectors above 16 deg latitude, and (2) approximately 70% of the post 1926 DEs in the expanded survey show a close association with corotated shorter-term polarity reversals. This result is difficult to reconcile with the predictions of the stream model, in which DEs should correlate more strongly with streams than with sector boundaries. It is concluded that the sector boundary model better describes the disconnection phenomenon and that the use of DEs as unique sector boundary markers is presently justified.

Niedner, M. B., Jr.↗
DOE OSTI2022

Burning plasma achieved in inertial fusion

Obtaining a burning plasma is a critical step towards self-sustaining fusion energy. A burning plasma is one in which the fusion reactions themselves are the primary source of heating in the plasma, which is necessary to sustain and propagate the burn, enabling high energy gain. After decades of fusion research, here we achieve a burning-plasma state in the laboratory. These experiments were conducted at the US National Ignition Facility, a laser facility delivering up to 1.9 megajoules of energy in pulses with peak powers up to 500 terawatts. We use the lasers to generate X-rays in a radiation cavity to indirectly drive a fuel-containing capsule via the X-ray ablation pressure, which results in the implosion process compressing and heating the fuel via mechanical work. The burning-plasma state was created using a strategy to increase the spatial scale of the capsule through two different implosion concepts. These experiments show fusion self-heating in excess of the mechanical work injected into the implosions, satisfying several burning-plasma metrics. Additionally, we describe a subset of experiments that appear to have crossed the static self-heating boundary, where fusion heating surpasses the energy losses from radiation and conduction. These results provide an opportunity to study α-particle-dominated plasmas and burning-plasma physics in the laboratory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2013

Radial Profiles of Plasma Electron Characteristics in a Low-Power Arcjet

Since 1960, the Arc Jet Complex at NASA Ames Research Center has been a source of long-duration, high-enthalpy flow for materials testing with application to the thermal protection of aerospace vehicle components. From their inception the facilities have played an integral role supporting many of NASA's space flight programs and numerous DoD projects. In recent years advancements in computational fluid dynamics (CFO) have made the resultant models a valuable tool for assessing and predicting performance, however, the inherent limitation of models to compensate for a dissociated, transitionally high temperature nonequilibrated flowfield have made further measurements necessary. The use of electrostatic probe diagnostics within similarly harsh plasma environments in previous studies have been met with much success. In this study, the use of a single Langmuir probe was implemented to characterize the plasma parameters of interest as they vary radially within a large volume of the plume. Classical Langmuir probe theory was applied to achieve first order estimates of the heavy particle temperature, the ratio T(sub i)/T(sub e), and the ionization fraction. As expected, both the electron temperature and electron density measurements show a dependence on radial distance from the plume centerline, with electron density profiles showing the largest dependence. This paper aims to validate and strengthen modeling work conducted at the NASA Ames Research Center by measuring the critical plasma parameters (electron temperature, electron density, and plasma potential) within the arc plume of a subscale arc jet. These parameters are intended to give physical insight into the flow characteristics while providing the necessary boundary conditions to validate full scale simulations.

Codron, Douglas↗
DOE OSTI2021

Constrained stellarator coil curvature optimization with FOCUS

Finding less complicated coils that have adequately low field errors is a crucial step in stellarator development. One coil metric that is of high importance is the maximum curvature of the coil centreline, or coil single filament. Conductors cannot be bent below some threshold minimum radius of curvature. High coil curvatures can cause strains to exceed acceptable levels, especially in superconducting coils. We investigate three ways to optimize coil curvature and find that applying penalty functions to the coil curvature solves for coils that have a constrained maximum curvature and low field error. Penalty functions are implemented in FOCUS and coil solutions optimized for an HSX-like ‘plasma boundary’ are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

X-ray induced grain boundary formation and grain rotation in Bi 2 Se 3

Optimizing grain boundary characteristics in polycrystalline materials can improve their properties. Many processing methods have been developed for grain boundary manipulation, including the use of intense radiation in certain applications. Here, in this work, we used X-ray free electron laser pulses to irradiate single-crystalline bismuth selenide (Bi 2 Se 3 ) and observed grain boundary formation and subsequent grain rotation in response to the X-ray radiation. Our observations with simultaneous transmission X-ray microscopy and X-ray diffraction demonstrate how intense X-ray radiation can rapidly change size and texture of grains.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Energy-dependent implementation of secondary electron emission models in continuum kinetic sheath simulations

The plasma-material interactions present in multiple fusion and propulsion concepts between the flow of plasma through a channel and a material wall drive the emission of secondary electrons. This emission is capable of altering the fundamental structure of the sheath region, significantly changing the expected particle fluxes to the wall. The emission spectrum is separated into two major energy regimes, a peak of elastically backscattered primary electrons at the incoming energy, and cold secondary electrons inelastically emitted directly from the material. The ability of continuum kinetic simulations to accurately represent the secondary electron emission is limited by relevant models being formulated in terms of monoenergetic particle interactions which cannot be applied directly to the discrete distribution function. As a result, rigorous implementation of energy-dependent physics is often neglected in favor of simplified, constant models. We present here a novel implementation of semi-empirical models in the boundary of continuum kinetic simulations which allows the full range of this emission to be accurately captured in physically-relevant regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Boundary conditions which lead to excitation of instabilities in plasma simulations

Two examples of two-dimensional electrostatic particle-code simulations are shown in which one exhibits characteristics of a stable plasma while the other exhibits unstable, long wavelength plasma oscillations. The only difference between the two simulations is a change in the boundary condition on the electrostatic potential. An energy theorem is derived which shows that the rate of change of field and particle energy within a closed volume is related to a surface integral involving the electrostatic potential and the normal component of the electric current. An analytic theory is developed for a one-dimensional plasma to show how boundary effects can excite spurious plasma instabilities. The theory is tested with a series of one-dimensional plasma simulations. Finally, practical considerations on means of avoiding the non-physical instabilities in simulation plasmas are given.

Swift, D. W.↗