Search NASA⌕ Search

SEARCH · Search NASA

Results for “Transition to turbulence”

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 199 records · Page 11

On the violation of the zeroth law of turbulence in space plasmas

The zeroth law of turbulence states that, for fixed energy input into large-scale motions, the statistical steady state of a turbulent system is independent of microphysical dissipation properties. This behaviour, which is fundamental to nearly all fluid-like systems from industrial processes to galaxies, occurs because nonlinear processes generate smaller and smaller scales in the flow, until the dissipation – no matter how small – can thermalise the energy input. Using direct numerical simulations and theoretical arguments, we show that in strongly magnetised plasma turbulence such as that recently observed by the Parker Solar Probe spacecraft, the zeroth law is routinely violated. Namely, when such turbulence is ‘imbalanced’ – when the large-scale energy input is dominated by Alfvénic perturbations propagating in one direction (the most common situation in space plasmas) – nonlinear conservation laws imply the existence of a ‘barrier’ at scales near the ion gyroradius. This causes energy to build up over time at large scales. The resulting magnetic-energy spectra bear a strong resemblance to those observed in situ , exhibiting a sharp, steep kinetic transition range above and around the ion-Larmor scale, with flattening at yet smaller scales. The effect thus offers a possible solution to the decade-long puzzle of the position and variability of ion-kinetic spectral breaks in plasma turbulence. The existence of the ‘barrier’ also suggests that, how a plasma is forced at large scales (the imbalance) may have a crucial influence on thermodynamic properties such as the ion-to-electron heating ratio.

Physics↗

Zonal magnetic fields regulate nonlinear edge-localized-mode dynamics via self-consistent force balance

Edge-localized modes (ELMs) eject intense bursts of heat and particles that threaten plasma-facing components in fusion reactors. Nonlinear full-torus BOUT++ simulations show that turbulence-driven zonal magnetic fields (ZMFs) play an essential role in nonlinear ELM evolution by maintaining self-consistent force balance. Zonal flows mitigate the initial crash through shear but do not prevent continued radial transport. When ZMFs are self-consistently included, turbulence-driven zonal currents modify the parallel current distribution and magnetic tension and are associated with a reduction of the axisymmetric (𝑛 = 0) perturbed radial force imbalance. This coincides with a transition from convective, bursty propagation to more localized, diffusive transport. Similar behavior is observed across the regimes considered, including both resistive-ballooning and peeling-ballooning cases. Finally, associated signatures, including radial electric field shear and parallel current redistribution, provide experimentally accessible diagnostics for present devices and ITER-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$ -pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al. , Phys. Plasmas , vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al. , J. Plasma Phys. , vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

Physics↗

Budgets for Decadal Variability in Pacific Ocean Heat Content

A slowdown in the rate of surface warming in the early 2000s led to renewed interest in the redistribution of ocean heat content (OHC) and its relationship with internal climate variability. In this work we use the Community Earth System Model version 1 to study the relationship between OHC and the interdecadal Pacific oscillation (IPO), a major mode of decadal sea surface temperature variability in the Pacific Ocean. By comparing the relative contributions of surface heat flux and ocean dynamics to changes in OHC for different phases of the IPO, we try to identify the underlying physical processes involved. Our results suggest that during IPO phase transitions, changes of 0–300-m OHC across the northern extratropical Pacific are positively contributed by both surface heat flux and oceanic heat transport. By contrast, oceanic heat transport appears to drive the OHC changes in equatorial Pacific whereas surface heat flux acts as a damping term. During a positive IPO phase, weakened wind-driven circulation acts to increase the OHC in the equatorial Pacific while the enhanced evaporation acts to damp OHC anomalies. In the Kuroshio–Oyashio Extension region, a dipole anomaly of zonal heat advection amplifies an OHC dipole anomaly that moves eastward, while strong turbulent heat fluxes act to dampen this OHC anomaly. In the northern subtropical Pacific, both the wind-driven evaporation change and the change of zonal heat advection along Kuroshio Extension contribute to the OHC change during phase transition. For the northern subpolar Pacific, both surface heat flux and enhanced meridional advection contribute to the positive OHC anomalies during the positive IPO phase.

