2025-04-07 | | Total: 10
Griffin plot relates the peak amplitudes of vortex-induced vibration to structrual mass-damping parameter, known as the Scruton number. Griffin plot serves as a fundamental tool in many engineering fields. This study confirms a general self-similarity in Griffin plots, where plots derived from transient responses at any Scruton number converge to a single, consisten curve. This self-similarity arises from weak aeroelastic nonlinearity in vortex-induced vibration, manifasting as amplitude-dependent aerodynamic damping. Based on this self-similarity property, we propose a numerical method to estimate Griffin plots from transient displacement responses at any Scruton number. The resulting plots align closely with experimental data for both cross-flow and torsional vortex-induced vibrations, highlighting robust self-similar behavior across different Scruton numbers. Furthermore, we observe a consistent trend in Griffin plots for a rectangular cylinder, closed-box, and double-girder bridge deck: the reciprocal of peak amplitudes shows an approximately linear relationship with the Scruton number, especially in torsional vortex-induced vibration. To generate this linearity, we develop a simple empirical model of vortex-induced forces. This model accurately reproduces the Griffin plot for a rectangular cylinder using aeroelastic parameters derived from a single Scruton number, significantly reducing the need for extensive experimental measurements.
We investigate horizontal dispersion of a passive scalar in a porous stratum with Rayleigh-Darcy convection initiated by a geothermal gradient. While increasing Rayleigh number (Ra) keeps enhancing convection, the horizontal dispersion coefficient (ˆD) only scales with Ra in a narrow range, and saturates at Ra≳. We rationalize this two-stage behavior: at low Ra, passive scalar migrates through bulk circumflex; at high Ra, boundary micro-plume layer becomes the dominant scalar pathway. Theory gives saturate \hat{D} value proportional to the Lewis number Le and molecular diffusivity D_0.
Patterning solid surfaces with varying wettability is important to manage droplets in microfluidics, heat transfer and printing. Solid surface roughness poses fundamental limitations including contact-line pinning and solid friction. Here, we report an experimental strategy and theoretical design principles for Patterned Liquid Surfaces (PaLS) that combines the controlled wettability from patterning with the ultra-smoothness of a lubricant-infused surface. In contrast to a solid, on PaLS a droplet can be in 10 different wetting states. This richness arises from the adaptation of the liquid lubricants, and can be harnessed to control the apparent contact angle of the droplet over the full range of wettability whilst removing contact-line pinning effects induced by the solid surface. In the limit of thin liquid films, we derive surface-averaged laws for the apparent contact angle for each wetting state, which capture both experimental and simulation data. Our results provide a distinct approach to surface patterning that exploits the interaction of fluids with lubricant-impregnated surfaces.
This study employs three-dimensional particle-resolved simulations of planar shocks passing through a suspension of stationary solid particles to study wake-induced gas-phase velocity fluctuations, termed pseudo-turbulence. Strong coupling through interphase momentum and energy exchange generates unsteady wakes and shocklets in the interstitial space between particles. A Helmholtz decomposition of the velocity field shows that the majority of pseudo-turbulence is contained in the solenoidal component from particle wakes, whereas the dilatational component corresponds to the downstream edge of the particle curtain where the flow chokes. One-dimensional phase-averaged statistics of pseudo-turbulent kinetic energy (PTKE) are quantified at various stages of flow development. Reduction in PTKE is observed with increasing shock Mach number due to decreased production, consistent with single-phase compressible turbulence. The anisotropy in Reynolds stresses is found to be relatively constant through the curtain and consistent over all the conditions simulated. Analysis of the budget of PTKE shows that the majority of turbulence is produced through drag and balanced by viscous dissipation. The energy spectra of the streamwise gas-phase velocity fluctuations reveal an inertial subrange that begins at the mean interparticle spacing and decays with a power law of -5/3 and steepens to -3 at scales much smaller than the particle diameter. A two-equation model is proposed for PTKE and its dissipation. The model is implemented within a hyperbolic Eulerian-based two-fluid model and shows excellent agreement with the particle-resolved simulations.
Hydrodynamic density functional theory (DFT) is applied to analyse dynamic contact angles of droplets in order to assess its predictive capability regarding wetting phenomena at the microscopic scale and to evaluate its feasibility for multiscale modelling. Hydrodynamic DFT incorporates the influence of fluid-fluid and solid-fluid interfaces into a hydrodynamic theory by including a thermodynamic model based on classical DFT for the chemical potential of inhomogeneous fluids. It simplifies to the isothermal Navier-Stokes equations far away from interfaces, thus connecting microscopic molecular modelling and continuum fluid dynamics. In this work we use a Helmholtz energy functional based on the perturbed-chain statistical associating fluid theory (PC-SAFT) and the viscosity is obtained from generalised entropy scaling, a one-parameter model which takes microscopic information of the fluid and solid phase into account. Deterministic (noise-free) density and velocity profiles reveal wetting phenomena including different advancing and receding contact angles, the transition from equilibrium to steady state and the rolling motion of droplets. Compared to a viscosity model based on bulk values, generalised entropy scaling provides more accurate results, which stresses the importance of including microscopic information in the local viscosity model. Hydrodynamic DFT is transferable as it captures the influence of different external forces, wetting strengths and (molecular) solid roughness. For all results good quantitative agreement with non-equilibrium molecular dynamics simulations is found, which emphasises that hydrodynamic DFT is able to predict wetting phenomena at the microscopic scale.
The extent to which statistical equilibrium theory is applicable to driven dissipative dynamics remains an important open question in many systems. We use extensive direct numerical simulations of the incompressible two-dimensional (2D) Navier-Stokes equation to examine the steady state of large-scale condensates in 2D turbulence at finite Reynolds number Re in the absence of bottom drag. Large-scale condensates appear above a critical Reynolds number Re_c\approx 4.19. Close to this onset, we find a power-law scaling of the energy with Re-Re_c, with the energy spectrum at large scales following the absolute equilibrium form proposed by Kraichnan. At larger Re, the energy spectrum deviates from this form, displaying a steep power-law range at low wave numbers with exponent -5, with most of the energy dissipation occurring within the condensate at large scales. We show that this spectral exponent is consistent with the logarithmic radial vorticity profile of the condensate vortices predicted by quasi-linear theory for a viscously saturated condensate. Our findings shed new light on the classical problem of large-scale turbulent condensation in forced dissipative 2D flows in finite domains, showing that the large scales are close to equilibrium dynamics in weakly turbulent flows but not in the strong condensate regime with Re\gg1.
In this work, Galerkin projection is used to build Reduced Order Models (ROM) for two-dimensional Rayleigh-Bénard (RB) convection with no-slip walls. We compare an uncoupled projection approach that uses separate orthonormal bases for velocity and temperature with a coupled formalism where the equations are projected onto a single basis combining velocity and temperature components. Orthonormal bases for modal projection are obtained as the eigenvalues of the controllability Gramian of the linearized RB equations. Various coupled and uncoupled ROMs with different number of modes are generated and validated against Direct Numerical Simulations (DNS) over a wide range of Rayleigh numbers, Ra. DNS and ROM results are compared in terms of mean vertical profiles, heat flux, flow structures, bifurcation diagrams and energy spectra. Coupled ROMs are found to be unstable at high Ra numbers with a stability limit that depends on the basis Ra. Uncoupled models show an increasing agreement with DNS as a function of the system dimension. It is found that for the system truncations investigated here, a quantitative agreement with DNS can be obtained up to Ra\simeq 4\times 10^5. ROMs are used to perform a bifurcation analysis for Pr=10 and the results compared to DNS. They qualitatively predict the transitions between periodic, quasiperiodic and chaotic states as well as the spectral characteristics over a wide range of Ra numbers. Overall, these results show that these ROMs reproduce the main flow features of RB convection and could be used as DNS-surrogates for the development of active control strategies and state estimation applications.
Wormlike micelle (WLM) solutions are abundant in energy, environmental, and industrial applications, which often rely on their flow through tortuous channels. How does the interplay between fluid rheology and channel geometry influence the flow behavior? Here, we address this question by experimentally visualizing and quantifying the flow of a semi-dilute WLM solution in millifluidic serpentine channels. At low flow rates, the base flow is steady and laminar, with strong asymmetry and wall slip. When the flow rate exceeds a critical threshold, the flow exhibits an elastic instability, producing spatially-heterogeneous, unsteady three-dimensional (3D) flow characterized by two notable features: (i) the formation and persistence of stagnant but strongly-fluctuating and multistable "dead zones" in channel bends, and (ii) intermittent 3D "twists" throughout the bulk flow. The geometry of these dead zones and twisting events can be rationalized by considering the minimization of local streamline curvature to reduce flow-generated elastic stresses. Altogether, our results shed new light into how the interplay between solution rheology and tortuous boundary geometry influences WLM flow behavior, with implications for predicting and controlling WLM flows in a broad range of complex environments.
We consider a simple thermodynamically consistent model that captures the interplay between autocatalytically reacting surfactants, the Marangoni effect and wetting dynamics. An ambient bath of surfactant acts as a chemostat and provides the system with chemical fuel, thereby driving it away from thermodynamic equilibrium. We find that a positive feedback loop between the local reactions and the Marangoni effect induces surface tension gradients that allow for self-propelled droplets. Besides simple directional motion, we find crawling and shuttling droplets as well as droplets performing random walks, thus exploring the entire substrate. We study the occurring dynamics and show how the observed states emerge from local and global bifurcations. Due to the underlying generic thermodynamic structure, we expect that our results are relevant not only to directly related biomimetic droplet systems but also to structurally similar systems like chemically active phase separating mixtures.
Physical systems with complex unsteady dynamics, such as fluid flows, are often poorly represented by a single mean solution. For many practical applications, it is crucial to access the full distribution of possible states, from which relevant statistics (e.g., RMS and two-point correlations) can be derived. Here, we propose a graph-based latent diffusion (or alternatively, flow-matching) model that enables direct sampling of states from their equilibrium distribution, given a mesh discretization of the system and its physical parameters. This allows for the efficient computation of flow statistics without running long and expensive numerical simulations. The graph-based structure enables operations on unstructured meshes, which is critical for representing complex geometries with spatially localized high gradients, while latent-space diffusion modeling with a multi-scale GNN allows for efficient learning and inference of entire distributions of solutions. A key finding is that the proposed networks can accurately learn full distributions even when trained on incomplete data from relatively short simulations. We apply this method to a range of fluid dynamics tasks, such as predicting pressure distributions on 3D wing models in turbulent flow, demonstrating both accuracy and computational efficiency in challenging scenarios. The ability to directly sample accurate solutions, and capturing their diversity from short ground-truth simulations, is highly promising for complex scientific modeling tasks.