2026-05-29 | | Total: 10
Many experimental studies have reported variations in interfacial tension. Isolating all the geometric and fluid material parameters and varying the interfacial tension can be useful to check their influence. Numerical investigations using Free Surface 3D (FS3D), have been conducted to compare varying values of interfacial tension and evaluate the sensitivity. A grid independence study compared the compound crown height of a splash to determine the required resolution for validation. A qualitative validation showed FS3D could correctly capture the impact morphology while varying the viscosity ratio of the drop and film liquid when compared to the experimental results. A quantitative validation for a water drop impacting onto an oil film shows a good match for the crown heights of the numerical and experimental data. The same setup was then extended to study the variation of interfacial tension, where the deviation of the overall compound crown height and spreading diameter of the internal crowns was compared. Results revealed minor changes in the compound crown height and spreading diameter of the drop liquid, but the internal crown composition showed significant differences. In order to run FS3D efficiently on the new supercomputer Hunter, which has a new APU architecture-based system, extensive work had to be done. To adapt to the new hardware architecture, large parts of FS3D have been ported to utilise the AMD Instinct MI300A accelerated processing units (APUs) at HLRS using OpenMP. Implementation of Umpire memory pools improved performance for larger workloads per APU. The GPU-accelerated code achieves a 4 times speedup compared to CPU-only execution on the same hardware. Strong and weak scaling tests have been conducted, showing good strong scaling for up to 4 APUs, and linear weak scaling for up to 512 APUs, resulting in a total of 4096**3 cells for the first time.
This paper develops a reduced-order framework for modelling the two-way coupling between gravity waves and turbulent wakes in large-scale wind farms. Linearising the non-hydrostatic Boussinesq equations and introducing simplifications appropriate to the boundary layer and the overlying stratified free atmosphere yield separate governing equations for the two regions. These are coupled through a dynamic boundary condition at the capping inversion, which directly captures the feedback of gravity waves on the boundary-layer flow. A mixed spectral-finite-difference discretisation yields a computationally efficient model while retaining vertical boundary-layer structure. Comparisons with large-eddy simulations (LES) confirm the model successfully reproduces both internal wind-farm flow and large-scale gravity-wave effects. It captures the upstream blockage induced by adverse pressure gradients, as well as the accelerated wake recovery within and downwind of the farm, driven by favourable pressure gradients.
The simplified lattice Boltzmann method (SLBM) is a recent development in the lattice Boltzmann method (LBM) community, addressing the intrinsic limitations of the traditional LBM by directly evolving macroscopic quantities and maintaining numerical stability in high Reynolds number simulations. However, fundamental understanding of the numerical dissipation and dispersion of SLBM is still lacking, and the origin of its good numerical stability remains unknown. In this work, a generalized formulation is developed, revealing that the SLBM recovers modified macroscopic equations containing both intrinsic physical deviations and numerical truncation errors. To remove these deviations, the macroscopic equation derived from the standard BGK-LBM is adopted as a reference model and solved by the predictor corrector strategy, which constitutes the reformulated SLBM. The proposed method uses the generalized SLBM formulation in the predictor step with tunable high-order parameters, while the corrector step is realized by the finite-difference discretization. Linear wave analysis clarifies the roles of these parameters in controlling numerical dissipation and dispersion, which is then validated in more complicated numerical examples. It is demonstrated that the reformulated SLBM preserves the second order accuracy, improves the dispersion and dissipation performance, enhances numerical stability, and resolves fine vortex structures on relatively coarse grids. Thus, the proposed method combines improved numerical properties with the simplicity of SLBM, offering a high fidelity and stable scheme for incompressible flow simulations.
Transport at small scales is classically understood within an equilibrium framework, where dispersion theory successfully describes shear-enhanced diffusion for passive particles in the continuum limit. However, as most bacteria can move on their own, their motility in flows, inherently out of thermal equilibrium, fundamentally challenges this framework. A minimal, predictive unified theory of bacterial transport in low-Reynolds-number flows remains lacking. Here, from first principles, we develop an analytical hydrodynamic model that enforces consistent no-flux boundary conditions and uses the method of images to characterize the flow-wall coupling. The model quantitatively reproduces measured bacterial distributions and reveals a hydrodynamic locking mechanism accompanied by mean-drift invariance -- an active counterpart to Taylor dispersion. We clarify that shear-induced depletion and alignment are dual manifestations of a single active phase-space topology, ruling out explanations based solely on the local shear magnitude. The theory is validated against microfluidic experiments spanning multiple bacterial species and shear geometries, from one-dimensional to fully three-dimensional flows. Our findings establish a unified phase-space framework for bacterial hydrodynamics, advancing the fundamental understanding of active matter.
This paper is associated with a poster winner of a 2025 American Physical Society's Division of Fluid Dynamics (DFD) Gallery of Fluid Motion Award for work presented at the DFD Gallery of Fluid Motion. The original poster is available online at the Gallery of Fluid Motion, https://gfm.aps.org/meetings/dfd-2025/692e39e5a7f805227b16fe9e
Closure-level accuracy in neural kinetic shock solvers is not guaranteed by accurate density, velocity and temperature profiles, because the relevant observables are velocity-weighted projections of the nonequilibrium distribution. We study this observability problem for one-dimensional Bhatnagar--Gross--Krook (BGK) shock waves using a positive macro--micro physics-informed neural network (PINN) in which the distribution is represented as a local Maxwellian multiplied by a bounded exponential correction. Independent discrete-velocity method (DVM) references are used for validation. Shock-tube tests show that sparse joint anchoring of heat flux and normal stress stabilises the primary nonequilibrium layer, whereas residual-only, macro-only and single-moment variants fail in distinct ways. In a stationary Mach-2 normal shock, a flux-locked compact model recovers $ρ$, $u_x$, $T$, $q_x$, $σ_{xx}$ and $m_{xxx}^{cl}$, but leaves $R_{xx}^{cl}$ with order-unity error. DVM diagnostics show that $R_{xx}^{cl}$ is controlled by a sign-changing, tail-weighted cancellation weakly observed by lower moments. A shock-local closure correction aligned with this missing projection reduces the relative $R_{xx}^{cl}$ error to $1.12\times10^{-1}$ while preserving the lower moments. A common-initialisation ablation shows that optional distribution-function probe losses are diagnostic rather than constitutive. A supplementary DVM--PINN comparison for the scalar fourth-order excess $Δ$ shows that the obstruction is anisotropic, sign-changing tail weighting rather than fourth-order polynomial degree alone.
Wave steepness is a key geometric variable for describing breaking occurrence and its consequences, including energy dissipation and air entrainment. Using three laboratory campaigns under varying spectral conditions and co-flowing wind forcing, we contrast two types of steepness commonly used for unsteady breaking waves: spectrally-informed wave-group steepness (prognostic), obtained from fixed-point surface-elevation records, and locally-measured crest steepness (diagnostic), obtained from spatial surface profiles extracted using the SDBW-I image-processing method developed herein. For the former, the long-adopted $\mathcal{S}_n$ (linear sum of Fourier-component steepness) increases appreciably within about two dominant wavelengths upstream of breaking because of its sensitivity to evolving high-frequency content. When measured sufficiently far upstream, however, wave-group steepness remains approximately linearly related to the local zero-crossing steepness $\mathcal{S}_b$ across bulk unforced conditions. Notwithstanding this, we argue that the crest-front steepness, $\mathcal{S}_{\mathrm{front}}(t_b)$, which delineates the front-face slope at incipient breaking, is the most physically meaningful metric examined here. It exhibits a consistent breaking-onset lower-bound threshold of $\mathcal{S}_{\mathrm{front}}(t_b)\approx0.2$, while values above this threshold decrease with wind speed as crests become less forward leaning. This may be attributed to wind-modified dispersion, enhanced high-frequency spectral content and aerodynamic sheltering, suggesting that wind--wave and wave--wave interactions act as competing mechanisms in triggering breaking through kinematic and energetic processes beyond what geometry alone can explain. Even so, $\mathcal{S}_{\mathrm{front}}(t_b)$ has strong potential as a controlling variable for future studies of breaking energetics and crest-scale dynamics.
Real-time thermal-hydraulic simulation is essential for digital twin (DT) technology that supports the safe and efficient operation of small modular reactors (SMRs). Computational fluid dynamics (CFD) provides high-fidelity flow analysis, but its computational cost prevents direct use in DT applications. AI-based surrogate modeling has been actively investigated to address this limitation, yet neural operator--based surrogates for CFD-level transient analysis of SMR-specific geometries have not been reported. This study presents an integrated framework that combines a reduced-order model (ROM) with neural operators, applied to the helical coil steam generator (HCSG) of the System-integrated Modular Advanced Reactor (SMART). Two ROM strategies tailored to each CFD data type were compared, an MLP-based autoencoder (AE) for unstructured mesh data and a convolutional autoencoder (CAE) for structured mesh data, and each was coupled with the deep operator network (DeepONet) to construct the latent DeepONet (L-DeepONet). The Fourier neural operator (FNO) was additionally adopted for comparison. A multi-scale technique was incorporated into both frameworks to mitigate spectral bias and improve the prediction of Kármán vortex streets developing inside the HCSG. The multi-scale L-DeepONet captured the instantaneous periodic vortex dynamics in both velocity and pressure fields, while the FNO and its multi-scale variant predicted the time-averaged mean flow and provided reliable pressure drop estimates. These complementary characteristics provide a practical model-selection guideline that links each architecture to specific DT objectives based on CFD data type and the required level of flow resolution.
Microfluidic channels have emerged as useful tools to control dynamic forcing on transported microscale objects, as encountered in emulsions, biological flows, and other soft matter systems. Tailored channel designs enable precise interfacial and bulk rheological measurements of complex materials over a wide range of forcing timescales. After a brief overview of recent experiments illustrating these techniques, we discuss perspectives for future research in this direction, including the study of lubrication films in highly confined droplets, the measurement of fast relaxation dynamics of complex interfaces, and the high-throughput rheological characterization of microscopic soft matter systems ranging from single macromolecules to cells.
In a low Reynolds number fluid environment that microswimmers encounter, back-and-forth motion cannot lead to net displacement. In mammalian sperm, the mechanical wave propagating along their single flagellum breaks the cancellation between back-and-forth motion and, therefore, assumed to define the movement direction. Here, we show experimentally that the movement direction deviates from the opposite of the wave propagation direction when sperm move at the interface of a viscoelastic fluid and a solid substrate. In fact, the oscillation of the movement direction is out of phase with the oscillation of the tail wave direction, in phase with the head, and the movement has a larger amplitude than the wave direction. When we tried to reconstruct the movement direction as a linear combination of the head orientation and the wave direction, we found that the contributions from these two varied dynamically in time. Further, the last bend of the flagellum does not move in the lab frame (as observed under the microscope). We characterized this as an approximate semi-holonomic constraint on the speed of wave propagation, flagellar sliding, and cell forward movement. Overall, our results highlight the appearance of head and tail taking up roles in directing sperm motility.