Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence Stabilization”

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 253 records · Page 14

Implementation and Evaluation of Physics-Driven Dynamic Entrainment-Mixing Parameterization in a Climate Model and Its Impact on Low-Cloud Simulation

The turbulent entrainment-mixing process in the Community Earth System Model version 1.2 (CESM1.2) is assumed to follow the extremely inhomogeneous entrainment-mixing. However, different entrainment-mixing scenarios can occur in real clouds. To address this deficiency, a unifying parameterization that represents different entrainment-mixing processes is implemented and evaluated in CESM1.2. The results indicate that the homogeneous mixing degree values simulated by the new parameterization in CESM1.2 are predominantly greater than 50%, suggesting a tendency toward homogeneous mixing. Compared to the extremely inhomogeneous mixing mechanism, the new parameterization increases the cloud droplet number concentration (Nc). More importantly, the new parameterization improves low-cloud fraction (CLDLOW) simulation in Northwest Pacific (NWP) and Southeast Pacific (SEP) regions, with relative improvements of 2.95% and 4.17%, respectively. Furthermore, the improvements reach up to 44.6% and 16.2% in the NWP and SEP regions, respectively, when considering the relationship between N c and CLDLOW. Further analysis reveals that the new parameterization enhances cloud optical depth, longwave radiative cooling effect, net condensation rate, cloud water mixing ratio, lower-troposphere stability, and CLDLOW by increasing N c . Additionally, these results underscore the importance of improving entrainment-mixing parameterization in climate models.

54 ENVIRONMENTAL SCIENCES↗

Tools Assessing Performance (FY2022 Quarter 4 Report)

We were finally successful in acquiring TTU’s 200-meter tower data for the duration of the Isolated Building Wake Experiment (IBWE). This data will provide an independent measurement of the undisturbed wind profile including heights far above those that are typical of meteorological sensor networks. These measurements will characterize the mean and turbulence throughout the atmospheric surface layer, which will be particularly useful for validation of highfidelity computational fluid dynamics models. This dataset can be shared freely among the TAP team but not outside of it as TTU sells this data to the wind power industry. The journal article overview has been submitted to LANL’s publication review and will be submitted to the Journal of Wind Engineering and Industrial Aerodynamics once the review is complete. We are nearing completion on another two articles. The first discusses analyses of the evolution of the measurements of mean wind speed along the centerline of the building as a function of atmospheric stability. Rather than using individual time steps, all of the wind speed measurements are normalized using the upwind measurement at 2 building heights (H) and filtered by wind direction and atmospheric stability. The second discusses the development of the diffusive wake model and uses the averaged profiles from the first article to validate the behavior of the model. The portion of the section below on the diffusive wake model discussing comparisons with the averaged profiles, will be included in the validation section of this paper.

17 WIND ENERGY↗

Development of an Ultrahigh-bandwidth Phase Contrast Imaging System for detection of electron scale turbulence and Gigahertz Radio-Frequency Waves

