SEARCH · Engineering Papers

Results for “Agreements”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Revisiting proton–proton fusion in chiral effective field theory

We calculate the S-factor for proton–proton fusion using chiral effective field theory interactions and currents. By performing order-by-order calculations with a variety of chiral interactions that are regularized and calibrated in different ways, we assess the uncertainty in the S-factor from the truncation of the effective field theory expansion and from the sensitivity of the S-factor to the short-distance axial current determined from three- and four-nucleon observables. We find that S(0) = (4.100 ± 0.024(syst) ± 0.013(stat) ± 0.008(g$_{A}$)) × 10$^{−23}$ MeV fm$^{2}$, where the three uncertainties arise, respectively, from the truncation of the effective field theory expansion, use of the two-nucleon axial current fit to few-nucleon observables and variation of the axial coupling constant within the recommended range. The increased value of S(0) compared to previous calculations is mainly driven by an increase in the recommended value for the axial coupling constant and is in agreement with a recent analysis based on pionless effective field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

The role of intermediate ΔΔ states in nucleon–nucleon scattering in the large-Nc and unitary limits, and ΔΔ and ΩΩ scattering

We explore potential explanations for why using large-Nc (Nc is the number of colors) scaling to determine the relative size of few-nucleon low-energy operators agrees with experiment even when dynamical Δ’s are not explicitly included. Given that the large-Nc analysis is predicated on the nucleons and Δ’s being degenerate, this is a curious result. We show that for purely S-wave interactions the relationships dictated by large-Nc scaling are unaffected whether the Δ is included or not. In the case of higher partial waves that do not mix with S-waves, the impact of the Δ is perturbative, which makes the agreement with naive (Δ-less) large-Nc ordering unsurprising. For higher partial waves that mix with S-waves, the nucleon and Δ would need to decouple to get agreement with naive large-Nc ordering. We find all NN, ΔN, and ΔΔ low energy coefficients for leading-order baryon–baryon scattering in Δ-full pionless effective field theory in terms of the two independent parameters dictated by the SU(2F) spin-flavor symmetry that arises in the Nc → ∞ limit. Because of recent lattice quantum chromodynamics results and experimental interest, we extend our analysis to the three-flavor case to study ΩΩ scattering. We show that in the unitary limit (where scattering lengths become infinite) one of the two SU(2F) parameters is driven to zero, resulting in enhanced symmetries, which agree with those found in spin-1/2 entanglement studies.

Richardson, Thomas R↗

DOE OSTI2021

Assessment of filtered cameras for quantitative 2D analysis of divertor conditions during detachment in JET L-mode plasmas

We report that estimates for 2D distributions of electron temperature, T e , electron density, n e , and atomic deuterium density, n o , in the JET divertor volume have been inferred from deuterium Balmer line intensity ratios obtained from tomographic reconstructions of divertor camera measurements. This enables also investigation of ionization, S ion , and recombination, S rec , rates. The analysis shows a decrease of T e to 0.5–1.0 eV throughout the outer divertor during detachment in low-confinement (L-mode) plasmas. Simultaneously, the high-n e region and the n 0 distribution in the outer divertor are observed to elongate and shift from the outer strike point towards the X-point. The observations are in qualitative agreement and follow the same sequence with modelling predictions of EDGE2D-EIRENE simulations of a density scan. While the method was found to provide good representation of the evolution of volumetric recombination during detachment, in agreement with the simulations, the movement of the ionization front upstream could not be followed due to lack of spatial overlap between the ionization region and the necessary emission distributions. Consequently, the representation of the ionization conditions and the particle balance in the detached outer divertor are compromised.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Transport and losses of fusion-born alpha particles in the presence of tearing modes using the new Toroidal Accelerated Particle Simulator (TAPaS)

Here, the transport and losses of fusion-born alpha particles are studied in the presence of a single-helicity tearing mode, characterized by (m = 2, n = 1). The analysis is performed by means of the recently developed Toroidal Accelerated Particle Simulator (TAPaS). Although such modes have usually been believed to result only in a local flattening of the radial profiles, it is shown that the density profile can exhibit a global modification leading to significant losses of alpha particles. This is due to the fact that, although the magnetic field does not exhibit any chaotic behaviour, the trajectories of alpha particles do, as revealed by their Poincaré maps. Such results are in qualitative agreement with past observations and simulations of energetic particles generated by neutral beam injection in TFTR, DIII-D and AUG tokamaks. In-depth analysis is carried out to characterize the impact of the tearing mode on the transport and losses of fusion-born alpha-particles with a realistic density profile. The impact of the amplitude is evidenced. Moreover, the effect of the island rotation frequency is assessed based on a detailed analysis of the linear resonances in phase-space, in agreement with the simulation results. Finally, the probability density function of the exit time has been computed and the transport of alpha particles has been found to be anomalous.

