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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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↗
DOE OSTI2020

H-mode pedestal improvements with neon injection in DIII-D

Improvements in the H-mode pedestal pressure and global energy confinement were determined on DIII-D with both 3 and 6 MW of neutral beam (NBI) heating with neon injection. The pedestal pressure increased, primarily associated with a significant increase in the pedestal density with a relatively smaller decrease in the pedestal temperature. At the same neon injection rate, a rapid decrease in edge localized mode (ELM) frequency and ~35% increase in plasma stored energy were observed in the 3 MW discharge, while in the 6 MW discharge, the ELM frequency was gradually decreased and the plasma stored energy gradually increased by up to 20%. The measured pedestal widths from both 3 and 6 MW discharges matched the EPED1.0 width scaling, ${\Delta _{{\psi _N}}} = 0.089\beta _{p,ped}^{1/2}$ ,data-, within 20% deviations before and after neon injection. The peeling-ballooning mode (PBM) stability analysis showed that after neon injection, the ballooning stability boundary was increased, while the peeling stability boundary only dropped once the neon level was relatively high. The ballooning stability boundary increase with neon injection is consistent with increase of the ion diamagnetic stabilization (i.e. including the contribution from carbon, neon as well as deuterium ions) which is used in the criterion to determine when the calculated linear growth rates of PBM are unstable. Consequently, the new ELITE computed PBM boundaries were more consistent with the operating points than when only deuterium was included. This in turn suggests that the impurity may help to improve pedestal pressure through affecting ion diamagnetic frequency in the pedestal region, which could positively lead to the ballooning stability boundary extension.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Boundary condition effects on runaway electron mitigation coil modeling for the SPARC and DIII-D tokamaks

Abstract Extended-MHD modeling of planned Runaway Electron Mitigation Coils (REMC) for SPARC and DIII-D is performed with the NIMROD code. A coil has been designed for each machine, with the two differing in shape and location, but both having n = 1 symmetry (with n the toroidal mode number). Compared to previous modeling efforts, three improvements are made to the simulations boundary conditions. First a resistive wall model is used in place of an ideal wall. Second, the ThinCurr code is used to compute the time-dependent 3D fields used as magnetic boundary conditions for the simulations. Third, the simulation boundary is moved from the first-wall location to the Vacuum Vessel (VV), which extends the boundary past the location of the internal REMC. To remove the 3D coil from the simulation domain, an equivalent set of 3D fields is calculated at the VV boundary that produce approximately the same field distribution at the last closed flux surface assuming vacuum between the two. Each of these three boundary condition improvements leads to an improvement in the predicted performance of the REMC for both machines. The resistive wall alone primarily effects the resonance of the coil with the plasma after the TQ, affecting the q-profile evolution in the SPARC modeling, and allowing the applied spectrum to be modified in response to the plasma in the DIII-D modeling. The movement of the simulation boundary has the most significant effect on the RE confinement overall, including in the early stages, particularly for a DIII-D inner wall limited equilibrium, where the RE loss fraction increases from 90% to > 99%, with SPARC RE losses also occurring much earlier when the boundary is placed at the VV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Space platforms and meteorological applications

The use of Shuttle launched and tended space platforms (SP) for research and earth observations in LEO are examined. SPs will accommodate instruments for extended time periods, requiring a Shuttle visit for maintenance and package exchange once every six months. Designs have included off-the-shelf hardware to reduce cost, and comprise instruments for studies of cosmic rays, solar physics, materials processing, space plasma physics, and environmental observations. Solar cell power arrays are detailed, along with thermal control, data management, and attitude control subsystems. Meteorological applications are discussed in terms of the Lower Atmospheric Research Satellite (LARS), for research into coupling of radiative, chemical, and dynamical processes, interactions between the earth's surface, boundary layer, and troposphere, and dynamic and energetic relationships between lower atmospheric processes. Flight configurations and add-on concepts for the SP to lead to a manned space station are indicated.