The study of waves and turbulence is vital to the development of future reactor-grade plasma devices developed in the quest for fusion energy. These fluctuations are responsible for moving heat and particles across the magnetic field, and a predictive understanding of them is needed to achieve the density and temperature required to sustain a plasma fusion reaction. While many techniques have been developed for measuring waves and fluctuations, every measurement method has limitations. There are relatively few techniques for measuring very high frequency fluctuations, such as radio frequency waves injected to heat the plasma, unstable waves driven by suprathermal particles, or short wavelength electrostatic waves driven by electron temperature or density gradients. The present project builds upon the proven phase contrast imaging (PCI) technique to extend the response of the diagnostic by orders of magnitude in frequency and almost a factor of ten in spatial resolution. PCI provides a measurement of electron density based on small angle scattering of a CO-2 laser beam by using optical techniques to render a phase shift as an intensity change on a detector. Due to the telecommunications revolution, technological development by manufacturers has focused on components in the near -infrared, so that high- quality lasers and detectors at 1.55 µm are readily available. Shifting PCI design to a new, shorter wavelength has numerous advantages and challenges. Similar detector performance is available with room-temperature arrays with GHz bandwidth, while the detector arrays for 10.6 µm required liquid nitrogen cooling and were therefore limited to a bandwidth of about 1 MHz. Shifting to a shorter wavelength reduces the angle at which the laser beam scatters off of plasma waves, allowing more such scattered components to pass through the aperture of the vacuum vessel port and be collected by the PCI, which increases the spatial resolution. Concomitant with these benefits, various questions of performance arise. At shorter laser wave- length, the sensitivity to mirror and lens quality is increased, the contribution of the laser to the overall system noise is increased, and the sensitivity to vibrations is increased. The custom optical components at the heart of the PCI technique were required to be properly scaled for the shorter wavelength, so fabrication technologies needed to be explored. This project was designed to show that a low noise, high response PCI system at 1.55 µm could be constructed and operated, and then to quantify potential issues to allow extrapolation to a full-size production system providing physics measurements on a large plasma device. The first stage, producing the custom optical component called a Phase Plate, was successfully achieved using two methods. First, an easily reproduced masking and coating technique was able to produce good phase plates with the required parameters. Second, a nanofabrication technique was found to produce extremely high quality phase plates at a competitive price. The PCI constructed with the new phase plates and 1.55 µm laser was found to provide excellent wavelength measurements with the theoretically expected response. The sensitivity to optical surface quality was found to be in line with previous measurements at 10.6 µm. The observed signal-to-noise ratio was similar to the theoretically expected value. The effect of vibrations on PCI was studied with the first measurement of the effect of beam motion on PCI response and comparison to theory, allowing for a quantitative prediction of the effect of vibrations on a production PCI system and the requirements for improved beam stabilization. PCI is an extremely cost-effective method to provide a low noise, absolutely calibrated measurement of plasma fluctuations across a wide spatial scale. This project has shown that a 1.55 µm PCI using modern techniques and components is less expensive than the 10.6 µm PCI of a few years ago, with the largest savings in phase plate fabrication and the infrared detector array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Navigating the Noise: Bringing Clarity to ML Parameterization Design With O $\boldsymbol{\mathcal{O}}$(100) Ensembles

Abstract Machine‐learning (ML) parameterizations of subgrid processes (here of turbulence, convection, and radiation) may one day replace conventional parameterizations by emulating high‐resolution physics without the cost of explicit simulation. However, uncertainty about the relationship between offline and online performance (i.e., when integrated with a large‐scale general circulation model) hinders their development. Much of this uncertainty stems from limited sampling of the noisy, emergent effects of upstream ML design decisions on downstream online hybrid simulation. Our work rectifies the sampling issue via the construction of a semi‐automated, end‐to‐end pipeline for size ensembles of hybrid simulations, revealing important nuances in how systematic reductions in offline error manifest in changes to online error and online stability. For example, removing dropout and switching from a Mean Squared Error to a Mean Absolute Error loss both reduce offline error, but they have opposite effects on online error and online stability. Other design decisions, like incorporating memory, converting moisture input from specific humidity to relative humidity, using batch normalization, and training on multiple climates do not come with any such compromises. Finally, we show that ensemble sizes of may be necessary to reliably detect causally relevant differences online. By enabling rapid online experimentation at scale, we can empirically settle debates regarding subgrid ML parameterization design that would have otherwise remained unresolved in the noise.

