EUROPEAN JOURNAL OF MECHANICS B-FLUIDS

Validations of the microchannel flow model for characterizing vascularized tissues
Poul SS, Ormachea J, Hollenbach SJ and Parker KJ
The microchannel flow model postulates that stress-strain behavior in soft tissues is influenced by the time constants of fluid-filled vessels related to Poiseuille's law. A consequence of this framework is that changes in fluid viscosity and changes in vessel diameter (through vasoconstriction) have a measurable effect on tissue stiffness. These influences are examined through the theory of the microchannel flow model. Then, the effects of viscosity and vasoconstriction are demonstrated in gelatin phantoms and in perfused tissues, respectively. We find good agreement between theory and experiments using both a simple model made from gelatin and from living, perfused, placental tissue .
Computational Modeling of Bubbles Growth Using the Coupled Level Set-Volume of Fluid Method
Taqieddin A, Liu Y, Alshawabkeh AN and Allshouse MR
Understanding the generation, growth, and dynamics of bubbles as they absorb or release dissolved gas in reactive flows is crucial for optimizing the efficiency of electrochemically gas-evolving systems like alkaline water electrolysis or hydrogen production. To better model these bubbly flow systems, we use a coupled level set and volume of fluid approach integrated with a one-fluid transport of species model to study the dynamics of stationary and rising bubbles in reactive two-phase flows. To accomplish this, source terms are incorporated into the continuity and phase conservation equations to allow the bubble to grow or shrink as the species moves through the interface. Verification of the hydrodynamics of the solver for non-reactive systems demonstrates the requisite high fidelity interface capturing and mass conservation necessary to incorporate transport of species. In reactive systems where the species impacts the bubble volume, the model reproduces the theoretically predicted and experimentally measured diffusion-controlled growth rate (i.e., () ∝ ). The model is then applied to rising bubbles to demonstrate the impact of transport of species on both the bubble velocity and shape as well as the concentration field in its wake. This improved model enables the incorporation of electric fields and chemical reactions that are essential for studying the physicochemical hydrodynamics in multiphysics systems.
Shock-induced collapse of a bubble inside a deformable vessel
Coralic V and Colonius T
Shockwave lithotripsy repeatedly focuses shockwaves on kidney stones to induce their fracture, partially through cavitation erosion. A typical side effect of the procedure is hemorrhage, which is potentially the result of the growth and collapse of bubbles inside blood vessels. To identify the mechanisms by which shock-induced collapse could lead to the onset of injury, we study an idealized problem involving a preexisting bubble in a deformable vessel. We utilize a high-order accurate, shock- and interface-capturing, finite-volume scheme and simulate the three-dimensional shock-induced collapse of an air bubble immersed in a cylindrical water column which is embedded in a gelatin/water mixture. The mixture is a soft tissue simulant, 10% gelatin by weight, and is modeled by the stiffened gas equation of state. The bubble dynamics of this model configuration are characterized by the collapse of the bubble and its subsequent jetting in the direction of the propagation of the shockwave. The vessel wall, which is defined by the material interface between the water and gelatin/water mixture, is invaginated by the collapse and distended by the impact of the jet. The present results show that the highest measured pressures and deformations occur when the volumetric confinement of the bubble is strongest, the bubble is nearest the vessel wall and/or the angle of incidence of the shockwave reduces the distance between the jet tip and the nearest vessel surface. For a particular case considered, the 40 MPa shockwave utilized in this study to collapse the bubble generated a vessel wall pressure of almost 450 MPa and produced both an invagination and distention of nearly 50% of the initial vessel radius on a 𝒪(10) ns timescale. These results are indicative of the significant potential of shock-induced collapse to contribute to the injury of blood vessels in shockwave lithotripsy.
On the three-dimensional vortical structure of early diastolic flow in a patient-specific left ventricle
Le TB and Sotiropoulos F
We study the formation of the mitral vortex ring during early diastolic filling in a patient-specific left ventricle using direct numerical simulation. The geometry of the left ventricle is reconstructed from Magnetic Resonance Imaging (MRI). The heart wall motion is modeled by a cell-based activation methodology, which yields physiologic kinematics with heart rate equal to 52 beats per minute. We show that the structure of the mitral vortex ring consists of the main vortex ring and trailing vortex tubes, which originate at the heart wall. The trailing vortex tubes play an important role in exciting twisting circumferential instability modes of the mitral vortex ring. At the end of diastole, the vortex ring impinges on the wall and the intraventricular flow transitions to a weak turbulent state. Our results can be used to help interprete and analyze three-dimensional in-vivo flow measurements obtained with MRI.
Alternating Flow in a Moving Corner
Laine-Pearson FE and Hydon PE
Recent experiments have shown that rapid kinematic mixing occurs in the pulmonary alveoli. Here the Reynolds number is very small, there is recirculation in the alveolar cavity and the alveolar walls move periodically. We have recently shown that non-diffusing particles move chaotically in a two-dimensional model flow with the above features. In parts of the lung, however, there is asynchrony between the wall motion and the ductal flow immediately outside the alveolus. The extent to which this asynchrony affects kinematic mixing in real alveoli is not yet known. The purpose of this paper is to describe the effect of asynchrony on chaotic advection in our two-dimensional model, in order to understand the circumstances in which this becomes significant.
Experimental Investigation of the Effect of Non-Newtonian Behavior of Blood Flow in the Fontan Circulation
Cheng AL, Pahlevan NM, Rinderknecht DG, Wood JC and Gharib M
The Fontan procedure for univentricular heart defects creates a unique circulation where all pulmonary blood flow is passively supplied directly from systemic veins. Computational simulations, aimed at optimizing the surgery, have assumed blood to be a Newtonian fluid without evaluating the potential error introduced by this assumption. We compared flow behavior between a non-Newtonian blood analog (0.04% xanthan gum) and a control Newtonian fluid (45% glycerol) in a simplified model of the Fontan circulation. Particle image velocimetry was used to examine flow behavior at two different cardiac outputs and two caval blood flow distributions. Pressure and flow rates were measured at each inlet and outlet. Velocity, shear strain, and shear stress maps were derived from velocity data. Power loss was calculated from pressure, flow, and velocity data. Power loss was increased in all test conditions with xanthan gum vs. glycerol (mean 10±2.9% vs. 5.6±1.3%, p=0.032). Pulmonary blood flow distribution differed in all conditions, more so at low cardiac output. Caval blood flow mixing patterns and shear stress were also qualitatively different between the solutions in all conditions. We conclude that assuming blood to be a Newtonian fluid introduces considerable error into simulations of the Fontan circulation, where low-shear flow predominates.
A note on the stability characteristics of the respiratory events
Vadivukkarasan M
The present outbreak enables the researchers from fluid mechanics to widen the understanding of expelling respiratory liquids from a unique perspective to diminish the persistence of COVID-19. This article focuses on uncovering the instability mechanism responsible for forming droplets and aerosols during respiratory events such as breathing, talking, coughing and sneezing. We illustrate a mathematical framework by revisiting the model (Vadivukkarasan and Panchagnula, 2017) and show the associated instabilities during respiratory events. We envisage the combined Rayleigh-Taylor-Kelvin-Helmholtz (R-T-K-H) model as a robust tool for respiratory events. This study highlights the distinct possibility of respiratory droplet formation over multiple instabilities and provides a fundamental understanding. We present the different dominant modes through a ternary phase diagram for three-dimensional numbers (Bond number and Weber numbers). Furthermore, this model can be extended phenomenologically to viscous fluids to satisfy mucus and saliva in the respiratory liquids.
Computational design and experimental analysis of a novel visor for COVID-19 patients receiving high-flow nasal oxygen therapy
Ijaz M, Fhrighil SN, Brett R, Connolly J, Conneely A, O'Connor G, O'Halloran M and Yousefian S
The Covid-19 global pandemic has reshaped the requirements of healthcare sectors worldwide. Following the exposure risks associated with Covid-19, this paper aims to design, optimise, and validate a wearable medical device that reduces the risk of transmission of contagious droplets from infected patients in a hospital setting. This study specifically focuses on those receiving high-flow nasal oxygen therapy. The design process consisted of optimising the geometry of the visor to ensure that the maximum possible percentage of harmful droplets exhaled by the patient can be successfully captured by a vacuum tube attached to the visor. This has been completed by deriving a number of concept designs and assessing their effectiveness, based on numerical analysis, computational fluid dynamics (CFD) simulations and experimental testing. The CFD results are validated using various experimental methods such as Schlieren imaging, particle measurement testing and laser sheet visualisation. Droplet capturing efficiency of the visor was measured through CFD and validated through experimental particle measurement testing. The results presented a 5% deviation between CFD and experimental results. Also, the modifications based on the validated CFD results improved the visor effectiveness by 47% and 38% for breathing and coughing events, respectively.
An experimental investigation of oscillatory flow in the cerebral aqueduct
Sincomb S, Moral-Pulido F, Campos O, Martínez-Bazán C, Haughton V and Sánchez AL
This study aims at clarifying the relation between the oscillatory flow of cerebrospinal fluid (CSF) in the cerebral aqueduct, a narrow conduit connecting the third and fourth ventricles, and the corresponding interventricular pressure difference. Dimensional analysis is used in designing an anatomically correct scaled model of the aqueduct flow, with physical similarity maintained by adjusting the flow frequency and the properties of the working fluid. The time-varying pressure difference across the aqueduct corresponding to a given oscillatory flow rate is measured in parametric ranges covering the range of flow conditions commonly encountered in healthy subjects. Parametric dependences are delineated for the time-averaged pressure fluctuations and for the phase lag between the transaqueductal pressure difference and the flow rate, both having clinical relevance. The results are validated through comparisons with predictions obtained with a previously derived computational model. The parametric quantification in this study enables the derivation of a simple formula for the relation between the transaqueductal pressure and the stroke volume. This relationship can be useful in the quantification of transmantle pressure differences based on non-invasive magnetic-resonance-velocimetry measurements of aqueduct flow for investigation of CSF-related disorders.
Transverse flow under oscillating stimulation in helical square ducts with cochlea-like geometrical curvature and torsion
Harte NC, Obrist D, Caversaccio M, Lajoinie GPR and Wimmer W
The cochlea, situated within the inner ear, is a spiral-shaped, liquid-filled organ responsible for hearing. The physiological significance of its shape remains uncertain. Previous research has scarcely addressed the occurrence of transverse flow within the cochlea, particularly in relation to its unique shape. This study aims to investigate the impact of the geometric features of the cochlea on fluid dynamics by characterizing transverse flow induced by harmonically oscillating axial flow in square ducts with curvature and torsion resembling human cochlear anatomy. We examined four geometries to investigate curvature and torsion effects on axial and transverse flow components. Twelve frequencies from 0.125 Hz to 256 Hz were studied, covering infrasound and low-frequency hearing, with mean inlet velocity amplitudes representing levels expected for normal conversation or louder situations. Our simulations show that torsion contributes significantly to transverse flow in unsteady conditions, and that its contribution increases with increasing oscillation frequency. Curvature alone has a small effect on transverse flow strength, which decreases rapidly with increasing frequency. Strikingly, the combined effect of curvature and torsion on transverse flow is greater than expected from a simple superposition of the two effects, especially when the relative contribution of curvature alone becomes negligible. These findings may be relevant to understanding physiological processes in the cochlea, including metabolite transport and wall shear stress. Further studies are needed to investigate possible implications for cochlear mechanics.