Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma models”

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 379 records · Page 21

Large gyro-orbit model of ion velocity distribution in plasma near a wall in a grazing-angle magnetic field

A model is presented for the ion distribution function in a plasma at a solid target with a magnetic field $\boldsymbol {B}$ inclined at a small angle, $\alpha \ll 1$ (in radians), to the target. Adiabatic electrons are assumed, requiring $\alpha \gg \sqrt {Zm_{e}/m_{i}}$ , where $m_{e}$ and $m_{i}$ are the electron and ion mass, respectively, and $Z$ is the charge state of the ion. An electric field $\boldsymbol {E}$ is present to repel electrons, and so the characteristic size of the electrostatic potential $\phi$ is set by the electron temperature $T_{e}$ , $e\phi \sim T_{e}$ , where $e$ is the proton charge. An asymptotic scale separation between the Debye length $\lambda _{D} = \sqrt {\epsilon _0 T_{{e}} / e^{2} n_{{e}} }$ , the ion sound gyro-radius $\rho _{s} = \sqrt { m_{i} ( ZT_{e} + T_{i} ) } / (ZeB)$ and the size of the collisional region $d_{c} = \alpha \lambda _{\textrm {mfp}}$ is assumed, $\lambda _{D} \ll \rho _{s} \ll d_{c}$ . Here $\epsilon _0$ is the permittivity of free space, $n_{e}$ is the electron density, $T_{i}$ is the ion temperature, $B= |\boldsymbol {B}|$ and $\lambda _{\textrm {mfp}}$ is the collisional mean free path of an ion. The form of the ion distribution function is assumed at distances $x$ from the wall such that $\rho _{s} \ll x \ll d_{c}$ , that is, collisions are not treated. A self-consistent solution of the electrostatic potential for $x \sim \rho _{s}$ is required to solve for the quasi-periodic ion trajectories and for the ion distribution function at the target. The large gyro-orbit model presented here allows to bypass the numerical solution of $\phi (x)$ and results in an analytical expression for the ion distribution function at the target. It assumes that $\tau =T_{i}/(ZT_{e})\gg 1$ , and ignores the electric force on the quasi-periodic ion trajectory until close to the target. For $\tau \gtrsim 1$ , the model provides an extremely fast approximation to energy–angle distributions of ions at the target. These can be used to make sputtering predictions.

Physics↗

Interpretive modelling of boron transport in the boundary plasma of WEST experiments with the impurity powder dropper

Abstract Boron (B) powder injection is a potential alternative to glow discharge boronization as a wall conditioning method for tokamaks. This technique is currently being studied in WEST experiments, during which B powder is injected by an Impurity Powder Dropper developed by PPPL. In order to interpret and analyse experimental trends, and to help develop future experiments, a modelling workflow using a boundary plasma simulation (SOLEDGE-EIRENE) and powder ablation simulation (Dust Injection Simulator) was developed and tested. The effect of adding a B neutral source to simulated deuterium + oxygen (D + O) plasmas was compared to experimental data from the WEST C5 campaign, where B powder was injected in a dedicated experiment. While the impact of B injection on radiated power P rad measurements at the upper divertor was similar, there were significant differences in measurements of P rad , outer strike point electron temperature T e OSP and O-II line intensity at the lower divertor between experiment and simulation. This discrepancy suggests that those parameters were affected by phenomena not present in the simulations, with the most likely candidates being reduced D recycling and a reduced O sourcing from the divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

RF sheath induced sputtering on Proto-MPEX part 2: Impurity transport modeling and experimental comparison

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a pulsed linear plasma device that uses a radio frequency (RF) helicon antenna and a helicon ceramic window as the plasma ionization source. Modeling performed within this work on the helicon operations predicts hot spots on the helicon ceramic window due to an RF-induced sheath. The RF sheath potential leads to physical sputtering of the ceramic window predominantly where the helicon antenna sits around the helicon window. The role of oxygen in sputtering within Proto-MPEX is also investigated and is needed to match the experimental data. Experiments on Proto-MPEX show reasonable agreement with the modeled erosion/redeposition pattern on the helicon window, and the modeled impurity flux to the target shows radial transport within the device as expected from classical collisional impurity transport. Temperature screening of impurities is shown to be needed for Proto-MPEX to reduce impurity flux to the target. The modeling approach used in this work can be applied to other plasma RF applications for the determination of impurity production and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dust and powder in fusion plasmas: recent developments in theory, modeling, and experiments

