We report results on rotating stratified turbulence in the absence of forcing, with large-scale isotropic initial conditions, using direct numerical simulations computed on grids of up to 4096^3 points. The Reynolds and Froude numbers are respectively equal to Re=5.4 x 10^4 and Fr=0.0242. The ratio of the Brunt-V\"ais\"al\"a to the inertial wave frequency, N/f, is taken to be equal to 4.95, a choice appropriate to model the dynamics of the southern abyssal ocean at mid latitudes. This gives a global buoyancy Reynolds number R_B=ReFr^2=32, a value sufficient for some isotropy to be recovered in the small scales beyond the Ozmidov scale, but still moderate enough that the intermediate scales where waves are prevalent are well resolved. We concentrate on the large-scale dynamics, for which we find a spectrum compatible with the Bolgiano-Obukhov scaling, and confirm that the Froude number based on a typical vertical length scale is of order unity, with strong gradients in the vertical. Two characteristic scales emerge from this computation, and are identified from sharp variations in the spectral distribution of either total energy or helicity. A spectral break is also observed at a scale at which the partition of energy between the kinetic and potential modes changes abruptly, and beyond which a Kolmogorov-like spectrum recovers. Large slanted layers are ubiquitous in the flow in the velocity and temperature fields, with local overturning events indicated by small Richardson numbers, and a small large-scale enhancement of energy directly attributable to the effect of rotation is also observed.

D. Rosenberg A. Pouquet R. Marino P. D. Mininni

We are investigating the SN 1987A remnant by modeling the circumstellar environment of the progenitor star. Interacting stellar winds models have been reasonably successful at reproducing the gross features of the nebula, but some details, such as the early rise of the radio and x-ray emission from the supernova blast and the presence of the outer rings, are not explained in a pure wind model. In this paper we describe preliminary results from 2-D models that incorporate the effects of photoionization due to radiation from the central star. These models have successfully produced a thick HII region, as needed to produce the early radio and X-ray emission. The HII region is thickest away from the equatorial plane. The models have also produced a feature resembling the outer rings, but we suspect this to be an artifact of the 2-D calculations that would not persist in 3-D.

Robert Link Duane L. Rosenberg Roger A. Chevalier

We present an object-oriented geophysical and astrophysical spectral-element adaptive refinement (GASpAR) code for application to turbulent flows. Like most spectral-element codes, GASpAR combines finite-element efficiency with spectral-method accuracy. It is also designed to be flexible enough for a range of geophysics and astrophysics applications where turbulence or other complex multiscale problems arise. For extensibility and flexibilty the code is designed in an object-oriented manner. The computational core is based on spectral-element operators, which are represented as objects. The formalism accommodates both conforming and nonconforming elements and their associated data structures for handling interelement communications in a parallel environment. Many aspects of this code are a synthesis of existing methods; however, we focus on a new formulation of dynamic adaptive refinement (DARe) of nonconforming h-type. This paper presents the code and its algorithms; we do not consider parallel efficiency metrics or performance. As a demonstration of the code we offer several two-dimensional test cases that we propose as standard test problems for comparable DARe codes. The suitability of these test problems for turbulent flow simulation is considered.

Duane Rosenberg Aime' Fournier Paul Fischer Annick Pouquet

The interplay between waves and eddies in stably stratified rotating flows is investigated by means of world-class direct numerical simulations using up to $3072^3$ grid points. Strikingly, we find that the shift from vortex to wave dominated dynamics occurs at a wavenumber $k_R$ which does not depend on Reynolds number, suggesting that partition of energy between wave and vortical modes is not sensitive to the development of turbulence at the smaller scales. We also show that $k_R$ is comparable to the wavenumber at which exchanges between kinetic and potential modes stabilize at close to equipartition, emphasizing the role of potential energy, as conjectured in the atmosphere and the oceans. Moreover, $k_R$ varies as the inverse of the Froude number as explained by the scaling prediction proposed, consistent with recent observations and modeling of the Mesosphere-Lower Thermosphere and of the ocean.

Raffaele Marino Duane Rosenberg Corentin Herbert Annick Pouquet

We observe the emergence of strong vertical drafts in direct numerical simulations of the Boussinesq equations in a range of parameters of geophysical interest. These structures, which appear intermittently in space and time, generate turbulence and enhance kinetic and potential energy dissipation, providing a possible explanation for the observed variability of the local energy dissipation in the bulk of oceanic flows, and the modulation of its probability distribution function. We show how, due to the extreme drafts, in runs with Froude numbers observable in geophysical scenarios, roughly 10% of the domain flow can account for up to 50% of the global volume dissipation, reminiscent of estimates based on oceanic models.

Raffaele Marino Fabio Feraco Leonardo Primavera Alain Pumir Annick Pouquet Duane Rosenberg Pablo D. Mininni

