The placement of a scaled-down Savonius (drag) vertical-axis wind turbine on model buildings is analysed experimentally by the use of turbine performance and flow field measurements in a wind tunnel. The set-up consists of two surface mounted cubes aligned in the flow direction. The turbine is tested at six different streamwise positions – three on each cube. Velocity field measurements are performed with particle image velocimetry along the centreline of the cubes with and without the turbine. The performance at each position is evaluated based on measurements of the produced torque and the rotational speed of the turbine. It is demonstrated that the common practice of estimating wind resources based on the urban flow field without the turbine present is insufficient. The turbine has a substantial influence on the flow field and thus also on the available power. The performance is found to be optimal in the front and centre of the first building with a significant drop-off to the back. This trend is reversed for the downstream building. Holistically, for more generic geometries and varying wind directions, the results suggest the central position on a building is a good compromise.

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.

Effects of gaps (rectangular surface cavities) on boundary-layer transition are investigated using a combination of linear stability theory and experiments, for boundary layers where the smooth-surface transition results from Tollmien–Schlichting (TS) instability. Results are presented for a wide range of gap characteristics, with the associated transition locations ranging from the smooth-surface location all the way forward to the gap location. The transition movement is well described by a variable $N$-factor, which links the gap characteristics to the level of instability amplification $e^N$ leading to transition. The gap effects on TS-wave transition are characterized by two limiting behaviours. For shallow gaps $d/w < 0.017$, the reduction in $N$-factor is a function of the gap depth $d$ and is independent of the gap width $w$. For deep gaps $d/w > 0.028$, the reduction in $N$-factor is a function of the gap width and is independent of the gap depth. When both the gap width and depth are sufficiently large relative to the displacement thickness $\delta ^*$, the TS-wave transition is bypassed, resulting in transition at the gap location. These behaviours are mapped out in terms of ($w/ \delta ^*$, $d/ \delta ^*$), providing a predictive model for gap effects on transition.

To mitigate the risk of wake vortex encounters during final approach, so-called plate lines have been developed. Data collected during a six-month measurement campaign at Vienna International Airport are used to assess the potential for reducing minimum aircraft separations facilitated by plate lines during approach and landing following the re-categorisation (RECAT-EU) methodology for revised wake turbulence categorisation. To ensure that no other parameters controlling wake vortex decay bias the analysis, it is verified that wind speed, atmospheric turbulence, thermal stratification and flight altitude reside in similar ranges with and without the plates. The analysis follows the steps of the RECAT-EU method to generate non-dimensional so-called reasonable worst-case circulation decay curves; one as a reference for nominal operations without plates and one representing the accelerated wake vortex decay brought about by the plate lines. The difference between the two circulation decay curves corresponds to the non-dimensional time-based aircraft separation reduction potential that can be translated into distance-based separation gains. Depending on the particular RECAT-EU category combination, the attained aircraft separation reduction potential ranges between 12% and 15%. Constricting the analysis to wake vortices generated by one aircraft type representing the Upper Heavy RECAT-EU category, the separation reduction potential amounts up to 24%.

When manufacturing a turbine engine, the combustor annulus and the turbine annulus are created as separate parts and assembled. This leads to an inter-platform gap between the two components, which must be supplied with leakage air to prevent ingestion of the extremely hot combustion gases into the interior of the engine. The combustor and the turbine are likely to misalign because of differential thermal expansion or assembly tolerances. This paper presents a direct numerical simulation study of inter-platform misalignment with leakage flow supplied at the junction of the platforms. The geometry is two misaligned plates with a cross-flow and a leakage flow simulated as a slot jet. The misalignment of the two plates gives rise to a forward misalignment configuration and a backward misalignment configuration, and the jet/cross-flow gives rise to a windward mixing layer and a leeward mixing layer. Compared with the aligned configuration, the cooling effectiveness immediately downstream of the gap decreases in the forward misalignment configuration and increases in the backward misalignment configuration; this response amplifies as the flow rate through the gap increases. In addition to the cooling effectiveness, we report flow statistics, including the velocity, the temperature, the turbulent kinetic energy and the relevant turbulent fluxes. We find strong turbulence generation in the leeward mixing layer and high turbulence level as a result. Mixing of the thermal energy, on the other hand, occurs predominantly in the windward mixing layer. The eddy viscosity and the eddy conductivity that are critical to turbulence modelling are also reported. We find negative eddy viscosity at regions where the incoming boundary layer starts to mix with the leakage jet. The analysis shows that the negative eddy viscosity is a result of flow hysteresis: it takes time, or travel distance, before the eddies in the incoming boundary layer and the eddies in the leakage jet come to an equilibrium, thereby favouring a transport Reynolds stress model over a local eddy viscosity type model. The novelty of this paper lies in the direct numerical simulations, which provide direct access to the near-wall flow field and clarify the effects of blowing ratio and platform misalignment on heat transfer. The novelty also lies in the data analysis, which sheds light on how this flow should be modelled.

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.

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.