The motion of small non-spherical particles is often studied using the unsteady Stokes equations. Zhang & Stone (J. Fluid Mech., vol. 367, 1998, pp. 329–358) reported an asymptotic treatment for nearly spherical particles, to first order in particle non-sphericity, i.e. $O(\epsilon )$, where $\epsilon$ quantifies the shape deviation from a sphere. Importantly, key physical phenomena are absent at $O(\epsilon )$, including (1) coupling between the torque experienced by the particle and its linear translation, (2) coupling between the force the particle experiences and its rotation and (3) the effect of non-sphericity on the orientation averages of these forces and torques. We present an explicit asymptotic theory to second order in particle non-sphericity, i.e. $O(\epsilon ^2)$, for the force and torque acting on a particle in a general unsteady Stokes flow. The derived analytical formulae apply to particles of arbitrary shape, providing the leading-order asymptotic theory for the three above-mentioned phenomena. The theory is demonstrated for several example nearly spherical particles including a spheroid, a ‘pear-shaped’ particle and a simple model for a SARS-CoV-2 virion. This includes formulae for force and torque as a function of particle orientation and their corresponding orientation averages. Our study reveals that the orientation-averaged forces and torques experienced by a nearly spherical particle cannot be generally represented by a perfect sphere. The reported formulae are validated using finite-amplitude three-dimensional direct numerical simulations of the Navier–Stokes equations. A Mathematica notebook is also provided, facilitating implementation of the theory for particle shapes of the user's choosing.

We investigate the effect of three-dimensionality on the synchronisation characteristics of the wake behind an oscillating circular cylinder at ${\textit {Re}} = 300$. Cylinder oscillations in rotation, transverse translation and streamwise translation are considered. We utilise phase-reduction analysis, which quantifies the phase-sensitivity function of periodic flows, to examine the synchronisation properties. Here, we present an ensemble-based framework for phase-reduction analysis to handle three-dimensional wakes that are not perfectly time-periodic. Based on the phase-sensitivity functions, synchronisability to three types of cylinder oscillations is evaluated. In spite of similar trends, we find that phase-sensitivity functions involving three-dimensional wakes are lower in magnitude compared with those of two-dimensional wakes, which leads to narrower conditions for synchronisation to weak cylinder oscillations. We unveil that the difference between the phase-sensitivity functions of two- and three-dimensional flows is strongly correlated to the amplitude variation of the three-dimensional flow by the cylinder motions. This finding reveals that the cylinder motion modifies the three-dimensionality of the wake as well as the phase of vortex shedding, which leads to reduced phase modulation. The synchronisation conditions of three-dimensional wakes, predicted by phase-reduction analysis, agree with the identification by parametric studies using direct numerical simulations for forced oscillations with small amplitudes. This study presents the potential capability of phase-reduction to study synchronisation characteristics of complex flows.

To aid in prediction of turbulent boundary layer flows over rough surfaces, a new model is proposed to estimate hydrodynamic roughness based solely on geometric surface information. The model is based on a fluid-mechanics motivated geometric parameter called the wind-shade factor. Sheltering is included using a rapid algorithm adapted from the landscape shadow literature, while local pressure drag is estimated using a piecewise potential flow approximation. Similarly to evaluating traditional surface parameters such as skewness or average slope magnitude, the wind-shade factor is purely geometric and can be evaluated efficiently from knowing the surface elevation map and the mean flow direction. The wind-shade roughness model is applied to over 100 different surfaces available in a public roughness database and some others, and the predicted sandgrain-roughness heights are compared with measured values. Effects of various model ingredients are analysed, and transitionally rough surfaces are treated by adding a term representing the viscous stress component.

We establish a theoretical framework for predicting friction and heat transfer coefficients in variable-property forced air convection. Drawing from concepts in high-speed wall turbulence, which also involves significant temperature, viscosity and density variations, we utilize the mean momentum balance and mean thermal balance equations to develop integral transformations that account for the impact of variable fluid properties. These transformations are then applied inversely to predict the friction and heat transfer coefficients, leveraging the universality of passive scalars transport theory. Our proposed approach is validated using a comprehensive dataset from direct numerical simulations (DNS), covering both heating and cooling conditions up to a friction Reynolds number $\textit {Re}_\tau \approx 3200$. The predicted friction and heat transfer coefficients closely match the DNS data with accuracy margin 1–2 %, representing a significant improvement over the current state of the art.

Understanding the linear growth of disturbances due to external forcing is crucial for flow stability analysis, flow control, and uncertainty quantification. These applications typically require a large number of forward simulations of the forced linearized dynamics, often in a brute-force fashion. When dealing with simple steady-state or periodic base flows, there exist powerful and cost-effective solution operator techniques. Once these solution operators are constructed, they can be used to determine the response to various forcings with negligible computational cost. However, these methods do not apply to problems with arbitrarily time-dependent base flows. This paper develops and investigates reduced-order modelling with time-dependent bases to build low-rank solution operators for forced linearized dynamics with arbitrarily time-dependent base flows. In particular, we use forced optimally time-dependent decomposition (f-OTD), which extracts the time-dependent correlated structures of the flow response to various excitations. Several demonstrations are included to illustrate the utility of the f-OTD low-rank approximation for performing global transient stability analysis. Additionally, we demonstrate the application of f-OTD in computing the post-transient response of linearized Navier–Stokes equations to a large number of impulses, which has applications in flow control.

