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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRSDate not supplied

Thermomechanical Relaxation of Vibrationally Excited NO in a Hypersonic Turbulent Boundary Layer

The effects of thermal nonequilibrium are studied using a DC glow discharge plasma to selectively modify the vibrational temperature of an otherwise cold flow. These experiments were performed at 𝑴=5.7and 𝑹𝒆=6×10 6 /m in the Actively Controlled Expansion tunnel at Texas A&M University on a canonical 2.75° half-angle wedge test article. A∼47 W plasma was generated across spanwise electrodes, and removable trips produced either a laminar or turbulent boundary layer. NO planar laser-induced fluorescence (PLIF) quantified the extent of the NO 𝑿 2 𝚷 (𝒗=1) vibrational perturbation due to the plasma and its subsequent interplay with the boundary layer, the first such data of which the authors are aware. NOPLIF measurements were taken at three locations along the test article, and quasi-temperature measurements are reported at each. Despite not quantifiably perturbing the boundary layer, the plasma was shown to alter the NO vibrational temperature profile therein. Independent analysis of the thermodynamics and fluid mechanics of the system suggests each set of physics could explain the results observed near and far from the wall, respectively.

Casey J. Broslawski↗
NASA NTRSDate not supplied

Thermomechanical Relaxation of Vibrationally Excited NO in a Hypersonic Turbulent Boundary Layer

The effects of thermal nonequilibrium are studied using a DC glow discharge plasma to selectively modify the vibrational temperature of an otherwise cold flow. These experiments were performed at 𝑴=5.7and 𝑹𝒆=6×10 6 /m in the Actively Controlled Expansion tunnel at Texas A&M University on a canonical 2.75° half-angle wedge test article. A∼47 W plasma was generated across spanwise electrodes, and removable trips produced either a laminar or turbulent boundary layer. NO planar laser-induced fluorescence (PLIF) quantified the extent of the NO 𝑿 2 𝚷 (𝒗=1) vibrational perturbation due to the plasma and its subsequent interplay with the boundary layer, the first such data of which the authors are aware. NOPLIF measurements were taken at three locations along the test article, and quasi-temperature measurements are reported at each. Despite not quantifiably perturbing the boundary layer, the plasma was shown to alter the NO vibrational temperature profile therein. Independent analysis of the thermodynamics and fluid mechanics of the system suggests each set of physics could explain the results observed near and far from the wall, respectively.

Casey J. Broslawski↗
DOE OSTI2025

Multiscale modeling and laser diagnostics to reveal non-equilibrium reaction chemistry at a plasma-liquid interface

Low-temperature, atmospheric-pressure plasmas in contact with liquid are at the core of a wide range of applications including water treatment, medicine, materials synthesis, and chemical transformation. In general, plasma-liquid processes are scientifically compelling because reactivity can be produced without a catalyst, in relatively inert molecules, such as air or nitrogen and liquid water, in their native states at ambient conditions. However, reactions at a plasma-liquid interface are extremely complex, occurring at a multiphase, gas-liquid boundary where excited or dissociate gaseous species dissolve and react with solution-phase species, and unique reaction pathways are induced by non-equilibrium chemistry. In particular, detailed knowledge of the physical and chemical processes, including what species are produced in the gas phase and how these species are subsequently transported across and react near the interfacial region, remains largely unanswered. In this project, modeling of the non-equilibrium reaction chemistry at the interface of low-temperature, atmospheric-pressure plasmas and liquid water is developed, supported by advanced laser diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Microphysics of the auroral acceleration region of other regions of the magnetosphere

Electric field and plasma measurements pertinent to the microphysics of the auroral acceleration region are reviewed. Observations in other regions of the magnetosphere, which appear to involve similar physical processes, are summarized in order to emphasize the requirement on Cluster instrumentation for wide dynamic range, high time resolution, time domain field, and particle measurements. It is suggested that the microphysics of current driven instabilities in the auroral acceleration region is duplicated in the high latitude boundary of the plasma sheet and that the macroscale structure of the tail is determined by this microphysics.

Mozer, F. S.↗
DOE OSTI2021

