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.

In several problems involving fluid flows, computational fluid dynamics (CFD) provides detailed quantitative information and allows the designer to successfully optimize the system by minimizing a cost function. Sometimes, however, one cannot improve the system with CFD alone, because a suitable cost function is not readily available; one notable example is diagnosis in medicine. The application considered here belongs to the field of rhinology; a correct air flow is key for the functioning of the human nose, yet the notion of a functionally normal nose is not available and a cost function cannot be written. An alternative and attractive pathway to diagnosis and surgery planning is offered by data-driven methods. In this work, we consider the machine learning study of nasal impairment caused by anatomic malformations, with the aim of understanding whether fluid dynamic features, available after a CFD analysis, are more effective than purely geometric features at the training of a neural network for regression. Our experiments are carried out on an extremely simplified anatomic model and a correspondingly simple CFD approach; nevertheless, they show that flow-based features perform better than geometry-based ones and allow the training of a neural network with fewer inputs, a crucial advantage in fields like medicine.

Constriction in the flow passage of the physiological circulatory system is central to the occurrence of several diseased conditions such as thrombosis and is also pivotal towards the understanding of several regulatory processes in the human microvasculature. It is, therefore, imperative to advance a mechanistic insight into the dynamics of lipid vesicles, cell mimicking, fluid-filled compartments, through a physiologically relevant microconfinement, with particular focus on deciphering the role of its mechanophysical properties. Here we bring out the role of membrane bending rigidity and the initial deflation (deviation in shape from sphericity) on the transient shape evolution of a lipid vesicle as it migrates through a microfluidic constriction, a paradigm that is unexplored thus far. Based on our experimental observations as well as theoretical insights, we construct a regime map to elucidate the range of the key dimensionless parameters orchestrating this dynamical transition. Furthermore, our observations on the vesicle's stretching dynamics emerging from selective mapping with viscosity contrast between the encapsulated and the suspending fluid medium offer potential physiologically relevant cues on the impact of cell aging on its deformability across a constricted path. Such mechanistic insights may help in establishing quantitative correlations between the dynamical transition of a lipid vesicle and its membrane mechanics, which may in turn have decisive implications in health and disease while circulating across microvascular fluidic pathways.

This paper investigates the impact of large-scale turbulence on the aerodynamic characteristics of a pitching wind turbine blade at Reynolds number 135 000, whose cross-section is a NACA0012 airfoil with constant chord length. Large-eddy simulations at reduced frequencies, $k_{red} = 0.05$ and 0.1, were validated against reference data from the literature. An efficient method capable of generating synthetic large-scale turbulence at the inlet was applied by using two streamwise integral length scales $L_{x} = 1c$ and 1.5$c$, which represent energetic turbulence eddies at the height where the wind turbine operates. For $k_{red}= 0.1$, the change in the maximum lift coefficient at the dynamic stall angle near the maximum angle of attack is on average 20 % lower and during the downstroke it is on average 22 % lower, compared with the smooth inflow. A higher reduced frequency ($k_{red} = 0.2$) apparently does not further change the lift, drag and moment coefficients, and the inflow turbulence disordered leading-edge vortices. The turbulent shear stress and the phase-averaged dispersive shear stress in the wake are of the same magnitude, but with negative and positive signs, respectively, suggesting that the large-scale phase-averaged fluctuations transfer momentum in the opposite direction compared with the turbulent fluctuations, reducing the drag on the suction-side flow, and subsequently increasing the averaged lift coefficient. This demonstrates the critical importance of the large-scale non-turbulent unsteadiness in the wake of a pitching wing.

In the field of experimental fluid dynamics, the direct measurement of vorticity remains a challenge, even though it plays a crucial role in understanding turbulent flows. The present study explores the influence of the rotation of nanoparticles on their luminescence anisotropy as a potential novel measurement method. This relation opens a new field of flow diagnostics, based on the measurement of polarized intensity components. Potentially, the method allows for the direct measurement of the vorticity. For this, the canonical flow in this study is a turbulent round jet at ${{Re}} = {12\,000}$ and 14 400. It is confirmed that the flow regime has an influence on the luminescence anisotropy. Using a model of such deterministic rotations according to another work by the authors (Schmidt & Rösgen, Phys. Rev. Res., vol. 5, no. 3, 2023, 033006), the magnitude of the vorticity components is computed, since the presented set-up is limited to sensing the magnitude of these quantities. The computed components indicate the self-similarity of the vorticity magnitude. A large-eddy simulation is conducted for comparison with the experiments, demonstrating good agreement.

