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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Discovering the 3D Structure and Dynamics of the Sun-Interstellar Medium System on a Global Scale Los Alamos LDRD Report

We discover the 3D plasma structure of the Sun-interstellar medium interaction region, called the heliosheath. This is enabled by (1) the completion of a full solar cycle of observations of hydrogen energetic neutral atoms (ENAs) emanating from this region and detected by the Los Alamos-led IBEX-Hi ENA imager on NASA’s Interstellar Boundary Explorer (IBEX) mission, and (2) the development of tools to fully exploit all data from the IBEX-Hi imager, leading to a significant increase in the effective instrument sensitivity. This effort prepares us for optimizing the design and data exploitation of the IMAP-Hi imager, the Los Alamos-led next generation imager that will fly on the Interstellar Mapping and Acceleration Probe (IMAP) mission, currently under production. This project realized its goals as evidenced by two main outcomes: (1) The publication of the first-ever empirical 3D map of the interstellar boundary, and (2) the development and release of a new validated ENA data set that has tripled the amount of data available for scientific analysis of the outer heliosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Magnetospheric energetic ions from the earth's ionosphere

In the decade and a half since the initial discovery that the earth's own ionosphere could at times contribute measurably to the hot plasma in the magnetosphere, significant progress has been made in both current knowledge and understanding of this connection. It is now known that ions of ionospheric origin are found in all major regions of the magnetosphere and at its boundaries. The source region in the ionosphere and the acceleration and transport processes involved in coupling the cold ionospheric plasma to the hot magnetospheric plasma are complex and variable. The large scale morphology of the ionospheric outflow and its distribution throughout the magnetosphere is now understood and progress is being made in the understanding of the fundamental physical processes involved. In this paper attention is given to the large scale morphology and current understanding of the sources for ionospheric ions found in various regions of the magnetosphere and their transport.

Shelley, E. G.↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An experimental platform for investigating astrophysically relevant magnetized plasma jets on the COBRA facility

A new platform has been developed for the 1-MA COBRA generator to investigate the physical processes affecting the formation, collimation, and stability of high-speed outflows in magnetically driven laboratory plasma jets. Such experiments serve as diagnostically accessible surrogates for astrophysical jets under the assumption that the underlying dynamics are scale invariant. In contrast to previous current driven high energy density laboratory jet experiments that use radial/conical wire arrays or foils, the platform described here uses azimuthally symmetric gas-puff injection. This avoids the ablation phase from a solid target, allowing the jets to develop earlier and be driven longer without depleting their mass source and disrupting. A permanent magnet provides an initial poloidal magnetic field, which links the two concentric electrodes and mimics the boundary conditions of a star-accretion disk system. Extended magnetohydrodynamic effects can be assessed using a polarity convolute, which allows for reversal of the electrode bias. The resulting plasma jets exhibit remarkable stability, persisting for hundreds of nanoseconds and achieving aspect ratios ≳30:1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Machine Learning Characterization of Alfvénic and Sub-Alfvénic Chirping and Correlation With Fast-Ion Loss at NSTX

Abrupt large events in the Alfvenic and sub-Alfvenic frequency bands in tokamaks are typically correlated with increased fast ion loss. Here, machine learning is used to speed up the laborious process of characterizing the behaviour of magnetic perturbations from corresponding frequency spectrograms that are typically identified by humans. Analysis allows for comparison between different mode character (such as quiescent, fixedfrequency, chirping, avalanching) and plasma parameters obtained from the TRANSP code such as the ratio of the neutral beam injection (NBI) velocity and the Alfven velocity (v inj. /v A ), the q-profile, and the ratio of the neutral beam beta and the total plasma beta (β beam,i /β). In agreement with previous work by Fredrickson et al., in this work, we find correlation between β beam,i and mode character. In addition, previously unknown correlations are found between moments of the spectrograms and mode character. Character transition from quiescent to non-quiescent behaviour for magnetic fluctuations in the 50 - 200 kHz frequency band is observed along the boundary v Φ ≲ (1/4) (v inj. - 3v A ) where v Φ is the rotation

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The influence of boundary and edge-plasma modeling in computations of axisymmetric vertical displacement