Effects of divertor electrical drifts on particle distribution and detachment near the divertor target plate in DIII-D

Strong impacts of drifts on the divertor plasma in–out asymmetry and detachment are demonstrated in DIII-D with an open divertor configuration. For forward toroidal field, B T , i.e., with the ion B × ∇B drift toward the divertor, the particle flux to the inner divertor, as represented by the Langmuir probe measured ion saturation current (J sat ), exhibits a double peak structure, with electron temperature, lower at the inner target. Reversing the B T direction reverses both the radial and poloidal E × B flows, leading to a broad particle flux profile in the outboard scrape-off layer (SOL) with a similar double-peak structure to that observed at the inner target with forward B T . The correlation of a double peak structure with divertor temperature profiles confirms physical coupling between the drift flow and sheath boundary condition and their strong impact on divertor profiles. In addition, under reversed B T conditions, increasing the density flattens the target temperature profile. However, J sat remains high away from the strike point, rendering it difficult to achieve an “effective” detached plasma, i.e., with effective reduction in both peak heat flux and peak temperature (in the far SOL). In contrast, divertor detachment with a cold and flat temperature profile can be achieved at both target plates with the forward B T .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Nonlinear bubble competition of the multimode ablative Rayleigh–Taylor instability and applications to inertial confinement fusion

The self-similar nonlinear evolution of the multimode ablative Rayleigh–Taylor instability (RTI) and the ablation-generated vorticity effect are studied for a range of initial conditions. We show that, unlike classical RTI, the nonlinear multimode bubble-front evolution remains in the bubble competition regime due to ablation-generated vorticity, which accelerates the bubbles, thereby preventing a transition into the bubble-merger regime. We develop an analytical bubble competition model to describe the linear and nonlinear stages of ablative RTI. Here, we show that vorticity inside the multimode bubbles is most significant at small scales with large initial perturbation. Since these small scales persist in the bubble competition regime, the self-similar growth coefficient ab can be enhanced by up to 30% relative to ablative bubble competition without vorticity effects. We use the ablative bubble competition model to explain the hydrodynamic stability boundary observed in OMEGA low-adiabat implosion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electromagnetic Waves in Cold Plasmas

A collisionless fluid model of a plasma is presented and solved to first order in the cold, high frequency limit. From this model, a homogeneous system of linear equations is obtained for a harmonic electric field perturbation from which the Appleton-Hartree equation is derived. Moreover, the polarization states of electromagnetic waves in cold, collisionless plasmas are examined. Complete sets of equations governing the reflection and transmission of electromagnetic waves in both non magnetized and magnetized plasmas are obtained. In solving the boundary problem, the Booker quartic is derived and the nature of its roots discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Some contributions to knowledge of the magnetospheric plasma by ISEE investigators

The ISEE project has made substantial contributions to the knowledge of the magnetosphere during the period of the IMS, especially in the discipline of Space Plasma Physics. Results obtained during approximately the first two years of the operation of ISEE-1 and -2, and touches on relevant results of ISEE-3 are reviewed. The ability to control the separation between ISEE-1 and -2, which are in nearly identical orbits, has permitted study of the motion and structure of the bow shock and magnetopause, the boundary layers, and the plasma sheet. Much evidence was obtained favoring the existence of reconnection and its relevance to the transfer of magnetic flux from the frontside to the rear of the magnetosphere, although not everyone agrees that it is the only important process. The presence of both reflected and accelerated particles was shown to lead to the development of a foreshock region between the bow shock and the interplanetary magnetic field line tangential to it.

Ogilvie, K. W.↗
DOE OSTI2024

Grad–Shafranov equilibria via data-free physics informed neural networks