Gravure printing is a type of printing method that uses metal cylinders with engraved cells that hold ink. The ink is transferred directly to the paper or other material by pressing it against the cylinder. The flow associated with gravure printing includes a flow in a liquid bridge formed in the contact region of the cylinders and a thin-film coating flow of the ink. The flow is governed by viscous and capillary forces. In many cases, the flow is unstable, which leads to the formation of instability patterns on the printed surfaces. The analysis of these instabilities is a very challenging problem, especially since industrial inks are usually rheologically complex. In this experimental and theoretical study, the flow of inks on a rotating cylinder is analysed, accounting for the shear-rate-dependent liquid viscosity. A theoretical solution for the film flow allows us to predict the width of the liquid bridge between two cylinders. Moreover, it is shown that the measured characteristic size of the printed pattern is of the same order as the predicted liquid bridge width. We observe a nearly linear dependence of pattern size and liquid bridge width.

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.

Large eddy simulations are performed to investigate the impact of a solid obstacle on the flow around a multiperforated plate typical of aeronautical combustion chambers. The reference configuration is a perforated plate with approximately 200 holes immersed between a cold vein and hot vein at a typical operating point of helicopter combustors. The micro-jet Reynolds number is of the order of 4000, while the blowing and momentum ratios are close to 4 and 8, respectively. A variant configuration is considered that features an additional cylindrical obstacle located in the cold vein and mimicking a spark plug. The study reveals that, downstream of the obstacle, the cooling effectiveness of the plate is reduced by approximately 40 % compared with the reference case, mainly due to the absence of perforation at the obstacle location. The mass flow rate within the holes in the wake produced by the obstacle is reduced by 7 %, which is likely to locally influence the plate cooling. The reduction is attributed to the wake's pressure loss and its impact on the discharge coefficient. Additionally, the cooling effectiveness outside the wake shows a 5 % increase that can be linked to the mass flow rate increase within corresponding holes.

The receptivity of a laminar boundary-layer flow to small-amplitude wall roughness is investigated on an ONERA-D swept aerofoil by introducing a dedicated transfer function from small-amplitude wall displacements to full-state velocity perturbations. The singular value decomposition of this operator for a given spanwise wavenumber provides optimal wall roughness and flow responses that maximise an input–output gain. At the most receptive spanwise wavenumber, the optimal response is a cross-flow mode associated with an optimal roughness located close to the attachment line and presenting a wavy shape with a wavevector nearly orthogonal to the external streamlines. The method therefore allows direct identification of the location and structure (chordwise and spanwise wavenumbers) of the most receptive roughness. For various given wall roughness shapes and locations (periodic or compact in the chordwise and/or spanwise directions), an approximation of the response based on the dominant optimal response is shown to accurately match the total response downstream of the roughness. The method therefore allows a straightforward computation of the response of the flow to any given small-amplitude roughness.

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.

Supersonic rectangular twin jets (SRTJ) are of interest for current and future generations of tactical aircraft. However, the adverse effects of screech-loop coupling have previously been documented to significantly increase the near-field (NF) pressure fluctuations. These high-amplitude pressure fluctuations have the potential to fatigue and damage nearby aft components of the aircraft. Previous studies have documented that the NF pressure fluctuation level depends on the coupling of the jets: in-phase coupling along the twin jets’ minor axes produces stronger NF pressure fluctuations than that of out-of-phase coupling. The objective of this work is to further investigate the effects of coupling modes on NF pressure fluctuations in SRTJ and to mitigate the adverse effects of coupling using active flow control. Localized arc filament plasma actuators are employed to alter the SRTJs’ coupling mode by leveraging natural flow instabilities with minimal power input. A NF microphone array is used for simultaneous coupling and NF pressure fluctuation measurements. Schlieren imaging and spectral proper orthogonal decomposition are used to assess the effects of control on the flow field. The effect of excitation at various frequencies and actuation patterns on coupling and NF pressure fluctuations in different flow regimes are explored and discussed.

Shallowly submerged oscillating structures may be found in wave energy devices or semi-submersible vessels. Predicting the force on such structures is critical for design purposes, but complicated due to nonlinear phenomena which can occur in shallow water, including wave breaking and bore formation. Such effects are particularly important around the first ‘resonance’ frequency of the fluid on top of the device, where linear theory predicts large flows on/off the cylinder and corresponding surface elevations and forces. In an effort to create a reliable and efficient model to predict the hydrodynamic force on a shallowly submerged truncated vertical cylinder, an axisymmetric nonlinear hybrid model is developed for forced heave oscillations. The flow above the cylinder is modelled using the nonlinear shallow water equations, and linear potential flow theory is used in the surrounding fluid. The model is compared with experimental results for forced heave oscillations and performs well for predicting the heave force. It is then used to examine linearised heave force for increasing amplitudes of (prescribed) harmonic heave motion. There is a significant reduction in the peaks of radiation damping and added mass coefficients with increasing amplitude, and associated shifts in the frequencies of the peaks.

