Search NASA⌕ Search

SEARCH · Search NASA

Results for “global instability”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A subgrid model for electron-scale turbulent transport in global ion-scale gyrokinetic simulations of tokamak plasmas

A subgrid electron–temperature–gradient (ETG) model is demonstrated here which averages local electron-scale turbulence over intermediate scales in space and time to include in global ion–temperature–gradient (ITG) simulations. This approach results in ion-scale equations, which incorporate the electron heat transport from ETG turbulence and the effects of electron-scale turbulence on the ion scale. Flux-tube ETG Cyclone Base Case simulations are performed at different radial locations and a kinetic form of the flux is added to global ion-scale (IS) simulations as a source term. Analytic radial profiles of ETG heat flux are constructed and compared to flux-tube results at multiple radial locations. Different ratios of ITG to ETG heat flux levels are considered and the results of capturing ETG heat transport in global ion-scale simulations are discussed. The potential coupling of the ETG streamer potential and intermediate-scale zonal flows to the ion scale is further addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Gyrokinetic particle simulations of interactions between energetic particles and magnetic islands induced by neoclassical tearing modes

Interactions between energetic particles (EPs) and neoclassical tearing mode (NTM) islands in the DIII-D tokamak are studied using the global gyrokinetic toroidal code (GTC). GTC simulations find that the EP radial profile is partially flattened within the magnetic island regions and that there are stochastic regions in the particle phase space. Radial particle flux is induced mainly around the magnetic island regions and decreases with time to almost zero when the initial EP distribution achieves a new steady-state in the absence of EP sources. Stochastic regions of magnetic field lines induced by the superposition of multiple islands have weak effects on the particle flux when the width of stochastic regions is smaller than the EP drift orbit width. In conclusion, the perturbed parallel EP current induced by the magnetic islands has weak stabilizing effects on the linear growth rate of the NTM instability in this DIII-D experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dynamics of multilayer Rayleigh–Taylor instability at moderately high Atwood numbers

This paper investigates the multilayer Rayleigh–Taylor instability (RTI) using statistically stationary experiments conducted in a gas tunnel. Employing diagnostics such as particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF), we make simultaneous velocity–density measurements to study how dynamics and mixing are linked in this variable density flow. Experiments are conducted in a newly built, blow-down three-layer gas tunnel facility. Mixing between three gas streams is studied, where the top and bottom streams are comprised of air, and the middle stream is an air–helium mixture. Shear is minimized between these streams by matching their inlet velocities. The four experimental conditions investigated here consist of two different density ratios (Atwood numbers 0.3 and 0.6), each investigated at two instability development times (or equivalently, two streamwise locations), and all experiments are with the same middle stream thickness of 3 cm. The growth of the middle layer is measured using laser-based planar Mie scattering visualization. The mixing width is found to grow linearly with time at late times. Various quantitative measures of molecular mixing indicate a very high degree of molecular mixing at late times in the multilayer RTI flow. The vertical turbulent mass flux $a_y$ is calculated. In addition to mostly negative values of $a_y$ , typical of buoyancy-dominated flows due to negative correlation between velocity and density fluctuations, positive regions are also observed in profiles of $a_y$ due to entrainment and erosion at the lower edge of the mixing region. Global energy budgets are calculated for the multilayer RTI flow at late times and it is found that the majority of potential energy released has been dissipated due to viscous effects, and a large value of mixing efficiency ( $\sim$ 60 %) is observed.

Mechanics↗
DOE OSTI2023

Nonlinear second order electromagnetic gyrokinetic theory for a tokamak plasma

The steep plasma pressure gradient that forms at the edge of the high confinement, H-mode regime of tokamak operation provides free energy to drive electromagnetic micro-instabilities that are widely believed to influence the transport processes in this so-called pedestal region. This high pressure gradient also provides a high current density (bootstrap current), known to influence ballooning mode stability and to be important for driving kink modes in the ideal magneto-hydrodynamic plasma model (so-called peeling-ballooning modes). Furthermore, efficient, steady state future tokamak power plants must operate with a large bootstrap current in the core and especially concerning spherical tokamaks, confinement will be influenced by electromagnetic turbulence. To accommodate these important situations, conventional electromagnetic gyrokinetic theory is extended to incorporate neoclassical effects in the equilibrium drives, allowing $B_{\vartheta} ~ B_0$ (B0 is the confining magnetic field, and $B_{\vartheta}$ is its poloidal component). This provides a global gyrokinetic model that self-consistently captures the consequences of large bootstrap current fractions on the equilibrium distribution functions.

