2026-09-30 | | Total: 16
Turbulent flows are highly chaotic, which makes their instantaneous states difficult to predict. Data assimilation aims to reduce this predictive uncertainty by synthesizing the predictions of a physical solver with partial observations of the system. Traditional data assimilation methods such as ensemble and variational approaches can incur substantial computational cost through evaluations of high-fidelity solvers and, for variational methods, adjoint calculations. Recently, generative diffusion models have been explored in the solution of inverse problems in fluid mechanics, demonstrating their potential for data assimilation. However, existing diffusion-based methods often use free parameters instead of the prescribed observation-error covariance to control the balance between the prior and likelihood, departing from the Bayesian formulation of data assimilation. We propose a framework for principled generative data assimilation in which observations are assimilated in the clean-state space and the prescribed observation-error covariance enters explicitly through an ensemble Kalman update. In the Lorenz-63 system, the proposed method yields a lower trajectory error and equation residual than the other diffusion-based posterior sampling methods considered, although its ensemble underestimates the posterior uncertainty. In two-dimensional Rayleigh-Bénard convection, the generated fields recover the statistics of the direct numerical simulation reference at large scales, with deviations confined to the smaller scales. These results suggest that the framework can produce physically plausible turbulent flow fields, indicating its potential for more challenging turbulent flow applications.
A droplet breaks up in tens of microseconds, and a recording captures perhaps a dozen frames. The states in between cannot be recovered without repeating the experiment, and simulating them is too costly to sweep an operating envelope. Yet they are present in the corpus as a whole: a campaign spanning a device's actuation range produces morphologies resembling those any single recording missed. Exploiting that requires a representation that is low-dimensional, invertible, and faithful to shape rather than to sampling. Proper orthogonal decomposition supplies the first two but measures distance in sampled coordinates; manifold learning supplies the third but no map back to a shape; elastic shape analysis supplies a shape metric but no reduced coordinates. SHROM composes all three. Interfaces are represented by their square-root velocity functions, a neighbour graph is built over that shape space, and an autoencoder is trained to reconstruct while penalising latents in which graph neighbours are not latent neighbours. The demonstration uses 301,539 inkjet droplet contours and four filmed break-up sequences. The graph term does not improve reconstruction. It determines whether position in the latent carries meaning: clustering the latent of an otherwise identical model recovers the shape-space partition at chance level (ARI = 0.063 +/- 0.059), and at 0.781 +/- 0.049 with the term active. Waveform parameters predict the full contour at R^2 = 0.878 +/- 0.026. Negative results are reported in the same terms. The graph metric proved immaterial across five choices, and interpolation error saturates at the reconstruction limit, so a plain autoencoder leads that task. The main limitation is interpolation across a topology change: no component of the regularised latent holds both a single-component and a post-break-up shape, whereas an unregularised autoencoder mixes them freely.
Adsorption and desorption redistribute soluble surfactant between a liquid film and its free surface, but cannot change the total inventory in a closed periodic domain. We audit this requirement in nonlinear weighted-residual models for soluble-surfactant-laden falling films. Although a previous correction restores conservation of the nonlinear surface Marangoni flux, a separate inconsistency remains in the depth-integrated bulk diffusion. We derive the resulting mass imbalance and construct explicit positive periodic concentration fields that exhibit it. The missing contribution arises from the streamwise concentration gradient appearing in the normal diffusive flux at a sloping interface. Including this contribution restores local conservation of the coupled bulk-surface transport system while leaving the hydrodynamic equations, the leading-order adsorption-desorption closure, and the linear stability properties of the uniform equilibrium unchanged. In an independent no-slip calculation using the parameters of an earlier travelling-wave study, the reference model loses $0.0228832\%$ of its initial surfactant mass by $t=4500$, whereas the corrected run has an absolute relative mass change of $1.2\times10^{-15}$. The main hump and capillary ripples remain closely comparable at $t\approx4500$. This benchmark therefore shows that agreement with travelling-wave calculations does not, by itself, guarantee that a reduced model respects the conservation laws inherited from the underlying transport problem. Although the drift is small in this case, the corrected equations conserve the reduced total inventory for every smooth periodic solution. We also distinguish this exact conservation property from higher-order compatibility with the full curved-interface boundary conditions.
Molecular mixing models for the conditional scalar diffusion (CSD) term in the filtered mass density function (FMDF) transport equation used in Large Eddy Simulation (LES-FMDF) approach are evaluated against direct numerical simulation (DNS) of spatially evolving jets. A database of DNS of turbulent mixing slot jets of $C_{7}H_{16} / O_2$ and $C_{7}H_{16} / Air$ is developed. The formulation includes the compressible form of the governing equations, a generalized multicomponent diffusion model, a cubic real gas equation of state and realistic property models. Simulations are conducted over a wide range of initial pressures ($1 atm < P_0 < 100 atm$) and jet width based Reynolds numbers of 850 and 1300. FMDF of mixture fraction at various filter widths are obtained from the simulation at several spatial locations within the flow. The CSD term in the exact transport equation of FMDF is calculated from the DNS. An \textit{a priori} analysis of Interaction by Exchange of Mean (IEM), Modified Curl (MC) and Mapping Closure (MAPPING) mixing models for CSD is conducted. At high pressures, significant deviations are observed between the CSD measured from the DNS and those predicted by the models. The source of the deviations is determined to be the mixing frequency used in the models. The significance of the mixing frequency and its variation with change in pressure and diffusion models is studied. The primary objective of this work is to assess the applicability of aforementioned mixing models to flows at large pressures and determine the causes for deviation in prediction.
The usage of chemical dispersant to promote the breakup of oil droplets remains one of the crucial response strategy towards controlling the rapidly rising droplets during an underwater oil spillage disaster like the Deepwater horizon. However, the detrimental effects of dispersant on the marine life are now well recognized. Here, we demonstrate that vortex ring presents an environmental friendly alternative to the debatable usage of chemical dispersant that primarily aids droplet break-up. We propose two mitigating strategies through the application of vortex rings: (i) atomization of large oil droplets into smaller droplets, which enhances natural dispersion and biodegradation (ii) targeted relocation of spilled oil droplets from ecologically sensitive regions to safer remote locations during minor spill events. This proof-of-concept study demonstrates that vortex rings can efficiently atomize underwater oil droplets while simultaneously transporting them over long distances. The efficiency strongly depends on the droplet size distribution and the strength of the vortex ring. Using time scale analysis and force balances arguments, we propose droplet response and atomization criteria by a traversing vortex ring. Beyond atomization and droplet transport which are center to many engineering processes, the work finds high relevance in the area of bubble turbulence, agricultural industry, multiphase reactors, combustion process, atmospheric fluid dynamics and oil extraction from rocks, to name a few.
Droplet collision governs droplet population dynamics in many chemical engineering processes, such as spray drying, spray cooling, agricultural spraying, and combustion. Existing analytical models impose deterministic, pairwise boundaries between collision outcomes, whereas machine-learning classifiers lack the explicit functional form required of analytical collision submodels. In this study, we develop a probabilistic symbolic-distillation model using nearly forty thousand experimental events spanning eight regimes and five dimensionless parameters, including over five thousand data for ambient pressure up to 50 atm. A machine-learning teacher learns the joint outcome-probability landscape from these data, and symbolic regression subsequently distils it into eight class-specific expressions that jointly define a coupled analytical model. The resulting analytical field replaces abrupt regime switching with finite-width fuzzy boundaries. It outperforms the evaluated conventional analytical boundary models and reveals that their main limitation is the inability of zero-width boundaries to represent gradual probability transitions. The "biased-dice" sampling scheme provides a statistically consistent and practically convenient model implementation for Eulerian-Lagrangian spray simulation.
We investigate the motion of a spherical microswimmer in a planar Poiseuille flow using analytical calculations and numerical simulations. Our results show that the microswimmer's orientation undergoes periodic oscillations governed by pendulum-like dynamics. The oscillation period is proportional to the channel width and to the complete elliptic integral of the first kind evaluated at the maximum orientation angle, and inversely proportional to the square root of the product of the maximum flow speed and the self-propulsion speed. Based on the net displacement along the flow direction within one oscillation period and the signs of the maximum and minimum velocities, we identify five motion states: upstream, oscillatory upstream, zero-drift oscillatory, oscillatory downstream, and downstream motion. The net displacement is jointly determined by the ratio of the self-propulsion speed to the maximum flow speed and the complete elliptic integrals of the first and second kinds. The signs of the maximum and minimum velocities are determined by the speed ratio and the cosine of the orientation angle. Flow nonuniformity promotes upstream migration, whereas inertial lift produces negligible radial displacement during one oscillation period. These findings identify the key parameters governing the motion of spherical microswimmers and provide significant implications for understanding microorganism motility and designing microrobots with prescribed motion states.
Pressure-driven radial flows of viscoelastic fluids are common in various industrial applications, such as injection molding and extrusion. It is well known that viscoelastic rheology can significantly impact the hydrodynamic features of non-Newtonian flows. However, these hydrodynamic features remain not fully understood compared to those observed in Newtonian flows. We analyze the pressure-driven radial flow of an Oldroyd-B fluid between parallel plates and present a theoretical framework together with finite-element numerical simulations for determining the flow rate-pressure drop relation. Unlike previous theoretical studies restricted to the weakly viscoelastic limit of low Deborah ($De$) numbers, we apply lubrication theory and consider the ultra-dilute limit, which allows us to study viscoelastic radial flows at order-one Deborah numbers. Using the one-way coupling between the Newtonian velocity profile and elastic stresses, we derive closed-form expressions for the conformation tensor and pressure drop in the ultra-dilute limit. We show that the pressure drop monotonically increases with $De$, identify the physical mechanisms governing this increase, and delineate the range of validity of the ultra-dilute approximation. We further reveal that the low-$De$ asymptotic analysis of the pressure drop may fail to accurately capture finite-element simulation results, even at low Deborah numbers, owing to its inability to satisfy the inlet elastic stresses. In contrast, our theoretical predictions based on the ultra-dilute limit are in excellent agreement with the finite-element simulation results, enabling us to elucidate the pressure drop behavior for order-one Deborah numbers.
Porous and perforated leading edges reduce turbulence-interaction noise, but predictions of the reduction have relied on material impedances adjusted to the acoustic data. This paper removes that adjustment. The finite-chord scattering problem is solved in compressible subsonic flow for an arbitrary chordwise admittance, with the shed wake, the Kutta condition and permeability-dependent edge singularities, and closed to the far field with the finite span and the array aperture represented. The admittance of a perforated insert is derived from the Rayleigh conductivity of one aperture with three computed corrections: plate thickness, interaction between neighbouring apertures, and Howe's grazing-flow blockage evaluated at the aperture-averaged velocity of the boundary layer. A bulk porous insert is represented by its measured permeability and pore-fluid inertia. For three perforates at four flow conditions the prediction agrees with measurements to 2.5 dB rms over 84 comparisons, against 4.4 dB with the unscaled impedance used previously. Although the Mach number is 0.05, the chord is acoustically non-compact, and an incompressible solution rises where the measurement falls. Two of three bulk porous inserts are predicted to within 1.7 dB rms where interaction noise dominates. Narrow, equally spaced peaks measured on perforated plates are reproduced by none of the admittances tested; their spacing implies a source convecting at 0.62 of the free-stream speed, absent from a frozen-gust model. A parameter study shows that the perforate geometry acts almost entirely through one low-frequency admittance group, in which the hole radius cancels for holes buried in the laminar boundary layer.
Transporting droplets along cylindrical surfaces is an emerging area of research with a wide range of practical applications such as fog collection and filtration. Recent experimental studies by Feng et al. (2024) have revealed a wealth of new droplet coalescence dynamics on a pre-wetted cylindrical fibre, illustrating the need for more advanced theory. In this paper, we study both the early-stage and late-stage dynamics of droplet coalescence on a sloped cylindrical fibre. In the early stage, when the minimum thickness of the liquid bridge connecting two droplets is small, we propose a lubrication-based model to analyze the capillary-driven self-similar dynamics of the liquid bridge. Our theory predicts the bridge growth rate and demonstrates that asymmetry in the adjoining contact angles contributes to the migration of the merged droplet, showing good agreement with experimental observations. In the late stage, as droplet merging progresses, inertia effects become significant and lead to rapid mass transfer. To capture this regime, we develop a weighted residual model for liquid films flowing on a sloped pre-wetted cylinder, incorporating key physical effects such as gravity, low-to-moderate inertia, surface tension, and intermolecular forces. This model explains the experimentally observed directional self-propelled transport during droplet coalescence, where the larger droplet migrates toward the smaller one. Numerical simulations based on this model further support these findings, showing good agreement between the predicted coalescence-induced droplet migration and experimental observations.
We formulate a scale dependent tensorial extension of a fractional dissipation operator motivated by a macroscopic Navier Stokes correspondence. The starting point is a fractional model in which the order of dissipation decreases continuously from the Laplacian value toward the Kolmogorov value as a local Reynolds number increases. At the critical value of the local Reynolds number, the model reduces to an isotropic operator. We then ask whether the same macroscopic description can be continued into a regime in which velocity fields have Onsager, critical or subcritical regularity without introducing an unrelated closure or an additional isotropic dissipation scale. In addition, we extend the original AFNS free, energy, construction by one nonnegative concentration order parameter c. The new term is quadratic in c and linearly coupled to the defect beta function, so that minimization leaves the original AFNS free energy exactly unchanged when the concentration branch is inactive. The resulting three regime mosaic is regular, inertial, and concentrated/anomalous. The construction uses the coarse-grained stress and its critical normalization. The normalized anisotropy tensor provides the lowest order objective tensorial information available at the critical scale, while the scale derivative of the Duchon Robert energy defect, supplies a scalar measure of departure from the Onsager threshold. We propose the minimal constitutive deformation of the critical operator: in Weyl quantization, the isotropic fractional symbol is dressed by this exponential tensorial factor, giving a pseudodifferential operator that reduces to the scalar critical operator. The exponential parametrization guarantees positive eigenvalues and unit determinant for symmetric traceless, so the deformation redistributes dissipation among directions without introducing a new isotropic amplitude.
The wave-particle decomposition of kinetic relaxation equations is extended to the full Boltzmann equation. An identity transformation writes the collision term as a relaxation toward the local Maxwellian plus a conservative remainder. Along characteristics, the equilibrium integral over a local kinetic horizon defines the wave, and its exact complement, the particle, carries the remainder. The moments of the wave and particle equations are extended Navier-Stokes equations, and the total conservation law and the particle equation form the wave-particle multiscale equations, a closed system, exact for every horizon, that spans a continuous spectrum from the Boltzmann equation to Navier-Stokes hydrodynamics, parametrized by the horizon-to-relaxation ratio. A coupled implicit wave-particle iteration solves this system, updating the conservative variables with the particle fixed and the particle with the predicted wave and an endpoint traction from the macroscopic prediction. Under discrete compatibility assumptions the iteration preserves the discrete Boltzmann solution as its fixed point, and its near-continuum convergence factor is bounded by a particle transport share decaying exponentially with this ratio. Tests on a normal shock, Couette and cavity flows, and hypersonic cylinder and Apollo flows agree with the conventional iterative scheme in rarefied and transition regimes and with Navier-Stokes solutions near the continuum limit. There the wave-particle iteration needs about three to more than four orders of magnitude fewer iterations than this scheme or, where it did not converge, its Shakhov version, and the Boltzmann solution costs one to four times the wall time of a Navier-Stokes solver on the same mesh.
Time-variable stellar heating may create new climate states, such as permanently hot or cold hemispheres, on asynchronous planets. Planets on eccentric orbits, circumbinary planets and planets on polar orbits around gravity-darkened stars all receive variable irradiation, yet variably irradiated planets (VIPs) have not been studied from a fundamental dynamical perspective. We present the first generalized study of the atmospheric response of variably irradiated tidally locked planets, using a 1.5-layer shallow-water model forced towards a time-dependent radiative equilibrium. We derive linearized analytic expressions for the mean, amplitude and phase lag of the response, which is that of a forced damped harmonic oscillator, and recover the established trends of tidally locked planets in the limit of steady forcing. When the instellation variation and wave timescales match and the damping is weak, a resonance occurs because waves propagate and interfere freely. The resonance produces a global standing wave mode and consequently a temporarily warmer nightside. Linear and nonlinear simulations reproduce these results and show that the hotspot oscillates longitudinally at the forcing frequency and reaches the antistellar point at resonance. The oscillation is eastward or westward of the substellar longitude depending on the rotation rate and the nonlinearity of the system. Synthetic phase curves carry signatures of this variability and invert in shape at resonance because of nightside warming. Among observed systems, we find 11 planets with at least 10% phase curve modulation from variable irradiation. Eccentric, pseudo-synchronous giants are the most promising targets, with HAT-P-2b having the largest predicted modulation.
We present Suêtes, a fully differentiable, non-hydrostatic, limited-area atmospheric dynamical core implemented in JAX. Suêtes includes both a semi-implicit semi-Lagrangian scheme and a Split-Explicit Runge-Kutta scheme with acoustic substepping, enabling applications from coarse regional simulations to convection-permitting dynamics. Its terrain-following geometry, spatial discretization, time integration, initial conditions, lateral boundary treatment, and included physical parameterizations are all compatible with reverse-mode automatic differentiation. A single reverse-mode pass can therefore compute sensitivities of forecast diagnostics to selected initial and boundary fields, physical parameters, topography, and computational geometry. We provide a series of numerical experiments---including a tracer-source inverse problem, a sensitivity analysis of a three-dimensional squall-line, adversarial perturbation construction an ERA5-driven downslope hurricane force wind event, topography optimization, and terrain-following coordinate optimization---demonstrating the versatility of Suêtes as a framework for gradient-based diagnosis and model development at regional and convection-permitting scale.
Odd pressure is a parity-odd transport coefficient that gives an isotropic pressure response to vorticity in compressible two-dimensional fluids. By asymmetrically coupling vortical and compressional modes, it provides a direct mechanism for a nonreciprocal hydrodynamic response reminiscent of several biological systems. Yet microscopic models realizing odd pressure as a transport coefficient are relatively scarce. Here, we construct one such model using a lattice-gas cellular automaton based on the Frisch-Hasslacher-Pomeau II model with rest particles and chiral collisions, whose Chapman-Enskog coarse-graining produces a nonzero odd pressure coefficient. The key ingredient is a local parity-breaking collision that weakly rotates moving particles in the presence of a rest particle. At hydrodynamic scales, this microscopic rule acts as an effective magnetic field, converting the model's bulk viscosity into odd pressure. More generally, we derive expressions relating Hall viscosity, odd pressure, and odd torque to their parity-even counterparts in our model.
Divergence-free vector fields are fundamental state variables in incompressible flows and many PDE systems. Redundant parameterizations, including Neural Conservation Law (NCL) potentials, map multiple auxiliary representations to the same physical field. Our experiments show that this redundancy can reduce static representation-fitting error by enlarging the set of equivalent solutions, but the resulting many-to-one mapping does not provide a unique state for operator learning. We introduce a minimal representation that encodes a real \(D\)-component divergence-free vector field on a \(D\)-dimensional domain as a real \((D-1)\)-component field on the same domain. Exploiting the transverse structure imposed by incompressibility in Fourier space, we use a Householder orthogonal transformation to construct the reduced coordinates directly. For periodic and closed impermeable fields, the transform is invertible, isometric, and angle-preserving. For open nonperiodic flows, Fourier extension constructs a compatible periodic field, and a minimum-energy rule selects a unique reduced representation. Neural operators then learn temporal evolution entirely in this reduced space. At inference, the predicted \((D-1)\)-component field is decoded directly into a physical divergence-free \(D\)-component field, without predicting an ambient field or applying post-hoc projection. Experiments on static fitting and temporal prediction reveal a task-dependent trade-off: redundancy facilitates static optimization, whereas unique invertible coordinates provide a well-defined state for temporal dynamics. By removing unconstrained longitudinal or null directions from the learned state space, the proposed formulation achieves lower prediction error and greater robustness while preserving divergence freedom by construction.