SEARCH · Engineering Papers
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.
The data in this repository correspond to the figures of the paper "A study of resistive peeling-ballooning modes across spherical tokamaks" published in the journal Plasma Physics and Controlled Fusion. The figures are created based on simulations with the M3D-C1 and NIMROD codes using Python post-processing tools.
Atmospheric particles play critical roles in climate. However, significant knowledge gaps remain regarding the vertically resolved organic molecular-level composition of atmospheric particles due to aloft sampling challenges. To address this, we use a tethered balloon system at the Southern Great Plains Observatory and high-resolution mass spectrometry to, respectively, collect and characterize organic molecular formulas (MF) in the ground level and aloft (up to 750 m) samples. We show that organic MF uniquely detected aloft were dominated by organonitrates (139 MF; 54% of all uniquely detected aloft MF). Organonitrates that were uniquely detected aloft featured elevated O/C ratios (0.73 ± 0.23) compared to aloft organonitrates that were commonly observed at the ground level (0.63 ± 0.22). Unique aloft organic molecular composition was positively associated with increased cloud coverage, increased aloft relative humidity (~40% increase compared to ground level), and decreased vertical wind variance. Furthermore, 29% of extremely low volatility organic compounds in the aloft sample were truly unique to the aloft sample compared to the ground level, emphasizing potential oligomer formation at higher altitudes. Overall, this study highlights the importance of considering vertically resolved organic molecular composition (particularly for organonitrates) and hypothesizes that aqueous phase transformations and vertical wind variance may be key variables affecting the molecular composition of aloft organic aerosol.
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.
We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain with a scale-separated four-field reduced MHD (RMHD) model using the BOUT++ framework. In the early nonlinear phase, the small first crash is triggered by linearly unstable RDBMs, and m/n=2/1 magnetic islands are nonlinearly excited by nonlinear coupling of RDBMs as well as their higher harmonics. Here, m is the poloidal mode number, n is the toroidal mode number, the q = 2 rational surface exists near the pressure gradient peak, and q is the safety factor. Simultaneously, middle-n RDBM turbulence develops but is poloidally localized around X-points of the magnetic islands, leading to the small energy loss. The second large crash occurs in the late nonlinear phase. Higher harmonics of m/n=2/1 magnetic islands well develop around the q = 2 surface via nonlinear coupling and make the magnetic field stochastic by magnetic island overlapping. Turbulence heat transport develops at X-points of higher harmonics of m/n=2/1 magnetic islands, resulting in the turbulence spreading in the poloidal direction. The large second crash is triggered when the turbulence covers the whole poloidal region so that the magnetic island generation and magnetic field stochastization before the large crash can be interpreted as ELM precursors. It is concluded that the ELM trigger is attributed to the turbulent spreading in the poloidal direction in synchronization with the magnetic field stochastization and the crash is driven by E × B convection rather than the conventional Rechester–Rosenbluth anomalous electron heat transport.
We use a new gyrokinetic threshold model to predict a bifurcation in tokamak pedestal width-height scalings that depends strongly on plasma shaping and aspect-ratio. The bifurcation arises from the first and second stability properties of kinetic-ballooning-modes that yields wide and narrow pedestal branches, expanding the space of accessible pedestal widths and heights. The wide branch offers potential for edge-localized-mode-free pedestals with high core pressure. For negative triangularity, low-aspect-ratio configurations are predicted to give steeper pedestals than conventional-aspect-ratio. Both wide and narrow branches have been attained in tokamak experiments.
We apply the global gyrokinetic code GENE-3D to perform kinetic-ballooning-mode (KBM) turbulence simulations in the Wendelstein 7-X (W7-X) stellarator using experimentally relevant parameters. In contrast to the flux-tube approach, we observe saturation of the turbulent fluctuations at finite amplitude. Analyzing the turbulent characteristics, we identify the presence of weakly driven KBM turbulence hybridizing with ion-temperature-gradient and trapped-electron-mode turbulence in the scenario under consideration, showing that it is possible that KBM activity can be driven in W7-X, opening up the path for studying this type of turbulence, which is highly relevant for future fusion power plants.
The electron parallel nonlinearity (EPN) is implemented in the gyrokinetic particle-in-cell turbulence code GEM [Y. Chen and S. E. Parker, J. Comp. Phys. 220, 839 (2007)]. Application to the Cyclone Base Case reveals a strong effect of EPN on the saturated heat transport above the kinetic ballooning mode (KBM) threshold. Evidence is provided to show that the strong effect is associated with the electron radial motion due to magnetic fluttering, which turns fine structures of the KBM eigenmode in radius into fine structures in velocity and increases the magnitude of the EPN term in the kinetic equation.
A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.
We present nonlinear magnetohydrodynamic (MHD) simulations of high- β Wendelstein 7-X plasmas using the stellarator extension of the M3D- C 1 code, building on the recent work that shows benign saturation of ideal ballooning modes above the designed β limit in the standard configuration [Zhou et al., Phys. Rev. Lett. 133, 135102 (2024)]. First, we examine the results' sensitivity to the parallel thermal conductivity. It is found that while an increased parallel conductivity reduces the linear growth rate, the saturated pressure profile is barely affected. Second, we consider the dependence on the profile shape. It is shown that an equilibrium with a peaked pressure profile and lower β is subject to more significant change than a broad profile with higher β and a larger growth rate, suggesting that benign saturation, or nonlinear stability, is not guaranteed and not dictated by linear growth. Third, we study the influence of the magnetic configuration, with the equilibrium rotational transform varied by adjusting the planar coil current. With similar growth rates, similar magnitudes of profile change are found regardless of the presence of a low-order resonance, which implies that the saturation mechanism is not specific to a resonant or non-resonant mode. These results indicate that MHD stability should still be treated seriously in stellarator operation and design, for which nonlinear modeling using tools like M3D- C 1 can play an instrumental role.
This study investigates the impact of various impurity species on peeling–ballooning (PB) modes and microturbulence in tokamak plasmas through the extension of traditional two-fluid and gyro-landau-fluid (GLF) models. By incorporating finite Larmor radius (FLR) effects, the analysis provides a comprehensive understanding of impurity-driven impact and its interaction with plasma turbulence. Depending on charge state and local plasma conditions, heavy impurities may exhibit gyro-radii larger than those of main ions, which are captured in the extended GLF model presented. Following the presentation of modified two-fluid equations incorporating impurity effects, we systematically analyze the distinctions between impurity and main ion dynamics and their resultant feedback mechanisms on plasma behavior. Derivation of the linear dispersion relation enables quantification of impurity-mediated modifications to: plasma vorticity, diamagnetic drift and gyroviscous effects, electron Hall physics, and FLR effects. BOUT++ – based linear simulations corroborate this formalism, demonstrating systematic stabilization of PB modes upon impurity seeding. And then operational implications for practical impurity control strategies in tokamak devices are proposed. The results underscore the necessity of impurity management to maintain stability and optimize plasma confinement, with specific focus on how FLR effects contribute to transport dynamics. This work paves the way for enhanced modeling and simulation efforts, supporting the development of strategies to control impurity-induced turbulence and improve overall reactor performance.
Abstract Plasmoid-mediated reconnection is investigated for nonlinear peeling–ballooning edge-localized modes (P-B ELMs). The formation of current sheets and the transition to 3D current sheet instability is demonstrated through fully nonlinear resistive MHD simulations of P-B ELMs in DIII-D discharges. Large-scale axisymmetric current sheets, as well as small-scale poloidally extended current sheets, are formed as the coherent P-B ELM filaments nonlinearly evolve. It is observed that, at high Lundquist numbers, these current sheets break during a reconnection burst, i.e. a secondary exponential growth of intermediate modes followed by relaxation due to the suppression of P-B drive.
Abstract The edge localized mode (ELM) frequency ( f ELM ) decreased by 63% when electron cyclotron heating (ECH) deposition location is shifted from ρ = 0.4 to ρ = 0.8 in DIII-D discharges where the power ratio between neutral beam injection (NBI) and ECH ( P NBI / P ECH ) is kept at ∼1. The performance of the pedestal in the ECH heated case is compared with a pure NBI reference discharge while keeping the total input power constant. All these discharges are performed at balanced input torque conditions. Furthermore, in the pure NBI discharge a strong decoupling of the peeling–ballooning (PB) thresholds is observed. The PB decoupling is preserved when the ECH is deposited at ρ = 0.8 and P NBI / P ECH ∼ 1, while the thresholds manifest a closed stability boundary when the ECH is deposited at ρ = 0.4. The inter-ELM pedestal recovery time is considerably larger for the ECH at ρ = 0.8 case. Increased pedestal turbulence is observed in beam emission spectroscopy (BES), Doppler backscattering and magnetic diagnostics for the ECH at the ρ = 0.8 case. Strong growth of a TEM-like mode is observed in BES and the mode growth is correlated with the decrease in f ELM . In view of these observations, the increased pedestal turbulence seems to be the plausible reason behind the delayed pedestal recovery following an ELM event in the ECH at ρ = 0.8 case, and the preservation of PB decoupling through temperature pedestal profile widening. TRANSP interpretative simulations show that the ECH at the ρ = 0.8 case is more susceptible to ITG/TEM turbulence.
We report the observation of a set of coherent high frequency electromagnetic fluctuations that leads to a turbulence induced self-regulating phenomenon in the DIII-D high bootstrap current fraction plasma. The fluctuations have frequency of 130~220kHz, the poloidal wave length and phase velocity are 16~30 m -1 and ~30 km/s, respectively in the outboard midplane with the estimated toroidal mode number n~5- 9. The fluctuations are located in the internal transport barrier (ITB) region at large radius and are experimentally validated to be kinetic ballooning modes (KBM). Furthermore, quasilinear estimation predicts the KBM to be able to drive experimental particle flux and non-negligible thermal flux, suggesting its significant role in regulating the ITB saturation.
Core collapse supernovae are thought to be one of the main sources in the galaxy of elements heavier than iron. Understanding the origin of the elements is thus tightly linked to our understanding of the explosion mechanism of supernovae and supernova nucleosynthesis. X-ray and gamma-ray observations of young supernova remnants, combined with improved theoretical modeling, have resulted in enormous improvements in our knowledge of these events. Here, the isotope Ti44 is one of the most sensitive probes of the innermost regions of the core collapse engine, and its spatial and velocity distribution are key observables. Hard x-ray imaging spectroscopy with the Nuclear Spectroscopic Telescope Array (NuSTAR) has provided new insights into the structure of the supernova remnant Cassiopeia A (Cas A), establishing the convective nature of the supernova engine. However, many questions about the details of this engine remain. We present here the concept for a balloon-borne follow-up mission called A SuperConducting ENergetic x-ray Telescope (ASCENT). ASCENT uses transition edge sensor gamma-ray microcalorimeter detectors with a demonstrated 55-eV full-width half maximum energy resolution at 97 keV. This 8- to 16-fold improvement in energy resolution over NuSTAR will allow for high-resolution imaging and spectroscopy of the Ti44 emission. This will allow for a detailed reconstruction of gamma-ray line redshifts, widths, and shapes, allowing us to address questions such as, What is the source of the neutron star kicks? What is the dominant production pathway for Ti44? Is the engine of Cas A unique?
During long-distance sound propagation in planetary atmospheres, acoustic waves may reflect off the air/surface interface one or more times. For low sound frequencies and flat interfaces, the incident and reflected wave tend to be nearly identical. However, this may not be the case when the downgoing acoustic wave encounters topography. Here, we describe a set of direct and reflected acoustic signals recorded on free-flying balloons during the hypersonic entry of the OSIRIS-REx sample return capsule (SRC). In two of the three cases presented here, an impulsive reflected arrival similar in form to the direct sonic boom of the SRC is observed, followed by a diffuse coda. In contrast, one of the floating stations lacked an impulsive reflection entirely, with only the coda present. We use full waveform modeling to show how reflection in the presence of complex topography can explain coda seen in all three examples as well as the lack of impulsive arrival on the third. Our results indicate that the complex signals often observed in long range acoustic propagation could be due, in part, to interactions with topography during transmission.
Infrasound data and balloon trajectories. These data were collected during the OSIRIS-REx sample return capsule re-entry.
Natural events and human activity often generate acoustic waves capable of traveling tens to tens of thousands of kilometers across the globe. Ground-based acoustic sensors are limited to dry land and often suffer from wind noise. In contrast, balloon borne acoustic sensors can cross oceans, polar ice caps, and other inhospitable areas, greatly expanding sensor coverage. Since they move with the mean wind speed, their background noise levels are exceptionally low. In the last six years, such sensors have recorded sounds from colliding ocean waves, surface and buried chemical explosions, thunder, wind/mountain interactions, wind turbines, aircraft, and possibly meteors and the aurora. These results have led to new insights on acoustic heating of the upper atmosphere, the detectability of underground explosions, and directional sound fields generated by ocean waves.
The Mesa Photonics' cloud droplet measurement system (CDMS) performs in situ measurement of droplet size distribution and droplet number density in clouds. These characteristics are important cloud microphysical properties that are critical input parameters for atmospheric models and are also useful for proper calibration and validation of performance of other atmospheric measurement instrumentation. This small campaign was the first project (out of two) that involved integration of the CDMS into the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility’s tethered balloon system (TBS) and its initial testing at the ARM Southern Great Plains (SGP) atmospheric observatory. The follow-up campaign (AFC07016) involved testing of the CDMS under relevant conditions at the third ARM Mobile Facility (AMF3) at Oliktok Point, Alaska and making in situ measurements in clouds. The goal of this small field campaign was to deploy the CDMS on the ARM TBS, demonstrate its integrity and compatibility with the TBS, and test the wireless telecommunication system in preparation for subsequent in-cloud measurements at ARM’s North Slope of Alaska (NSA) observatory at Utqiagvik (formerly Barrow). The specific technical objectives included: (1) Integration of the Mesa Photonics' CDMS into the ARM TBS, (2) testing of the wireless telecommunication system of the CDMS, (3) evaluating the data acquisition software and image processing algorithms especially under high-background-illumination conditions, and (4) testing the ruggedness of the instrument's optical alignment and other performance characteristics. All objectives have been successfully met. The mounting hardware was designed and built mostly at Sandia National Laboratories (SNL), which allowed reliable mounting of the CDMS on the TBS and co-locating the CDMS with other instruments. The field tests were performed by the principal investigator (PI) in collaboration with the TBS operational crew led by D. Dexheimer. The field testing demonstrated good compatibility of the CDMS with the TBS. During the TBS flights, the crew did not experience any operational issues with the CDMS. The instrument demonstrated reliable operational characteristics, including sufficient battery life and good thermal management. The wireless telecommunication system has been successfully tested at TBS flight altitudes of up to 1000 m. During all flights, the CDMS raw data were wirelessly transmitted to the ground station and processed in real time. The ruggedness of the instrument's optical alignment was evaluated by performing calibration of the CDMS (using Mesa Photonics' fixed-size monodisperse droplet generator) before and after the campaign. The calibration was stable within the instrument's measurement precision. In summary, the campaign demonstrated good compatibility of the Mesa Photonics' CDMS with the ARM TBS. The flight tests confirmed the mechanical integrity of the CDMS and stability of its optical layout, as well as reliability of the wireless telecommunication system. Successful completion of this campaign provided a solid basis for the next campaign (AFC07016, November 2020), that involved testing of the CDMS during in-cloud TBS flights at AMF3 at Oliktok Point.