A large number of magnetohydrodynamic (MHD) equilibrium calculations are often required for uncertainty quantification, optimization, and real-time diagnostic information, making MHD equilibrium codes vital to the field of plasma physics. In this paper, we explore a method for solving the Grad–Shafranov equation by using physics-informed neural networks (PINNs). For PINNs, we optimize neural networks by directly minimizing the residual of the partial differential equation as a loss function. We show that PINNs can accurately and effectively solve the Grad–Shafranov equation with several different boundary conditions, making it more flexible than traditional solvers. This method is flexible as it does not require any mesh and basis choice, thereby streamlining the computational process. We also explore the parameter space by varying the size of the model, the learning rate, and boundary conditions to map various tradeoffs such as between reconstruction error and computational speed. Additionally, we introduce a parameterized PINN framework, expanding the input space to include variables such as pressure, aspect ratio, elongation, and triangularity in order to handle a broader range of plasma scenarios within a single network. Parameterized PINNs could be used in future work to solve inverse problems such as shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Modeling rf sheath formation in turbulent tokamak boundary plasma

During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Observation of Nonaxisymmetric Standard Magnetorotational Instability Induced by a Free-Shear Layer

The standard magnetorotational instability (SMRI) with a magnetic field component parallel to the rotation axis is widely believed to be responsible for the fast accretion in astronomical disks. In conventional base flows with a Keplerian profile or an ideal Couette profile, most studies focus on axisymmetric SMRI, since excitation of nonaxisymmetric SMRI in such flows requires a magnetic Reynolds number (Rm) more than an order of magnitude larger. Here, we report that, in a magnetized Taylor-Couette flow, nonaxisymmetric SMRI with an azimuthal mode number m = 1 can be triggered by a free-shear layer in the base flow at Rm ≳1, the same threshold as for axisymmetric SMRI. Global linear analysis reveals that the free-shear layer reduces the required Rm, possibly by introducing an extremum in the vorticity of the base flow. Nonlinear simulations validate the results from linear analysis and confirm that a novel instability recently discovered experimentally [Wang et al., Nat. Commun. 13, 4679 (2022)] is the nonaxisymmetric m =1 SMRI. Further, our finding has astronomical implications as free-shear layers are ubiquitous in celestial systems, such as the disk-star boundary layer, the solar tachocline, and the edge of planet-opened gaps in protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Final report of activities for the LDRD-express project #223796 titled: “Fluid models of charged species transport: numerical methods with mathematically guaranteed properties”, PI: Ignacio Tomas, Co-PI: John Shadid

This report summarizes the findings and outcomes of the LDRD-express project with title “Fluid models of charged species transport: numerical methods with mathematically guaranteed properties”. The primary motivation of this project was the computational/mathematical exploration of the ideas advanced aiming to improve the state-of-the-art on numerical methods for the one-fluid Euler-Poisson models and gain some understanding on the Euler-Maxwell model. Euler-Poisson and Euler-Maxwell, by themselves are not the most technically relevant PDE plasma-models. However, both of them are elementary building blocks of PDE-models used in actual technical applications and include most (if not all) of their mathematical difficulties. Outside the classical ideal MHD models, rigorous mathematical and numerical understanding of one-fluid models is still a quite undeveloped research area, and the treatment/understanding of boundary conditions is minimal (borderline non-existent) at this point in time. This report focuses primarily on bulk-behaviour of Euler-Poisson’s model, touching boundary conditions only tangentially.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Vertical instability forecasting and controllability assessment of multi-device tokamak plasmas in DECAF with data-driven optimization

Abstract Reliable vertical position control will be an essential element of any future tokamak-based fusion power plant in order to reduce disruptions and maximize performance. We investigate methods to improve vertical controllability boundary determination in plasma operational space and demonstrate a data-driven approach based on direct pseudoinversion of operational space data that is rigorously quantitative, applicable in real-time plasma control systems, and physically intuitive to interpret. Applied to historical shot data from entire run campaigns on the MAST-U, KSTAR, and NSTX tokamaks, this approach, implemented in DECAF, improves vertical displacement event identification accuracy to 98.9%–100%. Further, we explore the application of a physics-based vertical stability metric as an early warning forecaster for vertical displacement events. The development of a linear surrogate model for the plasma current density profile, with a coefficient of determination of 0.992 on the training dataset, enables potential employment of this forecaster in real-time. The application of this approach on historical data from the MAST-U MU02 campaign yields a forecaster with 62.6% accuracy, indicating promise for this method when further refined and potentially coupled with other stability metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Cross-tail current - Resonant orbits