54 ENVIRONMENTAL SCIENCES↗

Hypersonic Fluid–Structure Interaction on a Cone–Slice–Ramp Geometry

In this study, fluid–structure interactions were measured between a representative control surface and the hypersonic flow deflected by it. The control surface is simplified as a spanwise finite ramp placed on a longitudinal slice of a cone. The front surface of the ramp contains a thin panel designed to respond to the unsteady fluid loading arising from the shock-wave/boundary-layer interactions. Experiments were conducted at Mach 5 and Mach 8 with ramps of different angles. High-speed schlieren captured the unsteady flow dynamics and accelerometers behind the thin panel measured its structural response. Panel vibrations were dominated by natural modes that were excited by the broadband aerodynamic fluctuations arising in the flowfield. However, increased structural response was observed in two distinct flow regimes: 1) attached or small separation interactions, where the transitional regime induced the strongest panel fluctuations. This was in agreement with the observation of increased convective undulations or bulges in the separation shock generated by the passage of turbulent spots, and 2) large separated interactions, where shear layer flapping in the laminar regime produced strong panel response at the flapping frequency. In addition, panel heating during the experiment caused a downward shift in its natural mode frequencies.

33 ADVANCED PROPULSION SYSTEMS↗

Investigating the Impacts of Daytime Boundary Layer Clouds on Surface Energy Fluxes and Boundary Layer Structure During CHEESEHEAD19

Studies of land-atmosphere interactions under a clear sky and low cumulus cloud conditions are common from long-term observatories like at the southern great plains. How well the relationships and responses of surface radiative and turbulent heat fluxes determined from these investigations hold for more heterogeneous surfaces in other climate regimes, however, is uncertain. In this study, detailed observations of the surface energy budget and daytime boundary layer properties are analyzed using measurements from the Chequamegon Heterogenous Ecosystem Energy-Balance Study Enabled by a High-Density Extensive Array of Detectors 2019 (CHEESEHEAD19) field campaign, July-October 2019, across a heterogeneous forested landscape of northern Wisconsin. A cloud regime framework is employed to classify consecutive periods of clear skies from lower atmosphere stratiform and cumulus clouds. A seasonal transition from low cumulus to low stratiform periods occurred, together with a diurnal pattern in cloudy or clear sky period dominance. Radiative forcing was highly dependent on sky conditions, leading to changes in the redistribution efficiency of radiative energy by the surface turbulent heat fluxes. During CHEESEHEAD19, small Bowen ratios dominated with daytime latent heat fluxes three times as large as sensible heat fluxes for all sky conditions studied; the forested region, therefore, falls within an energy-limited regime. The depth of the daytime mixed layer depended upon the sky condition and thermodynamic setting; deeper mixed layers occurred during periods of low cumulus and not clear skies. Profiles of vertical velocity were found to have enhanced variance under low cumulus compared to clear sky periods, suggesting potential for cloud feedbacks on boundary layer structure and surface energy fluxes.

54 ENVIRONMENTAL SCIENCES↗

Scrape-off layer stability, turbulence and transport in MAST-U (Final Report)

U.S. Department of Energy (DOE) supported work at Lodestar Research Corporation under grant DE-SC0019270 is summarized in this final report. The research conducted under this grant explored the properties of the edge plasma in spherical torus magnetic fusion devices. Harnessing the energy that is released in fusion reactions would provide a safe and abundant source of power to meet the growing energy needs of the world population. The edge and scrape-off layer (SOL) region of fusion plasmas is the region where the hot core plasma transitions to a cooler surrounding boundary plasma. The properties of this boundary plasma determine how interactions with material surfaces will impact both the overall device performance and the lifetime of crucial components such as the walls of the device chamber and other hardware. The research carried out studied filamentary plasma instabilities and turbulent transport in the edge and SOL of the MAST U spherical torus.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Statistical analysis of magnetic divertor configuration influence on H-mode transitions

