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 91 records · Page 5

The thermal catastrophe model of substorms

Equations are presented for the physical mechanisms involved in the resonant absorption of Alfven waves in the plasma sheet boundary layer. It is shown that energy absorbed by the plasma-sheet particles is a function of the central plasma sheet temperature. The heating curve, when coupled with convective transport, yielded an equation of state for the steady state plasma sheet, whose solution has the form of a mathematical catastrophe. The master equation includes dynamic terms describing the transition across the thermal catastrophe, making it possible to evaluate the time scale for the catastrophe to occur.

Goertz, C. K.↗

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NSTX-U theory, modeling and analysis results

Here, the mission of the low aspect ratio spherical tokamak NSTX-U is to advance the physics basis and technical solutions required for optimizing the configuration of next-step steady-state tokamak fusion devices. NSTX-U will ultimately operate at up to 2 MA of plasma current and 1 T toroidal field on axis for 5 s, and has available up to 15 MW of neutral beam injection power at different tangency radii and 6 MW of high harmonic fast wave heating. With these capabilities NSTX-U will develop the physics understanding and control tools to ramp-up and sustain high performance fully non-inductive plasmas with large bootstrap fraction and enhanced confinement enabled via the low aspect ratio, high beta configuration. With its unique capabilities, NSTX-U research also supports ITER and other critical fusion development needs. Super-Alfvénic ions in beam-heated NSTX-U plasmas access energetic particle (EP) parameter space that is relevant for both α-heated conventional and low aspect ratio burning plasmas. NSTX-U can also generate very large target heat fluxes to test conventional and innovative plasma exhaust and plasma facing component solutions. This paper summarizes recent analysis, theory and modelling progress to advance the tokamak physics basis in the areas of macrostability and 3D fields, EP stability and fast ion transport, thermal transport and pedestal structure, boundary and plasma material interaction, RF heating, scenario optimization and real-time control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observations of a plectonemic configuration in a stable magnetized plasma jet

Astrophysical jets are collimated high-speed outflows emerging from spinning and accreting matter around celestial objects and may spontaneously result from self-organized processes. Magnetic self-organization is commonly observed in laboratory plasma physics experiments; however, they require close-fitting flux conservers to constrain and stabilize the toroidal or cylindrical structures. Here we report the first observations of a long, stable, free-boundary plasma jet far from chamber walls, embedding a double-helix magnetic structure resembling a force-free plectonemic Taylor state. The jets arise from an experimental setup that mimics an accretion disk and has no close-fitting solid flux conserver. The results support the hypothesis that self-organization could be a universal, intrinsic explanation for jet formation, collimation, and stability and may help explain double-helix features in celestial observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reimagining full wave rf quasilinear theory in a tokamak

The velocity dependent resonant interaction of particles with applied radio frequency (rf) waves during heating and current drive in the presence of pitch angle scattering collisions gives rise to narrow collisional velocity space boundary layers that dramatically enhance the role of collisions as recently shown by Catto (J. Plasma Phys., vol. 86, 815860302, 2020). The behavior is a generalization of the narrow collisional boundary layer that forms during Landau damping as found by Johnston (Phys. Fluids, vol. 14, 1971, pp. 2719-2726) and Auerbach (Phys. Fluids, vol. 20, 1977, pp. 1836-1844). For a wave of parallel wave number k|| interacting with weakly collisional plasma species of collision frequency ν and thermal speed vth , the effective collision frequency becomes of order ν(k_||v_th /ν)^2/3>> ν . The narrow boundary layers that arise because of the diffusive nature of the collisions allows a physically meaningful wave-particle interaction time to be defined that is the inverse of this effective collision frequency. The collisionality implied by the narrow boundary layer results in changes in the standard quasilinear treatment of applied rf fields in tokamaks while remaining consistent with causality. These changes occur because successive poloidal interactions with the rf are correlated in tokamak geometry and because the resonant velocity space dependent interactions are controlled by the spatial and temporal behavior of the perturbed full wave fields rather than just the spatially local Landau and Doppler shifted cyclotron wave-particle resonance condition associated with unperturbed motion of the particles. The correlation of successive poloidal circuits of the tokamak leads to the appearance in the quasilinear operator of transit averaged resonance conditions localized in velocity space boundary layers that maintain negative definite entropy production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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

