2026-07-15 | | Total: 10
We develop a thermodynamically consistent framework for weakly nonlocal continua with tensor-valued internal variables. Let $\boldsymbol{L}=\operatorname{grad}\boldsymbol{v}$, with stretching tensor $\boldsymbol{D}=\operatorname{sym}\boldsymbol{L}$ and spin tensor $\boldsymbol{W}=\operatorname{skw}\boldsymbol{L}$. We introduce the generators $\boldsymbolΓ_α=\boldsymbol{W}+α\boldsymbol{D}$, $α\in\{0,1\}$, which unify corotational and upper-convected transport. The resulting kinematic structure induces canonical frame-indifferent evolutions of both the internal variable and its spatial gradient, thereby providing a closure for gradient-dependent theories. Starting from the balances of linear momentum and microforces, together with an internal power expenditure depending on the internal variable and its gradient, we derive a local free-energy imbalance for incompressible isothermal processes. Under isotropy and inherited symmetry assumptions, this imbalance admits a canonical decomposition into contributions associated with $\boldsymbol{D}$, $\operatorname{grad}\boldsymbol{L}$, the generator-induced rate $\mathfrak{D}_α\boldsymbol{J}$, and its gradient $\mathfrak{D}^{\nabla}_α(\operatorname{grad}\boldsymbol{J})$. This decomposition yields explicit constitutive restrictions ensuring thermodynamic consistency and identifies the induced higher-order stress contributions arising from gradient dependence. Finally, we construct a coupled gradient theory combining viscoelasticity and constrained orientational order, in which distinct internal variables evolve under different transport mechanisms. The framework extends classical theories with tensorial internal variables, including Oldroyd-B and Landau-de Gennes-type models.
Rogue wave formation and enhancement over coastal areas have been documented over the last decade. This seems to contradict the observed low rogue wave (RW) probability near the surf zone. Without considering wave breaking, RW amplification is expected in this regime. To address this gap, we consider fully nonlinear effects of wave breaking through the proxy of height-to-depth ratio, spatial changes on the wavenumber through the WKB approximation, and slope-corrected refraction mild-slope equations on the energetics of irregular wave fields travelling over a breakwater. By increasing the significant wave height towards the breaking limit, the kinetic energy grows faster than the variance of the surface elevation due to nonlinearity. Thus, the kurtosis decays, albeit not to the point of getting sub-Gaussian statistics. We thereby resolve the apparent paradox of the occurrence probability of RW increase at the beginning of shoaling but subsequently decrease when wave breaking becomes dominant. Motivated by these theoretical developments, we experimentally probe inhomogeneous wave fields nearing the wave-breaking regime. We conduct unidirectional irregular wave experiments in a 30 m long wave flume, generating broad-banded waves over a symmetric submerged breakwater featuring a bottom slope of 1/5, allowing detailed characterization of spectral evolution and the persistence of elevated excess kurtosis even near the breaking limit for this steep slope. We thereby confirm that the excess kurtosis can still be large if the bottom slope is steep, and its maximum value is at least 4 times larger than in other ocean processes, occurring atop the breakwater and about half of its deep water peak wavelength distance after the shoal. As these conditions are typical near shorelines, this understanding is key to coastal areas.
Understanding mixing at density interfaces is essential for predicting transport in stratified environmental flows. Laboratory studies have mostly relied on steady, spatially uniform forcing, whereas turbulence in nature is intermittent and heterogeneous. Here, we present experiments on a two-layer salt-stratified fluid forced by random turbulent bursts generated with a randomly actuated synthetic jet array (RASJA). Density fields are recorded with the light attenuation technique, allowing us to resolve the interface evolution. We measure that the upward velocity of the interface decreases with the density jump, in agreement with the power-law found in previous oscillating-grid studies. At large density differences, the interface sharpens during mixing, contrary to the smaller density jump case. Background potential energy analysis demonstrates irreversible mixing in both cases, with comparable energy changes. These results extend classical laboratory observations to a more isotropic forcing, offering new insights into the dynamics of mixing in geophysical settings.
Closed-loop wall control learnt by multi-agent reinforcement learning can lower skin-friction drag in turbulent channels, but these gradient-based policies are trained on small periodic boxes and exhibit reduced performance when carried over to a larger domain. We recently showed that such policies are also prone to saturated bang-bang actuations that collapse into standing streamwise waves whose scale is set by the computational box rather than by the near-wall cycle, and proposed architectural fixes that avoid these degeneracies. Here, we employ Evolution Strategy (ES) to optimise a recurrent closed-loop controller directly on a large turbulent channel at $\mathit{Re}_τ\simeq180$, evaluating policy performance over full flow episodes using an energy-aware criterion and processing candidate policies in parallel. To our knowledge, this is the first application of an evolution strategy to the control of a turbulent flow. The ES controller reduces the skin friction by about $26\%$, exceeding the gradient-based multi-agent controller of Cavallazzi et al. (2026), GRU-MARL, trained on a minimal box ($17\%$), and marginally exceeding classic opposition control (OC, $22\%$). A wall-normal decomposition of the friction, Reynolds-stress profiles and anisotropy invariants show that the ES and opposition-controlled flows follow separate trajectories through the buffer layer, reaching comparable drag reduction by different reorganisations of the near-wall turbulence. In particular, the ES actuation correlates predominantly with the streamwise velocity fluctuations rather than with the wall-normal velocity that classical OC targets.
Using theory and numerical simulations, here we describe the early stages of the impact with a velocity $V$ of a drop of radius $R_d$ of a low viscosity liquid such as water against a layer of generic thickness $H$ of the same liquid. Our predictions for the initial velocity $V_t\gg V$ and the diameter $H_t\ll R_d$ of the toroidal rim bordering the edge of the thin lamella which is ejected radially outwards after the impact, are in fair agreement with both the numerical results and also with the experimental measurements reported by Zhang et al [J. Fluid Mech, 690, 5, 2012]. Consequently, the present findings can be employed, for instance, to predict the initial diameters and velocities of the fastest tiny droplets which are ejected right after a rain drop falls on a liquid pool, an ubiquitous phenomenon with implications in the dispersal of contaminants and in the generation of aerosols and of ice condensation nuclei
This study extends variable-lattice-density optimization to high-Reynolds-number flows for the design of periodically arranged pin-fin heat sinks. Effective permeability and drag coefficient are identified from unit-cell Reynolds-averaged Navier--Stokes analyses of cylindrical pin-fin arrays and incorporated into a reduced model based on the Darcy--Forchheimer law for macroscopic design exploration. To enable stable optimization under high-Reynolds-number conditions, a dual-mesh framework is introduced, in which the flow field and sensitivities are evaluated on a fine mesh, whereas the design variables are updated on a coarse mesh corresponding to the unit-cell arrangement. For the base condition, the $L_2$ norm of the temperature deviation from the area-averaged temperature is decreased from 4.38 K to 1.00 K in the reduced model, and geometry-resolved analysis of the reconstructed design confirms a reduction from 5.98 K to 1.17 K. Additional calculations with higher inlet velocities and modified outlet locations show that the proposed method captures the dominant flow-redistribution trends under different operating and geometric conditions, although the thermal objective became less accurate when local solid-temperature variations became pronounced. These results indicate that the proposed approach is useful as practical design-exploration for identifying candidate pin-fin heat sink configurations under high-Reynolds-number conditions.
A longstanding challenge in biofluid dynamics research is a mechanistic understanding of the coordinated movement of motile cilia and its resulting ability to facilitate fluid transport. In this study, we develop numerical techniques to simultaneously compute the emergent coordination of and propulsion by filamentous model cilia covering the surface of a sphere. To accomplish this, we develop what we refer to as the filament oscillator model, in which each cilium has two dynamic degrees of freedom: a phase variable that maps to a specific shape in a prescribed sequence, and an angle that describes the overall orientation of the sequence. By varying a parameter related to cilium stiffness, we show that there is bistability between symplectic-like and diaplectic metachronal waves, provided that the stiffness is sufficiently low. Above the critical stiffness, only diaplectic waves emerge. Further, we analyse the propulsive capabilities and flow fields of the two emergent states, showing that diaplectic waves provide more efficient propulsion due to their shorter wavelengths. In addition, we examine how introducing beat-plane tilt leads to ciliate rotation while maintaining nearly identical emergent states and comparable swimming speeds.
We propose a log-Gaussian scale-space limiter for hybrid continuum-ballistic gas dynamics. The method defines complementary continuum and ballistic weights as Gaussian cumulative probabilities in logarithmic Knudsen-number space. This construction gives a smooth numerical transition between Navier-Stokes-Fourier and kinetic/free-molecular fluxes and suppresses asymptotically invalid correction branches in the continuum and free-molecular limits. The limiter is incorporated into a conservative finite-volume interface flux by blending a Navier-Stokes-Fourier flux with a half-range Maxwellian kinetic flux. Numerical tests using one-dimensional DVM/BGK Fourier and Couette benchmarks show that rarefaction corrections improve macroscopic profiles relative to NSF. A parameter scan over the transition center and width gives a DVM/BGK-calibrated pair K0=0.03 and sigma=2.5, reducing the combined mean profile error by about 40% for the tested planar wall-bounded flows. The paper is intended as a proof-of-concept numerical method and calibration study.
We investigate homothetic forward self-similar solutions of the incompressible Navier-Stokes equations: the solutions $\overline{U}$ for which both $\overline{U}$ and $β\overline{U}$ are self-similar profiles for some nontrivial $β$. Homothetic solutions are, in addition, the only solutions for which a singular limit argument can be used to prove non-uniqueness of Leray-Hopf solutions along the lines of the Jia-Šverák program. In three dimensions, and for sufficiently regular initial data, we prove a Liouville theorem that rules out the existence of non-trivial homothetic solutions. In two dimensions, our Liouville theorem proves that the only decaying homothetic solution is the Oseen vortex. In addition, we prove that the Euler operator linearized around the Oseen vortex is stable. On the other hand, we also discover a homothetic solution for which an unstable approximate eigenvalue of the linearized Euler operator exists.
The classical Gibbs-Duhem relation applies to quasi-static processes and neglects kinetic effects, leaving a fundamental gap between Gibbs thermodynamics and Newtonian mechanics. Here, we derive a generalized Gibbs-Duhem framework that incorporates kinetic contributions, thereby establishing a unified connection between classical thermodynamics and Newtonian mechanics. Based on this framework, we propose an alternative evolution equation governing density dynamics at fluid-fluid interfaces. In appropriate limiting cases, the resulting density evolution equation naturally recovers the definition of the speed of sound, Bernoulli's law, and the van der Waals equation of state (EOS).