Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space plasma”

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 757 records · Page 42

Numerical assessment of the impact of the guiding-centre approximation on fast ion simulations in NSTX

Guiding-centre (GC) and full-orbit (FO) simulations of the beam-injected fast ion distribution and the corresponding neutron emissivity have been carried out for magnetohydrodynamics-quiescent National Spherical Torus eXperiment (NSTX) plasmas, using a combination of ASCOT5 and DRESS, to assess the suitability of the GC approximation for fast ions in NSTX. It was found that GC and FO simulations predicted substantially different steady-state distributions in both position and velocity space and different neutron emissivity profiles, leading to a 15% reduction in the predicted global neutron rate for FO relative to GC. These changes accompany a higher magnetic moment in FO, and correspond to a change in particle orbits from co-passing to trapped and stagnation orbits. ASCOT5 was also benchmarked against TRANSP/NUBEAM with input loaded entirely from TRANSP/NUBEAM output files, with agreement found between the GC simulations when finite Larmor radius (FLR) corrections were omitted. ASCOT5 FO and TRANSP/NUBEAM with FLR produced fast ion distributions which differed in localised regions, but predicted global neutron rates which agree within 3%.

ASCOT↗

Dataset from: Improved multispecies Dougherty collisions

The Dougherty model Fokker-Planck operator is extended to describe nonlinear full-f collisions between multiple species in plasmas. Simple relations for cross-species primitive moments are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, avoids unphysical temperatures and satisfies the H-theorem unlike an analogous Bhatnagar-Gross-Krook operator. Formulas for both a Cartesian velocity-space and a gyroaveraged operator are provided for use in Vlasov as well as long-wavelength gyrokinetic models. We present an algorithm for the discontinuous Galerkin discretization of this operator, and provide results from relaxation and Landau damping benchmarks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electromagnetic Analysis of ITER Electron Cyclotron Emission Model

The ITER electron cyclotron emission (ECE) diagnostic system is located at diagnostics shielding module 2 (DSM2), equatorial port 9 (EP9), to measure electron temperature profile and electron temperature fluctuations, and also assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two hot sources for calibration and two couples of mirrors for two different optical views. It shall receive that the large electromagnetic (EM) loads up to 100-MN/m3 force density due to the eddy currents on the copper mirrors generated by the short transient plasma disruption. The simplified EM analysis model using magnetic field ( B ) data and flux variations (dB/dt) method based on the worst case of plasma disruption, MD_DW_EXP16MS_CATIII, has been developed to calculate the force and moments for each ECE component and support structure for EP9 DSM2 Bay2 and Bay3. The B and dB/dt methods use the constant B and dB/dt assumed for each local bay space, but the real B and dB/dt are basically varied along the radial direction of each bay. The complex global model with the local ECE components and support structure using Maxwell transient has been developed to do EM analysis to compare the results of the B and dB/dt models. The volumetric EM force density of the whole ECE components and support structure can be used for the subsequent structural analysis. In conclusion, combined with the thermal and nuclear loads and seismic and initial loads, the ECE integration analysis can be finalized for the operation case.