kink modes↗
DOE OSTI2021

Thermal confinement and transport in spherical tokamaks: a review

Here, we review the thermal plasma confinement and transport properties observed and predicted in low aspect ratio tokamaks, or spherical tokamaks (STs), which can depart significantly from those observed at higher aspect ratio. In particular, thermal energy confinement scalings show a strong, near linear dependence of energy confinement time on toroidal magnetic field, while the dependence on plasma current is more modest, the opposite of what is seen at higher aspect ratio. STs have revealed a very strong improvement in normalized confinement with decreasing collisionality, much stronger than at higher aspect ratio, which bodes well for an ST-based fusion pilot plant should this trend continue at an even lower collisionality than has already been accessed. These differences arise because of fundamental differences in transport in STs due to the more extreme toroidicity (i.e. reduced region of bad curvature), and to the relatively larger $E_r \times B$ shearing rates, both of which can suppress electrostatic drift wave instabilities at both ion and electron gyroradius scales. In addition, electromagnetic effects are much stronger in STs because they operate at high β T . Gyrokinetic (GK) studies, coupled with low- and high-k turbulence measurements, have shed light on the underlying physics controlling transport. At lower β T , both ion- and electron-scale electrostatic drift turbulence may be responsible for transport. At higher β T , microtearing, kinetic ballooning, and hybrid trapped electron/kinetic ballooning modes increasingly play a role, and they have a much stronger impact in the core of ST plasmas than at higher aspect ratio. Flow shear affects the balance between ion- and electron-scale modes. Non-linear GK simulations find regimes where the electron heat flux decreases with decreasing collisionality, consistent with the experimental global normalized confinement scaling. The ST is unique in that the relatively low toroidal magnetic field allows for localized measurements of electron-scale turbulence, and this coupled with turbulence measurements at ion-scales has facilitated detailed comparisons with GK simulations. These data have provided compelling evidence for the presence of ion temperature gradient and electron temperature gradient turbulence in some plasmas, and direct experimental support for the impact of experimental actuators like rotation shear, density gradient and magnetic shear on turbulence and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Isotopic mass effects of tritium-fueled high-performance TFTR supershots

An increase in total stored energy correlated with the addition of tritium fuel was observed in supershots during the TFTR DT campaign. This supershot regime had strikingly high, centrally peaked ion and electron temperatures, and the largest neutron emission rates observed in TFTR. This paper presents a study of the causes of this increase in stored energy in supershots. Twenty-six supershots have been recently reanalyzed with the TRANSP plasma analysis code. Early TRANSP simulations did not accurately match the measured magnitude and time evolutions of the neutron emission rates. This mismatch is attributed to neglecting apparent increases of trace amounts of heavy impurities during neutral beam injection. The new TRANSP runs were tuned to accommodate this and match the measured global neutron emission rates. These new runs also had improved fidelity in predicting the time histories and radial dependencies of measured DT neutron emission rates. That in turn adds confidence in the simulated thermal deuterium and tritium density profiles that are needed for calculating the average hydrogenic atomic mass profiles. Six subsets of these supershots had well matched toroidal field B tor , plasma current I p , flux geometry, and total injected neutral beam power. The mix of D and T beam ions was varied for different discharges. The magnitude of the increase of the thermal ion energy W i with added tritium was relatively small, and the total W tot increase is dominated by the increase in fast beam ions with T. Analyses at times before the occurrences of deleterious MHD instabilities yielded scaling of W tot with the volume-average isotopic mass consistent with previous publications. The relative fraction of fast energy ions is expected to be small in practical tokamak reactors. Thus the increase in stored energy W tot observed in TFTR supershots does not appear likely to be significantly helpful for producing useful fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Close coupling of the peeling ballooning and external kink modes in high poloidal-beta discharges with strong shape and large pedestal width

