Large-eddy simulation (LES) using an unstructured overset numerical method is performed to study the flow around a ducted marine propeller for the highly unsteady off-design condition called crashback. Known as one of the most challenging propeller states to analyse, the propeller rotates in the reverse direction to yield negative thrust while the vehicle is still in forward motion. The LES results for the marine propeller David Taylor Model Basin 4381 with a neutrally loaded duct are validated against experiments, showing good agreement. The simulations are performed at Reynolds number of 561 000 and an advance ratio $J=-0.82$. The flow field around the different components (duct, rotor blades and stator blades) and their impact on the unsteady loading are examined. The side-force coefficient $K_S$ is mostly generated from the duct surface, consistent with experiments. The majority of the thrust and torque coefficients $K_T$ and $K_Q$ arise from the rotor blades. A prominent contribution to $K_Q$ is also produced from the stator blades. Tip-leakage flow between the rotor blade tips and duct surface is shown to play a major role in the local unsteady loads on the rotor blades and duct. The physical mechanisms responsible for the overall unsteady loads and large side-force production are identified as globally, the vortex ring and locally, leading-edge separation as well as tip-leakage flow which forms blade-local recirculation zones.

Sash windows, or counter-balanced vertically sliding windows, are a common feature in domestic buildings. We describe a model for the ventilation flow rate through a sash window, when the indoor temperature is warmer than outdoors. Depending on the position of the neutral pressure level relative to the sash window, we identify three flow regimes and the critical height of the lower opening for the flow to transition from one regime to another. We perform laboratory experiments in a water tank to measure the flow rate for different sash window geometries, and compare our experimental results with the model. Using our results we assess the optimal sash window arrangement for different natural ventilation strategies. Our results have implications for optimal ventilation rates and control of ventilation in smart buildings.

We propose an improved adjoint-based method for the reconstruction and prediction of the nonlinear wave field from coarse-resolution measurement data. We adopt the data assimilation framework using an adjoint equation to search for the optimal initial wave field to match the wave field simulation result at later times with the given measurement data. Compared with the conventional approach where the optimised initial surface elevation and velocity potential are independent of each other, our method features an additional constraint to dynamically connect these two control variables based on the dispersion relation of waves. The performance of our new method and the conventional method is assessed with the nonlinear wave data generated from phase-resolved nonlinear wave simulations using the high-order spectral method. We consider a variety of wave steepness and noise levels for the nonlinear irregular waves. It is found that the conventional method tends to overestimate the surface elevation in the high-frequency region and underestimate the velocity potential. In comparison, our new method shows significantly improved performance in the reconstruction and prediction of instantaneous surface elevation, surface velocity potential and high-order wave statistics, including the skewness and kurtosis.

This work explores the relationship between wind speed and time-dependent structural motion response as a means of leveraging the rich information visible in flow–structure interactions for anemometry. We build on recent work by Cardona, Bouman and Dabiri (Flow, vol. 1, 2021, E4), which presented an approach using mean structural bending. Here, we present the amplitude of the dynamic structural sway as an alternative signal that can be used when mean bending is small or inconvenient to measure. A force balance relating the instantaneous loading and instantaneous deflection yields a relationship between the incident wind speed and the amplitude of structural sway. This physical model is applied to two field datasets comprising 13 trees of 4 different species exposed to ambient wind conditions. Model generalization to the diverse test structures is achieved through normalization with respect to a reference condition. The model agrees well with experimental measurements of the local wind speed, suggesting that tree sway amplitude can be used as an indirect measurement of mean wind speed, and is applicable to a broad variety of diverse trees.

While there have been numerous applications of large eddy simulations (LES) to complex flows, their application to practical engineering configurations, such as full aircraft models, have been limited to date. Recently, however, advances in rapid, high quality mesh generation, low-dissipation numerical schemes and physics-based subgrid-scale and wall models have led to, for the first time, accurate simulations of a realistic aircraft in landing configuration (the Japanese Aerospace Exploration Agency Standard Model) in less than a day of turnaround time with modest resource requirements. In this paper, a systematic study of the predictive capability of LES across a range of angles of attack (including maximum lift and post-stall regimes), the robustness of the predictions to grid resolution and the incorporation of wind tunnel effects is carried out. Integrated engineering quantities of interest, such as lift, drag and pitching moment will be compared with experimental data, while sectional pressure forces will be used to corroborate the accuracy of the integrated quantities. Good agreement with experimental $C_L$ data is obtained across the lift curve with the coefficient of lift at maximum lift, $C_{L,max}$, consistently being predicted to within five lift counts of the experimental value. The grid point requirements to achieve this level of accuracy are reduced compared with recent estimates (even for wall modelled LES), with the solutions showing systematic improvement upon grid refinement, with the exception of the solution at the lowest angles of attack, which will be discussed later in the text. Simulations that include the wind tunnel walls and aircraft body mounting system are able to replicate important features of the flow field noted in the experiment that are absent from free air calculations of the same geometry, namely, the onset of inboard flow separation in the post-stall regime. Turnaround times of the order of a day are made possible in part by algorithmic advances made to leverage graphical processing units. The results presented herein suggest that this combined approach (meshing, numerical algorithms, modelling, efficient computer implementation) is on the threshold of readiness for industrial use in aeronautical design.