Fang, J. [Princeton Plasma Physics Laboratory (PPP↗

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical simulation of a relativistic magnetron using a fluid electron model

An approach to numerically modeling relativistic magnetrons, in which the electrons are represented with a relativistic fluid, is described. A principal effect in the operation of a magnetron is space-charge-limited (SCL) emission of electrons from the cathode. We have developed an approximate SCL emission boundary condition for the fluid electron model. This boundary condition prescribes the flux of electrons as a function of the normal component of the electric field on the boundary. Here we show the results of a benchmarking activity that applies the fluid SCL boundary condition to the one-dimensional Child–Langmuir diode problem and a canonical two-dimensional diode problem. Simulation results for a two-dimensional A6 magnetron are then presented. Computed bunching of the electron cloud occurs and coincides with significant microwave power generation. Numerical convergence of the solution is considered. Sharp gradients in the solution quantities at the diocotron resonance, spanning an interval of three to four grid cells in the most well-resolved case, are present and likely affect convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spectroscopic diagnostics of high-temperature plasma in stellar coronae using Fe XXIV-XXIV K-shell lines with XRISM

The RS CVn type binary star GT Muscae was observed during its quiescence using the Resolve X-ray microcalorimeter spectrometer onboard XRISM. The main and satellite lines of the Fe XXIV-XXIV K-shell transitions were resolved for the first time from stellar sources. We conducted line ratio analysis to investigate any deviations from collisional ionization equilibrium and Maxwell electron energy distribution with a single temperature. By using five combinations of direct excitation lines and dielectronic recombination satellite lines in three line complexes (Fe He$\alpha$, Ly$\alpha$, and He$\beta$), we found that the plasma is well characterized by two-temperature thermal plasmas with temperatures of 1.7 and 4.3 keV, which is consistent with thermal broadening of Fe xxv, and the broad-band fitting results in the 1.7–10 keV band. Other forms of deviation from a single-temperature plasma, such as different ionization and electron temperatures or the $\kappa$ distribution for the electron energy distributions, are not favored, which is reasonable for stellar coronae at quiescence. This study demonstrates the utility of the Fe K-shell line ratio diagnostics to probe plasma conditions using X-ray microcalorimeters.

X-rays: stars↗

Benchmark of the Fe xvv 𝓡 ratio in photoionized plasma during eclipse of Centaurus X-3 with XRISM/Resolve

The $\mathcal {R}$ ratio is a useful diagnostic of the X-ray emitting astrophysical plasmas and is defined as the intensity ratio of the forbidden over the inter-combination lines in the K$\alpha$ line complex of He-like ions. The value is altered by excitation processes (electron impact or UV photoexcitation) from the metastable upper level of the forbidden line, thereby constraining the electron density or UV field intensity. The diagnostic has been applied mostly in electron density constraints in collisionally ionized plasmas using low-Z elements, as was originally proposed for the Sun (Gabriel & Jordan, 1969a, MNRAS, 145, 241), but it can also be used in photoionized plasmas. To make use of this diagnostic, we need to know its value in the limit of no excitation of metastables ($\mathcal {R}_{0}$), which depends on the element, how the plasmas are formed, how the lines are propagated, and the spectral resolution affecting line blending principally with satellite lines from Li-like ions. We benchmark $\mathcal {R}_0$ for photoionized plasmas by comparing calculations using radiative transfer codes and observation data taken with the Resolve X-ray microcalorimeter onboard XRISM. We use the Fe xxv He$\alpha$ line complex of the photo-ionized plasma in Centaurus X-3 observed during eclipse, in which the plasma is expected to be in the limit of no metastable excitation. The measured $\mathcal {R} = 0.65 \pm 0.08$ is consistent with the value calculated using xstar for the plasma parameters derived from other line ratios of the spectrum. We conclude that the $\mathcal {R}$ ratio diagnostic can be used for high-Z elements such as Fe in photoionized plasmas, which has wide applications in plasmas around compact objects at various scales.

X-rays: binaries↗

Measurement of the radius dependence of charged-particle jet suppression in Pb–Pb collisions at $\sqrt{s_{NN}}=5.02$ TeV

The ALICE Collaboration reports a differential measurement of inclusive jet suppression using pp and Pb–Pb collision data at a center-of-mass energy per nucleon–nucleon collision $\sqrt{s_{NN}}=5.02$ TeV. Charged-particle jets are reconstructed using the anti-$k_T$ algorithm with resolution parameters R=0.2, 0.3, 0.4, 0.5, and 0.6 in pp collisions and R=0.2, 0.4, 0.6 in central (0–10%), semi-central (30–50%), and peripheral (60–80%) Pb–Pb collisions. A novel approach based on machine learning is employed to mitigate the influence of jet background. This enables measurements of inclusive jet suppression in new regions of phase space, including down to the lowest jet $p_T$ ≥ 40 GeV/c at R=0.6 in central Pb–Pb collisions. This is an important step for discriminating different models of jet quenching in the quark–gluon plasma. The transverse momentum spectra, nuclear modification factors, derived cross section, and nuclear modification factor ratios for different jet resolution parameters of charged-particle jets are presented and compared to model predictions. A mild dependence of the nuclear modification factor ratios on collision centrality and resolution parameter is observed. The results are compared to a variety of jet-quenching models with varying levels of agreement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The radioactive nuclei 26 Al and 60 Fe in the Cosmos and in the solar system

The cosmic evolution of the chemical elements from the Big Bang to the present time is driven by nuclear fusion reactions inside stars and stellar explosions. A cycle of matter recurrently re-processes metal-enriched stellar ejecta into the next generation of stars. The study of cosmic nucleosynthesis and this matter cycle requires the understanding of the physics of nuclear reactions, of the conditions at which the nuclear reactions are activated inside the stars and stellar explosions, of the stellar ejection mechanisms through winds and explosions, and of the transport of the ejecta towards the next cycle, from hot plasma to cold, star-forming gas. Due to the long timescales of stellar evolution, and because of the infrequent occurrence of stellar explosions, observational studies are challenging, as they have biases in time and space as well as different sensitivities related to the various astronomical methods. Here, we describe in detail the astrophysical and nuclear-physical processes involved in creating two radioactive isotopes useful in such studies, 26 Al and 60 Fe. Due to their radioactive lifetime of the order of a million years, these isotopes are suitable to characterise simultaneously the processes of nuclear fusion reactions and of interstellar transport. We describe and discuss the nuclear reactions involved in the production and destruction of 26 Al and 60 Fe, the key characteristics of the stellar sites of their nucleosynthesis and their interstellar journey after ejection from the nucleosynthesis sites. This allows us to connect the theoretical astrophysical aspects to the variety of astronomical messengers presented here, from stardust and cosmic-ray composition measurements, through observation of γ rays produced by radioactivity, to material deposited in deep-sea ocean crusts and to the inferred composition of the first solids that have formed in the Solar System. We show that considering measurements of the isotopic ratio of 26 Al to 60 Fe eliminate some of the unknowns when interpreting astronomical results, and discuss the lessons learned from these two isotopes on cosmic chemical evolution. This review paper has emerged from an ISSI-BJ Team project in 2017–2019, bringing together nuclear physicists, astronomers, and astrophysicists in this inter-disciplinary discussion.

79 ASTRONOMY AND ASTROPHYSICS↗

New, improved analysis of correlation ECE data to accurately determine turbulent electron temperature spectra and magnitudes (invited)

Turbulent electron temperature fluctuation measurement using a correlation electron cyclotron emission (CECE) radiometer has become an important diagnostic for studying energy transport in fusion plasmas, and its use is widespread in tokamaks (DIII-D, ASDEX Upgrade, Alcator C-Mod, Tore Supra, EAST, TCV, HL-2A, etc.). The CECE diagnostic typically performs correlation analysis between two closely spaced (within the turbulent correlation length) ECE channels that are dominated by uncorrelated thermal noise emission. This allows electron temperature fluctuations embedded in the thermal noise to be revealed and fluctuation level and spectra determined. We have demonstrated a new, improved CECE coherency-based analysis for calculating the temperature fluctuation frequency spectrum and level, which has been verified both numerically through the simulation of synthetic ECE radiometer data and through analysis of experimental data from the CECE system on DIII-D. The new formulation places coherency-based analysis on a firm foundational footing and corrects some currently published methodologies. This new method accurately accounts for bias error in the coherence function and correctly calculates noise levels for a fixed data record length. It provides excellent accuracy in determining temperature fluctuation level (e.g., <10% error) even for a small realization number in the ensemble average. The method also has a smaller uncertainty (i.e., error bar) in the power spectrum when compared to the more standard cross-power method when evaluated at low coherency. Direct calculation of system noise level using correlation between randomized intermediate frequency signals is recommended.

Wang, G. (ORCID:0000000225739827)↗

Effects of a nuclear-disturbed environment on electromagnetic wave propagation through the atmosphere

This paper investigates the effects of a nuclear-disturbed environment on the transmission of electromagnetic (EM) waves through the atmosphere. An atmospheric nuclear detonation can produce heightened free electron densities in the surrounding atmosphere that can disrupt EM waves that propagate through the disturbed region. Radiation transport models simulated the ionization and free electron densities created in the atmosphere from a 1 MT detonation at heights of burst of 5 km, 25 km, and 75 km. Recombination rates for the free electrons in the atmosphere were applied, from previous work in the literature, to determine the nuclear-induced electron densities as a function of time and space after the detonation. A ray-tracing algorithm was applied to determine the refraction and reflection of waves propagating in the different nuclear-disturbed environments. The simulation results show that the free electron plasma created from an atmospheric nuclear detonation depend on the height of burst of the weapon, the weapon yield, and the time after detonation. Detonations at higher altitudes produce higher free electron densities for greater durations and over larger ranges. The larger the free electron densities, the greater the impact on EM wavelengths in regards to refraction, reflection, and absorption in the atmosphere. An analysis of modern infrastructure and the effects of nuclear-disturbed atmospheres on different signal wavelengths and systems is discussed.

42 ENGINEERING↗

Large-scale Linear Magnetic Holes with Magnetic Mirror Properties in Hybrid Simulations of Solar Wind Turbulence

Magnetic holes (MHs) are coherent magnetic field dips whose size ranges from fluid to kinetic scale, ubiquitously observed in the heliosphere and in planetary environments. Despite the long-standing effort in interpreting the abundance of observations, the origin and properties of MHs are still debated. In this Letter, we investigate the interplay between plasma turbulence and MHs, using a 2D hybrid simulation initialized with solar wind parameters. We show that fully developed turbulence exhibits localized elongated magnetic depressions, whose properties are consistent with linear MHs frequently encountered in space. The observed MHs develop self-consistently from the initial magnetic field perturbations by trapping hot ions with large pitch angles. Ion trapping produces an enhanced perpendicular temperature anisotropy that makes MHs stable for hundreds of ion gyroperiods, despite the surrounding turbulence. We introduce a new quantity, based on local magnetic field and ion temperature values, to measure the efficiency of ion trapping, with potential applications to the detection of MHs in satellite measurements. We complement this method by analyzing the ion velocity distribution functions inside MHs. Our diagnostics reveal the presence of trapped gyrotropic ion populations, whose velocity distribution is consistent with a loss cone, as expected for the motion of particles inside a magnetic mirror. Our results have potential implications for the theoretical and numerical modeling of MHs.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating the carbon footprint of the integrated DBD‐plasma bi‐reforming unit via laboratory scale experiments and scaled‐up process modeling

Catalytic dielectric barrier discharge (DBD) plasma reactor experiments were performed in a tubular glass reactor with a 2 mm gap at 550°C to facilitate the reaction kinetics of steam added dry reforming or bireforming. The best specific energy input obtained was 11.2 eV/molecule feed at CO 2 :CH 4 :H 2 O of 4.5:1:4.5 ratio and gas hour space velocity (GHSV) = 432 h −1 . This value was used to design a conceptual process and assess the environmental impact of methane steam reforming-based H 2 production 18.4 kmol/h CO 2 emission processing into H 2 :CO = 2 syngas, with an emphasis on the carbon footprint.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energy and particle balance during plasma detachment in a long-leg divertor configuration

Abstract Comprehensive studies of energy and particle balances in the transition to plasma detachment in an alternative divertor configuration with long outer legs are shown. Numerical simulations are performed with the 2D code suite SOLPS 4.3, using a disconnected double null grid with narrow, tightly baffled long poloidal leg divertors at the outer lower target and outer upper target. A particle count scan is performed using the ‘closed gas box’ model, where the tunable parameter in the simulations is the total number of deuterium particles in the simulation space and all other parameters are held fixed, including a constant input power and trace neon impurity radiation, to assess the physics of the transition to detachment in the system as the particle count increases. Three main aspects of the physics of divertor detachment are addressed: the criteria for the local onset of divertor detachment in each of the divertors, the distribution of heat flux and other plasma parameters between the four divertors as each divertor transitions to detachment, and the role of perpendicular transport in the transition to the detached regime. A synergistic mechanism by which the cross-field transport is reduced by factors associated with the onset of plasma recombination effects is identified. These results are compared to the existing understanding of the physics of the transition to plasma detachment in standard divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shafranov-shift destabilization of ballooning-type micro-instabilities

Contrary to common wisdom, we argue that the Shafranov shift is destabilizing for the ambient ballooning-type instabilities, which account for most of the transport flux in tokamak core plasmas. Higher Shafranov shift indeed reduces the magnetic drift frequency ω d , i.e. improves the bad curvature mildly, around the outboard midplane. However, this improvement is limited to a very narrow spatial region while ω d is increased over the remaining poloidal space. The eigenfunction averaged $\langle$ω d $\rangle$ is effectively enhanced due to the finite mode width of the eigenfunction and is thus destabilizing to the ballooning-type mode, as demonstrated by gyrokinetic simulations with the CGYRO code using local Miller equilibrium geometry. The predicted nonlinear flux also increases with Shafranov shift, consistent with linear simulations. The reduced transport model TGLF can capture the physics reasonably well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transfer learning driven design optimization for inertial confinement fusion

Transfer learning is a promising approach to create predictive models that incorporate simulation and experimental data into a common framework. In this technique, a neural network is first trained on a large database of simulations and then partially retrained on sparse sets of experimental data to adjust predictions to be more consistent with reality. Previously, this technique has been used to create predictive models of Omega [Humbird et al., IEEE Trans. Plasma Sci. 48, 61–70 (2019)] and NIF [Humbird et al., Phys. Plasmas 28, 042709 (2021); Kustowski et al., Mach. Learn. 3, 015035 (2022)] inertial confinement fusion (ICF) experiments that are more accurate than simulations alone. Here, in this work, we conduct a transfer learning driven hypothetical ICF campaign in which the goal is to maximize experimental neutron yield via Bayesian optimization. The transfer learning model achieves yields within 5% of the maximum achievable yield in a modest-sized design space in fewer than 20 experiments. Furthermore, we demonstrate that this method is more efficient at optimizing designs than traditional model calibration techniques commonly employed in ICF design. Such an approach to ICF design could enable robust optimization of experimental performance under uncertainty.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improved performance of stellarator coil design optimization

Following up on earlier work which demonstrated an improved numerical stellarator coil design optimization performance by the use of stochastic optimization (Lobsien et al., Nucl. Fusion, vol. 58 (10), 2018, 106013), it is demonstrated here that significant further improvements can be made – lower field errors and improved robustness – for a Wendelstein 7-X test case. This is done by increasing the sample size and applying fully three-dimensional perturbations, but most importantly, by changing the design sequence in which the optimization targets are applied: optimization for field error is conducted first, with coil shape penalties only added to the objective function at a later step in the design process. A robust, feasible coil configuration with a local maximum field error of 3.66 % and an average field error of 0.95 % is achieved here, as compared to a maximum local field error of 6.08 % and average field error of 1.56 % found in our earlier work. These new results are compared to those found without stochastic optimization using the FOCUS and ONSET suites. The relationship between local minima in the optimization space and coil shape penalties is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