This report documents the FY25 Theory and Simulation Performance Target (TSPT) of developing an integrated modeling framework for fusion reactor design and assessment (FREDA). Over Q1-Q4, new capabilities were developed across both plasma and engineering domains and demonstrated on an example representation of a Compact Advanced Tokamak with a Dual Cooled Lead Lithium blanket. This represents a first-of-a-kind demonstration of coupled core-to-wall-to-engineering for a reactor. Self-consistent CESOL workflows were applied to provide core, pedestal, and SOL prediction; new modules were developed for energetic particle stability (FAR3D) and transport (TGLF-EP) analysis; and boundary plasma modeling (SOLPS-ITER, BOUT++/Hermes-3) was expanded to evaluate wall and divertor heat fluxes and interface with engineering thermal analysis. A parameterized CAD tool, TRACER, was expanded to generate medium-fidelity divertor, blanket, and coil geometries; OpenFOAM and Diablo workflows were applied for first-wall and divertor thermal analyses with helium cooling; and reduced-order models were created for high-mass-flux divertor cooling. Magnet multiphysics capabilities were verified between Elmer, Diablo, and a new MFEM-based solver, and workflows enable stress, thermal, and neutron-fluence analysis of TF coils with neutronics-driven heating. Nuclear and blanket analysis workflows were demonstrated, including tritium breeding, transport, and CFD-informed thermo-mechanical assessment. Preliminary multi-fidelity uncertainty quantification workflows were applied to boundary modeling codes and shown to achieve variance reductions with fewer high-fidelity boundary simulations. Key findings highlight the challenges of resolving the ITEP gap to find suitable balance between wall and divertor loads, neutron heating, and practical limits of PFC cooling. Next step priorities are to develop automated workflows to check boundary code convergence and detachment, implement tighter physics-engineering CAD provenance tracking, and inclusion of plasma-material interface models for SLAG and tungsten cracking behavior. Collectively, these developments establish sophisticated capabilities for predictive, multi-fidelity, whole-device modeling that integrates plasma physics, materials, magnets, and nuclear engineering to guide pathways to viable Fusion Pilot Plant design points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The plasma environment of comets

U.S. research activities in the area of cometary plasma physics during 1987-1990 are reviewed. Consideration is given to mass loading and its consequences in the upstream region, the cometary shock, the cometosheath, the diamagnetic cavity boundary and the inner shock, and the plasma tail. Special attention is given to models and observations that have modified the pre-encounter understanding of cometary plasma environments.

Gombosi, Tamas I.↗

Effect of geometric and magnetic boundary conditions on magnetic islands in 3D force-free ideal MHD equilibria

Three-dimensional boundary features and magnetic perturbations affect the magnetic topology of the force-free ideal MHD equilibria formed within a given volume. The use of the PSI-Tet 3D finite-element equilibrium code reveals the existence of magnetic islands inside such equilibria within boundaries with three-dimensional features in configurations that directly depend on the natures of those features. Here, the nature and limits of the influence of such features are explored, as is the influence of flux boundary conditions to qualitatively identify the structure of magnetic islands in Taylor state equilibria with three-dimensional boundary conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics of magnetospheric boundary layers

The central ideas of this grant are that the magnetospheric boundary layers link disparate regions of the magnetosphere together, and the global behavior of the magnetosphere can be understood only by understanding the linking mechanisms. Accordingly the present grant includes simultaneous research on the global, meso-, and micro-scale physics of the magnetosphere and its boundary layers. These boundary layers include the bow shock, magnetosheath, the plasma sheet boundary layer, and the ionosphere. Analytic, numerical and simulation projects have been performed on these subjects, as well as comparison of theoretical results with observational data. Very good progress has been made, with four papers published or in press and two additional papers submitted for publication during the six month period 1 June - 30 November 1993. At least two projects are currently being written up. In addition, members of the group have given papers at scientific meetings. The further structure of this report is as follows: section two contains brief accounts of research completed during the last six months, while section three describes the research projects intended for the grant's final period.

Cairns, I. H.↗

Physics of magnetospheric boundary layers

This final report was concerned with the ideas that: (1) magnetospheric boundary layers link disparate regions of the magnetosphere-solar wind system together; and (2) global behavior of the magnetosphere can be understood only by understanding its internal linking mechanisms and those with the solar wind. The research project involved simultaneous research on the global-, meso-, and micro-scale physics of the magnetosphere and its boundary layers, which included the bow shock, the magnetosheath, the plasma sheet boundary layer, and the ionosphere. Analytic, numerical, and simulation projects were performed on these subjects, as well as comparisons of theoretical results with observational data. Other related activity included in the research included: (1) prediction of geomagnetic activity; (2) global MHD (magnetohydrodynamic) simulations; (3) Alfven resonance heating; and (4) Critical Ionization Velocity (CIV) effect. In the appendixes are list of personnel involved, list of papers published; and reprints or photocopies of papers produced for this report.

Cairns, Iver H.↗

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of CMAS viscosity on the infiltration depth in thermal barrier coatings of different microstructures