Abstract In this paper, we present a brief historic overview of the research on dust in fusion devices with carbon plasma-facing components and then highlight the most recent developments in the post-carbon era of the field. In particular, we consider how the metallic dust form, mobilize, and interact with fusion plasmas and plasma facing components. Achievements in wall conditioning and associated anomalous plasma transport modification, including ELM suppression, with the powder injection technique is another focus of the paper. Capabilities of the state-of-art simulation tools to describe different aspects of dust in fusion devices are exemplified and new directions for future dust studies are brought forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurement and Modeling of Electron Temperature in Laboratory Photoionized Plasmas Relevant to Astrophysics

The electron temperature of photoionized plasmas characterizes the thermalization of photoelectrons, impacts the charge-state distribution, emissivity, and opacity through atomic recombination processes, and is needed to perform detailed comparisons with theory predictions. We discuss temperature measurements in laboratory photoionized plasmas and a comparison with model calculations done with several theory approximations and codes. These include a radiation-hydrodynamics simulation and two nonequilibrium heating and ionization models that tracked the evolution of the internal energy of the electrons. Furthermore, for the same physics model and X-ray flux time history, calculations were performed assuming steady-state or time-dependent conditions. The time history of steady-state results correlates with that of the X-ray flux, while that of the time-dependent cases does not, and it is qualitatively and quantitatively different from the steady-state case. Steady-state results significantly overestimated temperature measurements, while time-dependent results produced better approximations, which suggests the importance of transient effects in the experiment and also the need for time-resolved measurements.

Laboratory astrophysics↗

Implementation of a 9-point stencil in SOLPS-ITER and implications for Alcator C-Mod divertor plasma simulations

The SOLPS-ITER code suite is used worldwide for plasma edge modeling, the interpretation of experiments, as well as for the design of the ITER divertor. The numerical scheme of the plasma solver of the code, B2.5, is based on the assumption of perfectly field-aligned grids, while in practice grids are often strongly distorted to match divertor target shapes. Neglecting these grid distortion leads to qualitatively and quantitatively incorrect results for fluid neutral simulations, and may affect results in cold (detached) divertors even when using kinetic neutral simulations. In this contribution, we present the first results of a newly implemented 9-point stencil in B2.5 to properly handle misaligned grids. The new scheme is then applied to fluid neutral simulations of a well-diagnosed and previously modeled Alcator C-Mod discharge. Results are compared with the original 5-point scheme neglecting grid distortion effects, as well as with simulations including a full kinetic neutral model. We conclude that the 9-point stencil is essential to correctly model the transport of fluid neutrals on distorted grids, and to capture the effects of divertor closure on the fluid neutral behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Model for the electrical conductivity in dense plasma mixtures

A new density functional theory, average atom based model for the electrical conductivity of dense plasmas with a mixture of ion species, containing no adjustable parameters, is presented herein. The model takes the temperature, mass density and relative abundances of the species as input. It takes into account partial ionization, ionic structure, and core-valence orthogonality, and uses quantum mechanical calculations of cross sections. Comparison to an existing high fidelity but computationally expensive method reveals good agreement. The new model is computationally efficient and can reach high temperatures. A new mixing rule is also presented that gives reasonably accurate conductivities for high temperature plasma mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atmospheric pressure plasmas interacting with wet and dry microchannels: reverse surface ionization waves

Atmospheric pressure plasma jets (APPJs) are increasingly being used to functionalize polymers and dielectric materials for biomedical and biotechnology applications. Once such application is microfluidic labs-on-a-chip consisting of dielectric slabs with microchannel grooves hundreds of microns in width and depth. The periodic channels, an example of a complex surface, present challenges in terms of directly and uniformly exposing the surface to the plasma. In this paper, we discuss results from computational and experimental investigations of negative APPJs sustained in Ar/N 2 mixtures flowing into ambient air and incident onto a series of microchannels. Results from two-dimensional plasma hydrodynamics modeling are compared to experimental measurements of electric field and fast-camera imaging. The propagation of the plasma across dry microchannels largely consists of a sequence of surface ionization waves (SIWs) on the top ridges of the channels and bulk ionization waves (IWs) crossing over the channels. The IWs are directed into electric field enhanced vertices of the next ridge. The charging of these ridges produce reverse IWs responsible for the majority of the ionization. The propagation of the plasma across water filled microchannels evolve into hopping SIWs between the leading edges of the water channels, regions of electric enhancement due to polarization of the water. Positive, reverse IWs follow the pre-ionized path of the initial negative waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of Injected Current Streams on MHD Equilibrium Reconstruction of Local Helicity Injection Plasmas in a Spherical Tokamak