DIII-D plasmas are compared for two upper divertor configurations: with the outer strike point on the small angle slot (SAS) divertor target and with the outer strike point on the horizontal divertor target (HT). Scanning the vertical distance between the magnetic null point and the divertor target over a range 0.10–0.16 m is shown to increase the threshold power, $P_\mathrm{th}$, and edge plasma power, $P_\mathrm{Loss},$ for the low-to-high confinement (L–H) and H–L transitions respectively, by up to a factor of 1.4. The X-point height scans were performed at three L-mode core plasma line average electron densities, $\bar{n}_\mathrm{e} = $ 1.2, 2.2 and 3.6 $\times 10^{19}\,\mathrm{m}^{-3}$, to investigate the density dependence of divertor magnetic configuration influence on $P_\mathrm{th}$. The X-point height, $Z_\textrm{x-pt}$, was further extended across the range 0.16–0.22 m with the more open HT divertor configuration, for which a clear decrease in $P_\mathrm{th}$ with increasing $Z_\textrm{x-pt}$ is observed. The dependence of $P_\textrm{th}$ on divertor magnetic geometry is further investigated using a time-dependent probability density function (PDF) model and information geometry to elucidate the roles played by pedestal plasma turbulence and perpendicular velocity flows. The degree of stochasticity of the plasma turbulence is observed to be sensitive to the plasma heating rate. The calculated square of the information rate shows changes in the relative density fluctuations and perpendicular velocity PDFs begin 2–5 ms prior to the L–H transition for three plasmas; providing a crucial measurement of the dynamic timescale of external transport barrier formation. Additionally, both information length and rate provide potential predictors of the L–H transition for these plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On FIRE mode in KSTAR

We report on the status of Fast Ion Regulated Enhancement (FIRE) mode experiments in the Korea Superconducting Tokamak Advanced Research. This regime is being developed for high-performance, steady-state operation which features a stationary ion internal transport barrier, enabling a central ion temperature approaching 10 keV to be sustained for up to 50 s, without the need for delicate profile control and with no significant impurity accumulation. As its key novelty lies in the significant contribution of fast ions that stabilize core turbulence, the regime has been named FIRE mode. To achieve this regime, neutral beam injection is applied at moderate power levels near the L–H power threshold, while maintaining low plasma density to avoid the L–H transition. The scenario is typically established in diverted magnetic configurations. The core features of FIRE mode were investigated through power balance analysis and fluctuation measurements, revealing a clear transport bifurcation in the ion channel. At the plasma edge, FIRE mode occasionally exhibits I-mode characteristics, particularly in unfavorable magnetic null configurations with q 95 ∼ 4, including the presence of weakly coherent modes. In terms of MHD activity, sawtooth oscillations are observed but appear to be stabilized during the high-performance phase. Fast ion-driven Alfvénic eigenmodes (AEs), indicated by strong frequency chirping near 200 kHz, are also observed. Additionally, lower-frequency MHD activities, distinct from the fast ion-driven AEs, are present and have some influence on plasma performance. The enhancement of core confinement is primarily attributed to fast ion effects, which were evaluated with respect to dilution, alpha stabilization, and resonant interactions with turbulence using gyrokinetic analyses. Among these, the dilution effect was found to be the most dominant. The characteristics of FIRE mode were compared with those of other hot ion plasma scenarios, such as supershot and hot ion mode. While they share many similarities, FIRE mode is distinguished by the accessibility to conditions with T i ≈ T e , presence of I-mode edge features, and its long-duration sustainment. Predictive simulations of FIRE mode were performed using integrated transport modeling with TRIASSIC, employing the TGLF anomalous transport model. These simulations confirmed the critical role of fast ions in achieving this regime. The future prospect of FIRE mode for application in fusion reactors is discussed, with an emphasis on possibly extending the regime to higher density operation.

FIRE mode↗

Turbulence-induced bias in time-averaged laser absorption tomography of correlated concentration and temperature fields with a first-order correction