Lin, Jerry [Department of Earth System Sciences Un↗

Time-dependent experimental identification of inter-ELM microtearing modes in the tokamak edge on DIII-D

In a series of discharges on the DIII-D tokamak, fast vertical plasma jogs are used to induce current perturbations in the steep gradient region of the H-mode edge. These current perturbations directly influence the edge q profile, decoupling the resonant location and instability drive of pedestal-localized microtearing modes (MTMs.) By exploiting this effect, we develop and apply a new experimental technique to track the dynamical frequency evolution of MTMs in the pedestal region, providing a compelling validation of the MTM model. The frequency of potential MTMs is calculated as the Doppler-shifted electron diamagnetic frequency at rational q=m/n surfaces, showing remarkable agreement with chirped frequency behavior of n=3, 4 and 5 modes detected with fast magnetics. Data is collected throughout multiple ELM cycles in order to build robust statistics describing the time-dependent frequency evolution of MTMs, which can be explained by examining the recovery of pedestal gradients after an ELM event. MTMs have a dominant transport contribution in the electron thermal channel, so the presented results indicate that reduced models of pedestal transport must be electromagnetic in nature and constructed with accurate calculations of MTM stability; inclusion of this physics is essential for accurate predictions of the electron temperature pedestal profile. Supporting measurements of mode saturation, propagation direction and transport fingerprints are made to support the dynamic frequency determination, unambiguously and experimentally identifying MTMs in the pedestal region of DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Basal Melting and Oceanic Observations Beneath Central Fimbulisen, East Antarctica

Abstract Basal melting of ice shelves is fundamental to Antarctic ice sheet mass loss, yet direct observations remain sparse. We present the first year‐round melt record (2017–2021) from a phase‐sensitive radar on Fimbulisen, one of the fastest flowing ice shelves in Dronning Maud Land, East Antarctica. The observed long‐term mean ablation rate at 350 m depth below the central ice shelf was 1.0 ± 0.5 m yr −1 , marked by substantial sub‐weekly variability ranging from 0.4 to 3.5 m yr −1 . 36‐h filtered basal melt rate fluctuations closely align with ocean velocity. On seasonal time scales, melt rates peak during austral spring to autumn (September–March), driven by both elevated ocean velocities and thermal driving near the base. The combined effect of thermal driving and current speed explains the majority of the melt rate variability ( r = 0.84), highlighting the dominant role of shear‐driven turbulence. This relationship enables parameterization of melt rates for the decade‐long ocean record (2010–2021), although deviations appear under low and high forcing conditions. Both observed and parameterized melt rates show similar yearly mean magnitudes compared to satellite‐derived melt rates but with a tenfold lower seasonal amplitude and a 3‐month delay in seasonality. These detailed concurrent ice–ocean observations provide essential validation data for remote sensing and numerical models that aim to quantify and project ice‐shelf response to a change in ocean forcing. In situ measurements and continued monitoring are crucial for accurately assessing and modeling future basal melt rates, and for understanding the complex dynamics driving ice‐shelf stability and sea‐level change.

54 ENVIRONMENTAL SCIENCES↗

Achieving designed texture and flows in bulk active nematics using optimal control theory

Being intrinsically nonequilibrium, active materials can potentially perform functions that would be thermodynamically forbidden in passive materials. However, active systems have diverse local attractors that correspond to distinct dynamical states, many of which exhibit chaotic turbulent-like dynamics and thus cannot perform work or useful functions. Designing such a system to choose a specific dynamical state is a formidable challenge. Motivated by recent advances enabling optogenetic control of experimental active materials, we describe an optimal control theory framework that identifies a spatiotemporal sequence of light-generated activity that drives an active nematic system toward a prescribed dynamical steady state. Active nematics are unstable to spontaneous defect proliferation and chaotic streaming dynamics in the absence of control. We demonstrate that optimal control theory can compute activity fields that redirect the dynamics into a variety of alternative dynamical programs and functions. This includes dynamically reconfiguring between states, selecting and stabilizing emergent behaviors that do not correspond to attractors, and are hence unstable in the uncontrolled system. Furthermore, our results provide a roadmap to leverage optical control methods to rationally design structure, dynamics, and function in a wide variety of active materials.

Complex systems theory↗

Nek5000 developments in support of industry and the NRC

This year, the Nuclear Energy Advanced Modeling Simulation program (NEAMS) thermal-hydraulics verification and validation (V&V) work has focused in three areas of Nek5000 V&V-driven development. First, in a close collaborative effort with the U. S. Nuclear Regulatory Commission (NRC) staff, we have continued V&V efforts for the HYMERES-2 project using the OECD/NEA sponsored testing in the PSI PANDA facility. This year’s focus of ANL-NRC collaboration involves Nek5000 setups and validation for a range of problems relevant to and including the HYMERES-2 benchmark from PSI. The primary outcome of this year efforts is a more efficient geometry and inlet modeling simplification after a careful sensitivity study of the inlet profiles and pipe geometries. The resulting modeling choice of a short recycling/fully-developed turbulent inlet is within the experimental uncertainty estimate. This finding simplifies the next step of the cross-V&V HYMERES-2 project. In addition, the ANL team continue to provide assistance to the NRC staff in the form of Nek5000 application support in general and on the use of the HPC platforms of ALCF and INL in particular. This supports the NRC’s assessment of Nek5000 for use with the NRC Blue CRAB code suite. Second, we have implemented and tested more robust model of URANS, namely the k – τ model, a variant of the k-ω model, along with other improvements to RANS Nek5000 modeling in general. Because of its demonstrated robustness and stability, the k – τ model is the only RANS model that has been implemented in the new GPU version of the Nek5000 code, nekRS. Lastly, we report the initial implementation of Jacobian-free Newton Krylov approach to the direct Newton method for steady fluid solvers aimed at acceleration of RANS modeling and at IC improvement for LES campaigns. Also leveraging the Exascale Computing Project (ECP) ANL/CEED & SMR team’s software development effort to support NEAMS problems at large scale of the advanced computing architectures, NekRS, a GPU variant of Nek5000, built on top of kernels from libParanumal using OCCA for portability, has been successfully run on the full system of Summit (4608 nodes, 27648 GPUs).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Flow-induced oscillation patterns for two tandem cylinders with turbulence stimulation and variable stiffness and damping

Herein, the oscillation patterns of two 1-DOF cylinder oscillators, undergoing VIV and galloping, are investigated in a free surface water channel for 3.2 × 104 ≤ Re ≤ 1.2 × 10 5 . The cylinders are arranged in tandem and supported by springs for a range of different spring stiffness and damping parameters. The efficiency of a current energy converter (CEC), based on flow-induced oscillation (FIO) of multiple cylinders in tandem, is critically related to the cylinder oscillation patterns. Due to the limited number of studies on the FIO for multiple cylindrical or prismatic bodies, their oscillation patterns have neither been identified nor classified. The surfaces of the cylinders are modified by turbulence stimulation to enhance FIO. Three different center-to-center spacing, five stiffness, and six damping ratios, for a total of 90 sets of experiments were conducted. The current velocity range is from 0.34 m/s to 1.32 m/s. From more than 2000 tests, five major patterns, nine sub-patterns are identified and classified. The patterns are defined based on the amplitudes, frequencies and the phase angle differences between the two cylinders. The characteristics and mechanics of each oscillation pattern are explored and explained from the perspective of fluid-structure interaction (FSI). By systematically varying the parameters, the underlying hydrodynamic mechanisms, including the coupling level between vortices and cylinders, the shielding effect, the wake effect, and the stability states are revealed. Few preliminary observations on the connection between oscillation pattern and harnessed power by the tandem cylinders are reported.

42 ENGINEERING↗

Dependence of Cloud Macrophysical Properties and Phase Distributions on Environmental Conditions Over the North Atlantic and Southern Ocean: Results From COMBLE and MARCUS

The accurate representation of Cold Air Outbreaks (CAOs) and affiliated mixed-phase boundary layer (BL) clouds in models is challenging. How BL cloud properties evolve during CAOs and their dependence on meteorological conditions is not well understood but is important for the simulation of Earth's energy budgets. Here the properties of polar BL clouds over the North Atlantic (NA) and Southern Ocean (SO) are compared using observations from the Measurements of Aerosol Radiation and CloUds over the SO (MARCUS) and CAOs in the Marine BL Experiment (COMBLE) conducted over the NA. MARCUS observations show a stronger BL inversion than COMBLE, with a higher mean EIS (estimated inversion strength)/LTS (lower tropospheric stability) of -0.03 K/13 K compared to COMBLE’s -3.2 K/9.3 K. 39% of CAOs observed during COMBLE were intense with M > 5 K, while MARCUS only had 1.3%. 78%/72% of clouds sampled in CAOs during COMBLE/MARCUS had cloud top heights <4 km. The mean BL cloud top height was over 400 m higher, and the BL was over 500 m deeper for M of 10 K compared to 0 K for both regions. MARCUS observed a 27% moister BL structure than COMBLE when M > 5 K due to stronger BL inversion trapping more moisture within the BL. Under the same LTS, EIS, and M conditions, MARCUS observed a 12% drier BL structure, and clouds were 46% more turbulent than COMBLE. During CAOs, 54% of single-layer BL clouds sampled during MARCUS had liquid-dominated bases compared to 39% during COMBLE.

54 ENVIRONMENTAL SCIENCES↗

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗

Similarity of the turbulent kinetic energy dissipation rate distribution in the upper mixed layer of the tropical Indian Ocean

Turbulence within the upper ocean mixed layer plays a key role in various physical, biological, and chemical processes. Between September and November 2011, a dataset of 570 vertical profiles of the turbulent kinetic energy (TKE) dissipation rate, as well as conventional hydrological and meteorological data, were collected in the upper layer of the tropical Indian Ocean. These data were used to statistically analyze the vertical distribution of the TKE dissipation rate in the mixed layer. The arithmetic-mean method made the statistical TKE dissipation rate profile more scattered than the median and geometric-mean methods. The statistical TKE dissipation rate were respectively scaled by the surface buoyancy flux and the TKE dissipation rate at the mixed-layer base. It was found that the TKE dissipation rate scaled by that at the mixed-layer base exhibited better similarity characteristics than that scaled by the surface buoyancy flux, whether the stability parameter D/|L MO | was greater or less than 10, indicating that the TKE dissipation rate at the mixed-layer base is a better characteristic scaling parameter for reflecting the intrinsic structure of the TKE dissipation rate in the mixed layer, where D and L MO are respectively the mixed-layer thickness and the Monin-Obukhov length scale. The parameterization of the TKE dissipation rate at the mixed-layer base on the shear-driven dissipation rate and the surface buoyancy flux was further explored. It was found that the TKE dissipation rate at the mixed-layer base could be well fitted by a linear combination of three terms: the wind-shear-driven dissipation rate, the surface buoyancy flux, and a simple nonlinear coupling term of these two.

54 ENVIRONMENTAL SCIENCES↗

A new discovery of edge localized modes suppression using ICRH

Here, the high-confinement mode (H-mode) is very important for the fusion performance of ITER and future fusion reactors. The H-mode pedestal is a region featuring strongly reduced turbulence and transport just inside the limiting flux surface with strong plasma gradients, which drive edge localized modes (ELMs) in tokamaks. The ELMs, however, are quasiperiodic relaxations of the pedestal, resulting in a series of hot plasma eruptions similar to magnetospheric substorms, solar and stellar flares, which could potentially damage the ITER divertor plates and first walls. To mitigate and suppress the ELMs in H-mode plasmas, the fusion community has tried to develop effective actuators over the past two decades. Two separate excellent inventions have been awarded to recognize their achievements. The 2018 APS excellence award in Plasma Physics Research was given to Drs. Todd E. Evans, Max E. Fenstermacher, and Richard Alan Moyer for their first experimental demonstration of the stabilization of ELMs in high-confinement diverted discharges by application of very small edge-resonant magnetic perturbations, which led to the adoption of suppression coils in the ITER design. The 2018 Fusion Technology Award from the IEEE’s Nuclear and Plasma Science Society (NPSS) was given to Dr. Larry Baylor for his work designing the fueling, pumping and disruption mitigation system for the US ITER Project and for the implementation of pellet ELM mitigation for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observations of airflow around a supertall curved building and its impact on temperature and humidity in Houston's urban center

Characterization of realistically shaped skyscrapers embedded in non-uniform neighbourhoods experiencing intricate weather patterns remains inadequately investigated. Aiming to close this gap, the Center for Multiscale Applied Sensing team deployed its mobile observatory in the street canyons around the curved Wells Fargo Plaza skyscraper in downtown Houston, TX. Three deployments allowed airflow observations under different inflow wind and thermodynamic stability conditions. Doppler lidar measurements reveal that when inflow hits the curved wall of the skyscraper, perpendicular canyons experience similar vortex configurations creating two windward and two leeward circulations. Windsond measurements support that buoyancy within the deep street canyons can generate thermal updrafts as strong as 2 m s -1 which is sufficient to overturn the mechanical downwash under gentle wind conditions. Canyons experiencing venting during the daytime were observed to be more thermodynamically stable at night while thermodynamically stable canyons during the day were observed to be more thermodynamically unstable at night owing to the accumulation of heat near street-level. Fourier decomposition of the vertical velocity measurements shows that in all cases flow exhibited high Reynolds numbers and was composed of turbulent eddies of predominantly 6 and 15 min periods. Here, this observational dataset provides insights to assess wind load, pedestrian comfort, urban air mobility, and natural ventilation and may be used as a benchmark for numerical model and wind tunnel studies attempting to represent realistically complex urban conditions.

54 ENVIRONMENTAL SCIENCES↗

A Statistical Evaluation of WRF-LES Trace Gas Dispersion Using Project Prairie Grass Measurements

In recent years, new measurement systems have been deployed to monitor and quantify methane emissions from the natural gas sector. Large-eddy simulation (LES) has complemented measurement campaigns by serving as a controlled environment in which to study plume dynamics and sampling strategies. However, with few comparisons with controlled-release experiments, the accuracy of LES for modeling natural gas emissions is poorly characterized. In this paper, we evaluate LES from the Weather Research and Forecasting (WRF) Model against Project Prairie Grass campaign measurements and surface layer similarity theory. Using WRF-LES, we simulate continuous emissions from 30 near-surface trace gas sources in two stability regimes: strong convection and weak convection. We examine the impact of grid resolutions ranging from 6.25 to 52 m in the horizontal dimension on model results. We evaluate performance in a statistical framework, calculating fractional bias and conducting Welch’s t tests. WRF-LES accurately simulates observed surface concentrations at 100 m and beyond under strong convection; simulated concentrations pass t tests in this region irrespective of grid resolution. However, in weakly convective conditions with strong winds, WRF-LES substantially overpredicts concentrations—the magnitude of fractional bias often exceeds 30%, and all but one t test fails. The good performance of WRF-LES under strong convection correlates with agreement with local free convection theory and a minimal amount of parameterized turbulent kinetic energy. The poor performance under weak convection corresponds to misalignment with Monin–Obukhov similarity theory and a significant amount of parameterized turbulent kinetic energy.

17 WIND ENERGY↗

A Nonlocal Magneto-curvature Instability in a Differentially Rotating Disk

A global mode is shown to be unstable to nonaxisymmetric perturbations in a differentially rotating Keplerian disk containing either vertical or azimuthal magnetic fields. In an unstratified cylindrical disk model, using both global eigenvalue stability analysis and linear global initial-value simulations, it is demonstrated that this instability dominates at strong magnetic fields where local standard magnetorotational instability (MRI) becomes stable. Unlike the standard MRI mode, which is concentrated in the high flow shear region, these distinct global modes (with low azimuthal mode numbers) are extended in the global domain and are Alfvén-continuum-driven unstable modes. As its mode structure and relative dominance over MRI are inherently determined by the global spatial curvature as well as the flow shear in the presence of a magnetic field, we call it the magneto-curvature (magneto-spatial-curvature) instability. Consistent with the linear analysis, as the field strength is increased in the nonlinear simulations, a transition from MRI-driven turbulence to a state dominated by global nonaxisymmetric modes is obtained. This global instability could therefore be a source of nonlinear transport in accretion disks at a higher magnetic field than predicted by local models.

79 ASTRONOMY AND ASTROPHYSICS↗

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

WEST L-mode record long pulses guided by predictions using Integrated Modeling

A new record was set on the WEST Tokamak, designed to operate long duration plasmas in a tungsten (W) environment, with an injected energy of 1.15 GJ and a plasma duration 364s. Scenario development was supported by integrated modeling using the High Fidelity Plasma Simulator (HFPS), the European IMAS-coupled version of JETTO/JINTRAC, which integrates physics-driven modules into a unified framework. In particular, a reduced model for Lower-Hybrid heating and Current-Drive (LHCD) and the quasi-linear turbulent transport model TGLF are crucial for long pulses predictions up to the Last Closed Flux Surface (LCFS). Using this workflow, a 100 s reference discharge was modeled and plasma kinetic profiles and loop voltage were quantitatively well matched. In preparation for the recent long duration experiments, non-inductive current-drive actuators (I P , n e , P LHCD ) were varied to determine the operational domain going towards fully non-inductive discharges. In particular, decreasing the plasma current is shown to ease the access to such conditions, with a careful monitoring of (n e , P LHCD ) to avoid machine limitations. In addition, post-prediction experiments conducted within the investigated parameter range validated the predicted dependencies and were shown to be in quantitative agreement. Exploratory work on the use of ECCD for MHD stability purpose is also introduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