25 Search Results

Abstract Ocean weather comprises vortical and straining mesoscale motions, which play fundamentally different roles in the ocean circulation and climate system. Vorticity determines the movement of major ocean currents and gyres. Strain contributes to frontogenesis and the deformation of water masses, driving much of the mixing and vertical transport in the upper ocean. While recent studies have shown that interactions with the atmosphere damp the ocean’s mesoscale vortices O (100) km in size, the effect of winds on straining motions remains unexplored. Here, we derive a theory for wind work on the ocean’s vorticity and strain. Using satellite and model data, we discover that wind damps strain and vorticity at an equal rate globally, and unveil striking asymmetries based on their polarity. Subtropical winds damp oceanic cyclones and energize anticyclones outside strong current regions, while subpolar winds have the opposite effect. A similar pattern emerges for oceanic strain, where subtropical convergent flow is damped along the west-equatorward east-poleward direction and energized along the east-equatorward west-poleward direction. These findings reveal energy pathways through which the atmosphere shapes ocean weather.

The mesoscale spectrum describes the distribution of kinetic energy in the Earth's atmosphere between length scales of 10 and 400 km. Since the first observations, the origins of this spectrum have been controversial. At synoptic scales, the spectrum follows a –3 spectral slope, consistent with two-dimensional turbulence theory, but a shallower –5/3 slope was observed at the shorter mesoscales. The cause of the shallower slope remains obscure, illustrating our lack of understanding. Through a novel coarse-graining methodology, we are able to present a spatio-temporal climatology of the spectral slope. We find convection and orography have a shallowing effect and can quantify this using “conditioned spectra.” These are typical spectra for a meteorological condition, obtained by aggregating spectra where the condition holds. This allows the investigation of new relationships, such as that between energy flux and spectral slope. Potential future applications of our methodology include predictability research and model validation.

The lunar surface contains a wide variety of topographic shapes and features, each with different distributions and scales, and any analysis technique to objectively measure roughness must respect these qualities. Coarse-graining is a naturally scale-dependent filtering technique that preserves scale-dependent symmetries and produces coarse elevation maps that gradually erase the smaller features from the original topography. In this study of the lunar surface, we present two surface variability metrics obtained from coarse-graining lunar topography: fine elevation and coarse curvature. Both metrics are isotropic, deterministic, slope-independent, and coordinate-agnostic. Fine (detrended) elevation is acquired by subtracting the coarse elevation from the original topography and contains features that are smaller than the coarse-graining length-scale. Coarse curvature is the Laplacian of coarsened topography, and naturally quantifies the curvature at any scale and indicates whether a location is elevated or depressed relative to its neighborhood at that scale. We find that highlands and maria have distinct roughness characteristics at all length-scales. Our topographic spectra reveal four scale-breaks that mark characteristic shifts in surface roughness: 100, 300, 1,000, and 4,000 km. Comparing fine elevation distributions between maria and highlands, we show that maria fine elevation is biased toward smaller-magnitude elevations and that the maria–highland discrepancies are more pronounced at larger length-scales. Here, we also provide local examples of selected regions to demonstrate that these metrics can successfully distinguish geological features of different length-scales.

The Rayleigh-Taylor Instability (RTI) is a ubiquitous flow phenomenon occurring in a multitude of natural and engineered systems in which buoyancy forces exist. RTI is primarily a limiting process in many HED applications and there are significant efforts within the NNSA-ICF program to account for its effect in large scale simulations of ICF implosions both in two- and three-dimensions. The goal of the proposed research is to develop a fundamental understanding of the evolution and role of magnetic fields in the ablative Rayleigh–Taylor instability (aRTI) in the deeply nonlinear stage.

