2026-04-17 | | Total: 15
Recent investigations of turbulent circulation fluctuations have uncovered substantial insights into the statistical organization of flow structures and revealed unexpected geometric features of turbulent intermittency. Of particular interest here is the observation that circulation probability distribution functions admit a superstatistical representation, namely a description based on "ensembles of Boltzmann-Gibbs ensembles". A fundamental phenomenological ingredient of this approach, which serves as a natural starting point for modeling, relies on the strong correlation between the dissipation field and the spatial distribution of elementary circulation-carrying structures, i.e., small-scale vortices. Within the language of superstatistics, this corresponds to characterizing circulation statistics through an appropriate choice of conditioned (Boltzmann-like) distributions and mixing distributions. We show that the superstatistical class of q-exponentials, known to have broad applicability in a wide range of multiscale and non-equilibrium systems, provides an accurate description of the observed circulation statistics in homogeneous and isotropic turbulence. This finding opens avenues for exploring the statistical structure of the turbulent cascade in the context of non-extensive statistical mechanics, rooted in the concept of non-additive entropies.
Interfacial deformation under electric fields is a common phenomenon in many industrial processes. Particularly, we are interested in the dynamics of sessile soap bubbles in a parallel-plate electric field which exhibits a stable deformation regime followed by conical instability. Using side-view imaging, we track the equilibrium shapes, the transition to the unstable regime, and the pre-jet apex dynamics within one experimental system. In the stable regime, the meridional profile is well described by a spheroidal fit, and the aspect ratio collapses across initial bubble sizes onto a single steady-state branch when plotted against the dimensionless field $E^\ast = \sqrt{\mathrm{Bo}_e}$ for data acquired within a fixed ambient session where the electric Bond number $\mathrm{Bo}_e$ is defined as $\varepsilon_0 E_0^2 R_0/(2γ)$. The endpoint of this branch marks the transition to the unstable regime. Above onset of instability, the apex sharpens into a cone with half-angle $30.0^{\circ}$ $\pm$ $0.6^{\circ}$, below the classical Taylor value. To quantify the late pre-jet stage, we define the axial distance $b(t)$ from the instantaneous apex to a fixed reference vertex determined from the terminal cone geometry and measure its evolution. The corresponding rate grows as jetting is approached, and a near-tip inertia-capillary model captures the observed logarithmic trend as an approximation. Together, these measurements establish a single-system experimental benchmark in which stable electrocapillary deformation is organized by a single steady-state branch that leads into conical instability and pre-jet dynamics.
We analyze a system of bead--spring polymers interacting with Navier--Stokes turbulence to investigate chain--stretching physics and finite-time Lyapunov exponents in ultra--dilute solutions with Weissenberg number \(Wi \approx 80\). They stretch predominantly as material line elements, yet finite deviations arising from elasticity and excluded--volume forces occur with measurable probability. The chain end--to--end distance exhibits a power--law scaling regime. Polymers preferentially sample regions of axisymmetric biaxial extension, where they reach their largest extensions and stretch most rapidly. The degree of stretching is directly correlated with strain intensity, while relaxation events are concentrated in high--enstrophy regions. The chains align strongly with the second strain--rate eigenvector and tend to anti--align with the third; consequently, the second eigenvalue contributes significantly to polymer compression, despite its magnitude typically being smaller than that of the third. Along polymer trajectories, vorticity tends to align with both the first and second eigenvectors, a behavior that differs from the corresponding Eulerian statistics and from Lagrangian vortex--stretching phenomenology. After approximately ten large--eddy turnover times, the Lagrangian Lyapunov exponent histories from different chains appear to synchronise, consistent with convergence of the mean logarithmic stretch rates. All Lyapunov--exponent probability density functions exhibit departures from Gaussianity, and the intermediate Lyapunov exponent is positive in all realizations. The ratio of the mean intermediate to largest exponents is \(E[λ_2]/E[λ_1] \approx 2/7\). The largest and intermediate exponents are positively correlated, whereas the intermediate and smallest exponents are anticorrelated.
An interpretable, physics-consistent turbulence model correction framework, termed FISR-Equation Learner (EQL), is proposed by embedding equation learning directly into a Partial Differential Equations (PDE)-constrained field inversion process based on the adjoint method. Unlike conventional two-stage approaches, the correction model is optimized end-to-end in parameter space using an EQL architecture, enabling the direct identification of compact analytical expressions while maintaining consistency with the governing equations. The method is applied to the shear-stress-transport (SST) model and trained on two canonical separated flows, the curved backward-facing step and the NASA hump. The resulting explicit expression significantly reduces separation bubble overprediction and improves reattachment prediction, achieving performance comparable to neural-network-based end-to-end methods while retaining full interpretability. Generalization is demonstrated on unseen configurations, including periodic hills, a surface-mounted cube, and the high-lift NLR7301 airfoil. The model improves separated-flow predictions and stall characteristics without degrading attached boundary-layer performance. Overall, FISR-EQL provides a practical pathway toward optimal yet transparent data-driven turbulence model correction.
This study presents an efficient and accurate discrete adjoint gas-kinetic scheme (GKS) for sensitivity analysis and aerodynamic shape optimization in continuum flow regimes. Developed using the backward mode of algorithmic differentiation (AD), the adjoint solver is rigorously verified against a duality-preserving linearized GKS solver generated via forward-mode AD. The robustness and practical effectiveness of the solver are evaluated through three benchmark cases: the inverse design of turbine blades, lift-to-drag ratio enhancement, and shock-strength reduction for a NACA 0012 airfoil. To capture realistic flow physics, fully turbulent optimizations are conducted using the one-equation Spalart--Allmaras (SA) model. Numerical results demonstrate excellent agreement between the discrete adjoint and linearized solvers, exhibiting matching sensitivity convergence behaviors, identical asymptotic residual decay rates, and negligible discrepancies in final sensitivity predictions. Furthermore, the optimization studies confirm that targeted design objectives are consistently achieved within a limited number of design cycles, highlighting the solver's computational efficiency, accuracy, and suitability for complex aerodynamic geometries.
We study the navigation of a self-propelled inertial particle in two-dimensional Rayleigh--Bénard convection at Prandtl number $Pr = 0.71$ and cell aspect ratio $Γ= 4$ for Rayleigh numbers $Ra$ ranging from $10^{7}$ to $10^{11}$. A reinforcement-learning (RL) controller selects the propulsive acceleration, subject to an upper bound $\mathcal{A}_{\max}$, to achieve a prescribed horizontal displacement. We find that the success rate increases abruptly with $\mathcal{A}_{\max}$ at moderate $Ra$, whereas at higher $Ra$ the transition becomes more gradual and shifts to larger $\mathcal{A}_{\max}$. Moreover, although the completion time increases with $Ra$, the propulsion energy required for successful traversal decreases. Proper orthogonal decomposition (POD) reveals that these performance differences arise from reorganisation of the carrier flow. At moderate $Ra$, the dominant large-scale circulation partitions the domain through robust transport barriers, requiring a finite thrust surplus to cross them; at higher $Ra$, energy is distributed across many modes, the barriers fragment, and transient plume-assisted pathways emerge. Compared with a constant-heading baseline, the learned policy aligns with local currents and consumes significantly less energy. Lagrangian coherent structure (LCS) analysis further shows that the RL agent inherently learns to cross repelling barriers and surf along attracting pathways. Finally, by mapping these behaviours onto the local Eulerian flow topology using Voronoi tessellation and the $Q$-criterion, we distil an interpretable, physics-based heuristic strategy that achieves robust navigability. These results connect turbulent-flow organisation with autonomous navigation under bounded actuation.
Scyphozoan jellyfish exhibit the highest locomotive efficiency in the animal kingdom making them of particular interest in fluid dynamics and bioinspired robotics. Despite this prevalent analytical models of jellyfish swimming have been based on the swimming traits of hydrozoan jellyfish which utilize jet propulsion, rather than scyphozoan jellyfish which utilize paddling propulsive methods. Additionally, while stroke frequency is a driving variable in speeds achieved by undulatory swimmers, a similar dependence has not been previously explored for jellyfish. This work investigates the relationship between stroke frequency and swimming speeds in two species of scyphozoan jellyfish, Aurelia aurita and Cassiopea xamachana. An experimental study was conducted using a biohybrid technique that controls the muscle contraction frequency of freely swimming, live jellyfish with portable, implanted microelectronics. Swimming speeds were measured from video recordings in a 2.4 m tall water tank. It was found that despite differences in their natural swimming frequencies, the Aurelia and Cassiopea displayed similar speed-frequency relationships with peak swimming speeds occurring at 0.55 +/- 0.05 Hz and 0.50 +/- 0.05 Hz respectively. The difference in natural stroke frequency displayed by scyphomedusea despite the shared relationship between swimming speed and stroke frequency in these two species, suggests that natural stroke frequency may be more related to other functions such as filter feeding, rather than locomotion. A new analytical model developed for scyphozoan, paddling jellyfish was shown to have closer agreement with the experimental results than existing models based on jet propulsion. The model demonstrated the driving factors in the relationship between swimming speed and stroke frequency to be the speed of the jellyfish bell margin and changes in body kinematics with stroke frequency.
Self-propelled particles can navigate complex environments, including viscous fluid interfaces with curved geometries. In this work, we study the emergent dynamics of a suspension of self-propelled particles confined to a stationary curved viscous interface. The evolution of the particle configurations is modeled using the Fokker-Planck equation on the curved surface, formulated using Cartan's moving frame method, and coupled to the bulk and surface Stokes equations with flows driven by an interfacial nematic active stress. Specifically, for a spherical vesicle, the flow field and the distribution of the particles are analyzed theoretically and numerically within the framework of spin-weighted functions and spin-weighted spherical harmonics, which provide a natural geometric description of the probability distribution function on the sphere. A linear stability analysis about the uniform, isotropic state is performed and predicts a finite-wavelength instability, with mode selection arising from the competition between the vesicle radius and the Saffman-Delbrück length. This instability and the associated mode-selection mechanism are also confirmed in nonlinear numerical simulations using a pseudo-spectral method based on spin-weighted spherical harmonics.
Rotating detonation engines (RDEs) are a promising propulsion concept that may offer higher thermodynamic efficiency and specific impulse than conventional systems, but nonlinear phenomena, including transitions to oscillatory or chaotic propagation modes, can hinder practical operation. Deep Reinforcement Learning (DRL) has emerged as a promising method for controlling complex nonlinear dynamics such as those observed in RDEs. However, the multi-timescale nature of the RDE system makes direct application of DRL challenging. We address this challenge by reformulating the DRL problem in a moving reference frame that follows the detonation-wave pattern, making the wave structure appear quasi-steady to the agent. This reformulation enables scale separation between fast detonation propagation and slower operating-mode dynamics. We train DRL controllers to modulate spatially segmented injection pressure in a one-dimensional reduced-order RDE model and induce rapid transitions between different mode-locked states. Across a range of actuation periods, initial states, and target modes, controllers trained in the moving frame learn more reliably than those trained in a stationary frame and remain effective over a broader range of actuation periods. These results suggest that symmetry-aware moving reference frame formulations may be useful for related multiscale flow-control problems and that scale separation should be exploited whenever possible to enable DRL control of multi-timescale systems.
We investigate the thermo solutal transport phenomena and deposition patterns during the evaporation of surfactant laden droplets experimentally and through theoretical scaling based analysis. Experiments were conducted using the sessile droplet configuration in the acrylic chamber for both hydrophilic and hydrophobic substrates. Infrared thermography and particle image velocimetry measurements were conducted during evaporation to illustrate the temperature and velocity distributions, respectively. Sodium dodecyl sulphate SDS surfactant molecules enhanced the evaporation rate with an increase in concentration for the hydrophobic surface. In contrast, the evaporation rate increased up to 0.5 CMC and then decreased for droplets on a hydrophilic substrate. The evaporation rates computed from the shadowgraphy imaging were explained using the average velocities obtained from the PIV analysis. It was found that advection within the droplet is strongly dependent on surfactant concentration and wettability. Further, the theoretically obtained Marangoni velocities were in close agreement with the experimental values. It was found that Marangoni solutal advection dominates other advection mechanisms, such as Marangoni thermal advection and buoyancy driven flow. However, surfactant crowding and viscous resistance with increasing surfactant concentration can dampen the increase in solutal advection. The surface tension and viscosity measurements were also conducted with variation in surfactant concentration to understand the suppression of advection by viscous forces. The computation of contact line velocities showed sudden fluctuations, illustrating stick slip behaviour during droplet drying, complementing microscopic visual observations.
We present a numerical investigation of kerosene droplet mist film cooling for the thermal protection of the rotating detonation combustor (RDC) and compare its performance with conventional air film cooling and combined mist/air cooling scheme. In the study, the cooling behavior of kerosene droplets injected through wall film holes is numerically examined and compared with air film cooling and a combined mist/air cooling strategy, building on a benchmark validation against flat-plate experimental data. The results show that air film cooling exhibits an optimal operating range, beyond which excessive injection degrades film stability due to strong interaction with the rotating detonation wave. In contrast, kerosene-based mist cooling forms a more persistent near-wall cooling layer, providing enhanced heat removal through phase change and exhibiting improved resistance to film separation. In mist cooling, the droplet size primarily affects the immediate downstream cooling performance, with intermediate-sized droplets offering the improved balance between evaporation rate and film continuity. A combined mist/air cooling scheme can further improve cooling efficiency and accelerate wall temperature recovery after detonation wave passage while maintaining moderate impacts on the mainstream flow. Additionally, although kerosene droplets partially participate in combustion under film hole injection, the associated thermal load does not offset the overall cooling benefit. These findings demonstrate the feasibility and advantages of kerosene-based cooling schemes for RDC thermal management.
Unstructured meshes are among the most versatile approaches for capturing non-canonical geometries in fluid dynamics simulations. Despite this, most high-fidelity first-principles phase-change models are developed and applied on structured meshes. We present a phase-change simulation method for unstructured meshes that combines the algebraic Volume-of-Fluid (VOF) technique with geometric interface reconstruction, implemented in an in-house open-source CFD code. Phase-change rates are computed from local temperature gradients evaluated at the reconstructed interface, without empirical closure models, using a reconstruction procedure that operates on arbitrary polyhedral cells. Because the method relies on the standard finite-volume framework, it can be integrated into other cell-centred codes supporting unstructured meshes. The approach is validated against the one-dimensional Stefan and Sucking problems and the three-dimensional Scriven bubble growth on both hexahedral and polyhedral meshes, showing good agreement with analytical solutions in all three cases. A detailed analysis of the Scriven problem reveals that the interface-modified least-squares gradient stencil on Cartesian meshes overestimates the interfacial temperature gradient, producing a persistent overshoot of the analytical bubble radius and a coherent four-fold anisotropy that elongates the bubble along grid diagonals. On polyhedral meshes, the irregular face orientations eliminate both effects, yielding isotropic growth and monotonic convergence. Finally, we demonstrate the framework on turbulent upward co-current annular boiling flow, where early transient results are qualitatively consistent with a previous LES study and experimental observations of wave-modulated evaporation.
A scale-by-scale analysis of energy flux in the turbulent cascade can be performed using the spatially filtered magnetohydrodynamic (MHD) equations, while the gradient tensor invariants are widely used to characterise the structure of velocity and magnetic fields. Physical mechanisms responsible for energy flux require specific field configurations whose strength is quantified by these tensor invariants. We explore this requirement, showing that the tensor invariants act as proxies for mechanistic energy fluxes under quantifiable conditions. As a special case, the purely hydrodynamic contributions to energy flux can be expressed exactly in terms of the invariants of the velocity gradient tensor. We also show that the invariants bound the available energy flux for distinct physical mechanisms, formalising the idea that each transfer mechanism requires field configurations with gradients of sufficient strength to support a given energy flux. Results are illustrated using 3D simulations of freely decaying MHD turbulence.
Most practical engineering design problems involve nonlinear spatio-temporal dynamical systems. Multi-physics simulations are often performed to capture the fine spatio-temporal scales which govern the evolution of these systems. However, these simulations are often high-fidelity in nature, and can be computationally very expensive. Hence, generating data from these expensive simulations becomes a bottleneck in an end-to-end engineering design process. Spatio-temporal surrogate modeling of these dynamical systems has been a popular data-driven solution to tackle this computational bottleneck. This is because accurate machine learning models emulating the dynamical systems can be orders of magnitude faster than the actual simulations. However, one key limitation of purely data-driven approaches is their lack of generalizability to inputs outside the training distribution. In this paper, we propose a physics-informed spatio-temporal surrogate modeling (PISTM) framework constrained by the physics of the underlying dynamical system. The framework leverages state-of-the-art advancements in the field of Koopman autoencoders to learn the underlying spatio-temporal dynamics in a non-intrusive manner, coupled with a spatio-temporal surrogate model which predicts the behavior of the Koopman operator in a specified time window for unknown operating conditions. We evaluate our framework on a prototypical fluid flow problem of interest: two-dimensional incompressible flow around a cylinder.
Steady-state electrothermal systems involve strongly coupled heat transfer, fluid flow, and electric-potential transport, creating severe numerical challenges for standard physics-informed neural networks (PINNs) due to stark disparities in gradient scales and residual stiffnesses across the physical fields. To resolve these multiphysics bottlenecks, we introduce a Long Short-Term Memory PINN (LSTM-PINN) framework that utilizes a depth-recursive memory mechanism to preserve long-range spatial feature dependencies and maintain strict cross-field consistency. The proposed architecture is rigorously evaluated against conventional and attention-based networks across a unified five-field formulation encompassing four complex convective and drag regimes: Boussinesq electrothermal flow, drift-potential gauge-constrained transport, strong buoyancy-coupled convection, and Brinkman--Forchheimer drift. Quantitative and visual analyses demonstrate that LSTM-PINN successfully suppresses non-physical artifacts and structural distortions, yielding the highest thermodynamic fidelity and consistently outperforming state-of-the-art baselines in global error metrics. Ultimately, this memory-enhanced approach provides a highly robust and accurate computational baseline for capturing localized boundary layers and complex energy-momentum feedback in advanced electrothermal energy systems.