A number of previously published studies apply extended-MHD computations to model macroscopic dynamics of tokamak vertical displacement events (VDEs). The characteristic plasma-wall contact and resistive-wall diffusion imply sensitivity to boundary conditions in a general sense, but which conditions have significant influence depends on how the plasma is modeled. This work examines these dependencies by comparing results of axisymmetric extended-MHD computations with different sets of boundary conditions, plasma thermal-conduction models, and wall resistivity values. The geometry and plasma parameters of the computations represent a modest-sized tokamak. The forced-displacement scenario and computational setup are similar to those used in previous 3D computations. The results show that for a given wall resistivity, the VDE time-scale is most sensitive to variations in the boundary and thermal-conduction parameters that affect electron thermal transport. The electrical conductivity depends on electron temperature, and the dependence on thermal transport stems from its influence on the electrical circuit that includes the open-field halo current. Lastly, conditions that lead to hotter, broader halo regions slow the evolution. Significant sensitivity to the boundary condition on plasma flow-velocity exists when electron thermal conduction is restricted and electron energy loss is convective, which is expected for conditions at the entrance of the magnetic presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

Extended Acceleration in Slot Gaps and Pulsar High-Energy Emission

We revise the physics of primary electron acceleration in the "slot gap" (SG) above the pulsar polar caps (PCs), a regime originally proposed by Arons and Scharlemann (1979) in their electrodynamic model of pulsar PCs. We employ the standard definition of the SG as a pair-free space between the last open field lines and the boundary of the pair plasma column which is expected to develop above the bulk of the PC. The rationale for our revision is that the proper treatment of primary acceleration within the pulsar SGs should take into account the effect of the narrow geometry of the gap on the electrodynamics within the gap and also to include the effect of inertial frame dragging on the particle acceleration. We show that the accelerating electric field within the gap, being significantly boosted by the effect of frame dragging, becomes reduced because of the gap geometry by a factor proportional to the square of the SG width. The combination of the effects of frame dragging and geometrical screening in the gap region naturally gives rise to a regime of extended acceleration, that is not limited to favorably curved field lines as in earlier models, and the possibility of multiple-pair production by curvature photons at very high altitudes, up to several stellar radii. We present our estimates of the characteristic SG thickness across the PC, energetics of primaries accelerated within the gap, high-energy bolometric luminosities emitted from the high altitudes in the gaps, and maximum heating luminosities produced by positrons returning from the elevated pair fronts. The estimated theoretical high-energy luminosities are in good agreement with the corresponding empirical relationships for gamma-ray pulsars. We illustrate the results of our modeling of the pair cascades and gamma-ray emission from the high altitudes in the SG for the Crab pulsar. The combination of the frame-dragging field and high-altitude SG emission enables both acceleration at the smaller inclination angles and a larger emission beam, both necessary to produce widely-spaced double-peaked profiles.

White, Nicholas E.↗
DOE OSTI2023

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1985

Investigation of the dayside boundary region of the magnetosphere

The region near the dayside boundary of the magnetosphere is of particular interest mainly because it is the most likely site for solar wind energy transfer. Exactly where and how this transfer occurs is one of the basic problems of magnetospheric physics. The global nature of the energy transfer is addressed in this research. In order to test whether energy transfer affects the large scale plasma flow in the magnetosheath, global properties of the magnetosheath were analyzed and compared to models with no energy transfer. Also, global models of merging sites on the magnetopause were modeled with the aid of the magnetosheath models. Since energy transfer and intrinsic magnetosheath properties depend upon the orientation of the interplanetary magnetic field (IMF), it was essential to have an IMF monitor for these studies. In the course of using ISEE 3 in this role, statistical studies were performed to test the accuracy of using data from an orbit so far upstream. Global properties of the magnetosheath energetic ion population, magnetic field strength and orientation, and flow pattern, determined with ISEE data are summarized.

Crooker, N.↗
DOE OSTI2021

Using extended MHD to explore lasers as a trigger for x-pinches