While kinetic helicity is not Galilean invariant locally, it is known that its spatial integral quantifies the degree of knottedness of vorticity field lines. Being a topological property of the flow, mean kinetic helicity is Galilean invariant. Furthermore, we provide a direct mathematical proof that kinetic helicity is Galilean invariant when spatially integrated over regions enclosed by vorticity surfaces, i.e., surfaces of zero vorticity flux. We also discuss so-called relative kinetic helicity, which is Galilean invariant when integrated over any region in the flow.

Abstract The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length‐scale of interest. To carry out a length‐scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse‐grained. At length‐scales greater than 1,000 km, the dynamics closely follow the Topographic‐Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length‐scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic‐Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of km in size, and the associated bottom pressure anomalies in local vorticity balances at length‐scales smaller than 1,000 km. Overall, there is a transition from the Topographic‐Nonlinear regime at scales smaller than 1,000 km to the Topographic‐Sverdrup regime at length‐scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea‐land boundaries, where bottom stress and horizontal viscous friction generally peak.

Here, we generalize the “filtering spectrum” to probe scales along different directions by spatial coarse-graining. This multidimensional filtering spectrum quantifies the spectral content of flows that are not necessarily homogeneous. From multidimensional spectral information, we propose a simple metric for shape anisotropy at various scales. The method is applied to simulations of two-dimensional (2D) and 3D Rayleigh-Taylor (RT) turbulence, which is inhomogeneous and anisotropic. We show that 3D RT has clear shape anisotropy at large scales with approximately 4:3 vertical to horizontal aspect ratio but tends toward isotropy at small scales as expected. In sharp contrast, we find that RT in 2D simulations, which are still the main modeling framework for many applications, is isotropic at large scales and its shape anisotropy increases at smaller scales where structures tend to be horizontally elongated. While this may be surprising, it is consistent with recent results that large-scale isotropy in 2D RT is due to the generation of a large-scale overturning circulation via an upscale cascade, while small-scale anisotropy is due to the stable stratification resultant from such overturning and the inefficient mixing in 2D.

Here, we present an estimate for the ocean's global scale transfer of kinetic energy (KE), across scales from 10 to 40,000 km. Oceanic KE transfer between gyre scales and mesoscales is induced by the atmosphere’s Hadley, Ferrel, and polar cells, and the intertropical convergence zone induces an intense downscale KE transfer. Upscale transfer peaks at 300 gigawatts across mesoscales of 120 km in size, roughly one-third the energy input by winds into the oceanic general circulation. Nearly three quarters of this “cascade” occurs south of 15°S and penetrates almost the entire water column. The mesoscale cascade has a self-similar seasonal cycle with characteristic lag time of ≈27 days per octave of length scales; transfer across 50 km peaks in spring, while transfer across 500 km peaks in summer. KE of those mesoscales follows the same cycle but peaks ≈40 days after the peak cascade, suggesting that energy transferred across a scale is primarily deposited at a scale four times larger.

We demonstrate a methodology for diagnosing the multiscale dynamics and energy transfer in complex HED flows with realistic driving and boundary conditions. The approach separates incompressible, compressible, and baropycnal contributions to energy scale-transfer and quantifies the direction of these transfers in (generalized) wavenumber space. We use this to compare the kinetic energy (KE) transfer across scales in simulations of 2D axisymmetric vs fully 3D laser-driven plasma jets. Using the FLASH code, we model a turbulent jet ablated from an aluminum cone target in the configuration outlined. We show that, in addition to its well known bias for underestimating hydrodynamic instability growth, 2D modeling suffers from significant spurious energization of the bulk flow by a turbulent upscale cascade. In 2D, this arises as vorticity and strain from instabilities near the jet's leading edge transfer KE upscale, sustaining a coherent circulation that helps propel the axisymmetric jet farther (⁠ ≈25% by 3.5 ns) and helps keep it collimated. In 3D, the coherent circulation and upscale KE transfer are absent. Here, the methodology presented here may also help with inter-model comparison and validation, including future modeling efforts to alleviate some of the 2D hydrodynamic artifacts highlighted in this study.

Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (~1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Lastly, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.