In plasmas with strong shape and large pedestal width, like DIII-D high poloidal-beta discharges, the eigenvalue spectrum of the peeling ballooning mode (PBM) is shown to shift to the very low-n regime near one, where n is the toroidal mode number. This peeling-type eigenvalue spectrum has a further shift to n = 1, thus being smoothly connected to the n = 1 external kink mode (EKM), as the normalized beta (β N ) increases. Once this connection occurs, the mode takes a mixture form of the PBM and EKM with its mode structure varying from the PBM-like to the EKM-like one as β N increases. The mode stability also becomes sensitive to both the local pedestal gradient and global β N , thus allowing an anti-correlation between the two driving forces. Further, these results appear to provide a qualitative explanation of the two unexpected features observed in the DIII-D high poloidal-beta discharges, that is, the dominance of the n = 1 mode in the edge-localized-modes and the negative correlation between the pedestal height and the internal transport barrier strength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Exploration of Sn 70 Li 30 alloy as possible material for flowing liquid metal plasma facing components

As an advanced alternative to solid materials, Liquid Metals (LM) may offer more resilient and feasible Plasma Facing Components (PFCs). Particularly, regarding the unavoidable material erosion/degradation produced by particle/heat fluxes in future fusion devices where much longer duty cycles are expected. Furthermore, configurations that propose a flowing LM surface can add the advantage of a continuously fresh and clean layer facing the plasma. Although lithium is the most widely tested option, tin-lithium (SnLi) alloys have been proposed to attempt to combine the positive characteristics of both pure elements and ameliorate the specific issues of lithium. In this work, the potential use of Sn70Li30 alloy in such flowing concepts has been explored by addressing several preliminary and mandatory aspects for its utilization. Key issues such as wettability and compatibility of the alloy with relevant substrates have been studied in a multidisciplinary approach. The data obtained from deposited liquid tin-lithium droplets indicates approximate wetting temperatures of 360 °C, 390 °C and 405 °C for the fresh alloy on 316 stainless steel, molybdenum, and tungsten, respectively. However, the alloy contamination appeared to strongly affect the wetting characteristics of materials, increasing their wetting temperature by ~130 °C in the worst observed cases. Interestingly, in some instances, the instability of the liquid alloy surface was observed in the form of sudden gaseous ejection. The deposited droplets were posteriorly characterized in terms of absolute composition and depth profile by Inductively Coupled PlasmaOptical Emission Spectroscopy (ICP-OES) and Secondary Ion Mass Spectrometry (SIMS-ToF). Additionally, the nature and composition of the boundaries between the substrates and alloy microparticles was investigated by Scanning Electron Microscopy (SEM), Energy Dispersive X-Ray Spectroscopy (EDS), and 3D Laser microscopy. The overall results of this post-mortem characterization revealed that first signs of corrosion induced by both alloy elements (lithium-chromium association and iron-tin intermetallic mixing) were present on 316 stainless steel after short exposures (≤3 h) at temperatures lower than 550 °C. Conversely, molybdenum and tungsten showed good compatibility with the alloy in equivalent conditions. The global implications of these results are finally addressed, focusing on the future perspectives and the more viable scenarios for the eventual utilization of these alloys in flowing liquid metal configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A Chip-Scale Oscillation-Mode Optomechanical Inertial Sensor Near the Thermodynamical Limits

Modern navigation systems integrate the global positioning system (GPS) with an inertial navigation system (INS), which complement each other for correct attitude and velocity determination. Here, the core of the INS integrates accelerometers and gyroscopes used to measure forces and angular rate in the vehicular inertial reference frame. With the help of gyroscopes and by integrating the acceleration to compute velocity and distance, precision and compact accelerometers with sufficient accuracy can provide small–error location determination. Solid–state implementations, through coherent readout, can provide a platform for high performance acceleration detection. In contrast to prior accelerometers using piezoelectric or capacitive readout techniques, optical readout provides narrow–linewidth high–sensitivity laser detection along with low–noise resonant optomechanical transduction near the thermodynamical limits. Here an optomechanical inertial sensor with an 8.2 µg Hz –1/2 velocity random walk (VRW) at an acquisition rate of 100 Hz and 50.9 µg bias instability is demonstrated, suitable for applications, such as, inertial navigation, inclination sensing, platform stabilization, and/or wearable device motion detection. Driven into optomechanical sustained–oscillation, the slot photonic crystal cavity provides radio–frequency readout of the optically–driven transduction with an enhanced 625 µg Hz –1 sensitivity. Measuring the optomechanically–stiffened oscillation shift, instead of the optical transmission shift, provides a 220× VRW enhancement over pre–oscillation mode detection.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