X-pinches have been shown to be a source of extremely intense x-ray emissions useful for diagnosing plasma dynamics and imaging biological objects. The most striking feature of an x-pinch is the hotspot, the point source from where all the x rays come from. Unfortunately, the exact timing and location of the hotspot are still unpredictable. Since an x-pinch hotspot forms from instabilities (like an m = 0 mode), we will computationally explore whether we can use a high-power laser to control the timing and location of these instabilities. Our goal is to reduce the temporal and spatial jitter associated with the x-ray burst. Using an extended MHD (XMHD) code, we explore the non-relativistic instability generation using a current profile of a 250 kA linear transformer driver system and laser characteristics of SLAC's Matter in Extreme Conditions laboratory. Our results include both laser-penetration results using a boundary-defined electromagnetic (EM) wave and instability results from a power-deposition method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Combined Normal and Disruption, Electromagnetic Transient, Thermal, and Structural Analysis of COMPASS Upgrade

We report that COMPASS Upgrade is a replacement for the COMPASS device at the Institute of Plasma Physics (IPP) of the Czech Academy of Sciences. It is a cryogenic copper machine with a major radius of 0.89 m, with 5 T at the plasma centerline and 2-MA plasma current. It is intended to develop ITER relevant plasma densities and high power fluxes in the divertor region. The entire vacuum vessel is planned to reach 500 °C. The COMPASS-U TF coil is a bitter plate-like design with radial conductor plates that extend the width of the TF inner leg. Currents redistribute radially based on resistive and inductive effects. This necessitates an electromagnetic (EM) transient analysis of the coil current distribution. For copper coils like COMPASS Upgrade, the solution must be coupled with a thermal solution to properly include the temperature-dependent resistive effects. The EM transient solution includes the TF magnetic loads, and a stress pass on the results produces in-plane loads. The EM solution chosen requires inclusion of all the conducting and nonconducting regions surrounding the TF coil. With the TF current distribution solved, the addition of the poloidal coils and resulting background fields allows the determination of the out-of-plane loads. Integration of the Joule heating produced from the EM solution produces temperature distributions throughout the shot, which can be read in to a structural pass to include thermal stresses in the evaluation. This is true of both TF and PF temperature prediction with appropriate packing fractions. Currents in the PF coils, with computed background fields, produce Lorentz loads. Time transients of the PF currents produce startup eddy currents in the structures included in the model; inclusion of the vessel and other passive structures allows these to be included in a stress pass. The EM transient solution includes all the elements needed for a disruption analysis with prescribed motions and quench of the plasma. From a model originally intended to simulate current diffusion in the TF coil, the analysis can be refined to be used to predict the performance of most of the components of the tokamak or provide boundary conditions for more detailed submodels. In this article, the use of the current diffusion model for bounding assessments of the COMPASS Upgrade TF, PF, vessel, external structure, and disruption analyses will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Validation of a comprehensive first-principles-based framework for predicting the performance of future stellarators

This paper presents the validation of the GENE-KNOSOS-Tango framework for recovering both the steady-state plasma profiles in the considered radial domain and selected turbulence trends in a stellarator. This framework couples the gyrokinetic turbulence code GENE, the neoclassical transport code KNOSOS, and the transport solver Tango in a multi-timescale simulation feedback loop. Ion-scale kinetic-electron and electron-scale adiabatic-ion flux-tube simulations were performed to evolve the density and temperature profiles for four OP1.2b W7-X scenarios. The simulated density and temperature profiles showed good agreement with the experimental data using a reasonable set of boundary conditions. Equally important was the reproduction of observed trends for several turbulence properties, such as density fluctuations and turbulent heat diffusivities. Key effects were also touched upon, such as electron-scale turbulence and the neoclassical radial electric field shear. The validation of the GENE-KNOSOS-Tango framework enables credible predictions of physical phenomena in stellarators and reactor performance based on a given set of edge parameters.

Gyrokinetic simulations↗
DOE OSTI2023

Boundary, SOL, and Divertor Physics Studies on TCV (Final Report for Research Funded by DoE Award DE-SC0020327)