Open field line currents are intrinsic to DC helicity injection plasma startup and pose a challenge for inferring the plasma equilibrium with standard reconstruction analysis. Local helicity injection (LHI) is a type of DC helicity injection which uses small, modular current sources to drive force-free current along helical field lines to produce tokamak plasmas. MHD modeling and magnetic measurements during LHI indicate the injected current streams remain coherent as helical structures on the outboard edge of a core toroidal plasma that is tokamak-like in a toroidally averaged sense. To extract core plasma equilibrium properties, external magnetic diagnostics corrected for contributions from the injected current streams are fitted by a standard Grad-Shafranov equilibrium code. An iterative approach for estimating and subtracting the stream contributions from the diagnostic signals is described and applied to a model equilibrium database to reduce systematic errors introduced by the streams. Convergence is usually attained with 2 to 4 iterations, with derived equilibrium parameters matching the prescribed axisymmetric core values to within estimated experimental uncertainties. Accurate recovery of core parameters occurs when the ratio of the net toroidal windup current from the streams to the core plasma current is less than 0.2, which is typically satisfied in most experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parasitic Modulation of Microwave Signals by a Hypersonic Plasma Layer

During hypersonic flight, compressional and viscous heating of the air can form a plasma layer which encases the aircraft. If the boundary layer becomes turbulent, then the electron density fluctuations can effect a parasitic modulation in microwave signals transmitted through the plasma. We developed an approach for studying the interaction of microwave signals with a turbulent, hypersonic plasma layer. The approach affords a great deal of flexibility in both the plasma layer model and the antenna configuration. We then analyzed a situation in which microwaves, transmitted from a rectangular aperture antenna, propagate through a turbulent plasma layer to a distant receiver. We characterized the first-and second-order statistics of the computed parasitic modulation and quantified the depolarization of the signal. The amplitude fluctuations are lognormally distributed at low frequencies and Rice-distributed at high frequencies. Fluctuations in the copolarized phase and amplitude of the far-field signal are strongly anticorrelated. Here, we used a multioutput Gaussian process (MOGP) to model these quantities. The efficacy of the MOGP model is demonstrated by recovering the time evolution of the copolarized phase given the copolarized amplitude and occasional measurements of the phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterizing W sources in the all-W wall, all-RF WEST tokamak environment

In this work, experimental data, together with interpretive modeling tools, are examined to study trends in the tungsten (W) source in the all-W environment of the WEST tokamak, both from the divertor and from the main chamber. In particular, a poloidal limiter protecting an ion cyclotron resonance heating (ICRH) antenna is used as proxy for main chamber sourcing. The key study is carried out by stepping up lower hybrid current drive (LHCD) power, as the only auxiliary power source. Limiter and divertor W sources exhibit a qualitatively similar proportionality to the total power crossing the separatrix, P SEP , although the main chamber source remains substantially lower than the divertor source, for the range of P SEP accessible in the experiments. Intepretive modeling of the limiter source is carried out with a particle-in-cell (PIC) sheath model coupled to a surface sputtering model. Oxygen is used as a proxy for all light impurity species allowing for characterization of the critical W erosion regions. To get a good quantitative match to the data, it is necessary to assume that the oxygen arrives at the surface mostly at high ionization stages (4+ and above). A separate simulation with SOLEDGE-EIRENE, constrained to measured upstream scrape-off-layer plasma profiles, gives oxygen fractional abundances that are compatible with the PIC simulation result. This is understood to arise from transport processes that dominate over recombination. Substituting the LHCD by ICRH, in an equivalent experiment, the local W source exhibits a 3× enhancement. This can be matched by the simulation, by assuming local RF electric field rectification, based on ~100 eV peak-to-peak, near-antennna electric field. This work has highlighted the particular importance of understanding the ion charge state balance of light impurities as these are most likely the dominant sputtering species in fusion devices with high-Z walls.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing atomic physics at ultrahigh pressure using laser-driven implosions

Abstract Spectroscopic measurements of dense plasmas at billions of atmospheres provide tests to our fundamental understanding of how matter behaves at extreme conditions. Developing reliable atomic physics models at these conditions, benchmarked by experimental data, is crucial to an improved understanding of radiation transport in both stars and inertial fusion targets. However, detailed spectroscopic measurements at these conditions are rare, and traditional collisional-radiative equilibrium models, based on isolated-atom calculations and ad hoc continuum lowering models, have proved questionable at and beyond solid density. Here we report time-integrated and time-resolved x-ray spectroscopy measurements at several billion atmospheres using laser-driven implosions of Cu-doped targets. We use the imploding shell and its hot core at stagnation to probe the spectral changes of Cu-doped witness layer. These measurements indicate the necessity and viability of modeling dense plasmas with self-consistent methods like density-functional theory, which impact the accuracy of radiation transport simulations used to describe stellar evolution and the design of inertial fusion targets.

