2025-12-05 | | Total: 11
The integration of artificial intelligence into experimental fluid mechanics promises to accelerate discovery, yet most AI applications remain narrowly focused on numerical studies. This work proposes an AI Fluid Scientist framework that autonomously executes the complete experimental workflow: hypothesis generation, experimental design, robotic execution, data analysis, and manuscript preparation. We validate this through investigation of vortex-induced vibration (VIV) and wake-induced vibration (WIV) in tandem cylinders. Our work has four key contributions: (1) A computer-controlled circulating water tunnel (CWT) with programmatic control of flow velocity, cylinder position, and forcing parameters (vibration frequency and amplitude) with data acquisition (displacement, force, and torque). (2) Automated experiments reproduce literature benchmarks (Khalak and Williamson [1999] and Assi et al. [2013, 2010]) with frequency lock-in within 4% and matching critical spacing trends. (3) The framework with Human-in-the-Loop (HIL) discovers more WIV amplitude response phenomena, and uses a neural network to fit physical laws from data, which is 31% higher than that of polynomial fitting. (4) The framework with multi-agent with virtual-real interaction system executes hundreds of experiments end-to-end, which automatically completes the entire process of scientific research from hypothesis generation, experimental design, experimental execution, data analysis, and manuscript preparation. It greatly liberates human researchers and improves study efficiency, providing new paradigm for the development and research of experimental fluid mechanics.
Recent advances in velocity and temperature transformations have enabled recovery of the law of the wall in compressible wall-bounded turbulent flows. Building on this foundation, a flux-controlled wall model (FCWM) for Large Eddy Simulation (LES) is proposed. Unlike conventional wall-stress models that solve the turbulent boundary layer equations, FCWM formulates the near-wall modeling as a control problem applied directly to the outer LES solution. It consists of three components: (1) the compressible law of the wall, (2) a feedback flux-control strategy, and (3) a shifted boundary condition. The model adjusts the wall shear stress and heat flux based on discrepancies between the computed and target transformed velocity and temperature, respectively, at the matching location. The proposed wall model is evaluated using LES of turbulent channel flows across a broad range of conditions, including quasi-incompressible cases with bulk Mach number \(M_b = 0.1\) and friction Reynolds number \(Re_τ= 180 \sim 10{,}000\), and compressible cases with \(M_b = 0.74 \sim 4.0\) and bulk Reynolds number \(Re_b = 7667 \sim 34{,}000\). The wall-modelled LES reproduce mean velocity and temperature profiles in agreement with direct numerical simulation data. For all tested cases with \(M_b \leq 3\), the wall model achieves relative errors of \(|ε_{C_f}| < 4.1\%\), \(|ε_{B_q}| < 2.7\%\), and \(|ε_{T_c}| < 2.7\%\) in friction coefficient, non-dimensional heat flux, and centerline temperature, respectively. In the quasi-incompressible regime, the wall model achieves \(|ε_{C_f}| < 1\%\). Compared to the conventional equilibrium wall model, the proposed FCWM achieves higher accuracy in compressible turbulent channel flows without solving the boundary layer equations, thereby reducing computational cost.
The incompressible three-dimensional Euler equations develop very thin pancake-like regions of exponentially increasing vorticity. The characteristic thickness of such regions decreases exponentially with time, while the other two dimensions do not change considerably, making the flow near each pancake strongly anisotropic. The pancakes emerge in increasing number with time, which may enhance the anisotropy of the flow, especially if they orient similarly in space. In the present paper, we study numerically the local isotropy of the flow by analyzing the evolution of the so-called isotropy markers [Phys. Rev. Fluids 10, L022602 (2025)]. We show that these functions drift slowly towards unity, indicating the process of slow isotropization, which takes place without the viscous scales getting exited and despite the similar orientation of the emerging pancakes.
High-order Discontinuous Galerkin (DG) methods offer excellent accuracy for turbulent flow simulations, especially when implemented on GPU-oriented architectures that favor very high polynomial orders. On modern GPUs, high-order polynomial evaluations cost roughly the same as low-order ones, provided the DG degrees of freedom fit within device memory. However, even with high-order discretizations, simulations at high Reynolds numbers still require some level of under-resolution, leaving them sensitive to numerical dissipation and aliasing effects. Here, we investigate the interplay between intrinsic DG dissipation mechanisms (implicit dissipation) -- in particular split forms and Riemann solvers -- and explicit subgrid-scale models in Large Eddy Simulations (LES). Using the three-dimensional Taylor--Green vortex at $Re = 1600$ and an inviscid case ($Re \to \infty$), we evaluate kinetic energy dissipation, spectral accuracy, and numerical stability. Our results show that when stability for under-resolved turbulence is ensured through split-forms (energy- or entropy-stable) schemes, subgrid-scale (SGS) LES models are not strictly necessary. At moderate Reynolds numbers, when the spatial resolution is sufficient to capture the relevant turbulence scales (i.e., in well-resolved LES), adding SGS models does not improve accuracy because the wavenumber range where they act overlaps with the inherent numerical dissipation of the DG scheme. In contrast, when the resolution is insufficient, as is typical at high Reynolds numbers, explicit subgrid-scale models complement the numerical dissipation and enhance accuracy by removing the excess energy that numerical fluxes alone cannot dissipate. These findings provide practical guidance for choosing numerical strategies in high-order turbulence simulations.
Highly resolved simulations reveal the fundamental influence of a carrier fluid's flow dynamics on triboelectric powder charging. We found that particles transported through a square-shaped duct charge faster than in a channel flow caused by secondary flows that led to more severe particle-wall collisions. Specifically, particles with a Stokes number of 4.69 achieve 85 % of their equilibrium charge approximately 1.5 times faster in duct flow than in channel flow. Also, charge distribution is more uniform in a duct cross-section compared with a channel cross-section. In channel flow, particles are trapped near the walls and collide frequently due to limited movement in the wall-normal direction, causing localized charge buildup. In contrast, duct flow promotes better mixing through secondary flows, reducing repeating collisions and providing uniform charge distribution across the cross-section. Upon charging, electrostatic forces significantly reshape particle behaviour and distribution. Once the powder achieves half of its equilibrium charge, particles increasingly accumulate at the wall, leading to a reduced concentration in the central region. These changes in particle distribution have a noticeable impact on the surrounding fluid phase and alter the overall flow dynamics. These findings open the possibility for a new measure to control powder charging by imposing a specific pattern.
Fast, geometry-generalizing surrogates for unsteady flow remain challenging. We present a time-dependent, geometry-aware Deep Operator Network that predicts velocity fields for moderate-Re flows around parametric and non-parametric shapes. The model encodes geometry via a signed distance field (SDF) trunk and flow history via a CNN branch, trained on 841 high-fidelity simulations. On held-out shapes, it attains $\sim 5\%$ relative L2 single-step error and up to 1000X speedups over CFD. We provide physics-centric rollout diagnostics, including phase error at probes and divergence norms, to quantify long-horizon fidelity. These reveal accurate near-term transients but error accumulation in fine-scale wakes, most pronounced for sharp-cornered geometries. We analyze failure modes and outline practical mitigations. Code, splits, and scripts are openly released at: https://github.com/baskargroup/TimeDependent-DeepONet to support reproducibility and benchmarking.
Interactions between internal solitary waves and floating canopies of varying length and porosity are examined via laboratory experiments and complementary simulations for a miscible, two-layer system. In both approaches, internal solitary waves of varying amplitudes are generated by a jet-array mechanism that is driven by the nonlinear eKdV solution. Pycnocline displacements, phase speeds, and velocity fields are obtained using synchronized planar laser-induced fluorescence and particle imaging velocimetry systems in the experiment. In the simulations, the canopy is represented as a porous zone with prescribed porosity and hydraulic conductivity determined by the Kozeny-Carman model, which is validated by comparing simulated and measured horizontal velocity profiles. The higher-porosity (transitional) canopy produces a nearly monotonic, albeit minor, amplitude reduction and negligible wave energy dissipation after the interaction. However, the shear layer developed at the bottom edge of the lower-porosity (dense) canopy grows to a comparable strength as the shear sustained by the internal solitary wave profile at the pycnocline. The vortex pair generated by this shear accelerates the upper-layer fluid beneath the canopy, leading to complex nonlinear amplitude modulation and significant wave transformation. With an extended canopy length, the internal solitary waves settle to a quasi-steady state with a significant phase speed reduction. Upon the wave exiting the canopy, flow separation at the downstream edge of the canopy again pairs with the shear at the pycnocline, inducing an intensified jet. This complex interaction leads to energy transfer between kinetic and potential energy under the dense canopy.
This work develops a computational and theoretical framework for determining effective permeabilities in anisotropic microscopic geometries containing dense, fibre-like obstacles, motivated by the need to model flow in coiled aneurysm domains accurately. Building on homogenisation theory and fully resolved simulations in Representative Elementary Volumes (REVs), we validate the permeability model introduced in [C. Boutin, Study of permeability by periodic and self-consistent homogenisation. Eur. J. Mech. A Solids, 19(4):603-632, 2000] and propose a systematic methodology for capturing the directional variations induced by fibre orientation. The resulting permeability tensors are incorporated into macroscopic flow simulations based on the Darcy equation, enabling direct comparison of anisotropic and isotropic permeability models across several benchmark configurations. Our findings show that anisotropy has a significant impact on local flow direction and magnitude, generating directional permeability contrasts which cannot be reproduced by classical isotropic approximations. By integrating coil-induced microstructural effects into continuum-scale hemodynamic models, the proposed approach enables more realistic assessment of post-treatment aneurysm flow behaviour. Beyond this clinical application, the framework is broadly applicable to other biomedical and engineering systems involving fibrous or filamentous porous microstructures.
Turbulence of weakly interacting waves displays a great deal of universality: independence of the details of the interaction and of the pumping and dissipation scales. Here we study how inverse turbulent cascades (from small to large scales) transition from weak to strong. We find that while one-loop corrections can be dependent on excitation and dissipation scales, new types of universality appear in strong turbulence. We contrast turbulence of spin waves in ferromagnets with turbulent cascades in the Nonlinear Schrödinger Equation (NSE) and in an MMT-like model in higher dimensions having a multiplicative interaction vertex: vertex renormalization gives rise to dependence on the pumping (UV scale) in the former but not in the latter. As a result of this spectral nonlocality, spin-wave turbulence stops being weak if one is sufficiently far from the pumping scale, even when the interaction of waves with comparable wavenumbers is weak. We paraphrase this as: nonlocality enhances nonlinearity. We then describe strong turbulence in a multi-component version of these models with a large number of components. We argue that strong spin-wave turbulence is similar to turbulence of the focusing NSE, as it realizes a critical-balance state. However, UV nonlocality causes the level of spin-wave turbulence at large scales to decrease with increasing pumping level, culminating in a state that is independent of the level of pumping.
This work introduces RedEigCD, the first self-adaptive timestepping technique specifically tailored for reduced-order models (ROMs) of the incompressible Navier-Stokes equations. Building upon linear stability concepts, the method adapts the timestep by directly bounding the stability function of the employed time integration scheme using exact spectral information of matrices related to the reduced operators. Unlike traditional error-based adaptive methods, RedEigCD relies on the eigenbounds of the convective and diffusive ROM operators, whose computation is feasible at reduced scale and fully preserves the online efficiency of the ROM. A central theoretical contribution of this work is the proof, based on the combined theorems of Bendixson and Rao, that, under linearized assumptions, the maximum stable timestep for projection-based ROMs is shown to be larger than or equal to that of their corresponding full-order models (FOMs). Numerical experiments for both periodic and non-homogeneous boundary conditions demonstrate that RedEigCD yields stable timestep increases up to a factor 40 compared to the FOM, without compromising accuracy. The methodology thus establishes a new link between linear stability theory and reduced-order modeling, offering a systematic path towards efficient, self-regulating ROM integration in incompressible flow simulations.
NORi is a machine-learned (ML) parameterization of ocean boundary layer turbulence that is physics-based and augmented with neural networks. NORi stands for neural ordinary differential equations (NODEs) Richardson number (Ri) closure. The physical parameterization is controlled by a Richardson number-dependent diffusivity and viscosity. The NODEs are trained to capture the entrainment through the base of the boundary layer, which cannot be represented with a local diffusive closure. The parameterization is trained using large-eddy simulations in an "a posteriori" fashion, where parameters are calibrated with a loss function that explicitly depends on the actual time-integrated variables of interest rather than the instantaneous subgrid fluxes, which are inherently noisy. NORi is designed for the realistic nonlinear equation of state of seawater and demonstrates excellent prediction and generalization capabilities in capturing entrainment dynamics under different convective strengths, oceanic background stratifications, rotation strengths, and surface wind forcings. NORi is numerically stable for at least 100 years of integration time in large-scale simulations, despite only being trained on 2-day horizons, and can be run with time steps as long as one hour. The highly expressive neural networks, combined with a physically-rigorous base closure, prove to be a robust paradigm for designing parameterizations for climate models where data requirements are drastically reduced, inference performance can be directly targeted and optimized, and numerical stability is implicitly encouraged during training.