The influence of correlated scalar fluctuations on time-averaged laser absorption tomography measurements of temperature and species in a piloted turbulent premixed flame was examined using a coupled spectroscopic and fluid-dynamic analysis. To understand bias associated with turbulence, spatio-temporally resolved temperature and species mole fraction profiles predicted by large eddy simulations (LES) were used to synthetically generate time-resolved line-of-sight absorption measurements at short time scales (microsecond) to reflect the unsteady nature of a canonical jet burner across various transverse measurement planes. Inversion methods were employed on the time-averaged line-of-sight data to produce radially-resolved temperature and mole fraction profiles, analogous to those produced by laser absorption tomography performed on a time-averaged axisymmetric flowfield. It is shown that bias in the measurements compared to true time-averaged scalar fields is a function primarily of temperature dependence in absorptivity and non-zero correlation between temperature and species concentration scalars. Finally, a first-order correction to tomography measurements is proposed to account for the bias based on estimated correlations and the known spectroscopic parameters of the probed absorption transitions.

42 ENGINEERING↗

Effects of a magnetic field on hot electron generation from laser-plasma instabilities

We performed two-dimensional simulations of turbulent laser-plasma instabilities in the presence and absence of external magnetic fields using the Laser Plasma Simulation Environment (LPSE) code. The results demonstrate that, in the presence of a magnetic field, the transition from ballistic to gyrating electron motion enhances the energy transfer from electron plasma waves to the electron population. Although stronger magnetic fields produce a larger population of hot electrons, these electrons tend to remain confined near the quarter-critical density, where the instabilities also localize, thereby reducing the potential for hot electron transport deeper into the target. Additionally, we present a scaling analysis that quantifies hot electron generation as a function of plasma electron temperature, density scale length, and applied magnetic field strength. These findings may have important applications for controlling hot electron flux and mitigating preheat in inertial confinement fusion targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Long‐Term Observations of Turbulence Vertical Velocity Spectra in a Convective Mixed Layer: Dependence on Land‐Surface Forcing in the U.S. Southern Great Plains

Doppler lidar vertical velocity retrievals were analyzed for the scale and structure of mixed-layer turbulence over a 7-year period, on fair-weather warm season days (442 cases) in the U.S. Southern Great Plains. Based on the hypothesized influence of land-surface forcing and updraft size on convective boundary layer clouds, a spectral analysis was performed to quantify effects of surface forcing (surface buoyancy flux, friction velocity, and evaporative fraction) on turbulence scale. Significant (order-of-magnitude) variations in spectral density were found in the energy-production subrange and mesoscale regimes. Integral scale decreased with increasing buoyancy flux (and Monin-Obukhov stability parameter), while spectral density in the energy-production subrange increased, implying a transition to buoyancy-driven cellular structures with narrower updrafts. The influence of stability parameter was limited to the neutral to convective transition, and could not explain the wide variation in spectral density in the mesoscale regime. However, high friction velocity (u* > 0.5 m s –1 ) was associated with larger integral scale (updraft width), and a greater portion of spectral density in the mesoscale regime. Lidar profiles and radar reflectivity for individual cases revealed evidence of roll and wave-like structures contributing to mesoscale variability on high friction velocity days, suggesting shear instabilities on days with wetter land surface conditions and smaller buoyancy flux. The role of friction velocity emphasizes the multivariate nature of surface influences on turbulence, and implies that variance-based turbulence closures may not be adequate for capturing effects of surface forcing on updraft size.

54 ENVIRONMENTAL SCIENCES↗

On the stability of a pair of vortex rings

The growth of perturbations subject to the Crow instability along two vortex rings of equal and opposite circulation undergoing a head-on collision is examined. Unlike the planar case for semi-infinite line vortices, the zero-order geometry of the flow (i.e. the ring radius, core thickness and separation distance) and by extension the growth rates of perturbations vary in time. The governing equations are therefore temporally integrated to characterize the perturbation spectrum. The analysis, which considers the effects of ring curvature and the distribution of vorticity within the vortex cores, explains several key flow features observed in experiments. First, the zero-order motion of the rings is accurately reproduced. Next, the predicted emergent wavenumber, which sets the number of secondary vortex structures emerging after the cores come into contact, agrees with experiments, including the observed increase in the number of secondary structures with increasing Reynolds number. Finally, the analysis predicts an abrupt transition at a critical Reynolds number to a regime dominated by a higher-frequency, faster-growing instability mode that may be consistent with the experimentally observed rapid generation of a turbulent puff following the collision of rings at high Reynolds numbers.