We present a theoretical model and experimental demonstration of thin liquid film deformations due to a dielectric force distribution established by surface electrodes. We model the spatial electric field produced by a pair of parallel electrodes and use it to evaluate the stress on the liquid–air interface through Maxwell stresses. By coupling this force with the Young–Laplace equation, we obtain the deformation of the interface. To validate our theory, we design an experimental set-up which uses microfabricated electrodes to achieve spatial dielectrophoretic actuation of a thin liquid film, while providing measurements of microscale deformations through digital holographic microscopy. We characterize the deformation as a function of the electrode-pair geometry and film thickness, showing very good agreement with the model. Based on the insights from the characterization of the system, we pattern conductive lines of electrode pairs on the surface of a microfluidic chamber and demonstrate the ability to produce complex two-dimensional deformations. The films can remain in liquid form and be dynamically modulated between different configurations or polymerized to create solid structures with high surface quality.

It is demonstrated how aqueous droplets with volumes down to the sub-femtolitre range can be manipulated, including the withdrawal of minute samples from the droplets. The underlying principle is that of partial coalescence with a liquid reservoir in an applied electric field. Upon partial coalescence, a droplet merges with a reservoir and reappears with a smaller diameter. The droplets studied here perform a reciprocating motion between two reservoirs during which their volume gets reduced. Manipulation of droplets with diameters down to 400 nm is reported. A similarity relation is derived expressing the ratio of droplet diameters before and after partial coalescence as a function of the ratio between electric and interfacial-tension forces. The presented scheme allows the withdrawal of minute samples from small droplets and could prove helpful in various applications where droplets are used as tiny reaction spaces or when the goal is to tailor the size of individual droplets.

Water wave resonance between two side-by-side vessels is a multimode resonant hydrodynamic phenomenon with low damping. The potential flow damping and viscous damping inside the gap play a significant role, influencing the amplitudes of the gap resonances. The frequencies of the gap modes can be well predicted by linear potential flow theory, while much effort has been made to explore the nature of the viscous damping. A series of experiments is conducted to explore the temporal (Zhao et al., Journal of Fluid Mechanics, vol. 812, 2017, 905–939) and spatial structure (Zhao et al., Journal of Fluid Mechanics, vol. 883, 2020, A22) of the resonant responses along the gap. Ultimately, it is of practical interest to understand the response statistics along the gap in random seas, to facilitate decision making for safe offshore operations. Following our previous studies which focused on new physics, here we identify the design waves that produce the most probable maximum responses under unidirectional random linear wave excitation. This is achieved through an efficient prediction model within linear theory. Combining the experimental data and linear potential flow calculations, we provide the lower and upper bounds of gap responses, bracketing possible responses at field scale. The statistical model is expected to be of practical importance for offshore operations.

A new guideline for mitigating indoor airborne transmission of COVID-19 prescribes a limit on the time spent in a shared space with an infected individual (Bazant & Bush, Proceedings of the National Academy of Sciences of the United States of America, vol. 118, issue 17, 2021, e2018995118). Here, we rephrase this safety guideline in terms of occupancy time and mean exhaled carbon dioxide (${\rm CO}_{2}$) concentration in an indoor space, thereby enabling the use of ${\rm CO}_{2}$ monitors in the risk assessment of airborne transmission of respiratory diseases. While ${\rm CO}_{2}$ concentration is related to airborne pathogen concentration (Rudnick & Milton, Indoor Air, vol. 13, issue 3, 2003, pp. 237–245), the guideline developed here accounts for the different physical processes affecting their evolution, such as enhanced pathogen production from vocal activity and pathogen removal via face-mask use, filtration, sedimentation and deactivation. Critically, transmission risk depends on the total infectious dose, so necessarily depends on both the pathogen concentration and exposure time. The transmission risk is also modulated by the fractions of susceptible, infected and immune people within a population, which evolve as the pandemic runs its course. A mathematical model is developed that enables a prediction of airborne transmission risk from real-time ${\rm CO}_{2}$ measurements. Illustrative examples of implementing our guideline are presented using data from ${\rm CO}_{2}$ monitoring in university classrooms and office spaces.