We perform numerical simulations of decaying rotating stratified turbulence and show, in the Boussinesq framework, that helicity (velocity-vorticity correlation), as observed in super-cell storms and hurricanes, is spontaneously created due to an interplay between buoyancy and rotation common to large-scale atmospheric and oceanic flows. Helicity emerges from the joint action of eddies and of inertia-gravity waves (with inertia and gravity with respective associated frequencies $f$ and $N$), and it occurs when the waves are sufficiently strong. For $N/f < 3$ the amount of helicity produced is correctly predicted by a quasi-linear balance equation. Outside this regime, and up to the highest Reynolds number obtained in this study, namely $Re\approx 10000$, helicity production is found to be persistent for $N/f$ as large as $\approx 17$, and for $ReFr^2$ and $ReRo^2 $ respectively as large as $\approx 100$ and $\approx 24000$.

Raffaele Marino Pablo D. Mininni Duane Rosenberg Annick Pouquet

As part of the development of a Poisson solver for the spectral element discretization used in the GeoFluid Object Workbench (GeoFLOW) code, we propose a solver for the linear system arising from a Gauss-Legendre-Lobatto global spectral method. We precondition using a p-multigrid {\gamma}-cycle with highly-vectorizable smoothers, that we refer to as line smoothers. Our smoothers are restrictions of spectral and finite element discretizations to low-order one-dimensional problems along lines, that are solved by a reformulation of cyclic reduction as a direct multigrid method. We illustrate our method with numerical experiments showing the apparent boundedness of the iteration count for a fixed residual reduction over a range of moderately deformed domains, right hand sides and Dirichlet boundary conditions.

José Pablo Lucero Lorca Duane Rosenberg Isidora Jankov Conor McCoid Martin Jakob Gander

Interest in generative Electrocardiogram-Language Models (ELMs) is growing, as they can produce textual responses conditioned on ECG signals and textual queries. Unlike traditional classifiers that output label probabilities, ELMs are more versatile, supporting domain-specific tasks (e.g., waveform analysis, diagnosis, prognosis) as well as general tasks (e.g., open-ended questions, dialogue). Retrieval-Augmented Generation (RAG), widely used in Large Language Models (LLMs) to ground LLM outputs in retrieved knowledge, helps reduce hallucinations and improve natural language generation (NLG). However, despite its promise, no open-source implementation or systematic study of RAG pipeline design for ELMs currently exists. To address this gap, we present the first open-source RAG pipeline for ELMs, along with baselines and ablation studies for NLG. Experiments on three public datasets show that ELMs with RAG consistently improves performance over non-RAG baselines and highlights key ELM design considerations. Our code is available at: https://github.com/willxxy/ECG-Bench.

Xiaoyu Song William Han Tony Chen Chaojing Duan Michael A. Rosenberg Emerson Liu Ding Zhao

Using a code that employs a self-consistent method for computing the effects of photoionization on circumstellar gas dynamics, we model the formation of wind-driven nebulae around massive Wolf-Rayet (W-R) stars. Our algorithm incorporates a simplified model of the photo-ionization source, computes the fractional ionization of hydrogen due to the photoionizing flux and recombination, and determines self-consistently the energy balance due to ionization, photo-heating and radiative cooling. We take into account changes in stellar properties and mass-loss over the star's evolution. Our multi-dimensional simulations clearly reveal the presence of strong ionization front instabilities. Using various X-ray emission models, and abundances consistent with those derived for W-R nebulae, we compute the X-ray flux and spectra from our wind bubble models. We show the evolution of the X-ray spectral features with time over the evolution of the star, taking the absorption of the X-rays by the ionized bubble into account. Our simulated X-ray spectra compare reasonably well with observed spectra of Wolf-Rayet bubbles. They suggest that X-ray nebulae around massive stars may not be easily detectable, consistent with observations.

Vikram V. Dwarkadas Duane Rosenberg

We study the partition of energy between waves and vortices in stratified turbulence, with or without rotation, for a variety of parameters, focusing on the behavior of the waves and vortices in the inverse cascade of energy towards the large scales. To this end, we use direct numerical simulations in a cubic box at a Reynolds number Re=1000, with the ratio between the Brunt-V\"ais\"al\"a frequency N and the inertial frequency f varying from 1/4 to 20, together with a purely stratified run. The Froude number, measuring the strength of the stratification, varies within the range 0.02 < Fr < 0.32. We find that the inverse cascade is dominated by the slow quasi-geostrophic modes. Their energy spectra and fluxes exhibit characteristics of an inverse cascade, even though their energy is not conserved. Surprisingly, the slow vortices still dominate when the ratio N/f increases, also in the stratified case, although less and less so. However, when N/f increases, the inverse cascade of the slow modes becomes weaker and weaker, and it vanishes in the purely stratified case. We discuss how the disappearance of the inverse cascade of energy with increasing N/f can be interpreted in terms of the waves and vortices, and identify three major effects that can explain this transition based on inviscid invariants arguments.

Corentin Herbert Raffaele Marino Duane Rosenberg Annick Pouquet