2026-02-10 | | Total: 16
We study the dimensionality of two-dimensional Kolmogorov flows over a wide range of Reynolds numbers and forcing wavenumbers $k_f=\{2,4,8\}$ using two complementary approaches: convolutional autoencoders and a Kaplan-Yorke estimation based on Lyapunov analysis. As the Reynolds number increases, two distinct transitions are observed: the first corresponds to the destabilization of a periodic orbit, while the second marks the saturation of the large-scale motions. When expressed in terms of the forcing Reynolds number, these transitions occur at nearly the same value for all forcing wavenumbers, suggesting a universal scaling with respect to the forcing scale. By filtering the data to retain only the large-scale range ($k < k_f$), we show that the dimensionality estimated by the autoencoders also saturates at the second transition, implying that once the large scales are fully developed, the subsequent increase in dynamical activity occurs predominantly at smaller scales. At higher Reynolds numbers, the Kaplan-Yorke dimension ceases to grow, revealing its limited sensitivity to the nonlinear interactions that dominate in this regime. Both the Kaplan-Yorke saturation dimension and the filtered large-scale dimensionalities exhibit a linear dependence on $k_f$, indicating that the number of active degrees of freedom scales with the forcing scale rather than with the total number of available Fourier modes.
Water-repellent properties of superhydrophobic surfaces make them promising for anti-icing and deicing applications. Through experimental visualization of frozen sessile droplets undergoing melting on superhydrophobic surfaces, we identify a melting mode with the unmelted ice layer deposited at the bottom of the melting droplet, even though the density of ice is lower than that of water. In the deposited mode of the melting process, the time required for the frozen droplet to melt completely is much shorter than that in the floating mode. Force analysis shows that the melted fluid flows along the gas-liquid interface toward the top of the melting droplet, thereby exerting force and then suppressing the upward movement of the unmelted ice layer. Moreover, the flow within the liquid film formed between the unmelted ice layer and the heating wall is dominated by the viscous force, which has a lubrication effect and maintains the deposition of the unmelted ice layer. High heating temperature, large contact angle, and low particle concentration are helpful for the occurrence of the deposited mode.
We conduct numerical simulations of rocket plume impingement on a lunar landing surface using two complementary frameworks: the Boltzmann equation, which naturally captures rarefied gas dynamics, and the Navier Stokes (NS) equations, the conventional workhorse for turbulent flow simulations. We show that subtle rarefaction effects, long considered negligible in turbulent regimes, can become locally dominant within shear layers where viscous stresses predicted by the NS constitutive relation undergo sign reversals. This phenomenon, which we term constitutive degeneracy, produces order-one relative errors in predicted surface shear stress and heat flux. Our results demonstrate that turbulence can expose hidden limits of NS equations with broad implications for high-speed aerodynamics and planetary exploration.
Pumped storage type hydropower plants play an important role in mitigating real-time energy flexibility. Reversible pump-turbines undergo extreme operating conditions such as runaway and speed-no-load. Very limited studies are undertaken to understand the stochastic flow under these conditions in the reversible pump-turbine. The present study investigates the unsteady vortical flow, its transportation, and interaction with the blades at speed-no-load. Large eddy simulations are conducted in both turbine and pump modes. The computational domain contains 120 million nodes. Numerical results provided evidence of a large longitudinal vortex that develops on the high-pressure side of the blade, and transports into the blade passage and develops the unsteady "string of swirls". The results also showed another "string of swirls" in the draft tube, where flow in the center is reversible (pumping). The resulting flow instability is very high, and it has the potential to induce fatigue damage to the blades.
A fully coupled implicit finite-volume algorithm for incompressible viscoelastic interfacial flows is proposed, whereby the viscoelasticity of the flow is described by an upper-convected Maxwell constitutive model, including limited extensibility and shear-thinning behaviour. The governing equations describing the conservation of continuity and momentum, as well as the constitutive model are discretized using standard finite-volume methods and are solved for pressure, velocity and the polymer stress tensor in a single linear system of equations. Treating all terms of the linearized and discretized governing equations implicit in velocity, pressure and/or the components of the polymer stress tensor, a tightly coupled system of equations is obtained. The interface separating the interacting bulk phases and the surface tension acting at the fluid interface are modelled using a state-of-the-art front-tracking method. We demonstrate the capabilities of the proposed numerical framework with four representative test cases, including the deformation of a viscoelastic droplet in shear flow at large Weissenberg numbers of up to Wi=10^4, and the jump discontinuity of the rise velocity of a bubble rising in a viscoelastic liquid as a result of a "negative wake". Contrary to previous studies using segregated algorithms, the proposed fully coupled implicit algorithm does not apply or require a log-conformation approach to predict these flows. Overall, the fully implicit coupled front-tracking formulation provides a robust framework to reliable numerical predictions of strongly elastic interfacial flows at large Weissenberg numbers.
Building upon our first paper [Hu, K.X, & Hu, Y. J. (2025). Hydroacoustic Absorption and Amplification by Turbulence, arXiv:2512.07920], the present work conducts a more in-depth investigation into the impact of turbulence on hydroacoustic waves. The study includes the influences of temperature, unsteady laminar flow, standing wave effects, the alteration of phase, the variation of amplification factor with frequency, and the temporal evolution of the acoustic wave. Experiments indicate that the temperature rise caused by friction between the turbulent flow and the pipe wall does not significantly affect the acoustic wave, and therefore is not the primary cause of changes in wave amplitude. When two transducers are placed opposite each other, the acoustic waves can be regarded as a superposition of traveling waves and standing waves. When the pump is shut down after the pipe flow has stabilized, the temporal evolution of the acoustic waves during the subsequent turbulence decay process can be classified into six types. In addition to altering the amplitude, turbulence also changes the phase of waves. The total phase shift of the acoustic wave along the entire pipe equals the sum of the phase shifts in each segment. Both the amplification factor and the phase shift due to turbulence vary periodically with frequency. Experiments for low-frequency waves show that acoustic waves with frequencies below a certain threshold are not affected by turbulence. These findings suggest that the interaction between hydroacoustic waves and turbulence constitutes stimulated emission in water.
A green transition in aviation requires a drastic upscaling of Sustainable Aviation Fuel (SAF). The power-to-liquid process for the production of CO2-neutral jet fuel via electricity, called e-SAF, directly replaces fossil jet fuel without having to change infrastructure, aeroplanes, or jet-engines. The process combines green hydrogen with industrial exhaust gas, or captured carbon dioxide, in a circular economy concept. A key element of the e-SAF production plant is the reactor where syngas is produced. Traditional reactors use catalytic technology, which faces severe challenges due to the reduced performance over time because of catalyst degradation, clogging, and breakup due to embrittlement. A high-potential alternative is the catalyst-free reverse water-gas-shift (RWGS) reactor concept. The primary aim of this paper is to investigate the fundamental aspects of the catalyst-free RWGS process, such as reaction kinetics and the interactions between turbulence and chemistry. The secondary aim is to identify how a typical combustion subgrid scale models for Large Eddy Simulations (LES) perform when the chemical reactions are endothermic, in contrast to the strong endothermicity associated with classical combustion. It is found that even small traces of O2 in the CO2 stream can significantly increase the production rate of CO. This is attributed to the increased pool of OH. The effect is strongest at atmospheric pressure and less pronounced at higher pressure. By using the temporal jet framework to study turbulence-chemistry interactions, an algebraic equation for the prediction of the CO conversion time in a turbulent flow as a function of Damkohler number and chemical timescale is employed. Finally, it is concluded that the PaSR LES subgrid model designed for combustion reactions perform well also for the endothermic reverse water-gas-shift reaction.
Turbulence is a complex, irregular flow phenomenon ubiquitous in natural processes and engineering applications. The Reynolds-averaged Navier-Stokes (RANS) method, owing to its low computational cost, has become the primary approach for rapid simulation of engineering turbulence problems. However, the inaccuracy of classical turbulence models constitutes the main drawback of the RANS framework. With the rapid development of data-driven approaches, many data-driven turbulence models have been proposed, yet they still suffer from issues of generalizability and accuracy. In this work, we propose a few-shot, physically restorable, symbolic regression turbulence model based on the normalized general effective-viscosity hypothesis. Few-shot indicates that our model is trained on limited flow configurations spanning only a narrow subset of turbulent flow physics, yet can still outperform the baseline model in substantially different turbulent flows. Physically restorable means our model can nearly revert to the baseline model in regimes satisfying specific physical conditions, using only the symbolic regression training results. The normalized general effective-viscosity hypothesis was proposed in our previous study. Specifically, we first formalize the concept of few-shot data-driven turbulence models. Second, we train our symbolic regression turbulence models using only direct numerical simulation (DNS) data for three-dimensional periodic hill flow slices. Third, we evaluate our models on periodic hill flows, zero pressure gradient flat plate flow, NACA0012 airfoil flows, and NASA Rotor 37 transonic axial compressor flows. One of our symbolic regression turbulence models consistently outperforms the baseline model, and we further demonstrate that this model can nearly revert to baseline behavior in certain flow regimes.
It is shown that the upper limit on the fraction of power that can be extracted from an airstream approaching a wind turbine is about 78% rather than the historical value of 59%. This higher limit is based on the assumption that the wind turbine cannot accelerate the flow at its exit plane above the freestream velocity. The derivation of the historical limit of 59% must either violate the angular momentum equation or the continuity equation.
Accurate prediction of aerothermal loads in hypersonic flows is critical yet challenging due to the coupling of Shock-Wave/Boundary-Layer Interactions (SBLI) and thermal non-equilibrium. This work presents the development of a three-dimensional two-temperature Gas-Kinetic Scheme (3D 2T-GKS) on unstructured meshes. The scheme resolves translational-rotational and vibrational energy modes within a unified kinetic framework. A key innovation is the integration of a Generalized Kinetic Boundary Condition (GKBC), which physically decouples the thermal accommodation of vibrational energy from the translational-rotational mode, thereby offering a more accurate model for gas-surface interactions. Additionally, a Discontinuity Feedback Factor (DFF) is employed to capture strong shock waves with reduced numerical dissipation compared to classical limiters. The method is rigorously validated against standard experimental benchmarks, including the sharp double-cone and hollow cylinder-flare configurations. Numerical results demonstrate that the proposed solver, augmented by the GKBC, accurately captures complex wave structures, separation topologies, and surface heat flux distributions. These findings confirm the robustness and fidelity of the 3D 2T-GKS for simulating complex hypersonic non-equilibrium flows.
Spatial filters have played a central role in large eddy simulation and, more recently, in reduced order model (ROM) stabilization for convection-dominated flows. Nevertheless, important open questions remain: in under-resolved regimes, which filter is most suitable for a given stabilization or closure model? Moreover, once a filter is selected, how should its parameters, such as the filter radius, be determined? Addressing these questions is essential for the reliable design and performance of filter-based stabilization strategies. To answer these questions, we propose a novel strategy that differs fundamentally from current filter-based approaches: we replace traditional spatial filters with a data-driven stabilization operator (StabOp) that yields accurate results for a given resolution, quantity of interest, and stabilization strategy. Although the new StabOp can be used for both classical discretizations and ROMs, and for different types of filter-based stabilization or closure, for clarity we focus on ROMs and the Leray stabilization. To build the new StabOp, we postulate its model form as a linear, quadratic, or nonlinear mapping, and then solve a PDE-constrained optimization problem to minimize a given loss function. Using the resulting StabOp in the Leray ROM (L-ROM) yields a new stabilized ROM, StabOp-L-ROM. To assess the StabOp-L-ROM, we compare it with the L-ROM and the standard ROM in numerical simulations of four flows: 2D flow past a cylinder at Re=500, lid-driven cavity at Re=10000, 3D flow past a hemisphere at Re=2200, and minimal channel flow at Re=5000. Our numerical results show that the StabOp-L-ROM can be orders of magnitude more accurate than the classical L-ROM tuned with an optimal filter radius in the predictive regime. Furthermore, while the new StabOp smooths the input flow fields, its smoothing mechanism differs from that of classical spatial filters.
In hypersonic rarefied flows, insufficient intermolecular collisions cause significant deviations between translational and rotational temperatures, leading to strong thermal nonequilibrium. For diatomic gases such as nitrogen and oxygen, the direct simulation Monte Carlo (DSMC) method commonly employs the Borgnakke-Larsen (BL) model to simulate translational-rotational energy exchange (relaxation) processes. Although widely used, the BL model lacks a rigorous theoretical foundation and assumes that only a fraction of collisions lead to rotational relaxation. To address these shortcomings, Pullin introduced a kinetically consistent relaxation model into the gas kinetic theory. By employing the Beta function for energy partitioning, a concrete collision cross section that satisfies the detailed balance condition is constructed. In this study, a comparative investigation of the BL and Pullin models is performed within the DSMC framework, where both original and simplified equations are considered and parameterized by physical accommodated coefficient in the Beta function. A series of test cases--including zero-dimensional rotational relaxation of nitrogen, one-dimensional planar Couette flow and normal shock wave, two-dimensional hypersonic flow past a cylinder, and three-dimensional hypersonic flow around an X38-like vehicle--are performed to assess the accuracy and efficiency of these models. The results confirm the consistency between the Pullin and BL models. Owing to its rigorous theoretical foundation and accurate physical representation, the Pullin model is expected to provide substantial support for the extension of subsequent theoretical studies and numerical simulations. Moreover, in the highly rarefied flow regime (Knudsen number greater than 1, or altitudes above 100 km), the simplified Pullin model exhibits performance comparable to that of the BL model.
We perform axisymmetric numerical simulations to investigate the coalescence dynamics of a liquid drop in a deep liquid pool. This study aims to generalize the mechanisms of partial coalescence across a range of drop shapes, elucidate the underlying mechanism of neck oscillations, and examine the roles of inertial, viscous and gravitational forces, quantified by the Weber, Ohnesorge, and Bond numbers, in governing the coalescence behavior. A phase diagram is constructed to delineate the boundaries between partial and complete coalescence regimes based on these dimensionless parameters. Our analysis of the height-to-neck ratio shows that, upon contact with the pool, the primary drop forms an upward liquid column that ultimately pinches off due to inwardly directed horizontal momentum. Additionally, the study suggests that as the dimensionless numbers increase, the effect of the vertical collapse rate plays a significant role in the outcome of the coalescence process. Notably, the Rayleigh-Plateau instability is found to be insignificant in driving partial coalescence within the explored parameter space. We identified a transition regime between partial and complete coalescence, characterized by multiple neck oscillations that delay the pinch-off of secondary droplets. The formation of secondary droplets is most prominent for prolate drops, followed by spherical and oblate drops of comparable volume. Furthermore, we observe that the tendency to form multiple droplets from elongated liquid columns diminishes with an increase in the impact velocity of the primary drop.
The Kelvin-Helmholtz (KH) instability occurs in multiple heliospheric (solar-wind stream interfaces, planetary magnetospheres, cometary tails, heliopause flanks) and interstellar (protoplanetary disks, relativistic jets, neutron star accretion disks) environments. While the KH instability has been well-studied in the magnetohydrodynamic (MHD) limit, only limited studies were performed in the collisionless regime, which is conducive to development of anisotropic pressures. Collisionless plasmas are often described using the Chew Goldberger and Low (CGL) equations which feature an anisotropic pressure tensor. This paper presents a comprehensive analysis of the CGL version of the KH instability using linearised and numerical techniques. We find that the largest growth rates and the greatest incidence of magnetic effects occur in the MHD limit. In the large relaxation time CGL limit, part of the energy goes into the formation of pressure anisotropies, resulting in smaller amounts of energy being available for bending the field lines. Consequently, when we cross-compare CGL and MHD simulations that are otherwise identical, the current densities are largest in the MHD limit, and the largest magnetic islands also form in that limit. Early and late time formation of pressure anisotropies have also been studied. We also find that the strongest trend for forming intermittencies in the flow also occurs in the MHD limit. The paper also discusses possible consequences of our results for turbulence and reconnection in the heliosheath (the layer between the solar wind termination shock and the heliopause).
Spatially varying flow laws have been identified in dry granular flow, yet their applicability to unjammed suspensions remains unclear. This study demonstrates that the quasistatic suspension flow combines dry granular rheology with nonlocal effects in the shear band and damped viscous flow in the outer creep region. Through rotary shear experiments on a hydrogel granular raft, we observe that the flow decays from the interface in the quasistatic regime, where the particles remain mobile even below the yield stress. These findings suggest the universal flow law across the transition between jammed/brittle granular behavior and unjammed/ductile viscous flow.
Space-time materials are obtained by modulating a physical medium with a traveling-wave perturbation of one or several of its constitutive parameters, such as the density or the bulk modulus in the case of acoustic materials. When this modulation has the form of a moving and abrupt (subwavelength) transition between two parameter values, we refer to a spatio-temporal interface, which may be considered as a building block for more complex space-time materials. This work considers the problem interaction and scattering of acoustic waves with a single spatio-temporal interface, and a sequence of two interfaces forming a slab. Several regimes defined by the relation between the sound propagation velocities and the interface velocity (namely subsonic, intersonic, and supersonic regimes) are discussed. Analytical expressions for the frequency conversions and scattering coefficients are obtained, and compared with numerical simulations based on an equivalent FTFD squeme.