42 ENGINEERING↗

A crisis for the verification and validation of turbulence simulations

Three algorithms have been proposed for solution of the Rayleigh–Taylor turbulent mixing problem. They are based upon three different physical principles governing the Euler equations for fluid flow. The principles serve to select the physically relevant solution from among many nonunique solutions. The admissibility principle is in dispute. The three different algorithms, expressing the three physical admissibility principles can be formulated in terms of the three energy dissipation rates or the entropy production rates, as selected by the size of the sub grid scale coefficients. These have maximal values or less than maximal values. The resulting solutions are markedly different. We find strong validation evidence that supports the maximum rate principle, based on a review of prior results and on new results presented here. We review experimental data used for validation and sufficient to discriminate among the three. We present a new analysis of this data. We show that the hypothesized long wave length perturbations in the initial conditions are not significant, so that validation can be based on this data in a straight forward manner. One of the algorithms is labeled direct numerical simulation, but is not, and as a consequence, the two algorithms with less than maximal SGS coefficients are variants of one another. Finally recommendations for the numerical modeling of the deflagration to detonation transition in type Ia supernova are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

The Turbulent Stress Spectrum in the Inertial and Subinertial Ranges

For velocity and magnetic fields, the turbulent pressure and, more generally, the squared fields such as the components of the turbulent stress tensor, play important roles in astrophysics. For both one and three dimensions, we derive the equations relating the energy spectra of the fields to the spectra of their squares. We solve the resulting integrals numerically and show that for turbulent energy spectra of Kolmogorov type, the spectral slope of the stress spectrum is also of Kolmogorov type. For shallower turbulence spectra, the slope of the stress spectrum quickly approaches that of white noise, regardless of how blue the spectrum of the field is. For fully helical fields, the stress spectrum is elevated by about a factor of two in the subinertial range, while that in the inertial range remains unchanged. We discuss possible implications for understanding the spectrum of primordial gravitational waves from causally generated magnetic fields during cosmological phase transitions in the early universe. We also discuss potential diagnostic applications to the interstellar medium, where polarization and scintillation measurements characterize the square of the magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗

Relative alignment between magnetic fields and molecular gas structure in molecular clouds

Here, we compare the structure of synthetic dust polarization with synthetic molecular line emission from radiative transfer calculations using a three-dimensional, turbulent collapsing-cloud magnetohydrodynamics simulation. The histogram of relative orientation (HRO) technique and the projected Rayleigh statistic (PRS) are considered. In our trans-Alfvénic (more strongly magnetized) simulation, there is a transition to perpendicular alignment at densities above ~4 × 10 3 cm –3 . This transition is recovered in most of our synthetic observations of optically thin molecular tracers; however, for 12 CO it does not occur and the PRS remains in parallel alignment across the whole observer space. We calculate the physical depth of the optical depth τ = 1 surface and find that for 12 CO it is largely located in front of the cloud mid-plane, suggesting that 12 CO is too optically thick and instead mainly probes low-volume density gas. In our super-Alfvénic simulation, the magnetic field becomes significantly more tangled, and all observed tracers tend towards no preference for perpendicular or parallel alignment. An observable difference in alignment between optically thin and optically thick tracers may indicate the presence of a dynamically important magnetic field, though there is some degeneracy with viewing angle. We convolve our data with a Gaussian beam and compare it with HRO results of the Vela C molecular cloud. We find good agreement between these results and our sub-Alfvénic simulations when viewed with the magnetic field in the plane of the sky (especially when sensitivity limitations are considered), though the observations are also consistent with an intermediately inclined magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Hyperdiffusion of dust particles in a turbulent tokamak plasma

In this work, the effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number St∝a d /n 0 , with a d the dust radius and n 0 the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories Δr 2 is shown to depend on the ratio Ku/St, with Ku∝u rms the Kubo number and urms the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting Δr 2 ∝(Ku/St) 2 t 3 , or Δr 2 ∝(u rms n 0 /a d ) 2 t 3 . As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible and drives the transition from a hyperdiffusive regime toward a ballistic one.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