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.

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.

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.

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.

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.

Atmospheric turbulent velocity fluctuations are known to increase wind turbine structural loading and accelerate wake recovery, but the impact of vortical coherent structures in the atmosphere on wind turbines has not yet been evaluated. The current study uses flow imaging with natural snowfall with a field of view spanning the inflow and near wake. Vortical coherent structures with diameters of the order of 1 m are identified and characterized in the flow approaching a 2.5 MW wind turbine in the region spanning the bottom blade tip elevation to hub height. Their impact on turbine structural loading, power generation and wake behaviour are evaluated. Long coherent structure packets $(\mathrm{\ \mathbin{\lower.3ex\hbox{$\buildrel> \over {\smash{\scriptstyle\sim}\vphantom{_x}}$}}\ }200\;\textrm{m)}$ are shown to increase fluctuating stresses on the turbine support tower. Large inflow vortices interact with the turbine blades, leading to deviations from the expected power generation. The sign of these deviations is related to the rotation direction of the vortices, with rotation in the same direction as the circulation on the blades leading to periods of power surplus, and the opposite rotation causing power deficit. Periods of power deficit coincide with wake contraction events. These findings highlight the importance of considering coherent structure properties when making turbine design and siting decisions.

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.

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.

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.

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.

We conduct a well-controlled model experiment for a wide variety of canopy flows. Examples of these include engineering flows such as wind flow, dispersion of scalars through and over urban areas, and the convective heat transfer in many heat exchangers, as well as natural canopies such as flows through terrestrial or aquatic vegetation. We aim to shed the light on fundamental flow and transport phenomena common to these applications. Specifically, the characteristics of mean flow and scalar concentration characteristics of a turbulent boundary layer flow impinging on a canopy, which comprises a cluster of tall obstacles (this can also be interpreted as a porous obstruction). The cluster is created with a group of cylinders of diameter $d$ and height $h$ arranged in a circular patch of diameter $D$. The solidity of the patch/obstruction is defined by $\phi$ (the total planar area covered by cylinders), which is systematically varied ($0.098 \leq \phi \leq 1$) by increasing the number of cylinders in a patch ($N_c$). A point source is placed at ground level upstream of the patch and its transport over and around the patch is examined. Time-averaged velocity and scalar fields, obtained from simultaneous planar particle image velocimetry-planar laser-induced fluorescence (PIV-PLIF) measurements, reveal that the characteristics of wake and flow above porous patches are heavily influenced by $\phi$. In particular, we observe that the horizontal and vertical extent of the wake and scalar concentration downstream of the patches decreases and increases with $\phi$, respectively. Here, the recirculation bubble is shifted closer to the trailing edge (TE) of the patches as $\phi$ increases, limiting the flow from convecting downstream, decreasing the scalar concentration and virtually ‘extending’ the patch in the streamwise direction. As the bubble forms in the TE, vertical bleeding increases and hence the concentration increases above the patch where the cylinders appear to ‘extend’ vertically towards the freestream.

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.

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.

Redox flow batteries (RFBs) are an emerging electrochemical technology envisioned towards storage of renewable energy. A promising sub-class of RFBs utilizes single-flow membraneless architectures in an effort to minimize system cost and complexity. To support multiple functions, including reactant separation and fast reactant transport to electrode surfaces, electrolyte flow must be carefully designed and optimized. In this work, we propose adding a secondary channel adjacent to a permeable battery electrode, solving for the flow field and analysing the effects on the reactant concentration boundary layer at the electrode. We find that an adjacent channel with gradually changing thickness leads to a desired nearly uniform flow through the electrode to the adjacent channel. Consequently, the thickness of the concentration boundary layer is significantly reduced, increasing reactant transport to the electrode surface to 140% of the rate of a battery with a constant width adjacent channel, and 350% of the rate with no adjacent channel. Overall, this theory provides insight into the important role of flow physics for this promising sub-class of flow batteries, and can pave the way to improved energy efficiency of such flow batteries.

The flow patterns in the heart, in health and disease, have been of great interest for several years. Modern fluid dynamics analyses elucidate how underlying inefficient energetic or mixing characteristics of these flow patterns correlate with adverse effects. Unfortunately, translation of such modern analyses to the clinical stage remains a challenge. In this experimental work, we propose and demonstrate that braids of random and sparse particle trajectories provide an intuitive, global and practical description of cardiovascular flows. Moreover, we expose the flow pattern in an experimental healthy left ventricle model as a highly effective blood mixer at the topological level. Flow topologies that deviate from this pattern are accompanied by a reduction in energetic efficiency, as shown through comparisons with diseased flow models. These results suggest an ideal clinical approach to patient follow-up and the evaluation of the performance of medical devices.

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.

The SARS-CoV-2 is transmitted not only through coughing, but also through breathing, speaking or singing. We perform direct numerical simulations of the turbulent transport of potentially infectious aerosols in short conversations, involving repetitive phrases separated by quiescent intervals. We estimate that buoyancy effects due to droplet evaporation are small, and neglect them. A two-way conversation is shown to significantly reduce the aerosol exposure compared with a relative monologue by one person and relative silence of the other. This is because of the ‘cancelling’ effect produced by the two interacting speech jets. Unequal conversation is shown to significantly increase the infection risk to the person who talks less. Interestingly, a small height difference is worse for infection spread, due to reduced interference between the speech jets, than two faces at the same level. For small axial separation, speech jets show large oscillations and reach the other person intermittently. We suggest a range of lateral separations between two people to minimize transmission risk. A realistic estimate of the infection probability is provided by including exposure through the eyes and mouth, in addition to the more common method of using inhaled virions alone. We expect that our results will provide useful inputs to epidemiological models and to disease management.

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.

Leidenfrost drops were recently found to host strong dynamics. In the present study, we investigate both experimentally and theoretically the flow structures and stability inside a Leidenfrost water drop as it evaporates, starting with a large puddle. As revealed by infrared mapping, the drop base is warmer than its apex by typically 10 $^{\circ }$C, which is likely to trigger bulk thermobuoyant flows and Marangoni surface flows. Tracer particles unveil complex and strong flows that undergo successive symmetry breakings as the drop evaporates. We investigate the linear stability of the base flows in a non-deformable, quasi-static, levitating drop induced by thermobuoyancy and the effective thermocapillary surface stress, using only one adjustable parameter. The stability analysis of nominally axisymmetric thermoconvective flows, parametrized by the drop radius $R$, yields the most unstable, thus, dominant, azimuthal modes (of wavenumber $m$). Our theory predicts well the radii $R$ for the mode transitions and cascade with decreasing wavenumber from $m=3,\, m=2$, down to $m=1$ (the eventual rolling mode that entails propulsion) as the drop shrinks in size. The effect of the escaping vapour is not taken into account here, which may further destabilize the inner flow and couple to the liquid/vapour interface to give rise to motion (Bouillant et al. Nat. Phys., vol. 14 (12), 2018, pp. 1188–1192; Brandão & Schnitzer Physical Review Fluids, vol. 5 (9), 2020, 091601).