This article presents a broad overview of the fluid mechanics of the human cardiovascular system. Beginning in the heart, we travel through the main features of our circulation to highlight important functions and diseases where fluid mechanics plays a central role. Of particular focus is the role of computational modelling in uncovering the dynamic flow phenomenon throughout our body, its association with cardiovascular disease mechanisms and progression and its importance in clinical treatment planning. We also emphasize the multiscale nature of the cardiovascular system, and associated challenges. The main aim of this review is to highlight progress and ongoing challenges in our understanding of cardiovascular haemodynamics, as well as the future outlook for translating the current state-of-the-art to widespread clinical application and improved patient outcomes.

This paper reviews 70 years of cricket ball swing literature and reassesses the results in light of new measurements. A comparison of ball tracking data with published experimental results shows that current understanding does not explain the behaviour observed in professional cricket matches. Descriptions of cricket ball boundary layer aerodynamics are updated with new results which show that the seam acts more like a series of vortex generators than a trip and that there is a laminar separation bubble on the seam side of the ball. Previous results for reverse swing are consolidated and compared with new data, showing that different magnitudes of swing occur depending on the condition of both sides of the ball. Variation in pressure and temperature should be included when the Reynolds number, a non-dimensional ball speed, is calculated, while humidity can be neglected. Studies on the effect of wind speed and direction are summarised and results considering the effect of free-stream turbulence are compared with new measurements. Throughout the paper, recommendations for future work are suggested. These include quantitative study of the effect of surface defects and roughness, an assessment of whether atmospheric turbulence can affect swing and investigation into the effect of backspin on swing.

The axially symmetric, swirl-free gas dynamics and interlinked motion of a cork stopper provoked by the opening of a champagne bottle are modelled rigorously and studied numerically. The experimental study by Liger-Belair et al. (Science Advances, 5(9), 2019) animated the present investigation. Inspection analysis justifies the inviscid treatment of the expanding jet of air enriched with dissolved carbonic acid gas initially pressurised in the bottle. Solving of the resulting Euler equations is facilitated by the open-source software Clawpack. Specific enhancements allow for resolving of the emerging supersonic pockets, associated with surprisingly complex shock structures, as well as the gas–stopper interaction with due accuracy. Our experimental effort provided modelling of the frictional behaviour, constitutive law and reversible (de-)compression of the cork material. Initially, the gas expands inside the bottleneck yet sealed by the stopper, and is hence accelerated by the gas but decelerated by dry sliding friction. Once the stopper has passed the bottle opening, the jet rapidly assumes locally supersonic speed, where a complex shock pattern is detected. Special attention is paid to the formation and dissolution of one or even two Mach discs between the opening and the released stopper. This simulated dynamics is found to be in fairly good agreement with recent experimental findings. It also provides a first insight into the generation of the typical popping sound.

Over the years, many facilities have been developed to study turbulent flow in the laboratory. Homogeneous isotropic turbulence (HIT) with zero mean flow provides a unique environment for investigating fundamental aspects and specific applications of turbulent flow. We provide an extensive overview of laboratory facilities that generate incompressible zero-mean-flow HIT using different types of actuators and configurations. Reviewed facilities cover a variety of geometries and sizes, as well as forcing style (e.g. symmetric versus asymmetric and unsteady versus steady). We divide facilities into four categories, highlighting links between their geometries and the statistics of the flows they generate. We then compare published data to uncover similarities and differences among various turbulence-generation mechanisms. We also compare the decay of turbulence in zero-mean-flow facilities with that observed in wind and water tunnels, and we analyse the connections between flow characteristics and physical aspects of the facilities. Our results emphasize the importance of considering facility geometry and size together with the strength and type of actuators when studying zero-mean-flow HIT. Overall, we provide insight into how to optimally design and build laboratory facilities that generate zero-mean-flow HIT.

The flying wing is an aerodynamic configuration with high efficiency, but the lack of lateral-directional stability has always been an obstacle that limits its application. In this study, the wing rock motion of a 65° swept flying-wing aircraft is studied via wind tunnel experiments and numerical simulations at a low speed, and various unsteady motion phenomena are focused on. Both the experimental and numerical results show that the flying wing has a bicyclic ${C_l}$–$\phi $ hysteresis loop during its wing rock, different from the slender delta wing, rectangular wing, generic aircraft configuration, etc., which have a tricyclic hysteresis loop. This form of hysteresis loop implies a different energy exchange manner of the flying wing in the wing rock oscillation. Further analysis shows that the flying wing forms a unilateral leading-edge vortex (LEV) under a high roll angle, with its wing rock oscillation driven by the ‘vortex–shear-layer’ structure, which is different from that of slender and non-slender delta wings. Moreover, the quantitative dynamic hysteresis characteristics of the LEV's strength and location for the flying wing and the slender delta wing are also different. These results have proven the existence of a wing rock mode which is different from previous investigations, which enriches the understanding of self-induced oscillation. Present discoveries are also conducive to the aerodynamic shape design and flight manipulation of a flying-wing aircraft, which is significant for its wider application.

We review the progress on the applications of the vortex-surface field (VSF). The VSF isosurface is a vortex surface consisting of vortex lines. Based on the generalized Helmholtz theorem, the VSF isosurfaces of the same threshold at different times have strong coherence. As a general flow diagnostic tool for studying vortex evolution, the numerical VSF solution is first constructed in a given flow field by solving a pseudo-transport equation driven by the instantaneous frozen vorticity, and then the VSF evolution is calculated by the two-time method. From the database of numerical simulations or experiments, the VSF can elucidate mechanisms in the flows with essential vortex dynamics, such as isotropic turbulence, wall flow transition, flow past a flapping plate and turbulence–flame interaction. The characterization of VSFs reveals the correlation between robust statistical features and the critical quantities needed to be predicted in engineering applications, such as the friction coefficient in transition, thrust in bio-propulsion and growth rate in interface instability. Since the VSF evolution captures the essential Lagrangian-based dynamics of vortical flows, it inspires novel numerical methods on cutting-edge hardware, e.g. graphic and quantum processors.

