2026-04-15 | | Total: 17
Recent advances in data-driven modeling have shown that diffusion models can successfully generate synthetic Lagrangian trajectories in turbulent flows. Building on this progress, we extend the method to the joint generation of pairs of Lagrangian velocity trajectories, enabling a fully data-driven representation of turbulent pair dispersion, a long-standing fundamental problem with broad relevance in fluid dynamics. We demonstrate that diffusion models accurately reproduce the evolution of particle-pair separation, including deviations from Richardson's classical scaling law, while simultaneously preserving all key single-particle statistical properties reported in previous studies. These findings underscore the potential of diffusion-based generative models to emulate high-dimensional, multi-scale turbulent dynamics, further establishing them as a powerful tool for scientific modeling and for future geophysical and astrophysical applications.
Neural operators have been increasingly used as data-driven surrogates for time-marching predictions of turbulent flows. However, long-horizon autoregressive prediction is sensitive to error accumulation and the choice of prediction interval. Excessively small time increments may increase temporal redundancy and lengthen rollouts, which can degrade the stability of neural operators in turbulence forecasting. This work pursues a unified objective: stable long-horizon autoregressive prediction at fine temporal resolution for three-dimensional turbulence. We propose a multi-stepsize mixture-of-experts (Ms-MoE) neural operator built on an implicit factorized Transformer (IFactFormer) backbone. The model conditions on a requested relative stride and uses a time-step router to activate scale-specific routed experts together with a shared expert, yielding a single architecture that represents a family of stride-parameterized time-advancement operators. We evaluate the approach on forced homogeneous isotropic turbulence (HIT) and turbulent channel flow using filtered direct numerical simulation datasets. Relative to sampling intervals used in previous studies, we construct training datasets with up to 20 times finer temporal resolution and report long-horizon autoregressive rollouts using qualitative time-slice comparisons and long-time-averaged statistics. Ms-MoE-IFactFormer yields more stable long-horizon rollouts and improved agreement with long-time-averaged statistics on both HIT and turbulent channel flow, suggesting potential for stable time-marching at fine temporal resolution in more complex turbulent flows.
We report on the shape taken by the interface of a liquid bath when hit by a smooth oblique steady jet. When the angle between the jet and the bath decreases below $50^\circ$, a cavity is formed in front of the jet. In the inertial regime we explore, the jet boundary layer detaches in the impact region, thereby delimiting a core jet region outside of which the liquid is mainly in hydrostatic equilibrium. The shape of the outer meniscus is shown to be related to the one outside a tilted fiber piercing the fluid interface. In order to unravel the flow features and separation, we perform direct numerical simulations and show that the flow detachment displays an asymmetry, which results in the acceleration of the liquid below the surface, thereby creating a depression. With this observation, we propose a model balancing the suction force of this depression with the weight of the displaced water and the surface tension force to obtain a prediction for the typical width of the cavity.
This work proposes a statistically enhanced framework to address the instability and limited predictive capability of conventional Galerkin-Proper Orthogonal Decomposition (Galerkin-POD) models. The method reformulates the correction of the Galerkin-projected ODE system as a statistical inverse problem, in which the coefficients are inferred through Bayesian inference. By accounting for model uncertainty arising from POD mode truncation and data uncertainty introduced by data noise and numerical postprocessing, the framework systematically updates the ODE system coefficients using an analytical, sampling-free solution based on Gaussian likelihood and inverse-Gamma priors. The approach is first validated using a self-sustained oscillating flow over a dimpled surface at a moderate Reynolds number (Re=3000), demonstrating stable and accurate reproduction of the temporal dynamics and phase trajectories of coherent structures when compared with direct numerical simulation (DNS). It is then applied to a centrifugal compressor featuring strong tip-leakage vortex breakdown and impeller-diffuser interactions at Re=100000, where the model successfully captures dominant unsteady structures and frequency characteristics despite limited mode retention. Overall, the results show that Bayesian inference substantially enhances the robustness, stability, and predictive fidelity of Galerkin-POD models for compressible flow systems. The proposed methodology combines the physical interpretability of Galerkin projection with the statistical rigour of Bayesian inference, offering a general, computationally efficient, and uncertainty-aware reduced-order modelling framework for complex fluid dynamic applications.
The flow of a liquid metal (LM) in a rectangular duct segment, subject to a uniform transverse magnetic field and uniform heating at the side walls is explored in an ample parameter space using Direct Numerical Simulation (DNS). We modify electrical wall conductivity, (either highly conducting or perfectly insulating) and investigate the effects of the buoyancy force, both in horizontally and vertically orientated ducts. In the latter case, it may be directed either with the flow or against the flow, creating backflow regions. In this parameter space and with the presence of vortex promoters at the inlet of the duct we identify $4$ types of flow. We calculate the Nusselt number $Nu(t)$ for each of them and study the statistical properties to compare their heat transfer capabilities in future fusion reactor blankets.
As one of the cornerstones of fluid mechanics, stability analyses provide essential physical insights into the growth of perturbations and eventual transition to turbulence. However, classical \enquote{frozen-time} stability analyses implicitly assume a time-independence of their base flow and thus fail for \enquote{rapidly} diffusing shear layers. Here, we propose a self-similar ansatz to naturally incorporate the \enquote{diffusive} base-state expansion into the stability operator. Our approach reveals two competing physical mechanisms: an \enquote{expansion wind} delays the Kelvin-Helmholtz instability whereas a diminishing effective viscosity sustains this instability far beyond classical predictions. Direct numerical simulations confirm that our framework accurately captures the instability's extended lifespan, growth rate, and spectral topology, eventually revising the timeline of shear-induced mixing fundamentally.
Dynamic soaring enables sustained flight by extracting energy from wind shear, yet it is commonly understood as a cycle-level maneuver that assumes stable flow conditions. In realistic unsteady environments, however, such assumptions are often violated, raising the question of whether explicit cycle-level planning is necessary. Here, we show that dynamic soaring can emerge from step-level, state-feedback control using only local sensing, without explicit trajectory planning. Using deep reinforcement learning as a tool, we obtain policies that achieve robust omnidirectional navigation across diverse shear-flow conditions. The learned behavior organizes into a structured control law that coordinates turning and vertical motion, giving rise to a two-phase strategy governed by a trade-off between energy extraction and directional progress. The resulting policy generalizes across varying conditions and reproduces key features observed in biological flight and optimal-control solutions. These findings identify a feedback-based control structure underlying dynamic soaring, demonstrating that efficient energy-harvesting flight can emerge from local interactions with the flow without explicit planning, and providing insights for biological flight and autonomous systems in complex, flow-coupled environments.
We study the modulational stability problem for the traveling periodic waves (called Stokes waves) in an infinitely deep fluid by using pseudo-differential operators in conformal variables. We derive the criteria and the normal forms for four bifurcations which are repeated recurrently when the steepness of the Stokes wave is increased towards the highest wave with the peaked profile. The four bifurcations are observed in the following order: (a) new figure-8 bands appearing at each extremal point of speed, (b) degeneration of figure-8 bands resulting in vertical slopes, (c) new circular bands around the origin appearing at each period-doubling bifurcation, and (d) reconnection of figure-$\infty$ bands at each extremal point of energy. Our work uses the analytic theory of Stokes waves developed previously for Babenko's equation. The novelty of our work is the analytic extension of the modulational stability problem for singular pseudo-differential operators in terms of the Floquet parameter. The derivation of the normal form uses some structural assumptions which are known to be true for the Stokes waves. For the first and second bifurcation cycles, we compute numerically with a higher-order accuracy the actual values of wave steepness for which the structural assumptions are satisfied and the numerical coefficients of the normal forms to show the excellent agreement between the normal form theory and the numerical approximations of the spectral bands.
Multiscale Proper Orthogonal Decomposition (mPOD) decomposes fluid flows into energy-optimal modes within prescribed frequency bands by combining Proper Orthogonal Decomposition with a multiresolution analysis (MRA). In its classical formulation, mPOD relies on a filter bank of finite impulse response (FIR) filters, enabling lossless reconstruction while mitigating Gibbs oscillations and temporal ringing. However, the smooth transition bands required for this purpose introduce partial spectral overlap between adjacent scales and require, for each band, the solution of an eigenvalue problem spanning the full temporal dimension. This work introduces a fast spectral formulation of the mPOD that substantially reduces the computational cost. The proposed approach replaces time-domain FIR filters with compact spectral masks enforcing strictly disjoint frequency supports, thereby exactly decoupling the problem across scales. This leads to a block-diagonal correlation operator in spectral space, so that each band can be treated independently. The resulting eigenvalue problems reduce to small systems whose size depends on the number of active frequencies per band rather than the full time dimension. The approach is validated on a synthetic dataset highlighting spectral windowing effects and on experimental particle image velocimetry (PIV) data of a cylinder wake at Reynolds number \(\mbox{Re} \approx 5000\). In both cases, the proposed formulation accurately recovers the modal structures and singular values of the classical mPOD while reducing the computational cost by orders of magnitude.
Turbulent channel flow controlled by spanwise wall oscillations is studied using direct numerical simulations to improve how spanwise forcing reduces skin-friction drag. Harmonic wall oscillations generate a periodic transverse Stokes layer whose thickness $δ$ is determined by the forcing period $T$. Although an optimal $T$ that maximizes drag reduction is known to exist, its physical significance remains unclear. To elucidate it, we extend the spanwise Stokes layer by augmenting wall oscillation with an additional spanwise body force. In this formulation, $δ$ and $T$ become decoupled and can be varied independently. The oscillating wall thus appears as a special and suboptimal case of spanwise forcing. Optimal performance is obtained for substantially smaller $T$ and larger $δ$ than those of the classical Stokes layer. For the conditions examined, with Reynolds number and forcing amplitude held fixed, the maximum drag reduction increases by approximately one third, while the maximum net energy saving improves markedly from $-35\%$ to $+16\%$. These findings suggest that drag-reduction strategies based on spanwise forcing deserve renewed scrutiny: wall oscillation represents only one possible actuation method, and not necessarily the most effective one.
Recent study \cite{rajamanickam2025shear} of non-Boussinesq fluids in narrow channels identified a novel shear-induced horizontal buoyancy force that emerges upon depth-averaging the Navier--Stokes equations. This note demonstrates that this force is formally equivalent to the divergence of a Korteweg stress tensor. Unlike classical Korteweg stresses, which are typically attributed to molecular-scale cohesive potentials or implemented through assumed constitutive relations, we show that this emergent stress arises purely from self-coupled transport where the internal Ostroumov flow is "enslaved" to the local density gradient. We derive explicit expressions for the effective stress coefficients, revealing a fundamental dependence on the Prandtl number and Grashof number and identifying a transition in the effective internal pressure at $Pr=1/2$, which marks the crossover between the internal inertia of the shear flow and the hydrostatic tilting induced by the shear. This correspondence is contrasted with classical Taylor dispersion, where the absence of self-coupling yields only a uniaxial stress. Our results suggest that quadratic Korteweg-type stresses may be a universal manifestation of sub-scale transport in gradient-driven flows, providing a rigorous macro-scale origin for capillary-like stresses in miscible fluids.
Uncertainty in estimating the log-law parameters is arguably the greatest obstacle to establishing definitive conclusions regarding their numerical values and universality. This challenge is exacerbated by the limited number of studies that provide thorough uncertainty analyses of experimental data and fitting procedures, and those that do often adopt different approaches, undermining direct comparisons. The present study applies the generalised least squares (GLS) principle to the log-law velocity profile to establish a standardised, comprehensive framework for quantifying uncertainty in the log-law parameters across datasets. GLS contrasts with ordinary least squares (OLS) and weighted least squares (WLS), which do not account for correlation in errors across measured quantities, as well as with alternative heuristic methods that independently sample primitive variables. Instead, it incorporates a full covariance matrix of the residuals, propagated from the uncertainties in the primitive variables and consistent with the experimental methods employed. The study presents a systematic analysis of the response of the log-law regression model using synthetic data, emulating measurements from a hot-wire anemometer mounted on a linear traverse. This analysis serves as a predictive tool for experimental design, identifying a priori the dominant sources of uncertainty in the log-law parameters and potential mitigation strategies. The study also provides new insights into the correlation between the log-law parameters and proposes a new fitting procedure that eliminates the need to prescribe the location and extent of the log region. The open-source Python implementation of the log-law regression model is available for download on GitHub at https://github.com/ma2ferreira/gls_loglaw.git.
A new radiation solver, RAPRAL (RAdiation Prediction based on RAy tracing and Line-by-line) implemented in C++, is developed for simulating high-temperature thermochemical nonequilibrium radiative processes. RAPRAL integrates detailed line-by-line spectral modeling with a ray-tracing solution of the radiative transfer equation, enabling accurate resolution of both spectral features and spatial radiation transport. The adopted methods and their implementation are described in detail. To assess the overall capability and accuracy of RAPRAL, we first focus on the computation of atomic and molecular bulk spectral coefficients. Through comparison with the established code in the literature, RAPRAL demonstrates its ability to accurately capture key spectral features across a wide range of conditions. Moreover, RAPRAL is applied to predict afterbody radiative heating in the Fire II flight experiment, based on a two-temperature, 11-species air flowfield. The results demonstrate that the present approach provides reliable predictions of radiative heat flux and effectively captures the dominant radiation mechanisms. Overall, the presented results demonstrate that RAPRAL is a robust tool for simulating radiative processes in hypersonic air flows, and future versions will extend its capabilities to include species relevant to planetary atmospheres.
We propose a theory based on dynamical systems to explain and predict the occurrence of extreme events, of which critical transitions form a subset. In fast-slow nonlinear systems, we identify a cascade of events preceding extreme events: (i) a slow regime, in which the fast covariant Lyapunov vectors (CLVs) are both tangent to the fast eigenvectors and remain transversal to the slow subspace; (ii) a transition regime, in which the fast eigenvalues become neutrally stable while the fast CLVs are no longer tangent to the fast eigenvectors; and (iii) a critical regime, in which a strong spectral gap in the eigenvalues causes both fast and slow CLVs to become tangent along the dominant fast direction, breaking the transversality between fast and slow subspaces. Building on this cascade, we propose two precursors to forewarn the occurrence of extreme events. We numerically test the theory and precursors on low- and higher-dimensional systems. The proposed precursors predict extreme events and critical transitions with 100% precision and recall. This work opens opportunities for time-forecasting extreme events using theoretically grounded precursors.
Context. Be stars are presently viewed as B-type stars surrounded by a disc fueled by the star itself during episodicexcretion events. The origin of these events are poorly understood.Aims. This article aims to determine whether or not surface equatorial Kelvin waves can be unstable and therefore canplay a role in the triggering of the Be phenomenon.Methods. We first derive an analytical expression for gravito-inertial modes in the shallow-water framework. Then, weinvestigate numerically the evolution of equatorial Kelvin modes as system parameters vary. The study is extended tothick-layer configurations with a constant density fluid. We then analyze the stability of these modes under differentialrotation and viscous effects.Results. We show that equatorial Kelvin waves still exist in a spherical shell of finite thickness, but that their equatorialconfinement is weaker. At low azimuthal wavenumbers, Kelvin waves are in the inertial waves frequency band and thusget specificities of inertial waves like shear layers associated with singularities of the Poincaré equation. These shearlayers are new dissipative structures for Kelvin waves. When a radial (shellular) differential rotation is imposed, we showthat equatorial Kelvin waves can be destabilised provided that differential rotation and viscosity are in an appropriaterange. The non-monotonic behaviour of the growth rate of the instability is traced back to the rise of a critical layerwhere the fluid azimuthal velocity equals the phase speed of the surface waves.Conclusions. This study provides new insights into the behavior of equatorial Kelvin waves in astrophysics, particularlyin rapidly rotating stars. The results reinforce the idea that gravito-inertial waves, and more specifically the equatorialKelvin waves, can be unstable and thus be key parts in the mechanisms leading to the Be phenomenon.
The Euler ensemble was recently supported by large-scale ($4096^3$) direct numerical simulations as the universal statistical attractor of decaying fluid turbulence. Previous mathematical derivations of this ensemble relied on measure-theoretic limits of discrete polygonal loop equations. In this Letter, we present a continuous algebraic derivation. By reformulating the Navier-Stokes equation as a covariant derivative operator flow in the Lagrangian frame, we analytically eliminate advection. Applying Feynman's operational calculus, the 3D non-commutative operator algebra maps to ordering discontinuities (finite-difference jumps) on a 1D momentum loop. This continuous formalism reduces to the discrete, number-theoretic geometric quantization of the Euler ensemble via roots of unity without requiring spatial lattice approximations, demonstrating that macroscopic fluid chaos is a deterministic projection of the Farey sequence.
Turbulent mixing remains one of the primary uncertainties in the modeling of stellar interiors. In stellar evolution simulations, regions where mixing occurs are typically identified using instability criteria. A particularly interesting situation arises when nuclear reactions produce inversions in the mean molecular weight within stellar interiors. Under these conditions, the material can become unstable to either thermohaline or a Rayleigh-Taylor instabilities. We demonstrate that the standard criterion adopted in stellar evolution calculations does not accurately distinguish between these two regimes. We derive an alternative criterion and show that chemically driven convection in stellar interiors might be viable under much smaller mean molecular weight inversions than it is normally assumed. We investigate whether inversions in the mean molecular weight can trigger chemically driven convection above the red giant branch bump (RGBB) or during the helium core flash. We find that the inversion at the base of the convective envelope above the RGBB is too weak and short-lived to sustain steady-state convection. In contrast, rapid carbon production at the base of the He-flash-driven convective zone can maintain a steady chemically driven convective region. This process could significantly alter our understanding of the He-core flash and warrants further study.