Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “feature fusion”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods↗
DOE OSTI2020

Effect of Nuclear Data Covariances on Integral Experiment Design with Sensitivity and Uncertainty Analysis

Washington River Protection Solutions (WRPS) uses MCNP6.2 and the Whisper code for criticality safety analyses of the Hanford Tank Farm. Together the codes derive baseline upper subcritical limits (USLs) for the waste models using experimental benchmarks. Whisper returns higher USLs, i.e. , has less of a conservative penalty, when the neutronic similarity of the experimental benchmarks to the application is high. Unfortunately, few critical benchmarks have high similarity to the Hanford tanks. The waste in the tanks is highly dilute in plutonium and contains large masses of weakly neutron-absorbing elements like iron and manganese. Experimental benchmarks typically have low sensitivity to these absorbers because they are present as structural materials. Lacking similar benchmarks, new Thermal Epithermal eXperiment (TEX) configurations with high Pu content and interstitial iron absorbers have been designed for the criticality safety validation. The features of the design have been iterated upon to maximize the similarity between the experiment and different Hanford waste models. The similarity is quantified with sensitivity analysis and uncertainty quantification using the representativity coefficient, or c k . The representativity calculation requires nuclear data covariances, which may differ between nuclear data libraries and between library versions. Because of these variations, the optimal design may depend on the nuclear data covariances library. A scenario can be envisioned where an experiment is designed, and c k is maximized, with one set of covariance data. However, when the covariance data is changed, say from ENDF/B-VII.1 to ENDF/B-VIII.0, and the benchmark is used in a criticality safety evaluation, the experiment becomes suboptimal with respect to c k . In this paper, we present how the optimal design of the new TEX experiments varied depending on the nuclear data covariances used to calculate c k . We compare ENDF/B-VII.1 and ENDF/B-VIII.0, as if the library had been updated since the design of the experiment. Additionally, we use JEFF3.3 to simulate if the covariance data of a different library had been used. The results show that the covariances do have an important effect on the designs, less so for thermal systems (where the data are more consistent between evaluations) and more so for epithermal systems where more differences exist.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Multidimensional modeling of non-equilibrium plasma generated by a radio-frequency corona discharge

Low-temperature plasma (LTP) ignition concepts rely on the production of radical and charged species to speed up the onset of combustion in spark-ignition engines. These features are responsible for the superior performance of LTP igniters under extremely dilute combustion operation that is not achievable by conventional spark igniters. Additionally, LTP discharges extend the lifetime of the igniters, due to the avoidance of spark processes. For these reasons, the engine research community and the automotive industry have shown growing interest in this technology in the recent years. As of today, computational fluid-dynamics (CFD) codes typically used by the multi-dimensional engine modeling community do not have reliable models to describe LTP ignition processes. One key missing piece of information is the physical and chemical properties of the plasma and their effect on combustion ignition. Most non-equilibrium plasma simulations reported in literature are based on simplified, canonical geometries, with simple discharge excitation schemes. Here we conduct multi-dimensional modeling of the non-equilibrium plasma generated by an application-relevant radio-frequency (RF) corona discharge in air. Three test cases are simulated, characterized by different environmental pressure levels and peak electrode voltage values at room temperature. Streamer penetration, electron number density, atomic oxygen production, and bulk gas temperature distribution in the first 10 sinusoidal pulses are presented and discussed. This model can be used as a key tool for an in-depth understanding of RF-corona discharge for automotive applications and provides the basis for future implementations of dedicated LTP ignition models in CFD codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Role of fast-ion transport manipulating safety factor profile in KSTAR early diverting discharges