Superhydrophobic surfaces (SHSs) are considered to be a promising technology for achieving skin-friction drag reduction. Development of more efficient techniques for simulating the turbulent boundary layer on SHSs continues to be a subject of interest. In this study, numerical simulations were carried out to capture near-wall behaviours due to the effect of the SHS on wall-bounded flows. To achieve this, high- to intermediate-fidelity turbulence models including Reynolds-averaged Navier–Stokes, detached eddy simulation and large eddy simulation were utilized. With regard to slip conditions, the well-known Navier slip velocity method was used over the SHS. For validating the numerical solutions, the slip velocity and skin friction over the SHS were compared with the experimental output. Results showed that the velocity profile and Reynolds stresses on the SHS were comparable to the reported results. Then, the developed models were further extended to investigate the drag reduction effect of SHSs with rectangular grooves. The subsequent results showed that the combination of superhydrophobicity and rectangular grooves led to a better performance with a maximum drag reduction of 46.1%. This is due to the surface slip caused by the SHS and the secondary vortex effect created by the grooves. Our results revealed that Reynolds stresses of the slippery grooved surface were higher than those of the case in which a shear-free condition was employed for the grooved surface. More importantly, the numerical results indicate the previous assumption of the shear-free condition is inaccurate for the geometrically simplified grooved SHSs. Therefore, geometry modifications rather than an overly simplified shear-free boundary condition should be applied in computational fluid dynamics simulations for SHSs with grooves or other complex structures.

Cohesive particulate flows play an important role in environmental fluid dynamics, as well as in a wide variety of civil and process engineering applications. However, the scaling laws, constitutive equations and continuum field descriptions governing such flows are currently less well understood than for their non-cohesive counterparts. Grain-resolved simulations represent an essential tool for addressing this shortcoming, along with theoretical investigations, laboratory experiments and field studies. Here we provide a tutorial introduction to simulations of fine-grained sediments in viscous fluids, along with an overview of some representative insights that this approach has yielded to date. After a brief review of the key physical concepts governing van der Waals forces as the main cohesive effect for subaqueous sediment suspensions, we discuss their incorporation into particle-resolved simulations based on the immersed boundary method. We subsequently describe simulations of cohesive particles in several model turbulent flows, which demonstrate the emergence of a statistical equilibrium between the growth and break-up of aggregates. As a next step, we review the influence of cohesive forces on the settling behaviour of dense suspensions, before moving on to submerged granular collapses. Throughout the article, we highlight open research questions in the field of cohesive particulate flows whose investigation may benefit from grain-resolved simulations.

The present paper simplifies the naturally formed dunes (riverbeds) as large-scale three-dimensional staggered wavy walls to investigate the features of the accompanying secondary flows and streamwise vortices via large-eddy simulation. A comparison between the swirling strength and the mean velocities suggests where a secondary flow induces upwash or downwash motions. Moreover, we propose a pseudo-convex wall mechanism to interpret the directionality of the secondary flow. The centrifugal instability criterion is then used to reveal the generation of the streamwise vortices. Based on these analytical results, we found that the streamwise vortices are generated in the separation and reattachment points on both characteristic longitudinal–vertical and horizontal cross-sections, which is related to the curvature effect of the turbulent shear layer. Furthermore, the maximum Görtler number characterized by the ratio of centrifugal to viscous effects suggests that, for fixed ratio of spanwise- to streamwise-wavelength cases, the strongest centrifugal instability occurring on the longitudinal–vertical cross-section gradually dominates with the increases in amplitude. A similar trend for the cases with varied spanwise wavelength can also be found. It is also found that the streamwise vortices are generated more readily via transverse flow around the crest near the separation and reattachment points when the ratio of spanwise- to streamwise-wavelength equals 1.

The knuckleball is considered to be one of the hardest pitches to hit in baseball due to its seemingly unpredictable motion. It has gained popularity in cricket in recent times. It is shown that the delivery referred to as knuckleball in cricket, at present, does not exhibit a zigzag motion and is, therefore, a misnomer. We propose a delivery in cricket that is associated with an erratic trajectory similar to the knuckleball pitch in baseball. Force measurement experiments in a wind tunnel on a new cricket ball in various orientations of the seam to the incoming flow and at different Reynolds number are carried out. The results are utilized to estimate the trajectory of knuckleball deliveries. The key parameters are the seam angle, speed and spin rate of the ball at the time of its release. Their effect on the trajectory is studied in detail. The optimal combination of these parameters that result in a knuckleball, which is likely hard for the batter to play, is identified.

We performed large-eddy simulations in a neutral atmospheric boundary layer to study the interaction between two identical wind farms with 72 turbines each. We demonstrate that the wind farm wake created by the upstream farm affects the entire flow in and around the downstream farm. The vertical entrainment fluxes above the downstream wind farm are strengthened, resulting in a faster wind farm wake recovery behind the downstream farm. These findings illustrate that interaction between extended wind farms affects flow structures beyond the wind farm scale. Furthermore, we demonstrate that wind farm wakes can reduce the power production of turbines throughout the downstream wind farms. We additionally observe that a staggered wind farm extracts more energy from the flow and thus creates a stronger wind farm wake than an aligned wind farm.

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.

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.

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.