Research supported under this grant was organized around two diagnostic systems installed on the TCV tokamak located at EPFL in Lausanne, Switzerland: Multi-Spectral Imaging (MSI) and Gas-Puff Imaging (GPI). MSI monitors 2-d time-history images of up to 10 spectral lines, in the visible, using identical views for all wavelengths. The primary application of MSI is to infer 2-d time histories of electron temperature and density in the Scrape-Off-Layer (SOL) of diverted plasmas, utilizing ratios of emission from multiple neutral helium lines. The atomic models are validated against 1-d Thomson scattering measurements. SOL density and temperature maps are analyzed to infer divertor dynamics, including particle and power flows, for a variety of divertor configurations, and specifically as flux expansion along the long-legged outboard divertor is systematically varied. GPI is used to measure details of edge fluctuations in the SOL, particularly as a function of plasma triangularity (δ), including both positive and negative δ. A threshold is found, whereby for δ<-0.3, edge turbulence in the far SOL is completely suppressed. Dynamics of the filamentary turbulence was also analyzed using advance machine learning techniques. Research in support of four PhD theses was carried out within the scope of this Award, including 2 at MIT and 2 at EPFL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Electron Cyclotron Emission–Based Separatrix Identification in ITER with OMFIT Synthetic Modeling

Accurate determination of the separatrix location is essential for understanding edge plasma behavior and optimizing confinement in tokamaks, especially in next-generation devices such as ITER. In this study, a synthetic microwave diagnostics module was developed and implemented in the OMFIT framework to assess the feasibility of an electron cyclotron emission–based separatrix detection method in ITER plasmas. Simulations were carried out using ITER H-mode equilibrium scenarios with different plasma density profiles and different pedestal widths. Here, the results show that the electron emission temperature profiles consistently exhibit an inversion pattern near the edge, with a well-defined minimum point that could serve as a proxy for the separatrix location. However, unlike in the DIII-D, the minimum point in ITER is systematically offset by approximately 2 cm into the scrape-off layer, independent of density or pedestal width, which is within the radial resolution range (2 to 5 cm) determined by the 500-MHz channel spacing. While the method does not provide the exact separatrix location, it offers a reliable indicator of the boundary region and has potential applications for real-time boundary monitoring in ITER and other future fusion devices.

Electron cyclotron emission↗
DOE OSTI2023

Fluctuations and Self-Organization in Plasma Boundary Layers (Final Report)

Final technical report on an experimental study of the boundary of an ionized gas (plasma) where it interacts with solid surfaces. The plasma self-organizes into a sheath which experiences plasma flow and is naturally unstable to various fluctuations which can be observed with advanced laser-spectroscopy techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Plasma-assisted surface modification of biopolymers

The interaction of non-thermal plasmas with surfaces is an important problem of interest in various applications. This work focuses on two main aspects of the problem - simulation and experiments. The development of a low-temperature plasma modeling is described along with experiments focusing on the interaction of a low-temperature plasma with liquid water. The simulation framework is capable of handling arbitrary geometries, chemistry and a fully-coupled boundary condition for plasma-dielectric interface. Representative examples for the verification and validation of the framework are provided. The experimental work deals with the study of low-temperature plasmas with water with the nitrate injection rate compared for three different configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Verification and validation of the open-source plasma fluid code: Zapdos

Zapdos is an open-source finite element plasma fluid solver based on the MOOSE multiphysics framework. This paper outlines Zapdos verification, benchmarking, and validation efforts for 1D and 2D RF capacitively coupled plasma discharge models for mid-range pressures (0.1 - 1 Torr). The verification process involved using the method of manufactured solutions to assess Zapdos spatial and temporal error convergence. L2 errors ranged from 10 -2 to 10 -4 , while the convergence’s slope were in agreement with the predicted slopes for the tested variable and time integration orders. The benchmarking process involved comparisons to previously results from the validated finite element code, LSODI. These works included 1D and 2D simulations for a range of plasma parameters (densities, temperatures, voltage, etc.). For the 1D cases, Zapdos and LSODI results were in very good agreement. In the 2D cases, variable behaviors matched, with slight discrepancies in peak values. The validation process involved comparisons to experimental works including electron density measurements by microwave interferometry and metastable density measurements by planar laser-induced fluorescence imaging. Results shown reasonable agreement at higher pressure, with results starting to diverge at low pressures. Furthermore, probable causes for this diverges are the limitation of the fluid assumption for plasmas at low pressure, or the need for more robust boundary conditions. Overall, Zapdos shown reasonable results for the verification, benchmarking, and validation efforts, and Zapdos can be downloaded at https://github.com/shannon-lab/zapdos.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