Mitchell, K. L.↗
NASA NTRS2017

Formation of Heliospheric Arcs of Slow Solar Wind

A major challenge in solar and heliospheric physics is understanding the origin and nature of the so-called slow solar wind. The Sun's atmosphere is divided into magnetically open regions, known as coronal holes, where the plasma streams out freely and fills the solar system, and closed regions, where the plasma is confined to coronal loops. The boundary between these regions extends outward as the heliospheric current sheet (HCS). Measurements of plasma composition strongly imply that much of the slow wind consists of plasma from the closed corona that escapes onto open field lines, presumably by field-line opening or by interchange reconnection. Both of these processes are expected to release closed-field plasma into the solar wind within and immediately adjacent to the HCS. Mysteriously, however, slow wind with closed-field plasma composition is often observed in situ far from the HCS. We use high-resolution, three-dimensional, magnetohydrodynamic simulations to calculate the dynamics of a coronal hole with a geometry that includes a narrow corridor flanked by closed field and is driven by supergranule-like flows at the coronal-hole boundary. These dynamics produce giant arcs of closed-field plasma that originate at the open-closed boundary in the corona, but extend far from the HCS and span tens of degrees in latitude and longitude at Earth. We conclude that such structures can account for the long-puzzling slow-wind observations.

Higginson, A. K.↗
DOE OSTI2024

Artificial correlation heating in PIC simulations

The Particle-in-Cell (PIC) method, a cornerstone in plasma modeling, is widely employed for its ability to simulate kinetic phenomena in device-scale domains. Part of what makes this possible is that computational macroparticles represent many physical particles. It converges under certain constraints, including a grid spacing that resolves the Debye length and a time step small enough to respect the Courant–Friedrichs–Lewy condition and plasma frequency stability limit. Here, we introduce a new constraint necessary to avoid Artificial Correlation Heating (ACH). This requires that the macroparticle coupling strength be smaller than one, Γ ω < 1, where Γ ω ≡ Γω 2/3 ⁠, Γ = Z 2 e 2 /(4πε 0 ak B T) is the physical coupling strength, and w is the macroparticle weight. This is particularly relevant to 3D simulations of dense plasmas, which are becoming common with modern computing power. If this condition is violated, the finite macroparticle weight artificially enhances the coupling strength and causes the plasma to heat until the macroparticle coupling strength is near unity, depending on the grid resolution. A comprehensive model of ACH is developed that incorporates electron density, temperature, macroparticle weight, and grid resolution. It is then tested using PIC simulations, delineating the boundaries of the method's applicability and offering a predictive framework for ACH. Moreover, the research explores a runaway heating process induced by ACH in the presence of ionization, which can lead to numerical instability. A conclusion of this study is that the onset of ACH can impose a more stringent constraint on the macroparticle weight and average number of macroparticles per cell than what is typically expected, particularly in 3D simulations of dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Minnowbrook VI: 2009 Workshop on Flow Physics and Control for Internal and External Aerodynamics