MHD, disruptions and control physics: Chapter 4 of the special issue: on the path to tokamak burning plasma operation

In this chapter, we review the progress in MHD stability, disruptions and control in magnetic fusion research that has occurred over the past (more than) one and a half decades since the publication by Hender et al in 2007 on the same topic as part of the update of ITER Physics Basis. During this period, remarkable progress has been achieved in the understanding of the basic physics and overall control of MHD instabilities through a wide spectrum of dedicated experiments, theory and modeling. The sawtooth activities are probably today one of the best understood of MHD events and very robust control schemes have been developed for reliable operation of tokamaks through core heating. Similarly, significant improvements have been achieved in understanding and control of neoclassical tearing modes, resistive wall modes or locked modes and their control through ECCD or error field control. The field of disruption prediction through application of artificial intelligence, machine learning or deep learning methods, which had already started at the time of the 2007 review, has progressed significantly due to general progress in these fields and application of newer, more sophisticated algorithms. However, although remarkable progress has been achieved in the field of Disruptions, their understanding, prediction, possible avoidance and mitigation still remain probably the most active fields of R&D globally in this field. This is especially because reactor grade machines like ITER and DEMO will be much less tolerant in respect of disruptions and runaway currents, and their occurrences must be either avoided altogether or minimized to an acceptable value without causing any significant hindrance to robust machine operations. This review is intended to present a broad spectrum of the R&D that has occurred in this field in support of ITER, which will also be of immense significance for all future machines, especially reactors like DEMO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Machine Learning for Mapping Multipactor Susceptibility in RF Systems: Capabilities and Generalization Constraints

Multipactor is a surface-driven electron avalanche phenomenon that degrades the performance and reliability of radio-frequency (RF) systems in particle accelerator and vacuum electronics applications. Multipactor behavior in a given device structure is conventionally assessed through susceptibility charts, which provide a parameter-space characterization of the instability. In this work, we assess the capabilities of machine-learning (ML) models to learn and predict such susceptibility charts and analyze the constraints governing their generalization across materials. Using a simulation-derived dataset spanning six distinct secondary-electron-yield material profiles in a canonical two-surface planar geometry, we train supervised regression models and artificial neural networks to predict the time-averaged electron growth rate, δavg, across the relevant parameter space. Model performance is evaluated using metrics that explicitly probe the structure of susceptibility charts, including Intersection over Union, Structural Similarity Index, and correlation analysis. Tree-based ensemble models outperform neural-network models in reconstructing susceptibility regions and in generalizing across material domains. Principal-component analysis reveals disjoint material feature distributions, indicating that the piecewise mode structure of multipactor susceptibility is difficult to represent with a single global model and that generalization is constrained by data coverage rather than by model complexity. An exhaustive reduced-coverage study further shows that sparse material-space coverage can yield mean performance in the same general range but producing large variability in the susceptibility-region overlap. These results clarify the capabilities of ML-based surrogate models for parameter-space characterization of multipactor discharge. They also provide guidance for their appropriate use in RF system design.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Orthoimagery and Shapefiles Documenting Pre- and Post-August 2019 Slope Disturbances, Teller Road Site, Seward Peninsula, Alaska, 2018-2019