Magnetohydrodynamics (MHD) activity is analyzed for a particular KSTAR discharge category, the main features of which are a high-$q_\mathrm{min}$ and a broad q profile. An interesting finding in this scenario is that a steady-state high-$q_\mathrm{min}$ and broad q profile has been constructed even without a strong off-axis current drive scheme. Time traces of magnetics spectrograms and electron cyclotron emission imaging (ECEI) clearly show that Alfvénic type MHD activity appears as a high-$q_\mathrm{min}$ and broad q profile is formed. Conversely, this activity disappears as the current profiles evolve toward a low-$q_\mathrm{min}$ profile. In this study, the effect of MHD modes and accompanying fast-particle transport on the formation of various current profiles is explored. The frequency patterns of the magnetics spectrograms and ECEI coherence are compared using a NOVA analysis to determine which series of Alfvénic modes are active. Then, the kick-model is applied to potential modes, and used to estimate how much fast-ions are transported and beam driven current profiles are varied. Here, by comparing the reconstructed total current profile from the kinetic EFIT with the beam-driven current profile estimated from the kick model, it is determined how energetic particle transport is responsible for maintaining the safety factor profile and which q profile can provide high β N operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A mechanism for reduced compression in indirectly driven layered capsule implosions

High-yield implosions on the National Ignition Facility rely on maintaining low entropy in the deuterium–tritium fuel, quantified by its adiabat, in order to efficiently couple energy to the hot spot through high compression of the fuel layer. We present very-high-resolution xRAGE simulation results that study the impacts of interfacial mixing and the jetting of materials due to surface defects, defects on internal interfaces, voids, and engineering features on fuel layer compression. Defects and voids are typically neglected in implosion simulations due to their small size and three-dimensional geometry. Our results showed that supersonic jets of material arise through weak spots in the shell at peak implosion velocity that prevent uniform compression of the fuel layer even when they do not introduce contaminant into the hot spot. This occurs despite maintaining low fuel entropy, since the formation of the weak spots involves nonradial displacement of fuel mass. In contrast, simulations show that fuel–ablator mixing due to interfacial instabilities has a much smaller impact on compression. We show that defects on interior interfaces of plastic capsules decrease compression by 15% to 25% and interfacial mixing between the ablator and fuel decreases compression by less than 1% for implosions with plastic or high-density carbon (HDC) ablators. For low adiabat implosions, the impact of jetting seeded by the support tent can also decrease the compression by 25%. We demonstrate that the inclusion of interior defects in simulations can explain the inferred compression in two fielded plastic capsule implosions and that the inclusion of voids, for which available characterization has large uncertainties, in simulations of HDC capsule implosions has a qualitatively consistent impact. This mechanism offers a potential explanation for persistently overestimated fuel compression in design simulations of layered implosions on the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Runaway electron beam dynamics at low plasma density in DIII-D: energy distribution, current profile, and internal instability