Topics covered include: Flow Physics and control for Internal and External Aerodynamics (not in TOC...starts on pg13); Breaking CFD Bottlenecks in Gas-Turbine Flow-Path Design; Streamwise Vortices on the Convex Surfaces of Circular Cylinders and Turbomachinery Blading; DNS and Embedded DNS as Tools for Investigating Unsteady Heat Transfer Phenomena in Turbines; Cavitation, Flow Structure and Turbulence in the Tip Region of a Rotor Blade; Development and Application of Plasma Actuators for Active Control of High-Speed and High Reynolds Number Flows; Active Flow Control of Lifting Surface With Flap-Current Activities and Future Directions; Closed-Loop Control of Vortex Formation in Separated Flows; Global Instability on Laminar Separation Bubbles-Revisited; Very Large-Scale Motions in Smooth and Rough Wall Boundary Layers; Instability of a Supersonic Boundary-Layer With Localized Roughness; Active Control of Open Cavities; Amplitude Scaling of Active Separation Control; U.S. Air Force Research Laboratory's Need for Flow Physics and Control With Applications Involving Aero-Optics and Weapon Bay Cavities; Some Issues Related to Integrating Active Flow Control With Flight Control; Active Flow Control Strategies Using Surface Pressure Measurements; Reduction of Unsteady Forcing in a Vaned, Contra-Rotating Transonic Turbine Configuration; Active Flow Control Stator With Coanda Surface; Controlling Separation in Turbomachines; Flow Control on Low-Pressure Turbine Airfoils Using Vortex Generator Jets; Reduced Order Modeling Incompressible Flows; Study and Control of Flow Past Disk, and Circular and Rectangular Cylinders Aligned in the Flow; Periodic Forcing of a Turbulent Axisymmetric Wake; Control of Vortex Breakdown in Critical Swirl Regime Using Azimuthal Forcing; External and Turbomachinery Flow Control Working Group; Boundary Layers, Transitions and Separation; Efficiency Considerations in Low Pressure Turbines; Summary of Conference; and Final Plenary Session Transcript.

LaGraff, John E.↗
DOE OSTI2026

Permanent magnets for ELM suppression in tokamaks: feasibility and operational compatibility

Permanent magnets can provide static magnetic fields without power supplies or feed lines, offering a simpler alternative to conventional electromagnets in tokamaks. This work examines permanent magnet arrays (PMAs) for edge localized mode (ELM) control in double-null (DN) configurations, where traditional RMP coils exhibit limited effectiveness. Linear plasma response calculations using IPEC assess high-field-side (HFS) permanent magnet placement, benchmarking performance against DIII-D low-field-side (LFS) internal coils (I-coils) while evaluating engineering constraints. Modeling results demonstrate that HFS permanent magnet configurations generate HFS plasma response amplitudes more than 5 times larger than conventional I-coil systems. While other resonant response metrics show comparable or reduced response relative to I-coils, the combination of conventional RMP systems and permanent magnets is expected to provide enhanced performance. Operational impact assessments of persistent magnetic fields, including sideband field effects and startup/ramp-up compatibility, reveal acceptable performance within established operational boundaries. This analysis establishes PMAs as a technically viable approach for DN ELM control with reductions in system complexity, motivating further experimental validation on existing tokamak facilities.

3D magnetic perturbation↗
NASA NTRS1986

Destiny of earthward streaming plasma in the plasmasheet boundary layer

The dynamics of the earth's magnetotail have been investigated, and it has become clear that the plasmasheet boundary layer field lines map into the Region I Field-Aligned Currents (FAC) of the auroral zone. It is pointed out that the role of earthward streaming ions in the plasmasheet boundary layer may be of fundamental importance in the understanding of magnetotail dynamics, auroral zone physics, and especially for ionospheric-magnetospheric interactions. The present paper has the objective to evaluate propagation characteristics for the earthward streaming ions observed in the plasmasheet boundary layer. An investigation is conducted of the propagation characteristics of protons in the plasmasheet boundary layer using independent single particle dynamics, and conclusions are discussed. The density of earthward streaming ions found in the plasmasheet boundary layer should include the ring current as well as the auroral zone precipitaiton and inner plasmasheet regions of the magnetosphere.

Green, J. L.↗
DOE OSTI2025

Doubling fusion power with volumetric optimization in magnetic confinement fusion devices

A technique, volumetric power optimization, is presented for enhancing the power output of magnetic confinement fusion devices. Applied to a tokamak, this approach involves shifting the burning plasma region to a larger plasma volume while introducing minimal perturbations to the plasma boundary shape. This edge perturbation—squareness—is analogous to pinching and stretching the edge boundary. Stability calculations confirm that this edge alteration is compatible with maintaining plasma stability. This optimization method for optimizing fusion power output could improve the performance of magnetic confinement fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2016

Electron-Scale Measurements of Magnetic Reconnection in Space