anomalous transport↗

DOE OSTI2022

An assessment of full-wave effects on Maxwellian lower-hybrid wave damping

Lower-hybrid current drive (LHCD) actuators are important components of modern day fusion experiments as well as proposed fusion reactors. However, simulations of LHCD often differ substantially from experimental results, and from each other, especially in the inferred power deposition profile shape. Here we investigate some possible causes of this discrepancy; ‘full-wave’ effects such as interference and diffraction, which are omitted from standard raytracing simulations and the breakdown of the raytracing near reflections and caustics. We compare raytracing simulations to state-of-the-art full-wave simulations using matched hot-plasma dielectric tensors in realistic tokamak scenarios for the first time. We show that differences between full-wave simulations and raytracing in previous work were primarily due to numerical and physical inconsistencies in the simulations, and we demonstrate that quantitative agreement between raytracing and converged full-wave simulations can be obtained in reactor relevant-scenarios and qualitative agreement can be obtained in situations with weak damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Analysing the effects of heating and gas puffing in Proto-MPEX helicon and auxiliary heated plasmas

The Material Plasma Exposure Experiment (MPEX) is being constructed at Oak Ridge National Laboratory to investigate critical fusion reactor issues, such as plasma–material interactions (PMIs) under reactor-relevant conditions and time scales. The linear device Proto-MPEX was used as a test bed to address anticipated research and development issues associated with heating scenarios and establish the physics basis for MPEX. The SOLPS-ITER code suite has been applied to understand plasma and neutral transport in Proto-MPEX and to increase confidence in predictive simulations for MPEX. Coupling between COMSOL and SOLPS is performed to implement a 2D electron heating profile of the helicon source. The simulations show reasonable agreement with the experimental data for plasma with helicon and auxiliary electron cyclotron heating (ECH). Both Bohm and constant diffusion (D$_\perp$: 0.5 m 2 s –1 and $\chi_\perp$: 1 m 2 s –1 ) simulations show similar levels of agreement with respect to the sparse experimental data available, assuming a few per cent impurity concentration. ECH significantly increases the target electron temperature, however, the target electron density is reduced compared to helicon-only heated plasmas due to an increase in flow velocity and radial losses. The simulations show that further increasing the ECH power results in an increase in the target electron density due to increased recycling flux and ionization. The results indicate that ECH significantly enhances the target heat flux, with ECH power of 50 kW increasing the target heat flux from 0.4 to 17 MW m –2 . It is found that a small amount of gas puffing (GP) near the target plate can further increase the target heat fluxes at the higher ECH power cases, but the target heat flux is reduced at higher GP conditions due to a significant reduction in electron temperature via radiation. ECH and GP scenarios can generate a higher target flux, facilitating improved PMI studies with more reactor-relevant plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Modeling turbulent impurity transport in the SOL of DIII-D with a reduced model

A novel impurity transport model that approximates SOL turbulence as a fluctuating poloidal electric field is shown to be an acceptable replacement for the traditional approach of assigning an arbitrary radial diffusion coefficient to the impurity ions. The model is implemented in the DIVIMP impurity transport code and applied to an L-Mode tungsten divertor experiment on DIII-D. The poloidal electric field is represented as fluctuating between ±1000 V m –1 based on previous measurements. The resulting intermittent v r = E θ × B T transport causes ions to transport both into the core as well as into the far-SOL. Simultaneous agreement with estimates of the W density just inside the separatrix as well as in the far-SOL is obtained (n W ~ 10 14 m –3 and n W ~ 10 12 m –3 , respectively). Prompt re-deposition of the W ions was necessary to obtain agreement (f redep ~ 99%). Here, we conclude that simulating impurity transport using a physics-based approximation for turbulence in the SOL, versus arbitrarily assigning diffusion coefficients, may enable better reactor scale predictions of core impurity contamination.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗

DOE OSTI2025

Assessing time-dependent temperature profile predictions using reduced transport models for high performing NSTX plasmas