74 ATOMIC AND MOLECULAR PHYSICS↗

Influence of up-down asymmetry in plasma shape on RMP response

Shaping effect on the plasma response to the externally applied resonant magnetic perturbation field is investigated for both DIII-D and MAST experiments, utilizing toroidal modeling. The plasma boundary shape is systematically varied ranging from single-null (SN) to double-null (DN) configurations, while other equilibrium quantities are kept largely unchanged. The relative amplitude of the computed plasma surface displacement, between the top/bottom of the torus and the outboard mid-plane, is identified as the most reliable indicator that distinguishes the plasma response between the SN and DN configurations. Here, the underlying physics is the weakening of the edge-peeling component in the plasma response with increasing up-down symmetry of the plasma boundary shape.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Applicability of semiclassical theories in the strong-field plasma regime

For many purposes, classical plasma dynamics models can work surprisingly well, even for strong electromagnetic fields, approaching the Schwinger critical fields, and high frequencies, approaching the Compton frequency. However, the applicability of classical models tends to depend rather sensitively on the details of the problem. In the present paper, we study the specific case of plasma oscillations to draw a line between the classical and quantum relativistic regimes. Here, due to the field geometry of study, mechanisms like radiation reaction and Breit-Wheeler pair production, which tend to be important for electromagnetic fields, are rather effectively suppressed. Moreover, we find that the polarization current due to the electron spin is generally negligible for frequencies below the Compton frequency, compared with the free current, whose magnitude is well-approximated by the classical Vlasov theory. However, we show that pair creation due to the Schwinger mechanism can sometimes be important for surprisingly modest field strengths, of the order of 10% of the critical field or even smaller. A rough guideline for when the classical Vlasov theory can be applied is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Informing ARC divertor design and plasma facing material selection through integrated modeling

This INFUSE 2023 project between UCLA and Commonwealth Fusion Systems used computer modeling to test whether tungsten materials can survive in CFS's ARC fusion reactor divertor. The team simulated plasma conditions and material responses, finding that tungsten-rhenium alloys resist grain growth better than pure tungsten, and that hydrogen buildup depends more on particle flux than temperature. The work helps CFS design durable plasma-facing components for their fusion power plant.

36 MATERIALS SCIENCE↗

Dense plasma opacity via the multiple-scattering method

The calculation of the optical properties of hot dense plasmas with a model that has self-consistent plasma physics is a grand challenge for high energy density science. Here we exploit a recently developed electronic structure model that uses multiple scattering theory to solve the Kohn-Sham density functional theory equations for dense plasmas. Here we calculate opacities in this regime, validate the method, and apply it to recent experimental measurements of opacity for Cr, Ni, and Fe. Good agreement is found in the quasicontinuum region for Cr and Ni, while the self-consistent plasma physics of the approach cannot explain the observed difference between models and the experiment for Fe.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Model Membrane Systems Used to Study Plasma Membrane Lipid Asymmetry

It is well known that the lipid distribution in the bilayer leaflets of mammalian plasma membranes (PMs) is not symmetric. Despite this, model membrane studies have largely relied on chemically symmetric model membranes for the study of lipid–lipid and lipid–protein interactions. This is primarily due to the difficulty in preparing stable, asymmetric model membranes that are amenable to biophysical studies. However, in the last 20 years, efforts have been made in producing more biologically faithful model membranes. Here, we review several recently developed experimental and computational techniques for the robust generation of asymmetric model membranes and highlight a new and particularly promising technique to study membrane asymmetry.

59 BASIC BIOLOGICAL SCIENCES↗

On the relationship between the multi-region relaxed variational principle and resistive inner-layer theory

We show that the variational energy principle of the multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by 'ideal interfaces', and in each volume the magnetic field relaxes to a Taylor state, where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the Stepped-Pressure Equilibrium Code (SPEC) so that the equilibrium solution in each region is computed while preserving the force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low-beta plasma, we have been able to illustrate a relationship between the resistive singular-layer theory (Coppi et al 1966 Nucl. Fusion6 101; Glasser et al 1975 Phys. Fluids18 875–88) and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in a cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner-layer theories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