We study single bubble deformation statistics in an homogeneous and isotropic turbulent flow by means of direct numerical simulations. We consider bubbles at low Weber number ($We <3$) that have not been broken. We show that we can reproduce bubble deformations with a linear dynamics for each spherical harmonic mode. Inferring the coefficients of the linear model from the DNS data, we find that the natural frequency corresponds to the Rayleigh frequency, derived in a quiescent flow. However, the effective damping increases by a factor 7 compared with the quiescent case, at Taylor Reynolds number $\textit {Re}_\lambda = 55$. Looking at the flow structure around the bubble, we argue that the enhanced damping originates from a thick boundary layer surrounding the bubble. We demonstrate that the effective forcing, originating from the turbulent flow forcing on the bubble surface, is independent of bubble deformability. Therefore, the interface deformations are only one-way coupled to the flow. From this model we conclude that bubbles break rather from turbulent fluctuations than from a resonant mechanism. Eventually, we investigate the pressure modes’ statistics in the absence of bubbles and compare them with the effective forcing statistics. We show that both fields share the same probability distribution function, characterized by exponential tails, and a characteristic time scale corresponding to the eddy turnover time at the mode scale.

Floating fluid-filled membrane breakwater (FFMB) is a temporary structure that can attenuate waves in the deep sea. In this paper the hydrodynamic performance of the FFMB is analysed by using the eigenfunction expansion boundary element method (EEBEM) and physical model experiments. A general motion equation is derived that considers both the dynamic tension and curvature of the membrane. Moreover, an integral expression for the dynamic tension is provided. On this basis, a linear model for solving wave–membrane interaction is established through the EEBEM. Newly designed physical experiments are performed to verify the model and elucidate the nonlinear characteristics of the FFMB. Following verification of the model, this paper investigates the effects of various structural parameters of the FFMB on the wave transmission coefficient, reflection coefficient, horizontal wave force, vertical wave force and dynamic tension. Furthermore, the interrelationship between the structural resonant response and the hydrodynamic performance is elucidated, and the optimal density and filling ratio of the FFMB for engineering applications are proposed. The results demonstrate that the numerical and experimental results are in good agreement, indicating that the model and the motion equation are both practical and highly accurate. By optimizing the structural parameters, the FFMB is capable of effectively attenuating waves within a specific frequency band, while minimizing the wave force.

Liquid plug formation in thin channels due to the Plateau–Rayleigh instability of a liquid film is observed in a variety of fields. In this paper, complementarity between theoretical solutions and direct numerical simulations (DNS) based on a front-tracking algorithm is explored to evaluate the importance of inertia for the case of a cylindrical capillary. A linear stability analysis is first performed and DNS results are then used to investigate the spatial distributions of inertial, convective and viscous terms of the Navier–Stokes equation. The existence of both viscous and inertial regimes is evidenced with a threshold given by the film thickness. The presence of the core fluid slows down the instability. In the viscous regime, predictions of the lubrication theory are verified. An example of liquid water as the outer fluid film and water vapour as the inner core fluid is simulated with application to the fuel cells.

We propose a novel approach for non-Newtonian viscoelastic steady flows based on a decomposition of the rate-of-deformation tensor which here, in a simplified version, leads to an anisotropic generalised Newtonian fluid-like model with separated treatment of kinematics pertaining to shear and extensional flows. Care is taken to assure that the approach is objective and does not introduce spurious effects due to dependence on a specific reference frame. This is done by separating the rate-of-deformation tensor into shear and elongational components, with the local scalar shear and elongation rates being the second invariant of each tensor separately, and by modifying the method used to separate the rate-of-deformation tensor so that it becomes independent of superimposed rigid rotations, thus satisfying the principle of material indifference. We assess the model with two test cases: planar contraction flow, often employed to evaluate numerical methods or constitutive equations and for which experimentally observed corner vortex enhancement and pressure drop increase are seldom found in numerical simulations; and flow past confined or unconfined cylinders, for which experiments indicate a drag increase due to elasticity and most predictions give a drag decrease. With our anisotropic model, incorporating additional elongational-flow-related terms, large vortices and accentuated pressure drop coefficients can be predicted for the contraction problem and enhanced drag coefficient for the cylinder problems. This is the first work where problems of non-Newtonian fluid mechanics are solved numerically with separation of strain rate into shear and elongational components.

Dryland cities are important locations for human–environmental interactions and differ in many key characteristics from cities in wetter environments. Defined by chronic water deficit, these cities face challenges that include securing essential resources, reducing vulnerability to hazards and conserving threatened species. The resilience of dryland cities depends on interactions across the entire urban continuum, from urban cores and suburban areas to teleconnected zones and wildland–urban interfaces. Resilience solutions must enhance the well-being of residents and institutions while fostering adaptive capacity throughout the urban continuum. Key axes of solutions include hydrologic integration, including stormwater capture and reuse, nature-based solutions, including expanding urban tree cover for cooling and health benefits, and landscape sustainability, including the incorporation of spatial heterogeneity into planning and development. Addressing the large uncertainties in ensuring more resilient cities requires convergence research, the integration of theoretically driven science that brings researchers and stakeholders together to identify problems, solutions and opportunities for action. While convergence approaches look to address pressing scientific uncertainties, they also are inherently place-based and address compelling case studies to understand system dynamics and improve decision-making and land management. New research is needed to address the trade-offs resulting from decision-making and urban management activities, to meet the needs of diverse stakeholders and to ensure that policies do not marginalize underserved communities. By leveraging innovative technologies, sustainable practices and community involvement, dryland cities can overcome the challenges posed by chronic water limitations and thrive in their diverse environments.