Time-dependent, predictive simulations were performed with the 1.5D tokamak integrated modeling code TRANSP on a large set of well-analyzed, high performing discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate how well modern reduced transport models can reproduce experimentally observed temperature profiles in spherical tokamaks. Overall, it is found that simulations using the Multi-Mode Model (MMM) more consistently agree with the NSTX observations than those using the Trapped Gyro-Landau Fluid (TGLF) model, despite TGLF requiring orders of magnitude greater computational cost. When considering all examined discharges, MMM has median overpredictions of electron temperature (T e ) and ion temperature (T i ) profiles of 28% and 27%, respectively, relative to the experiment. TGLF overpredicts T e by 46%, with much larger variance than MMM, and underpredicts T i by 25%. As the ratio of kinetic to magnetic field pressure (β) is increased across NSTX discharges, TGLF predicts lower T e and significant flattening of the T i profile, conflicting with NSTX observations. When using an electrostatic version of TGLF, both T e and T i are substantially overpredicted, underscoring the importance of electromagnetic turbulence in the high β spherical tokamak regime. Additionally, calculations with neural net surrogate models for TGLF were performed outside of TRANSP with a time slice flux matching transport solver, finding better agreement with experiment than the TRANSP simulations, highlighting the impact of different transport solvers and simulation techniques. Furthermore, the reasonable agreement with experiment of temperature profiles predicted by MMM motivates a more detailed examination of the sensitivities of the TRANSP simulations with MMM to different NSTX plasma regimes in a companion paper (Lestz et al Plasma Phys. Control. Fusion 67 105030), in preparation for self-consistent, time-dependent predictive modeling of NSTX-U scenarios.