A technique to generate self-consistent 1D current sheets is described. Groups of monoenergetic protons were followed in a modified Harris magnetic field. This sample current sheet is characterized by resonant quasi-adiabatic orbits. The magnetic moment of a quasi-adiabatic ion which is injected from outside a current sheet changes substantially during the orbit but returns to almost its initial value by the time the ion leaves. Several ion and electron groups were combined to produce a plasma sheet in which the charged particles carry the currents needed to generate the magnetic field in which the orbits were traced. An electric field also is required to maintain charge neutrality. Three distinct orbit types, one involving untrapped ions and two composed of trapped ions, were identified. Limitations associated with the use of a 1D model also were investigated; it can provide a good physical picture of an important component of the cross-tail current, but cannot adequately describe any region of the magnetotail in which the principal current sheet is separated from the plasma sheet boundary layer by a nearly isotropic outer position of the central plasma sheet.

Kaufmann, Richard L.↗
NASA NTRS1994

Self-consistent steady state model of the low-latitude boundary layer

Results are presented from a steady state numerical model of the low-latitude boundary layer (LLBL) on closed field lines and its coupling to the dayside auroral ionosphere. In the model the boundary layer approximation is used, the result being that inertial forces are taken into account only in the main flow direction (-x) where they are balanced by pressure forces, j x B forces, and viscous forces. Motion in the transverse directions (y and z) is treated kinematically, the force balances in these two directions being purely static. Computationally, the model is two dimensional, describing the motion of plasma and frozen-in magnetic field in the equatorial (xy) plane but allowing for lowest-order polynomial variation of some quantities with the coordinate (z) perpendicular to that plane. The plasma expands and compresses isentropically; the magnetic field is calculated self-consistently, which leads to approximately parabolic field line shape in planes parallel to the magnetopause (the xz plane), with maximum field curvature near the magnetopause edge of the LLBL. Coupling to the ionosphere via region 1 field-aligned currents is included. The effects of the ionosphere are represented by two parallel resistive plates at fixed height above and below the equatorial plane. The model can be used to investigate the influence of various physical parameters, for example, viscous and magnetic Reynolds numbers, and of boundary conditions at the magnetopause and in the magnetosphere on the LLBL development in the -x direction. Special attention is given to viscous effects which, under suitable circumstances, lead to a region 1 current that first increases and then reduces with increasing longitude away from local noon. Asymptotic matching of the antisunward motion of the cool LLBL plasma to sunward convection of hot plasma in the magnetosphere is illustrated along with the entrainment of magnetospheric plasma by the antisunward LLBL flow.

Drakou, E.↗
DOE OSTI2021

Simulation of the SPARC plasma boundary with the UEDGE code

Here in this work, the UEDGE edge transport code is used to examine conditions in the SPARC divertor and edge plasma for various levels of carbon impurity and power from the core (P SOL ). A double-null magnetic configuration is simulated assuming up-down symmetry in geometry and physics. The anomalous heat and particle transport coefficients are tuned to match empirical predictions for SPARC's midplane density profiles, target plate heat flux profiles, and inner/outer divertor power sharing. Convective transport is included on the low-field side, while on the high-field side the transport is modeled as purely diffusive. Hydrogen neutrals are modeled as a fluid with inertial effects, and a carbon impurity is included using the fixed-fraction model. We find that detachment induced by impurity seeding could significantly reduce the heat flux to the divertor surfaces in the SPARC tokamak. At P SOL = 28 MW (the value predicted for SPARC's full-power H-mode scenario) cases with both divertor legs detached were obtained with a carbon impurity fraction between 0.3%–1.4%, far below Z eff limits for SPARC. When the plasma in the outer leg is detached, the peak heat flux density perpendicular to the target plate is below 1 MW m -2 , electron and ion temperatures are less than 1.5 eV, and momentum detachment is observed. However, the detachment state is found to be sensitive to the side-wall boundary conditions, the level of neutral pumping, and the target plate tilt. Finally, a broadly similar SOLPS simulation of SPARC is used to assess the appropriateness of the simpler impurity and neutral models used in UEDGE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