Parameters of the post-disruption runaway electron (RE) beam in the collisionless background plasma achieved after deuterium injection are investigated in DIII-D. The spatially resolved RE energy distribution function is measured for the first time during the RE plateau stage by inverting hard X-ray bremsstrahlung spectra. It has maximum energy up to 20 MeV and a non-monotonous feature at 5{6 MeV observed only in the core of the beam supporting the possibility of kinetic instabilities. The RE current profile is reconstructed for the first time using the spatially resolved RE energy distribution. It is found to be more peaked than the pre-disruption plasma current, with higher internal inductance, suggesting preferential formation of REs in the core plasma or potentially a radially inward motion of the REs. The accessed relatively low-current (180 kA) RE beam is found stable, likely due to its elevated q profile. From this base stable equilibrium, instability is accessed by ramping the solenoid current. Under a large applied accelerating loop voltage an internal kink mode is observed. The kink mode leads to sawtooth-like relaxation of the RE current profile, but drives no RE loss. Under conditions of a large decelerating voltage, chirping Alfvénic instabilities driven by REs during the RE plateau stage are observed for the first time. The low-frequency (0.1{7 MHz) Alfvénic instabilities, supposedly compressional Alfvén eigenmodes, are found to be correlated with a RE loss signal. Finally, this provides novel plasma conditions to study these instabilities and investigate opportunities to utilize them for RE control. Bifurcated electron cyclotron emission (ECE) spectrum is measured during the collisionless RE plateau, with a break point at ≈100 GHz suggesting resonant absorption of the ECE at low frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

On hybrid scenarios in KSTAR

We report the status of hybrid scenario experiments in Korea Superconducting Tokamak Advanced Research (KSTAR). The hybrid scenario is defined as stationary discharges with ${\beta _{\rm{N}}}$ $ \geqslant $ 2.4 and ${H_{89}}$ $ \geqslant $ 2.0 at ${q_{95}}$ < 6.5 without or with very mild sawtooth activities in KSTAR. It is being developed towards reactor-relevant conditions. High performance of ${\beta _{\rm{N}}} \lesssim 3.0$, ${H_{89}} \lesssim 2.4$ and G-factor ($ \equiv {\beta _{\rm{N}}}{H_{89}}/q_{95}^2$) $ \lesssim 0.4$6 has been achieved and sustained for $ \gtrsim 40{{\tau }_{\text{E}}}$ at ne/nGW ~ 0.7 with heating power of $ \lesssim $5 MW. Some KSTAR hybrid discharges exhibit a unique feature of a slow transition from conventional H-mode to hybrid mode after the third neutral beam injection. The reason for the confinement enhancement is extensively studied in this transition period of a representative discharge exhibiting a common feature of KSTAR hybrid scenarios. 0D performance analysis with magnetohydrodynamic activities, 1D kinetic profile dynamics, power balance analysis, linear gyro-kinetic analysis and edge pedestal stability analysis were conducted. The enhancement is thought to be from both the core and the pedestal. The improvement in the core region of the ion energy channel is observed from the linear gyro-kinetic analysis considering the electromagnetic, the fast ion, the Shafranov shift, ${\it{\omega }_{E \times B}}$, and the magnetic shear effect. The electromagnetic finite $\beta $ stabilisation plays a role in the inner core region at ${\rho _{\rm{tor}}}\sim0.35$ together with the fast ion effect. The alpha stabilisation effect is also found at ${\rho _{\rm{tor}}}\sim0.5$. ${\omega _{E \times B}}$, which could reduce the linear growth of the ion temperature gradient mode in the outer core region at ${\rho _{\rm{tor}}}\sim0.5 - 0.7$ with the highest contribution from the toroidal rotation. Furthermore, regarding the improvement in the pedestal, Shafranov shift broadens the stability boundary of the pedestal in support of the diamagnetic effect. The pedestal height and width could be reproduced by the EPED model, while a realistic current profile is used to calculate the internal inductance for Shafranov shift. Based on these findings, a comprehensive confinement enhancement mechanism has been proposed by considering the core-edge interplay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Development of an Enhanced Radiation Physics Toolset for Modeling Spectral and Imaging Signatures in the Warm Dense Matter Experiments

Radiative and atomic processes in plasmas play a critical role in a wide variety of high energy density laboratory plasma (HEDLP) experiments. The emission, absorption, and transport of radiation can strongly affect the overall energetics and evolution of such plasmas. In addition, radiation-based diagnostics – including imaging, spectroscopy, and absolute flux measurements – are widely used to determine key features of HEDLPs. To advance our understanding of HEDLP science, it is vital to have high-fidelity computational physics tools that have well-tested radiation physics modeling, and that are readily accessible to researchers in the HEDLP community. Simulations play an extremely important role for planning and designing the experiments, as well as for post-experiment data analysis. Prism Computational Sciences develops software that is used by National Laboratories and universities (including five members of LaserNetUS network). Prominent examples of such research efforts include z-pinch and short-pulse laser experiments designed to study the basic physics of photoionized plasmas and photoionization fronts, as well as their application to astrophysical plasmas. The main effort was dedicated to the development of non-equilibrium equation-of-state (EOS) models within the HELIOS-CR code, a hydrodynamics code with inline collisional-radiative atomic kinetics. Gas cell experiments on Z and Omega demonstrated the importance of non-equilibrium effects on atomic kinetics in photoionized plasmas. Recent proof-of-principle experiments on Omega EP confirmed the advantages of using a short-pulse laser to create an intense radiation drive, leading to additional experiments being proposed. Photoionization front experiments at LLE also emphasizes the importance of radiation and atomic physics. In both studies, HELIOS-CR simulations played a crucial role in computing non-equilibrium opacities and ionization distributions. A newly developed non-LTE EOS model will help addressing possible non-equilibrium effects, on for example specific heat, and their influence in plasma evolution. Prism also implemented support for open-source atomic data generated by the Flexible Atomic Code. This allows researchers to generate custom atomic tables and use them within the complex framework of simulation tools developed by Prism. The ability to use open-source atomic data would be extremely valuable for hydrodynamics and spectroscopic simulations that include high-Z materials, e.g., picosecond x-ray pulse generation experiments. Support for new atomic structures was fully implemented, and the data can be used by all simulation tools developed at Prism: radiation-hydrodynamics, imaging and spectroscopy, EOS and opacity. The development resulted in a significant fidelity enhancement to the simulations tools developed by Prism that are currently used in other cutting-edge experiments including: opacity measurement experiments performed to both understand the basic radiative and atomic properties of plasmas as well as provide data for more accurately modeling the internal structure of the Sun and other stars, high-intensity short-pulse laser experiments performed to develop short-wavelength light sources for use as backlighters and to investigate fast ignition concepts for inertial fusion energy; capsule implosion experiments designed to develop inertial fusion as an energy source, etc.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Measurements and modeling of type-I and type-II ELMs heat flux to the DIII-D divertor

Type-I and type-II edge-localized-modes (ELMs) heat flux profiles measured at the DIII-D divertor feature a peak in the vicinity of the strike-point and a plateau in the scrape-off-layer (SOL), which extends to the first wall. The plateau is present in attached and detached divertors and it is found to originate with plasma bursts upstream in the SOL. The integrated ELM heat flux is distributed at ~65% in the peak and ~35% in this plateau. The parallel loss model, currently used at ITER to predict power loads to the walls, is benchmarked using these results in the primary and secondary divertors with unprecedented constraints using experimental input data for ELM size, radial velocity, energy, electron temperature and density, heat flux footprints and number of filaments. The model can reproduce the experimental near-SOL peak within ~20%, but cannot match the SOL plateau. Employing a two-component approach for the ELM radial velocity, as guided by intermittent data, the full radial heat flux profile can be well matched. The ELM-averaged radial velocity at the separatrix, which explains profile widening, increases from ~0.2 km s –1 in attached to ~0.8 km s –1 in detached scenarios, as the ELM filaments' path becomes electrically disconnected from the sheath at the target. The results presented here indicate filaments fragmentation as a possible mechanism for ELM transport to the far-SOL and provide evidence on the beneficial role of detachment to mitigate ELM flux in the divertor far-SOL. However, these findings imply that wall regions far from the strike points in future machines should be designed to withstand significant heat flux, even for small-ELM regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

I-mode pedestal relaxation events in the Alcator C-Mod and ASDEX Upgrade tokamaks

Abstract In some conditions, I-mode plasmas can feature pedestal relaxation events (PREs) that transiently enhance the energy reaching the divertor target plates. To shed light on their appearance, characteristics and energy reaching the divertor targets, a comparative study between two tokamaks — Alcator C-Mod and ASDEX Upgrade — is carried out. It is found that PREs appear only in a subset of I-mode discharges, mainly when the plasma is close to the H-mode transition. Also, a growing oscillating precursor before the PRE onset is observed in the region close to the separatrix in both devices, and a discussion on a possible triggering mechanism is outlined. The PRE relative energy loss from the confined region is found to increase with decreasing pedestal top collisionality ν ped * . Similarly, also the relative electron temperature drop at the pedestal top, which is related to the conductive energy loss, rises with decreasing ν ped * . Based on these relations, the PRE relative energy loss in future devices such as DEMO and ARC is estimated. Finally, the divertor peak energy fluence due to the PRE is measured on each device. Those values are then compared to the model introduced in Eich et al (2017 Nucl. Mater. Energy 12 84–90) for type-I edge localized modes. The model is shown to provide an upper boundary for PRE energy fluence data, while a lower boundary is found by dividing the model by three. These two boundaries are used to make projections of the PRE divertor energy fluence to DEMO and ARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Linear simulation of magnetohydrodynamic plasma response to three-dimensional magnetic perturbations in high-β P plasmas

In this work, we report the numerical analyses of linear magnetohydrodynamics (MHD) plasma response to applied three-dimensional magnetic perturbations (MPs) in a joint DIII-D/EAST collaboration on high-β P (poloidal beta) plasmas, utilizing the extended-MHD code M3D-C1, with the purpose of realizing a better understanding of the existing experiment in which the n=3 MPs were applied to such high-β P plasmas attempting to control large amplitude type-I ELMs. Such high-β P plasmas obtained at the DIII-D tokamak feature an upper-biased double null configuration, a high edge safety factor q 95 ~6.4, and a stable internal transport barrier (ITB) leading to relatively high core pressures. Single-fluid simulations show that the plasma response to n=3 MPs, including both non-resonant/kinking and resonant components, is significantly weaker than that to n=1 or 2 MPs. To survey the impact of q 95 on plasma response to applied MPs, the SEGWAY (Self-consistent Equilibrium Generating Workflow for AnalYsis) module, developed in the OMFIT integrated modelling framework, is employed to generate a series of equilibria with a wide range of q 95 while other key parameters including the normalized beta, electron density at pedestal top, and plasma shape are kept fixed. Compared to the vacuum response, single-fluid M3D-C1 simulations predict a much more significant decrease of resonant plasma response to the applied n=3 MPs at the maximum penetration radii as q 95 increases. In contrast to single-fluid simulation results showing resonant penetration occurs only near the pedestal top where the E×B toroidal rotation frequency is zero, two-fluid simulations show two comparable resonant penetrations locating near the pedestal top and the ITB foot, where the perpendicular electron rotation frequency is zero. Such resonant field penetration near the ITB foot may be responsible for the observed formation of a staircase structure in both electron density and temperature profiles and thereby a considerable deterioration of global plasma performance when MPs are applied in high-β P plasmas. Motivated by this numerical work, we provide some ideas for the future research, with the purpose of realizing effective ELM control in such high-β P plasmas on the DIII-D and EAST devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

SOLEDGE3X full vessel plasma boundary simulations of ITER non-active phase plasmas

The onset of detachment in the ITER machine is analyzed in this work through the help of 2D-axisymmetric boundary plasma simulations with the SOLEDGE3X-EIRENE code, which features a numerical domain for the plasma solver extending up to the first wall. The plasma boundary is computed in scenarios from the first non-active phase of ITER, in pure H and at 20 MW. This set of simulations is used in two aspects: first, to study the plasma detachment in the divertor, and second, the plasma conditions, fluxes, and beryllium erosion at the first wall. Here, the code results are also compared to those obtained with the well-established SOLPS-ITER code, which includes a plasma numerical domain only covering the main SOL. Results show an increase in the SOL width λ q with increasing density, and a detailed analysis is carried out, for the first time, on each of the different plasma-neutral interactions in the code’s physics model in EIRENE. The gross beryllium erosion rates of first wall panels are estimated from 2D simulations, with the aim of assessing their sensitivity to two parameters: the divertor density regime, and the presence of density shoulders in the far-SOL formed by enhanced perpendicular transport at this location. The erosion contributions from neutrals and ions are considered in each case, and the charge-exchange atoms fluxes and energy distributions are provided, highlighting the two atom populations (cold and charge-exchange).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

LandScan HD: a high-resolution gridded ambient population methodology for the world

Unwarned population distributions accounting for routine human activities are needed to address many global human security challenges, including disasters, conflict, and infrastructure demand. LandScan High Definition (LSHD) supports this need through gridded ambient population estimates that measure average human presence between daytime and nighttime at a high spatial resolution of 3 arcseconds (approximately 90 m). Although LSHD has traditionally been produced on a country-specific basis, advances in global foundational data and computational resources now enable scaling its methodology to the world. Combining aspects of top-down and bottom-up gridded population methods, LSHD allocates subnational population totals from authoritative statistics to built-up areas based on occupancy estimates for multiple facility types (e.g., residential, commercial) and then reaggregates these estimates to a global population grid. We scale this approach by organizing the LSHD data stack into a 1° resolution tileset of vector analytic features, enabling an efficient and repeatable workflow for all countries worldwide. Examining the Philippines as an output of the global LSHD baseline dataset, we contrast unwarned and residential (WorldPop) population distributions by (1) exploring a practical application of flood risk assessment and (2) evaluating their congruence with outcomes of collective human activities (subnational CO 2 emissions). Finally, we discuss plans to address current LSHD limitations through data/modeling and uncertainty quantification improvements and provide outlook for workflow automation and extending the model to social, demographic and economic population characteristics.

Building morphology↗