The variety of configurations for vertical-axis wind turbines (VAWTs) make the development of universal scaling relationships for even basic performance parameters difficult. Rotor geometry changes can be characterized using the concept of solidity, defined as the ratio of solid rotor area to the swept area. However, few studies have explored the effect of this parameter at full-scale conditions due to the challenge of matching both the non-dimensional rotational rate (or tip speed ratio) and scale (or Reynolds number) in conventional wind tunnels. In this study, experiments were conducted on a VAWT model using a specialized compressed-air wind tunnel where the density can be increased to over 200 times atmospheric air. The number of blades on the model was altered to explore how solidity affects performance while keeping other geometric parameters, such as the ratio of blade chord to rotor radius, the same. These data were collected at conditions relevant to the field-scale VAWT but in the controlled environment of the lab. For the three highest solidity rotors (using the most blades), performance was found to depend similarly on the Reynolds number, despite changes in rotational effects. This result has direct implications for the modelling and design of high-solidity field-scale VAWTs.

We present a microfluidic method to measure the elastic properties of a population of microcapsules (liquid drops enclosed by a thin hyperelastic membrane). The method is based on the observation of flowing capsules in a cylindrical capillary tube and an automatic inverse analysis of the deformed profiles. The latter requires results from a full numerical model of the fluid–structure interaction accounting for nonlinear membrane elastic properties. For ease of use, we provide them under the form of databases, when the initially spherical capsule has a membrane governed by a neo-Hookean or a general Hooke's law with different surface Poisson ratios. Ultimately, the microfluidic method yields information on the type of elastic constitutive law that governs the capsule wall material together with the value of the elastic parameters. The method is applied to a population of ovalbumin microcapsules and is validated by means of independent experiments of the same capsules subjected to a different flow in a microrheological device. This is of great interest for quality control purposes, as small samples of capsule suspensions can be diverted to a measuring test section and mechanically tested with a 10 % precision using an automated process.

Mathematical models promise new insights into the mechanisms underlying the emergence of collective behaviour in fish. Here, we establish a mathematical model to examine collective behaviour of zebrafish, a popular animal species in preclinical research. The model accounts for social and hydrodynamic interactions between individuals, along with the burst-and-coast swimming style of zebrafish. Each fish is described as a system of coupled stochastic differential equations, which govern the time evolution of their speed and turn rate. Model parameters are calibrated using experimental observations of zebrafish pairs swimming in a shallow water tank. The model successfully captures the main features of the collective response of the animals, by predicting their preference to swim in-line, with one fish leading and the other trailing. During in-line swimming, the animals share the same orientation and keep a distance from each other, owing to hydrodynamic repulsion. Hydrodynamic interaction is also responsible for an increase in the speed of the pair swimming in-line. By linearizing the equations of motion, we demonstrate local stability of in-line swimming to small perturbations for a wide range of model parameters. Mathematically backed results presented herein support the application of dynamical systems theory to unveil the inner workings of fish collective behaviour.

An on-demand painting system with a simple structure device that ejects highly viscous liquids as microjets is introduced. An impulsive motion of the container results in the ejection of a viscous liquid jet from the nozzle. This system enabled us to paint letters on a section of a car body using commercial car paint with a zero-shear viscosity of 100 $\textrm {Pa} \cdot \textrm {s}$. To understand the jet velocity, we conducted systematic experiments. Experimental results showed that the jet velocity increases with the ratio between the liquid depths in the container and the nozzle, up to approximately 30 times faster than the initial velocity. However, a linear relation between the jet velocity and the ratio predicted by the previous model, which considers only the pressure impulse, does not hold for the high length ratios since the actual position of the stagnation point is different from the position predicted by the previous model. By solving the Laplace equation and using the model proposed by Gordillo et al. (J. Fluid Mech., vol. 894, 2020, pp. A3–11), we reproduce the non-monotonic behaviour of the jet velocity as a function of the length ratio. For practical use, we improve the jet-velocity model by considering mass conservation as well as the pressure impulse.