Calcium magnesium aluminum silicate (CMAS) is one of the leading concerns for the gas turbine industry. The effects of CMAS viscosity and the coating microstructure on CMAS infiltration depth were explored by conducting a time dependent interaction study. Three CMAS compositions were used from literature, and their viscosities predicted through FactSage viscosity module were drastically different. The interaction was carried out on three different TBCs synthesized using the solution precursor plasma spray process (SPPS): two of the TBCs were made of yttrium aluminum garnet (YAG) having different microstructures that promote different modes of CMAS infiltration, and one TBC was made of gadolinium zirconate (GZO). All samples had stress relieving vertical cracks and different intensities of horizontally banded porosity, (inter pass boundaries IPBs). A concentration of 100 mg/cm 2 of CMAS was applied on the TBCs which were then subjected to a 5-minute interaction at 1300 °C. Samples were analyzed using scanning electron microscopy (SEM), electron dispersive X-ray spectroscopy (EDXS), and transmission electron microscopy (TEM). Low viscosity CMAS readily penetrated the TBCs while more viscous CMAS showed less penetration. The depth of CMAS infiltration depended on the coating microstructure. In the YAG with IPBs, the CMAS spread horizontally in the IPBs before infiltrating deeper, resulting in reduced infiltration depth compared to other samples in spite of having wider vertical cracks. TEM and EDXS analysis were performed to investigate the phases present in the CMAS-TBC interaction region in YAG. Two regions were chosen, the top TBC surface in direct contact with the sea of CMAS, and the region at the CMAS penetration ended within the coating. Finally, the results showed that no secondary phases like apatite were observed in YAG, thus it can be concluded that the arrest of CMAS happened solely because of CMAS viscosity and the short infiltration time.

36 MATERIALS SCIENCE↗

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collisional effects on resonant particles in quasilinear theory

A careful examination of the effects of collisions on resonant wave-particle interactions leads to an alternate interpretation and deeper understanding of the quasilinear operator originally formulated by Kennel and Engelmann (Phys. Fluids vol. 9, 1966, pp. 2377- 2388) for collisionless, magnetized plasmas, and widely used to model radio frequency heating and current drive. The resonant and nearly resonant particles are particularly sensitive to collisions that pitch angle scatter them out of and into resonance. As a result, the resonant particle-wave interactions occur in the center of a narrow collisional boundary when the collision frequency nu is very small compared to the wave frequency omega. The diffusive nature of the pitch angle scattering combined with the wave-particle resonance condition enhances the collision frequency by (omega/nu)2/3 >>1, resulting in an effective resonant particle collision time of tau_int ~ (nu /omega)2/3 nu <<1/ nu . A rigorous collisional boundary layer analysis generalizes the standard quasilinear operator to a form that is fully consistent with Kennel-Englemann, but allows replacing the delta function appearing in the diffusivity with a simple integral (having the appropriate delta function limit) retaining the new physics associated with the narrow boundary layer, while preserving the entropy production principle. The limitations of the collisional boundary layer treatment are also estimated, and indicate that substantial departures from Maxwellian are not permitted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic Simulation of Inhomogeneous Plasma with a Variable Sized Grid System

Space plasmas are usually inhomogeneous and irregular. Particle in cell (PIC) plasma simulations, however, often consider the plasma to be homogeneous and periodic. Here non-periodic PIC simulations with inhomogeneities in the density, temperature and background magnetic field are considered. Boundary conditions and a variable sized grid system are discussed, along with an application to a magnetospheric plasma physics problem in the auroral zone.

Schriver, David↗

Bound states and electromagnetic radiation of relativistically rotating cylindrical wells

Here, we compute the effect of rigid rotation on the non-relativistic bound states. The energy levels of the bound states increase with the angular velocity of rotation until at certain value of the angular velocity they are completely pushed out into the continuum which corresponds to dissociation of the bound states. When the angular velocity exceeds the critical value at which the ground state disappears into the continuum, no bound state is possible. This effect should have important consequences for the phenomenology of the quark-gluon plasma. One of the ways to study it experimentally is to observe the electromagnetic radiation emitted by a rotating bound state. We compute the corresponding intensity of electromagnetic radiation and show that it strongly depends on the angular velocity of rotation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Turbulent field fluctuations in gyrokinetic and fluid plasmas

A key uncertainty in the design and development of magnetic confinement fusion energy reactors is predicting edge plasma turbulence. An essential step in overcoming this uncertainty is the validation in accuracy of reduced turbulent transport models. Drift-reduced Braginskii two-fluid theory is one such set of reduced equations that has for decades simulated boundary plasmas in experiment, but significant questions exist regarding its predictive ability. To this end, using a novel physics-informed deep learning framework, we demonstrate the first ever direct quantitative comparisons of turbulent field fluctuations between electrostatic two-fluid theory and electromagnetic gyrokinetic modeling with good overall agreement found in magnetized helical plasmas at low normalized pressure. This framework presents a new technique for the numerical validation and discovery of reduced global plasma turbulence models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