Magnetic reconnection is a fundamental physical process in plasmas whereby stored magnetic energy is converted into heat and kinetic energy of charged particles. Reconnection occurs in many astrophysical plasma environments and in laboratory plasmas. Using measurements with very high time resolution, NASA's Magnetospheric Multiscale (MMS) mission has found direct evidence for electron demagnetization and acceleration at sites along the sunward boundary of Earth's magnetosphere where the interplanetary magnetic field reconnects with the terrestrial magnetic field. We have (i) observed the conversion of magnetic energy to particle energy; (ii) measured the electric field and current, which together cause the dissipation of magnetic energy; and (iii) identified the electron population that carries the current as a result of demagnetization and acceleration within the reconnection diffusion/dissipation region.

Burch, J. L.↗
DOE OSTI2022

VPIC-Kokkos

VPIC is a particle-in-cell (PIC) code that solves the Maxwell-Boltzmann system of equations and can be used to simulate a wide variety of plasmas, including those that are fully kinetic and relativistic. Utilizing the Kokkos performance-portable framework, VPIC achieves high performance on multiple CPU and GPU architectures and is adaptable to future platforms with minimal developer effort. VPIC features very powerful input decks, allowing insertion of arbitrary C++ code for custom diagnostics, boundary conditions, and additional physics models.

Luedtke, Scott↗
NASA NTRS1974

The internal structure of the convective envelope of the sun

An axisymmetric hydrodynamic model of the sun consisting of a core at radiant equilibrium and a convective envelope is studied with consideration of gravitational forces, rotation, turbulence effects, and electromagnetic forces. The dimensionless parameters that control the equations of plasma motion within the sun and the boundary conditions are investigated; estimates are given for these parameters on the basis of experimental data and the existing theoretical models. Stability conditions are derived for convective equilibrium, developed cellular convection, and developed turbulent convection with consideration of rotation and turbulent Reynolds stresses. These conditions characterize the possible physical state of the matter in the convective envelope. Using estimates of the dimensionless parameters and models of some of them, a method of successive approximations is proposed for solution of the problem as a whole.

Yavorskaya, I. M.↗
NASA NTRS1992

A numerical code for a three-dimensional magnetospheric MHD equilibrium model

Two dimensional and three dimensional MHD equilibrium models were begun for Earth's magnetosphere. The original proposal was motivated by realizing that global, purely data based models of Earth's magnetosphere are inadequate for studying the underlying plasma physical principles according to which the magnetosphere evolves on the quasi-static convection time scale. Complex numerical grid generation schemes were established for a 3-D Poisson solver, and a robust Grad-Shafranov solver was coded for high beta MHD equilibria. Thus, the effects were calculated of both the magnetopause geometry and boundary conditions on the magnetotail current distribution.

Voigt, G.-H.↗
DOE OSTI2021

Collisional broadening of nonlinear resonant wave–particle interactions

A general procedure for understanding plasma behaviour when resonant wave–particle interactions are the sole destabilizing and transport mechanism or only heating and/or current drive source is highlighted without recourse to involved numerical or analytical treatments. These phenomena are characterized by transport that appears to be collisionless even though collisions play a central role in narrow collisional boundary layers. The order of magnitude estimates, which include nonlinear effects, are shown to provide expressions in agreement with the principal results of recent toroidal Alfvén eigenmode (TAE), toroidal magnetic field ripple, and heating and current drive treatments. More importantly, the retention of nonlinearities leads to new estimates of the alpha particle energy diffusivity at saturation for TAE modes, and the ripple threshold at which superbanana plateau evaluations of alpha particle transport are modified by nonlinear radial drift effects. In addition, the estimates indicate when quasilinear descriptions for heating and current drive will begin to fail. The phenomenological procedure demonstrates that in magnetic fusion relevant plasmas, narrow collisional boundary layers must be retained for resonant wave–particle interactions as they enhance the role of collisions, and make stochastic particle motion unlikely to be more important than other nonlinear processes.

Physics↗
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↗