In recent years, observations of the atmospheric surface layer have greatly promoted research on high-Reynolds-number wall-bounded turbulence, especially observations of wind-blown sand flows/sandstorms, which are typical sand-laden two-phase flows; these successes have advanced the science of gas–solid two-phase wall-bounded turbulence to very-high-Reynolds-number conditions. Based on a review of existing atmospheric surface layer observations and the development process, this paper summarizes the important promoting effect played by these observations in understanding the very-large-scale structure characteristics, turbulent kinetic energy fraction and amplitude modulation effect, and in reconstructing the spatial electric field under high-Reynolds-number wall turbulence. This review focuses on the main successes achieved by the observation of sand-laden two-phase flows and the three-dimensional turbulent flow field, especially in the streamwise direction. Finally, some suggestions and outlooks for further research on particle-laden two-phase wall-bounded turbulence under high-Reynolds-number conditions are presented.

This study aims to leverage the relationship between fluid dynamic loading and resulting structural deformation to infer the incident flow speed from measurements of time-dependent structure kinematics. Wind tunnel studies are performed on cantilevered cylinders and trees. Tip deflections of the wind-loaded structures are captured in time series data, and a physical model of the relationship between force and deflection is applied to calculate the instantaneous wind speed normalized with respect to a known reference wind speed. Wind speeds inferred from visual measurements showed consistent agreement with ground truth anemometer measurements for different cylinder and tree configurations. These results suggest an approach for non-intrusive, quantitative flow velocimetry that eliminates the need to directly visualize or instrument the flow itself.

Vertical axis turbine (VAT) arrays can achieve larger power generation per land area than their horizontal axis counterparts, due to the positive synergy from clustering VATs in close proximity. The VATs generate a three-dimensional wake that evolves unevenly over the vertical and transverse directions according to two governing length scales, namely the rotor's diameter and height. Theoretical wake models need to capture such a complex wake dynamics to enable reliable array design that maximises energy output. This paper presents two new theoretical VAT wake models based on super-Gaussian and Gaussian shape functions, which account for the three-dimensional velocity deficit distribution in the wake. The super-Gaussian model represents the initial elliptical shape with the superposition of vertical and lateral shape functions that progressively converge into an axisymmetric circular-shaped wake at a downstream distance that depends on the rotor's height-to-diameter aspect ratio. Our Gaussian model improves the initial wake width prediction taking into account the rectangular rotor's cross-section. Our models were well validated with large-eddy simulations (LES) of single VATs with varying aspect ratios and thrust coefficients operating in an atmospheric boundary layer. The super-Gaussian model attained a good agreement with LES in both near and far wake, whilst the Gaussian model represented well the far-wake region.

The study of fluid flow has enabled milestones in technology and a deeper understanding of the natural world. Existing journals cover fundamental aspects of fluid mechanics very well and other journals target specific application fields. We here introduce a new journal titled Flow that is dedicated to covering and highlighting the application of fluid mechanics to concrete problems across all fields. Flow benefits from the roster of reviewers of its sister publication, the Journal of Fluid Mechanics, and offers rapid reviews and an open-access format. For its readers, Flow offers accessible introductions to new application areas and an introduction to the varied tools of fluid mechanics. For authors, Flow offers a venue to reach a broad audience, a focus on enabling translational research, and a way to disseminate the tools of fluid mechanics.

Current methods for fabricating lenses rely on mechanical processing of the lens or mould, such as grinding, machining and polishing. The complexity of these fabrication processes and the required specialized equipment prohibit rapid prototyping of optical components. This work presents a simple method, based on free-energy minimization of liquid volumes, which allows us to quickly shape curable liquids into a wide range of spherical and aspherical optical components, without the need for any mechanical processing. After the desired shape is obtained, the liquid can be cured to produce a solid object with nanometric surface quality. We provide a theoretical model that accurately predicts the shape of the optical components, and demonstrate rapid fabrication of all types of spherical lenses (convex, concave, meniscus), cylindrical lenses, bifocal lenses, toroidal lenses, doublet lenses and aspheric lenses. The method is inexpensive and can be implemented using a variety of curable liquids with different optical and mechanical properties. In addition, the method is scale invariant and can be used to produce even very large optical components, without a significant increase in fabrication time. We believe that the ability to easily and rapidly create optical components, without the need for complex and expensive infrastructure, will provide researchers with new affordable tools for fabricating and testing optical designs.