This dataset was derived from UAS aerial photos and dGPS data collected at the Teller mile marker 47 site in 2018 and 2019 and used to support research quantifying the timing and rate of surface movements and analyze soil transport process in Arctic landscapes. In July 2018, a Phantom uncrewed aerial system (UAS) collected >6800 high resolution aerial photos of the lower portion of the Teller 47 watershed (see for raw photos; pending archive NGA281). During the UAS survey, 25 x 25 cm tile ground control points (GCPs) were laid out and secured with one rebar rod in the center. The four corners and rebar tops were surveyed with differential GPS, and these coordinates were used to construct and validate an orthomosaic of the site. These points were resurveyed in August 2019 and used to track annual movement. Changes in position between the two surveys are reported in (Lathrop et al. 2022; NGA254). Agisoft Metashape photogrammetry software was used to construct a georeferenced orthomosaic image (*.tif file) of a portion (0.68 km2) of the watershed with a final resolution of ~1 cm. Manual delineation of the perimeters of failures visible in the UAS imagery (*.tif files) was conducted to create shapefile polygons (two *.zip files). The failures were identified by the exposure of bare mineral soils, which were made visible by disruption of the overlying tundra vegetation. The shapefiles were then used to analyse the topographic distribution and sizes of the failures. In August 2019 an additional ~300 georeferenced UAS images of slope instability features were collected (data pending submission) and compared with the 2018 orthomosaic. Overview maps of the failure locations included as a *.pdf.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Detecting shearless phase-space transport barriers in global gyrokinetic turbulence simulations with test particle map models

In magnetically confined fusion plasmas, the role played by zonal 𝐸 x 𝐵 flow shear layers in the suppression of turbulent transport is relatively well understood. However, less is understood about the role played by the weak shear regions that arise in the non-monotonic radial electric field profiles often associated with these shear layers. In electrostatic simulations from the global total-𝑓 gyrokinetic particle-in-cell code XGC, we demonstrate how shearless regions with non-zero flow curvature form zonal ‘jets’ that, in conjunction with neighbouring regions of shear, can act as robust barriers to particle transport and turbulence spreading. By isolating quasi-coherent fluctuations radially localised to the zonal jets, we construct a map model for the Lagrangian dynamics of gyrokinetic test particles in the presence of drift waves. We identify the presence of shearless invariant tori in this model and verify that these tori act as partial phase-space transport barriers in the simulations. We also demonstrate how avalanches impinging on these shearless tori cause eddy detachment events that form ‘cold/warm core ring’ structures analogous to those found in oceanic jets, facilitating transport across the barriers without destroying them completely. We discuss how shearless tori may generically arise from tertiary instabilities or other types of discrete eigenmodes, suggesting their potential relevance to broader classes of turbulent fluctuations.

fusion plasma↗
DOE OSTI2022

Coupling a Lagrangian–Eulerian Spark-Ignition (LESI) model with LES combustion models for engine simulations

In the United States transportation sector, Light-Duty Vehicles (LDVs) are the largest energy consumers and CO 2 emitters. Electrification of LDVs is posed as a potential solution, but SI engines can still contribute to decarbonization. Car manufacturers have turned to unconventional engine operation to increase the efficiency of Spark-Ignition (SI) engines and reduce the carbon emissions of their fleets. Dilute, lean, and stratified-charge engine operation has the potential for engine efficiency improvements at the expense of increased cyclic variability and combustion instability. At such demanding engine conditions, the spark ignition event is key for flame initiation and propagation and for enhanced combustion stability. Reliable and accurate spark ignition models can help design ignition systems that reduce cyclic variability. Multiple computational spark-ignition models exist that perform well under conventional conditions, but the underlying physics needs to be expanded, for unconventional engine operation. In this paper, a hybrid Lagrangian–Eulerian Spark-Ignition (LESI) model is coupled with different turbulent flame propagation models for engine simulations. LESI relies on Lagrangian arc tracking and Eulerian energy deposition. The LESI model is coupled with the Well-Stirred Reactor (WSR), Thickened Flame Model (TFM), and g-equation model and used to simulate several cycles of a Direct-Injection Spark-Ignition (DISI) engine using a commercial Computational Fluid Dynamics (CFD) engine solver. The results showcase the successful coupling of LESI with the combustion models. Global engine metrics, such as pressure and Apparent Heat Release Rate (AHRR), for each simulation setup are compared to experimental engine results, for validation. In addition, results highlight the successful prediction of spark channel movement by comparing simulation images to experimental optical engine images. Finally, the successful coupling of LESI to combustion models, making it a usable model in the engine modeling community, is emphasized and future development details are discussed.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Subseasonal Tropical Convection Characteristics in the Energy Exascale Earth System Model Version 2

