2025-05-09 | | Total: 10
Subsurface gas storage is a critical technology in global efforts to mitigate climate change. In particular, underground hydrogen storage offers a promising solution for integrating large-scale renewable energy into the power grid. When injected into the subsurface, hydrogen's low viscosity compared to the resident brine causes it to spread rapidly, forming a thin gas layer above the brine, complicating recovery. In long aquifers, the large viscous pressure drop between the source and the outlet induces significant pressure variations, which may lead to substantial density changes in the injected gas. To examine the role of gas compressibility in the spreading dynamics, we use long-wave theory to derive coupled nonlinear evolution equations for the gas pressure and liquid/gas interface height, focusing on the limit of long domains, weak gas compressibility and low liquid/gas mobility ratio. Simulations using these equations are supplemented with a comprehensive asymptotic analysis of parameter regimes. Unlike the near-incompressible limit, in which gas spreading rates are dictated by the source strength and mobility ratio, and any compressive effects are transient, we show how, in general, compression of the main gas bubble can generate dynamic pressure changes that are coupled to those in the thin gas layer that spreads over the liquid, with compressive effects having a sustained influence along the layer. This allows compressibility to reduce spreading rates and lower gas pressures. We characterise this behaviour via a set of low-order models that reveal dominant scalings, highlighting the role of compressibility in mediating the evolution of the gas layer.
This paper introduces an inviscid Computational Fluid Dynamics (CFD) approach for the rapid aerodynamic assessment of Flettner rotor systems on ships. The method relies on the Eulerian flow equations, approximated utilizing a state-of-the-art Finite Volume Method with a dynamic momentum source term to enforce rotor circulation. The method offers substantial computational savings by avoiding near-wall refinement and easing time step constraints, making it ideal for early design phases such as design space exploration. Validation against potential flow theory and viscous reference simulations confirms that the method reliably predicts lift-induced forces despite its limitations in capturing parasitic drag. Three-dimensional simulations, including idealized wind tunnel setups and full-scale ship applications at high Reynolds numbers (up to ReL=1E08), demonstrate that results based on low-order convection featuring a solid numerical viscosity yield deviations with respect to viscous reference data of around O(10%). Accepting this reasonable loss of predictive accuracy provides a simulation framework with response times in the order of minutes compared to hours or even days.
We present a new characterization of small-scale intermittency of turbulence based on an information-theoretic measure, that inherently segregates turbulence intermittency from kinematic intermittency. Instead of the commonly studied higher-order moments, Kullback-Leibler (KL) divergence is used to quantify the deviation of turbulence pseudodissipation rate (or dissipation rate/enstrophy) from that of a Gaussian random velocity field as an accurate measure of the small-scale intermittency arising from turbulence dynamics. In addition, Shannon entropy is used to assess the uncertainty of these turbulence small-scale quantities. Analysis of direct numerical simulation data of forced homogeneous isotropic turbulent flow reveals the existence of two critical Taylor Reynolds numbers where the variation of uncertainty changes. The KL-divergence measure reveals two striking results about intermittency: (i) intermittency grows logarithmically with Reynolds number and (ii) dissipation rate and enstrophy are equally intermittent in a turbulent flow of any Reynolds number, if one considers only the intermittency arising from turbulence dynamics.
Accurate prediction of the hydrodynamic forces on particles is central to the fidelity of Euler-Lagrange (EL) simulations of particle-laden flows. Traditional EL methods typically rely on determining the hydrodynamic forces at the positions of the individual particles from the interpolated fluid velocity field, and feed these hydrodynamic forces back to the location of the particles. This approach can introduce significant errors in two-way coupled simulations, especially when the particle diameter is not much smaller than the computational grid spacing. In this study, we propose a novel force correlation framework that circumvents the need for undisturbed velocity estimation by leveraging volume-filtered quantities available directly from EL simulations. Through a rigorous analytical derivation in the Stokes regime and extensive particle-resolved direct numerical simulations (PR-DNS) at finite Reynolds numbers, we formulate force correlations that depend solely on the volume-filtered fluid velocity and local volume fraction, parametrized by the filter width. These correlations are shown to recover known drag laws in the appropriate asymptotic limits and exhibit a good agreement with analytical and high-fidelity numerical benchmarks for single particle cases, and, compared to existing correlations, an improved agreement for the drag force on particles in particle assemblies. The proposed framework significantly enhances the accuracy of hydrodynamic force predictions for both isolated particles and dense suspensions, without incurring the prohibitive computational costs associated with reconstructing undisturbed flow fields. This advancement lays the foundation for robust, scalable, and high-fidelity EL simulations of complex particulate flows across a wide range of industrial and environmental applications.
Pipe flow is a canonical example of a subcritical flow, the transition to turbulence requiring a finite perturbation. The Reynolds number (Re) serves as the control parameter for this transition, going from the ordered (laminar) to the chaotic (turbulent) phase with increasing Re. Just above the critical Re, where turbulence can be sustained indefinitely, turbulence spreads via the self-replication of localized turbulent structures called puffs. To reveal the workings behind this process, we consider transitions between one and two-puff states, which dynamically are transitions between two distinct chaotic saddles. We use direct numerical simulations to explore the phase space boundary between these saddles, adapting a bisection algorithm to identify an attracting state on the boundary, termed an edge state. At Re=2200, we also examine spontaneous transitions between the two saddles, demonstrating the relevance of the found edge state to puff self-replication. Our analysis reveals that the process of self-replication follows a previously proposed splitting mechanism, with the found edge state as its tipping point. Additionally, we report results for lower values of Re, where the bisection algorithm yields a different type of edge state. As we cannot directly observe splits at this Re, the self-replication mechanism here remains an open question. Our analysis suggests how this question could be addressed in future studies, and paves the way to probing the turbulence proliferation mechanism in other subcritical flows.
A model for the pseudo-turbulent Reynolds stress tensor in compressible flows through monodisperse particle clouds is developed based on data from particle resolved numerical simulations. This model extends previous models for the incompressible regime to subsonic bulk Mach numbers. A Lagrangian model for the local pseudo-turbulent Reynolds stress around a particle, based on the correlation of a local estimate of the Reynolds stress and the drag force on the particle, is also presented. This model can be employed in Euler-Lagrange type simulations. Additionally, corresponding models for the pseudo-turbulent transport of turbulence kinetic energy, the velocity triple-correlation in the volume averaged total energy equation, is also introduced as a first effort aimed at closing the energy equation.
This paper presents an adjoint-assisted, topology-optimization-inspired approach for analyzing topological sensitivities in fluid domains based on porous media formulations -- without directly utilizing the porosity field as a design variable. Instead, the method evaluates the sensitivity with respect to the porosity parameter via an indirect, adjoint formulation, enabling intuitive visualization of cost-functional improving or deteriorating regions. These sensitivity fields are not used in a classical optimization loop but rather to qualitatively assess spatial areas that indicate potential for flow improvement. A key benefit of the proposed approach is its non-intrusive nature: no modifications to existing adjoint solvers, e.g., for shape optimization purposes, are required. As such, the method can be interpreted as an alternative post-processing step that supports early-stage design insights. The approach is demonstrated on a Flettner rotor operating in atmospheric wind at a diameter-based Reynolds number and spinning ratio of ReD=2E+06 as well as k=3 to highlight its practical relevance and interpretability. Placing a container stack in cost-functional improving or deteriorating areas validates the method's prediction quality of topological flow sensitivities.
Surface bubbles in the ocean are critical in moderating several fluxes between the atmosphere and the ocean. In this paper, we experimentally investigate the drainage and lifetime of surface bubbles in solutions containing surfactants and salts, subjected to turbulence in the air surrounding them modelling the wind in the ocean. We carefully construct a setup allowing us to repeatably measure the mean lifetime of a series of surface bubbles, while varying the solution and the wind speed or humidity of the air. To that end, we show that renewing the surface layer is critical to avoid a change of the physical properties of the interface. We show that the drainage of the bubbles is well modeled by taking into account the outwards viscous flow and convective evaporation. The lifetime of surface bubbles in solutions containing no salt is given by the time when evaporation becomes the main draining mechanism and independent on surfactant concentration. When salt is added, the same scaling is valid only at high surfactant concentrations. At low concentrations, the lifetime is always smaller and independent of wind speed, owing to the presence of impurities triggering a thick bursting event.
We systematically investigated a reinforcement learning (RL)-based closed-loop active flow control strategy to enhance the lift-to-drag ratio of a wing section with an NLF(1)-0115 airfoil at an angle of attack 5 degree. The effects of key control parameters, including actuation location, observed state, reward function, and control update interval, are evaluated at a chord-based Reynolds number of Re=20,000. Results show that all parameters significantly influence control performance, with the update interval playing a particularly critical role. Properly chosen update intervals introduce a broader spectrum of actuation frequencies, enabling more effective interactions with a wider range of flow structures and contributing to improved control effectiveness. The optimally trained RL controller is further evaluated in a three-dimensional numerical setup at the same Reynolds number. Actuation is applied using both spanwise-uniform and spanwise-varying control profiles. The results demonstrate that the pretrained controller, combined with a physics-informed spanwise distribution, achieves substantial performance gains. These findings extend the feasibility and scalability of a pretrained RL-based control strategy to more complex airfoil flows.
The motion of an immersed granular material in a channel is characterised by complex interactions among the grains, between the grains and the permeating liquid, and between the grains and the channel walls. Here, we develop a reduced-order continuum model for the bulldozing of an immersed, sedimented granular material by a piston in a channel. In our continuum approach, the granular pile and the overlying fluid layer evolve as a system of coupled thin films. We model the granular phase as a dense, porous, visco-plastic material that experiences Coulomb-like friction with the walls. Conservation of mass and momentum under a thin-film approximation leads to an elliptic equation for the velocity of the grains that is coupled with an evolution equation for the height of the granular pile. We solve our model numerically for a variety of different scenarios to explore the interactions between wall friction, internal viscous-like stresses, and fluid flow above and through the pile. We complement our numerical results with a series of experiments that provide insight into the validity and limitations of the model.