Marine debris pollution is a growing problem impacting aquatic ecosystems, coastal recreation and human society. Beaches are known to be a sink for debris, and beaching needs to be accounted for in marine debris mass balances, but the process of buoyant debris beaching is not yet sufficiently well understood in order for it to be included in coastal models. We develop a simplified model for buoyant marine debris transport in bore-driven swash (where swash refers to the area that the water wets the beach with each incoming wave). We validate the model with laboratory experiments and use the combined results from the model and experiments to understand the parameters that are important for dictating particle beaching. The most relevant parameters are the particle inertia and the initial conditions with which debris particles enter the swash zone.

Martian dust storms in the planetary boundary layer share many qualitative similarities to terrestrial sandstorms. Both of these turbulence-driven, particle-laden boundary layer flows are known to generate electric fields due to the transport of differentially charged particles; this charge separation can be strong enough to lead to dielectric breakdown in the form of sparks or lightning. Using wall-modelled large-eddy simulations supplemented with conservation of equations for the charged particle transport, representative simulations of neutrally stable Martian and terrestrial particle-laden boundary layer flows are compared. The simulations, albeit canonical in nature, provide evidence to support previous observations of the less frequent occurrence of lightning on Mars but a higher occurrence of localised electric discharge events due to the much lower breakdown potential. The rarefied Martian atmosphere impedes charged particle transport, resulting in a weaker electric field than the equivalent terrestrial sandstorm. The lower drag force in the rarefied Martian atmosphere means that the electrostatic force plays a more significant role in the particle transport, which results in a self-regulation of the electric field. The strongest Martian dust storms show evidence of significant breakdown events and these discharge events only occur very close to the ground despite the very large boundary layer on Mars.

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.

Proper orthogonal decomposition (POD) is applied to three-dimensional (3-D) velocity fields collected from large-eddy simulations (LES) of a baffled stirred tank. In the LES, the tank operates with a Rushton-type impeller under turbulent conditions (at least in the near-impeller region) and the working fluid exhibits either Newtonian or shear-thinning rheology. The most energetic POD modes are analysed, and a POD reconstruction based on the higher modes is proposed to approximate the fluctuating component of the velocity field. Subsequently, the POD reconstruction is used to identify vortical structures and characterise them in terms of their shape. The structures are identified by considering a frame-invariant formulation of a popular, Eulerian, local-region-type method: the $Q$-criterion. Statistics of shape-related parameters are then investigated to address the morphology of the structures. It is found that: (i) regardless of the working fluid rheology, it seems feasible to decompose the 3-D field into its mean, most energetic periodic and fluctuating components using POD, allowing, for instance, reduced-order modelling of the energetic periodic motions for mixing enhancement purposes, and (ii) vortical structures related to turbulence are mostly tubular. Finding (ii) implies that, as starting point, phenomenological models for the interaction between fluid particles (drops and bubbles) and vortices should consider the latter as cylindrical structures rather than of spherical shape, as classically assumed in these models.

Fish often swim in crystallized group formations (schooling) and orient themselves against the incoming flow (rheotaxis). At the intersection of these two phenomena, we investigate the emergence of unique schooling patterns through passive hydrodynamic mechanisms in a fish pair, the simplest subsystem of a school. First, we develop a fluid dynamics-based mathematical model for the positions and orientations of two fish swimming against a flow in an infinite channel, modelling the effect of the self-propelling motion of each fish as a point-dipole. The resulting system of equations is studied to gain an understanding of the properties of the dynamical system, its equilibria and their stability. The system is found to have five types of equilibria, out of which only upstream swimming in in-line and staggered formations can be stable. A stable near-wall configuration is observed only in limiting cases. It is shown that the stability of these equilibria depends on the flow curvature and streamwise interfish distance, below critical values of which, the system may not have a stable equilibrium. The study reveals that simply through passive fluid dynamics, in the absence of any other feedback/sensing, we can justify rheotaxis and the existence of stable in-line and staggered schooling configurations.

The unsteady flow behaviour of two side-by-side rotors in ground proximity is experimentally investigated. The rotors induce a velocity distribution interacting with the ground causing the radial expansion of the rotor wakes. In between the rotors, an interaction of the two wakes takes place, resulting in an upward flow similar to a fountain. Two types of flow topologies are examined and correspond to two different stand-off heights between the rotors and the ground: the first one where the height of the fountain remains below the rotor disks, and a second one where it emerges above, being re-ingested. The fountain unsteadiness is shown to increase when re-ingestion takes place, determining a location switch from one rotor disk to the other, multiple times during acquisition. Consequently, variable inflow conditions are imposed on each of the two rotors. The fountain dynamics is observed at a frequency that is about two orders of magnitude lower than the blade passing frequency. The dominant characteristic time scale is linked to the flow recirculation path, relating this to system parameters of thrust and ground stand-off height. The flow field is analysed using proper orthogonal decomposition, in which coupled modes are identified. Results from the modal analysis are used to formulate a simple dynamic flow model of the re-ingestion switching cycle.