Accurate simulation of subseasonal tropical moist convection remains a key challenge for Earth system models. The difficulties stem from the reliance of cumulus cloud processes on model parameterizations and the need to represent the multiscale nature of interactions among clouds, radiation, moisture, circulation, and surface energy fluxes. Equatorial convection drives circulation anomalies that can affect weather patterns and extremes globally, motivating efforts to better understand and simulate these tropical disturbances. Here, a detailed review of subseasonal tropical convective behavior as simulated in the Energy Exascale Earth System Model version 2 (E3SMv2) is presented, with comparison to its predecessor version 1 (E3SMv1) and reference data sets. Model structural changes to the deep convective trigger function and surface fluxes, along with parametric tuning of the cloud and microphysics schemes, together result in an improved depiction of organized tropical convection across scales. In particular, E3SMv2 exhibits a more realistic Madden‐Julian oscillation (MJO) and low‐frequency Kelvin waves—owing to a sharper time mean equatorial meridional moisture gradient and improved convection‐circulation coupling —as well as a better depiction of MJO Northern Hemisphere teleconnections. Despite these improvements, subseasonal precipitation variance continues to be strongly underestimated in E3SMv2. Use of a cloud plume model also reveals that the coupling between daily averaged tropical precipitation and lower tropospheric instability in E3SM is inconsistent with observations, a bias that could potentially impact the simulation of intraseasonal disturbances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modulation of the Eastern Equatorial Pacific Seasonal Cycle by Tropical Instability Waves

Feedbacks from tropical instability waves (TIWs) on the seasonal cycle of the eastern Pacific Ocean are studied using two eddy-rich ocean simulations, with and without TIWs. By warming the equatorial waters by up to 0.4°C through nonlinear advection in boreal summer and fall, TIWs reduce the amplitude of the seasonal cycle in upper ocean temperatures. In addition, TIWs stabilize the upper part of the Equatorial Undercurrent (EUC) through enhanced barotropic energy conversion, leading to a year-round weakening by -0.15 m s -1 and preventing an unrealistic re-intensification in boreal fall usually found in non-eddy resolving models. A coarser simulation at 1-degree horizontal resolution fails to reproduce the TIW-induced nonlinear warming of equatorial waters, but succeeds in inhibiting the EUC re-intensification. This suggests a threshold effect in TIW strength, associated with the model's ability to simulate eddies, which may be responsible for long-standing biases displayed by global climate models in this region.

58 GEOSCIENCES↗
DOE OSTI2026

Density limit in peeling-limited pedestals at and above the Greenwald value in DIII-D high poloidal beta plasmas

High pressure, peeling limited pedestals with pedestal normalized beta β N,ped >2 and pedestal top density n ped at or above the Greenwald density nG have been achieved in DIII-D high poloidal beta plasmas, with high global normalized beta β N >3 and energy confinement H 98 ~1.2-1.7. Higher β N allows higher pedestal density above the Greenwald value and higher pedestal pressure, even with a low injected torque. MHD modeling confirms that the experimental profiles lie near the peeling-mode unstable boundary with high normalized pressure gradient βMHD and high edge current density. Experimental analysis and stability calculations indicate that the high poloidal beta with strong Shafranov shift, high βMHD and weak/negative magnetic shear improves the pedestal stability by decoupling the peeling and ballooning modes and stabilizing the ballooning modes, thus facilitating access to the second stable region of peeling-ballooning mode. The access to the second stable peeling-ballooning stability region opens Super-H-like channels without extremely strong shaping or strong torque injection. The high-pressure peeling pedestal allows the pedestal density to go beyond the Greenwald limit with strong ExB shear maintained: pedestal pressure increases with pedestal density even when n ped >n G , until reaching the ideal MHD instability boundary, where giant ELMs occur. The giant ELMs are dominated by a strong n=1 component and cause a large reduction of the edge pressure, but a negligible change of the core pressure, consistent with kink/peeling-mode induced instability. The pedestal recovers from the collapse and typically sustains a high baseline density, around the Greenwald limit, during the whole discharge duration. Experiments also found that internal transport barriers and n ped ~0.9nG, peeling limited pedestals could be simultaneously achieved in high β N plasmas, while an internal feedback between ITB strength and pedestal performance is found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