Lestz, J. B. [General Atomics, San Diego, CA (Unit↗

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 OSTI2026

Benchmark for two-dimensional large scale coherent structures in partially magnetized E × B plasmas—community collaboration & lessons learned

Low-temperature plasmas (LTPs) are essential to both fundamental scientific research and critical industrial applications. As in many areas of science, numerical simulations have become a vital tool for uncovering new physical phenomena and guiding technological development. Code benchmarking remains crucial for verifying implementations and evaluating performance. This work continues the Landmark benchmark initiative, a series specifically designed to support the verification of LTP codes. In this study, seventeen simulation codes from a collaborative community of nineteen international institutions modeled a partially magnetized E × B Penning discharge. The emergence of large scale coherent structures, or rotating plasma spokes, endows this configuration with an enormous range of time scales, making it particularly challenging to simulate. The codes showed excellent agreement on the rotation frequency of the spoke as well as key plasma properties, including time-averaged ion density, plasma potential, and electron temperature profiles. Achieving this level of agreement came with challenges, and we share lessons learned on how to conduct future benchmarking campaigns. Comparing code implementations, computational hardware, and simulation runtimes also revealed interesting trends, which are summarized with the aim of guiding future plasma simulation software development.

benchmarking↗

DOE OSTI2021

The impact of new d(p,γ)3 rates on Big Bang Nucleosynthesis

We consider the effect on Big Bang Nucleosynthesis (BBN) of new measurements of the d(p,γ) 3 He cross section by the LUNA Collaboration. These have an important effect on the primordial abundance of D/H which is also sensitive to the baryon density at the time of BBN. We have re-evaluated the thermal rate for this reaction, using a world average of cross section data, which we describe with model-independent polynomials; our results are in good agreement with a similar analysis by LUNA. We then perform a full likelihood analysis combining BBN and Planck cosmic microwave background (CMB) likelihood chains using the new rate combined with previous measurements and compare with the results using previous rates. Concordance between BBN and CMB measurements of the anisotropy spectrum using the old rates was excellent. The predicted deuterium abundance at the Planck value of the baryon density was (D/H) BBN+CMB old = (2.57 ± 0.13) × 10 –5 which can be compared with the value determined from quasar absorption systems (D/H) obs = (2.55 ± 0.03) × 10 –5 . Using the new rates we find (D/H) BBN+CMB = (2.51 ± 0.11) × 10 –5 . We thus find consistency among BBN theory, deuterium and 4He observations, and the CMB, when using reaction rates fit in our data-driven approach. We also find that the new reaction data tightens the constraints on the number of relativistic degrees of freedom during BBN, giving the effective number of light neutrino species N ν = 2.880 ± 0.144 in good agreement with the Standard Model of particle physics. Lastly, we note that the observed deuterium abundance continues to be more precise than the BBN+CMB prediction, whose error budget is now dominated by d(d,n) 3 He and d(d,p) 3 H.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Cosmological constraints from the cross-correlation of DESI Luminous Red Galaxies with CMB lensing from Planck PR4 and ACT DR6

Here, we infer the growth of large scale structure over the redshift range 0.4 ≲ z ≲ 1 from the cross-correlation of spectroscopically calibrated Luminous Red Galaxies (LRGs) selected from the Dark Energy Spectroscopic Instrument (DESI) legacy imaging survey with CMB lensing maps reconstructed from the latest Planck and ACT data. We adopt a hybrid effective field theory (HEFT) model that robustly regulates the cosmological information obtainable from smaller scales, such that our cosmological constraints are reliably derived from the (predominantly) linear regime. We perform an extensive set of bandpower- and parameter-level systematics checks to ensure the robustness of our results and to characterize the uniformity of the LRG sample. We demonstrate that our results are stable to a wide range of modeling assumptions, finding excellent agreement with a linear theory analysis performed on a restricted range of scales. From a tomographic analysis of the four LRG photometric redshift bins we find that the rate of structure growth is consistent with ΛCDM with an overall amplitude that is ≃ 5-7% lower than predicted by primary CMB measurements with modest (∼ 2σ) statistical significance. From the combined analysis of all four bins and their cross-correlations with Planck we obtain S 8 = 0.765 ± 0.023, which is less discrepant with primary CMB measurements than previous DESI LRG cross Planck CMB lensing results. From the cross-correlation with ACT we obtain S 8 = 0.790 +0.024 -0.027 , while when jointly analyzing Planck and ACT we find S 8 = 0.775 +0.019 -0.022 from our data alone and σ 8 = 0.772 +0.020 -0.023 with the addition of BAO data. These constraints are consistent with the latest Planck primary CMB analyses at the ≃ 1.6-2.2σ level, and are in excellent agreement with galaxy lensing surveys.

cosmological parameters from LSS↗

DOE OSTI2025

DESI 2024 VII: cosmological constraints from the full-shape modeling of clustering measurements

We present cosmological results from the measurement of clustering of galaxy, quasar and Lyman-α forest tracers from the first year of observations with the Dark Energy Spectroscopic Instrument (DESI Data Release 1). We adopt the full-shape (FS) modeling of the power spectrum, including the effects of redshift-space distortions, in an analysis which has been thoroughly validated in a series of supporting papers as summarised in [1]. We combine the full-shape information with DESI's DR1 constraints from the baryon acoustic oscillations (BAO) of these tracers. In the flat ΛCDM cosmological model, DESI (FS+BAO), combined with a baryon density prior from Big Bang Nucleosynthesis and a weak prior on the scalar spectral index, determines matter density to Ω m = 0.2962 ± 0.0095, and the amplitude of mass fluctuations to σ 8 = 0.842 ± 0.034. The addition of the cosmic microwave background (CMB) data tightens these constraints to Ω m = 0.3056 ± 0.0049 and σ 8 = 0.8121 ± 0.0053, while further addition of the joint clustering and lensing analysis from the Dark Energy Survey Year-3 (DESY3) data further improves these measurements, and leads to a 0.4% determination of the Hubble constant, H 0 = (68.40 ± 0.27) km s -1 Mpc -1 . In models with a time-varying dark energy equation of state parametrised by w 0 and w a , combinations of DESI (FS+BAO) with CMB and type Ia supernovae continue to show the preference, previously found in the DESI DR1 BAO analysis, for w 0 > -1 and w a < 0 with similar levels of significance. DESI data, in combination with the CMB, improve the upper limits on the sum of the neutrino masses relative to the case when only the DR1 BAO was available, giving ∑m ν < 0.071 eV at 95% confidence. We finally constrain deviations from general relativity represented by two modified gravity parameters. DESI (FS+BAO) data alone measure the parameter that controls the clustering of massive particles, μ 0 = 0.11 +0.45 -0.54 , in agreement with the zero value predicted by general relativity. The combination of DESI with the CMB and the clustering and lensing analysis from DESY3 constrains both modified-gravity parameters, giving μ 0 = 0.04 ± 0.22 and Σ 0 = 0.044 ± 0.047, again in agreement with general relativity.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Simulation of divertor heat flux width on EAST by BOUT++ transport code

The BOUT++ edge plasma transport code is applied to study the effects of neutral recycling, drifts and heating scheme on edge plasma profiles and the divertor heat flux width of two discharges of the experimental advanced superconducting tokamak (EAST) steady-state H-mode plasmas heated by low hybrid wave (LHW) and neutral beam (NB), respectively. Neutral recycling seems to have played an important role in the plasma density profile. The edge plasma density drops dramatically for the case w/o neutral recycling, while it can be sustained for the case w/ neutral recycling. Drifts are found to have significant influences on edge plasma profiles. Both the amplitude and width of the divertor heat flux are found to have increased a lot due to drifts. The simulated heat flux width w/ drifts for the two discharges shows reasonable agreement with the experiments. However, the width from the simulation and experiment for the LHW heated discharge is much larger than that of the NB heated discharge. Comparison with Goldston's drift-based model and more detailed analysis on the heat flux contributions from drifts versus turbulence show that drifts are the dominant factor in edge plasma transport for both LHW and NB heated discharges, while turbulence may have played a more important role in determining the heat flux width in LHW heated discharges than that in NB heated discharges. This may account for the larger heat flux width in the LHW heated discharge. The magnetic topology and the equilibrium change by the LHW power may also be a potential reason for the larger heat flux width in the LHW heated discharges, which still needs more evidence and studies to support it. Finally, additional SOLPS simulation w/o drifts for the two discharges show reasonable agreement with the counterpart from the BOUT++ simulation, suggesting the two are both suitable codes for EAST edge plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Turbulence driven widening of the near-SOL power width in ASDEX Upgrade H-Mode discharges

Divertor heat flux data from infra red (IR) from various tokamaks in H-Mode regime show that the power width, $\lambda_{q}$, scales approximately like 1.6 $(R/a) ρ_{s,pol}$. However, the IR based scaling comes with the restriction that only discharges were considered with attached divertor conditions. Here we set up a new data base covering low edge density data but also including high density plasmas reaching the H-mode density limit. We use Thomson-scattering to measure the electron temperature decay length which will set the near-SOL power width through parallel heat conduction. We present scaling results for the density, temperature and pressure separatrix fall-off lengths. As the principal result we derive a generalized power width scaling which reads as $\lambda_{q}\propto\,\rho_{s,pol}\cdot(1+2.1\alpha_{t}^{1.7})$ where $\alpha_{t}$ describes a normalized collisionality ($\alpha_{t} = 3\cdot\,10^{-18}\,R\,q^2\,n\,Z_{eff}T^{-2}$). The parameter $\alpha_{t}$ describes the relative importance of the interchange effect on drift-wave turbulence as proposed by Scott and is found to be closely linked to the diamagnetic parameter αd introduced by Rogers, Drake and Zeiler (RDZ). The edge operational space for ASDEX Upgrade displays a boundary qualitatively in agreement with the region labelled as density limit in the edge plasma phase space proposed by RDZ. The data base shows in the limit of low edge densities ($\alpha_{t}\simeq\,0.2$) accurate agreement to the IR based power width scaling. At elevated separatrix densities ($\alpha_{t}\simeq\,1$) the power width is broadened by a factor of up to about three accompanied by a reduction of the $\text{H}_{98,y2}$ confinement factor. Finally, plasmas with higher triangularity show higher values for $\text{H}_{98,y2}$ in the same range of normalized collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Collisionality driven turbulent particle transport changes in DIII-D H-mode plasmas

The results of the experimental dimensionless scan in this paper confirm that there is an increase in density peaking towards lower collisionality and that this can be partly linked to a shift in the turbulence regime from ITG towards TEM. However at the lowest collisionality, the changes in turbulence and transport are much more pronounced than expected from direct collisionality effect on the turbulence. In this paper, the collisionality, ν * is varied by a factor 5, while keeping ρ * , q , β , M , fixed. Additionally, a 3 Hz gas puff modulation is applied to modulate the electron density profile and extract the perturbed transport coefficients using two diagnostics. The transport analysis shows that the increase in density peaking at low ν * is linked to an increase in the inward particle pinch and not an increase in core fueling. These observations are not only in agreement with prior modeling scans of how turbulence changes as a function of collisionality and its impact upon the particle fluxes, but also with the multi-machine database (Fable E. et al 2010 Plasma Phys. Control. Fusion 52 015007) (Angioni C. et al 2003 Phys. Rev. Lett. 90 205003). The changes in turbulence across the collisionality scan were captured at large scale by the BES and at smaller scale by the DBS. A comparison with gradient-driven GENE simulations showed similar trends at both scales. Moreover, the changes observed in overall transport are in agreement with gradient-driven TGLF particle flux simulations. This indicates that TGLF/GENE when given the gradients as input, are able to reproduce the experimentally observed turbulence changes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