2026-03-02 | | Total: 7
Reconstructing fluid flows from sparse sensor measurements is a fundamental challenge in science and engineering. Widely separated measurements and complex, multiscale dynamics make accurate recovery of fine-scale structures difficult. In addition, existing methods face a persistent tradeoff: high-accuracy models are often computationally expensive, whereas faster approaches typically compromise fidelity. In this work, we introduce BLISSNet, a model that strikes a strong balance between reconstruction accuracy and computational efficiency for both flow reconstruction and nudging-based data assimilation. The model follows a DeepONet-like architecture, enabling zero-shot inference on domains of arbitrary size. After the first model call on a given domain, certain network components can be precomputed, leading to low inference cost for subsequent evaluations on large domains. Consequently, the model can achieve faster inference than classical interpolation methods such as radial basis function or bicubic interpolation. This combination of high accuracy, low cost, and zero-shot generalization makes BLISSNet well-suited for large-scale real-time flow reconstruction and data assimilation tasks.
The impulse response of a flame to acoustic velocity perturbations is a key quantity for predicting thermoacoustic stability, but its identification from sparse, noisy observations requires solving an ill-posed inverse convolution problem. This is typically achieved with system identification methods, which require hand-tuning of regularization, model order, and sampling parameters, and provide no principled mechanism for incorporating prior physical knowledge. In this paper, we reformulate the identification problem within a Bayesian framework. The impulse response is represented as a physically motivated distributed time delay model, whose parameters correspond to convective delays and dispersive broadening. For a given number of pulses, the model parameters are inferred from the data using Bayesian parameter inference. The number of pulses is then selected using Bayesian model comparison, which balances data fit against model complexity to identify the simplest model capable of explaining the data. The framework is demonstrated on broadband-forced large eddy simulation data from a turbulent swirl-stabilized burner. Bayesian model comparison selects a three-Gaussian impulse response for this flame, consistent with physical interpretations in previous work. Compared with system identification, the Bayesian approach produces impulse responses with fewer spurious features and enables straightforward enforcement of a known low-frequency gain. Finally, we show that the Bayesian approach is robust to significant reductions in recording length, making it appealing for impulse response identification from costly simulations, where there is an incentive to minimize computational cost.
We experimentally investigate the Lagrangian dynamics of finite-sized, neutrally buoyant droplets in homogeneous isotropic turbulence. The droplet size follows a log-normal distribution whose average value decreases with increasing Reynolds number, reflecting enhanced turbulent breakup. While size-conditioned velocity and acceleration statistics show only weak finite-size dependence, temporal measures reveal clear size-dependent dynamics: larger droplets exhibit longer Lagrangian velocity integral times and an extended ballistic regime in their mean squared displacement. These findings indicate that though droplets exhibit mild deformation and internal circulation, they behave similarly to finite-size rigid particles in terms of Lagrangian dynamics. Our study opens the way to study droplet-laden turbulence and droplet-flow interactions.
The transport of the Antarctic Circumpolar Current (ACC) has been shown to increase with friction. Previous studies explained this counter-intuitive relationship called frictional control based on the eddy geometric parametrizations. They focused on the eddy momentum transfer and eddy energetics. To maintain the balance between wind stress and eddy interfacial form stress, eddy energy must remain unchanged as friction increases; this requires enhanced baroclinicity to compensate for stronger eddy energy dissipation. However, the independence of eddy energy has not been fully verified, and this interpretation assumes negligible barotropic energy conversion. To address these gaps, we conduct sensitivity experiments in an idealized stratified reentrant channel with varying bottom drag. Numerical simulations show that eddy energy changes substantially with friction. Furthermore, in the high-drag regime, baroclinic energy conversion dominates eddy energy generation, whereas in the low-drag regime barotropic energy conversion contributes substantially. Despite these differences, baroclinicity increases with eddy energy dissipation across all regimes, although the relationship is weak in the low-drag regime owing to barotropic energy conversion. To explain this phenomenon, we extend the frictional control framework independent of any specific eddy parameterization. A simple scaling argument leads to a generalized frictional control, s~D(E)/τ_w, where s is baroclinicity, D(E) is eddy energy dissipation, and τ_w is wind stress. This framework provides a natural extension of the existing framework and successfully explains the numerical results. These results indicate that eddy dissipation controls the baroclinicity; therefore, properly parameterizing the eddy dissipation rate is essential for representing ACC dynamics in ocean models.
Clustering techniques offer a powerful framework for analyzing complex flow dynamics and reducing computational costs in large-scale simulations. In this work, we propose a novel clustering-based approach using Vector Quantization Principal Component Analysis (VQPCA) to identify structural sensitivity zones, namely the regions where the fluid flow is more receptive to changes. To the authors knowledge, this is the first application of VQPCA to a fluid dynamics problem for the identification of flow patterns and dynamically relevant regions. As a fully data-driven technique, it does not rely on adjoint methods; therefore, this approach has the advantage of having low computational cost, since it depends exclusively on data from the direct problem. The VQPCA technique demonstrates its ability to extract dominant flow features by clustering the flow field into regions characterized by their intrinsic dynamics. To assess the validity of this method, it is used to investigate the wake behind a circular cylinder, revealing similarities to previously established structural sensitivity regions. The robustness of the approach is further assessed through validation and calibration in different operating conditions in this flow scenario. As an extension of the analysis, we address the complex dynamics of two planar synthetic jets, where the clustering insights can lead to develop flow control strategies. These results highlight the potential of clustering-based methods as practical and effective tools to analyze and optimize fluid flows.
Based on everyday experience fluid flows tend to be ordered and quiescent if inertial forces are low and held in check by viscosity. This intuition spectacularly fails in the case of complex macromolecular fluids like polymer melts, paints and biofluids. In such cases elastic fluid properties can drive turbulent motions at moderate and even vanishing Reynolds numbers. By studying viscoelastic flows in curved pipes we demonstrate that this low inertia phenomenology results from the competition of two hydrodynamic instabilities and respectively from the co-existence and interdependence of two distinct turbulent states. Unexpectedly the established categories of elastic and elasto-inertial turbulence (ET and EIT) fail to demarcate the actual turbulent states, fundamentally changing the perception of this phenomenon a century after its discovery.
Data assimilation (DA) for compressible flows with shocks is challenging because many classical DA methods generate spurious oscillations and nonphysical features near uncertain shocks. We focus here on the ensemble Kalman filter (EnKF). We show that the poor performance of the standard EnKF may be attributed to the bimodal forecast distribution that can arise in the vicinity of an uncertain shock location; this violates the assumptions underpinning the EnKF, which assume a forecast which is close to Gaussian. To address this issue we introduce the new neural EnKF. The basic idea is to systematically embed neural function approximations within ensemble DA by mapping the forecast ensemble of shocked flows to the parameter space (weights and biases) of a deep neural network (NN) and to subsequently perform DA in that space. The nonlinear mapping encodes sharp and smooth flow features in an ensemble of NN parameters. Neural EnKF updates are therefore well-behaved only if the NN parameters vary smoothly within the neural representation of the forecast ensemble. We show that such a smooth variation of network parameters can be enforced via physics-informed transfer learning, and demonstrate that in so-doing the neural EnKF avoids the spurious oscillations and nonphysical features that plague the standard EnKF. The applicability of the neural EnKF is demonstrated through a series of systematic numerical experiments with an inviscid Burgers' equation, Sod's shock tube, and a two-dimensional blast wave.