Buoyancy-driven attraction of active droplets

For dissolving active oil droplets in an ambient liquid, it is generally assumed that the Marangoni effect results in repulsive interactions, while the buoyancy effects caused by the density difference between the droplets, diffusing product and the ambient fluid are usually neglected. However, it has been observed in recent experiments that active droplets can form clusters due to buoyancy-driven convection (Krüger et al. Eur. Phys. J. E, vol. 39, 2016, pp. 1-9). In this study, we numerically analyze the buoyancy effect, in addition to the propulsion caused by Marangoni flow (with its strength characterized by Péclet number Pe). The buoyancy effects have their origin in (i) the density difference between the droplet and the ambient liquid, which is characterized by Galileo number Ga, and (ii) the density difference between the diffusing product (i.e. filled micelles) and the ambient liquid, which can be quantified by a solutal Rayleigh number Ra. We analyze how the attracting and repulsing behaviour of neighbouring droplets depends on the control parameters Pe, Ga, and Ra. We find that while the Marangoni effect leads to the well-known repulsion between the interacting droplets, the buoyancy effect of the reaction product leads to buoyancy-driven attraction. At sufficiently large Ra, even collisions between the droplets can take place. Our study on the effect of Ga further shows that with increasing Ga, the collision becomes delayed. Moreover, we derive that the attracting velocity of the droplets, which is characterized by a Reynolds number Red, is proportional to Ra1/4/(ℓ/R), where ℓ/R is the distance between the neighbouring droplets normalized by the droplet radius. Finally, we numerically obtain the repulsive velocity of the droplets, characterized by a Reynolds number Rerep, which is proportional to PeRa-0.38. The balance of attractive and repulsive effect leads to Pe ~ Ra0.63, which agrees well with the transition curve between the regimes with and without collision.

Bistability in Radiatively Heated Melt Ponds

Melting and solidification processes, intertwined with convective flows, play a fundamental role in geophysical contexts. One of these processes is the formation of melt ponds on glaciers, ice shelves, and sea ice. It is driven by solar radiation and is of great significance for Earth's heat balance, as it significantly lowers the albedo. Through direct numerical simulations and theoretical analysis, we unveil a bistability phenomenon in the melt pond dynamics. As solar radiation intensity and the melt pond's initial depth vary, an abrupt transition occurs: this tipping point transforms the system from a stable fully frozen state to another stable equilibrium state, characterized by a distinct melt pond depth. The physics of this transition can be understood within a heat flux balance model, which exhibits excellent agreement with our numerical results. Together with the Grossmann-Lohse theory for internally heated convection, the model correctly predicts the bulk temperature and the flow strength within the melt ponds, offering insight into the coupling of phase transitions with adjacent turbulent flows and the interplay between convective melting and radiation-driven processes.

GPU accelerated digital twins of the human heart open new routes for cardiovascular research

The recruitment of patients for rare or complex cardiovascular diseases is a bottleneck for clinical trials and digital twins of the human heart have recently been proposed as a viable alternative. In this paper we present an unprecedented cardiovascular computer model which, relying on the latest GPU-acceleration technologies, replicates the full multi-physics dynamics of the human heart within a few hours per heartbeat. This opens the way to extensive simulation campaigns to study the response of synthetic cohorts of patients to cardiovascular disorders, novel prosthetic devices or surgical procedures. As a proof-of-concept we show the results obtained for left bundle branch block disorder and the subsequent cardiac resynchronization obtained by pacemaker implantation. The in-silico results closely match those obtained in clinical practice, confirming the reliability of the method. This innovative approach makes possible a systematic use of digital twins in cardiovascular research, thus reducing the need of real patients with their economical and ethical implications. This study is a major step towards in-silico clinical trials in the era of digital medicine.

Aspect Ratio Dependence of Heat Transfer in a Cylindrical Rayleigh-Bénard Cell

While the heat transfer and the flow dynamics in a cylindrical Rayleigh-Bénard (RB) cell are rather independent of the aspect ratio Γ (diameter/height) for large Γ, a small-Γ cell considerably stabilizes the flow and thus affects the heat transfer. Here, we first theoretically and numerically show that the critical Rayleigh number for the onset of convection at given Γ follows Ra_{c,Γ}∼Ra_{c,∞}(1+CΓ^{-2})^{2}, with C≲1.49 for Oberbeck-Boussinesq (OB) conditions. We then show that, in a broad aspect ratio range (1/32)≤Γ≤32, the rescaling Ra→Ra_{ℓ}≡Ra[Γ^{2}/(C+Γ^{2})]^{3/2} collapses various OB numerical and almost-OB experimental heat transport data Nu(Ra,Γ). Our findings predict the Γ dependence of the onset of the ultimate regime Ra_{u,Γ}∼[Γ^{2}/(C+Γ^{2})]^{-3/2} in the OB case. This prediction is consistent with almost-OB experimental results (which only exist for Γ=1, 1/2, and 1/3) for the transition in OB RB convection and explains why, in small-Γ cells, much larger Ra (namely, by a factor Γ^{-3}) must be achieved to observe the ultimate regime.

Diffusion-Free Scaling in Rotating Spherical Rayleigh-Bénard Convection

Direct numerical simulations are employed to reveal three distinctly different flow regions in rotating spherical Rayleigh-Bénard convection. In the high-latitude region I vertical (parallel to the axis of rotation) convective columns are generated between the hot inner and the cold outer sphere. The mid-latitude region I I is dominated by vertically aligned convective columns formed between the Northern and Southern hemispheres of the outer sphere. The diffusion-free scaling, which indicates bulk-dominated convection, originates from this mid-latitude region. In the equator region I I I , the vortices are affected by the outer spherical boundary and are much shorter than in region I I .

Droplets generated from toilets during urination as a possible vehicle of carbapenem-resistant Klebsiella pneumoniae

BACKGROUND

In the health care setting, infection control actions are fundamental for containing the dissemination of multidrug-resistant bacteria (MDR). Carbapenemase-producing Enterobacterales (CPE), especially Klebsiella pneumoniae (CR-KP), can spread among patients, although the dynamics of transmission are not fully known. Since CR-KP is present in wastewater and microorganisms are not completely removed from the toilet bowl by flushing, the risk of transmission in settings where toilets are shared should be addressed. We investigated whether urinating generates droplets that can be a vehicle for bacteria and explored the use of an innovative foam to control and eliminate this phenomenon.

METHODS

To study droplet formation during urination, we set up an experiment in which different geometrical configurations of toilets could be reproduced and customized. To demonstrate that droplets can mobilize bacteria from the toilet bowl, a standard ceramic toilet was contaminated with a KPC-producing Klebsiella pneumoniae ST101 isolate. Then, we reproduced urination and attached culture dishes to the bottom of the toilet lid for bacterial colony recovery with and without foam.

RESULTS

Rebound droplets invariably formed, irrespective of the geometrical configuration of the toilet. In microbiological experiments, we demonstrated that bacteria are always mobilized from the toilet bowl (mean value: 0.11 ± 0.05 CFU/cm2) and showed that a specific foam layer can completely suppress mobilization.

CONCLUSIONS

Our study demonstrated that droplets generated from toilets during urination can be a hidden source of CR-KP transmission in settings where toilets are shared among colonized and noncolonized patients.

Droplet plume emission during plasmonic bubble growth in ternary liquids

Plasmonic bubbles are of great relevance in numerous applications, including catalytic reactions, micro/nanomanipulation of molecules or particles dispersed in liquids, and cancer therapeutics. So far, studies have been focused on bubble nucleation in pure liquids. Here we investigate plasmonic bubble nucleation in ternary liquids consisting of ethanol, water, and trans-anethole oil, which can show the so-called ouzo effect. We find that oil (trans-anethole) droplet plumes are produced around the growing plasmonic bubbles. The nucleation of the microdroplets and their organization in droplet plumes is due to the symmetry breaking of the ethanol concentration field during the selective evaporation of ethanol from the surrounding ternary liquids into the growing plasmonic bubbles. Numerical simulations show the existence of a critical Marangoni number Ma (the ratio between solutal advection rate and the diffusion rate), above which the symmetry breaking of the ethanol concentration field occurs, leading to the emission of the droplet plumes. The numerical results agree with the experimental observation that more plumes are emitted with increasing ethanol-water relative weight ratios and hence Ma. Our findings on the droplet plume formation reveal the rich phenomena of plasmonic bubble nucleation in multicomponent liquids and help to pave the way to achieve enhanced mixing in multicomponent liquids in chemical, pharmaceutical, and cosmetic industries.

Anatomic features in SCAD assessed by CCT: A propensity score matching case control study

PURPOSE

Spontaneous coronary artery dissection (SCAD) may occur in middle age population without any cardiovascular risk factor. We retrospectively evaluated anatomic features of 11 patients with SCAD using a coronary arteries computed tomography (CCT), compared to age and sex balanced patients who underwent CCT.

MATERIAL AND METHODS

CCT was performed in 11 patients (7 females and 4 males) as follow-up in patients with SCAD (left anterior descending - LAD or circumflex artery - Cx) and compared, using the propensity score matching analysis, with 11 healthy patients. Several anatomic features were evaluated: Left main (LM) length, angle between descending coronary artery (LAD) and its first branch, angle between LAD and LM, distance from the annulus to RCA (a-RCA distance) and LM (a-LM distance) ostia and their ratio; ratio between LM length and length a-LM and tortuosity score of the vessel with SCAD. A fluid dynamic analysis has been performed to evaluate the effects on shear stress of vessels wall.

RESULTS

LM length was significantly shorter in patients with SCAD versus healthy subjects (P=0.01) as well as LM length/a-LM (P=0.03) and the angle between LAD and the first adjacent branch was sharper (P<0.01). Tortuosity score showed a statistically significant difference between groups (P<0.001). Fluid dynamic analysis demonstrates that, in SCAD group, an angle<90 degree is present at the first bifurcation and it can be a cause of increased strain on vessel wall in patients with high tortuosity of coronary artery.

CONCLUSION

Tortuosity and angle between the LAD and the adjacent arterial branch combined may determine increased shear stress on the vessel wall that increases the risk of SCAD.

Extended Lifetime of Respiratory Droplets in a Turbulent Vapor Puff and Its Implications on Airborne Disease Transmission

To quantify the fate of respiratory droplets under different ambient relative humidities, direct numerical simulations of a typical respiratory event are performed. We found that, because small droplets (with initial diameter of 10  μm) are swept by turbulent eddies in the expelled humid puff, their lifetime gets extended by a factor of more than 30 times as compared to what is suggested by the classical picture by Wells, for 50% relative humidity. With increasing ambient relative humidity the extension of the lifetimes of the small droplets further increases and goes up to around 150 times for 90% relative humidity, implying more than 2 m advection range of the respiratory droplets within 1 sec. Employing Lagrangian statistics, we demonstrate that the turbulent humid respiratory puff engulfs the small droplets, leading to many orders of magnitude increase in their lifetimes, implying that they can be transported much further during the respiratory events than the large ones. Our findings provide the starting points for larger parameter studies and may be instructive for developing strategies on optimizing ventilation and indoor humidity control. Such strategies are key in mitigating the COVID-19 pandemic in the present autumn and upcoming winter.

Extended Lifetime of Respiratory Droplets in a Turbulent Vapor Puff and Its Implications on Airborne Disease Transmission

To quantify the fate of respiratory droplets under different ambient relative humidities, direct numerical simulations of a typical respiratory event are performed. We found that, because small droplets (with initial diameter of 10  μm) are swept by turbulent eddies in the expelled humid puff, their lifetime gets extended by a factor of more than 30 times as compared to what is suggested by the classical picture by Wells, for 50% relative humidity. With increasing ambient relative humidity the extension of the lifetimes of the small droplets further increases and goes up to around 150 times for 90% relative humidity, implying more than 2 m advection range of the respiratory droplets within 1 sec. Employing Lagrangian statistics, we demonstrate that the turbulent humid respiratory puff engulfs the small droplets, leading to many orders of magnitude increase in their lifetimes, implying that they can be transported much further during the respiratory events than the large ones. Our findings provide the starting points for larger parameter studies and may be instructive for developing strategies on optimizing ventilation and indoor humidity control. Such strategies are key in mitigating the COVID-19 pandemic in the present autumn and upcoming winter.

Sensitivity analysis of an electrophysiology model for the left ventricle

Modelling the cardiac electrophysiology entails dealing with the uncertainties related to the input parameters such as the heart geometry and the electrical conductivities of the tissues, thus calling for an uncertainty quantification (UQ) of the results. Since the chambers of the heart have different shapes and tissues, in order to make the problem affordable, here we focus on the left ventricle with the aim of identifying which of the uncertain inputs mostly affect its electrophysiology. In a first phase, the uncertainty of the input parameters is evaluated using data available from the literature and the output quantities of interest (QoIs) of the problem are defined. According to the polynomial chaos expansion, a training dataset is then created by sampling the parameter space using a quasi-Monte Carlo method whereas a smaller independent dataset is used for the validation of the resulting metamodel. The latter is exploited to run a global sensitivity analysis with nonlinear variance-based indices and thus reduce the input parameter space accordingly. Thereafter, the uncertainty probability distribution of the QoIs are evaluated using a direct UQ strategy on a larger dataset and the results discussed in the light of the medical knowledge.

Periodically Modulated Thermal Convection

Many natural and industrial turbulent flows are subjected to time-dependent boundary conditions. Despite being ubiquitous, the influence of temporal modulations (with frequency f) on global transport properties has hardly been studied. Here, we perform numerical simulations of Rayleigh-Bénard convection with time periodic modulation in the temperature boundary condition and report how this modulation can lead to a significant heat flux (Nusselt number Nu) enhancement. Using the concept of Stokes thermal boundary layer, we can explain the onset frequency of the Nu enhancement and the optimal frequency at which Nu is maximal, and how they depend on the Rayleigh number Ra and Prandtl number Pr. From this, we construct a phase diagram in the 3D parameter space (f, Ra, Pr) and identify the following: (i) a regime where the modulation is too fast to affect Nu; (ii) a moderate modulation regime, where Nu increases with decreasing f, and (iii) slow modulation regime, where Nu decreases with further decreasing f. Our findings provide a framework to study other types of turbulent flows with time-dependent forcing.

Flow organization and heat transfer in turbulent wall sheared thermal convection

We perform direct numerical simulations of wall sheared Rayleigh–Bénard convection for Rayleigh numbers up to [Image: see text], Prandtl number unity and wall shear Reynolds numbers up to [Image: see text]. Using the Monin–Obukhov length [Image: see text] we observe the presence of three different flow states, a buoyancy dominated regime ([Image: see text]; with [Image: see text] the thermal boundary layer thickness), a transitional regime ([Image: see text]; with [Image: see text] the height of the domain) and a shear dominated regime ([Image: see text]). In the buoyancy dominated regime, the flow dynamics is similar to that of turbulent thermal convection. The transitional regime is characterized by rolls that are increasingly elongated with increasing shear. The flow in the shear dominated regime consists of very large-scale meandering rolls, similar to the ones found in conventional Couette flow. As a consequence of these different flow regimes, for fixed [Image: see text] and with increasing shear, the heat transfer first decreases, due to the breakup of the thermal rolls, and then increases at the beginning of the shear dominated regime. In the shear dominated regime the Nusselt number [Image: see text] effectively scales as [Image: see text] with [Image: see text], while we find [Image: see text] in the buoyancy dominated regime. In the transitional regime, the effective scaling exponent is [Image: see text], but the temperature and velocity profiles in this regime are not logarithmic yet, thus indicating transient dynamics and not the ultimate regime of thermal convection.

Multiple States in Turbulent Large-Aspect-Ratio Thermal Convection: What Determines the Number of Convection Rolls?

Wall-bounded turbulent flows can take different statistically stationary turbulent states, with different transport properties, even for the very same values of the control parameters. What state the system takes depends on the initial conditions. Here we analyze the multiple states in large-aspect ratio (Γ) two-dimensional turbulent Rayleigh-Bénard flow with no-slip plates and horizontally periodic boundary conditions as model system. We determine the number n of convection rolls, their mean aspect ratios Γ_{r}=Γ/n, and the corresponding transport properties of the flow (i.e., the Nusselt number Nu), as function of the control parameters Rayleigh (Ra) and Prandtl number. The effective scaling exponent β in Nu∼Ra^{β} is found to depend on the realized state and thus Γ_{r}, with a larger value for the smaller Γ_{r}. By making use of a generalized Friedrichs inequality, we show that the elliptical shape of the rolls and viscous damping determine the Γ_{r} window for the realizable turbulent states. The theoretical results are in excellent agreement with our numerical finding 2/3≤Γ_{r}≤4/3, where the lower threshold is approached for the larger Ra. Finally, we show that the theoretical approach to frame Γ_{r} also works for free-slip boundary conditions.

Biomechanical properties and histomorphometric features of aortic tissue in patients with or without bicuspid aortic valve

Background

We sought to investigate and compare biomechanical properties and histomorphometric findings of thoracic ascending aorta aneurysm (TAA) tissue from patients with bicuspid aortic valve (BAV) and tricuspid aortic valve (TAV) in order to clarify mechanisms underlying differences in the clinical course.

Methods

Circumferential sections of TAA tissue in patients with BAV (BAV-TAA) and TAV (TAV-TAA) were obtained during surgery and used for biomechanical tests and histomorphometrical analysis.

Results

In BAV-TAA, we observed biomechanical higher peak stress and lower Young modulus values compared with TAV-TAA wall. The right lateral longitudinal region seemed to be the most fragile zone of the TAA wall. Mechanical stress-induced rupture of BAV-TAA tissue was sudden and uniform in all aortic wall layers, whereas a gradual and progressive aortic wall breakage was described in TAV-TAA. Histomorphometric analysis revealed higher amount of collagen but not elastin in BAV-TAA tunica media.

Conclusions

The higher deformability of BAV-TAA tissue supports the hypothesis that increased wall shear stress doesn’t explain the increased risk of sudden onset of rupture and dissection; other mechanisms, likely related to alteration of specific genetic pathways and epigenetic signals, could be investigated to explain differences in aortic dissection and rupture in BAV patients.

Modeling mitral valve stenosis: A parametric study on the stenosis severity level

New computational techniques providing more accurate representation of human heart pathologies could help uncovering relevant physical phenomena and improve the outcome of medical therapies. In this framework, the present work describes an efficient computational model for the evaluation of the ventricular flow alteration in presence of mitral valve stenosis. The model is based on the direct numerical simulation of the Navier-Stokes equations two-way coupled with a structural solver for the left ventricle and mitral valve dynamics. The presence of mitral valve stenosis is mimicked by a single-parameter constraint acting on the kinematics of the mitral leaflets. Four different degrees of mitral valve stenosis are considered focusing on the hemodynamic alterations occurring in pathologic conditions. The mitral jet, generated during diastole, is seen to shrink and strengthen when the stenosis gets more severe. As a consequence, the kinetic energy of the flow, the tissues shear stresses, the transvalvular pressure drop and mitral regurgitation increase. It results that, as the stenosis severity level increases, the geometric and effective orifice areas decrease up to 50% with respect the normal case due to the reduced leaflets mobility and stronger blood acceleration during the diastolic phase. The modified intraventricular hemodynamics is also related to a stronger pressure gradient that, for severe stenosis, can be more than ten times larger than the healthy valve case. These computational results are fully consistent with the available clinical literature and open the way to the virtual assessment of surgical procedures and to the evaluation of prosthetic devices.

Rough-wall turbulent Taylor-Couette flow: The effect of the rib height

In this study, we combine experiments and direct numerical simulations to investigate the effects of the height of transverse ribs at the walls on both global and local flow properties in turbulent Taylor-Couette flow. We create rib roughness by attaching up to 6 axial obstacles to the surfaces of the cylinders over an extensive range of rib heights, up to blockages of 25% of the gap width. In the asymptotic ultimate regime, where the transport is independent of viscosity, we emperically find that the prefactor of the [Formula: see text] scaling (corresponding to the drag coefficient [Formula: see text] being constant) scales with the number of ribs [Formula: see text] and by the rib height [Formula: see text]. The physical mechanism behind this is that the dominant contribution to the torque originates from the pressure forces acting on the rib which scale with the rib height. The measured scaling relation of [Formula: see text] is slightly smaller than the expected [Formula: see text] scaling, presumably because the ribs cannot be regarded as completely isolated but interact. In the counter-rotating regime with smooth walls, the momentum transport is increased by turbulent Taylor vortices. We find that also in the presence of transverse ribs these vortices persist. In the counter-rotating regime, even for large roughness heights, the momentum transport is enhanced by these vortices.

Transition to the Ultimate Regime in Two-Dimensional Rayleigh-Bénard Convection

The possible transition to the so-called ultimate regime, wherein both the bulk and the boundary layers are turbulent, has been an outstanding issue in thermal convection, since the seminal work by Kraichnan [Phys. Fluids 5, 1374 (1962)PFLDAS0031-917110.1063/1.1706533]. Yet, when this transition takes place and how the local flow induces it is not fully understood. Here, by performing two-dimensional simulations of Rayleigh-Bénard turbulence covering six decades in Rayleigh number Ra up to 10^{14} for Prandtl number Pr=1, for the first time in numerical simulations we find the transition to the ultimate regime, namely, at Ra^{*}=10^{13}. We reveal how the emission of thermal plumes enhances the global heat transport, leading to a steeper increase of the Nusselt number than the classical Malkus scaling Nu∼Ra^{1/3} [Proc. R. Soc. A 225, 196 (1954)PRLAAZ1364-502110.1098/rspa.1954.0197]. Beyond the transition, the mean velocity profiles are logarithmic throughout, indicating turbulent boundary layers. In contrast, the temperature profiles are only locally logarithmic, namely, within the regions where plumes are emitted, and where the local Nusselt number has an effective scaling Nu∼Ra^{0.38}, corresponding to the effective scaling in the ultimate regime.

Flow-induced dissolution of femtoliter surface droplet arrays

The dissolution of liquid nanodroplets is a crucial step in many applied processes, such as separation and dispersion in the food industry, crystal formation of pharmaceutical products, concentrating and analysis in medical diagnosis, and drug delivery in aerosols. In this work, using both experiments and numerical simulations, we quantitatively study the dissolution dynamics of femtoliter surface droplets in a highly ordered array under a uniform flow. Our results show that the dissolution of femtoliter droplets strongly depends on their spatial positions relative to the flow direction, drop-to-drop spacing in the array, and the imposed flow rate. In some particular cases, the droplet at the edge of the array can dissolve about 30% faster than the ones located near the centre. The dissolution rate of the droplet increases by 60% as the inter-droplet spacing is increased from 2.5 μm to 20 μm. Moreover, the droplets close to the front of the flow commence to shrink earlier than those droplets in the center of the array. The average dissolution rate is faster for the faster flow. As a result, the dissolution time (Ti) decreases with the Reynolds number (Re) of the flow as Ti ∝ Re-3/4. The experimental results are in good agreement with the numerical simulations where the advection-diffusion equation for the concentration field is solved and the concentration gradient on the surface of the drop is computed. The findings suggest potential approaches to manipulate nanodroplet sizes in droplet arrays simply by dissolution controlled by an external flow. The obtained droplets with varying curvatures may serve as templates for generating multifocal microlenses in one array.

Diffusive interaction of multiple surface nanobubbles: shrinkage, growth, and coarsening

Surface nanobubbles are nanoscopic spherical-cap shaped gaseous domains on immersed substrates which are stable, even for days. After the stability of a single surface nanobubble has been theoretically explained, i.e. contact line pinning and gas oversaturation are required to stabilize it against diffusive dissolution [Lohse and Zhang, Phys. Rev. E, 2015, 91, 031003(R)], here we focus on the collective diffusive interaction of multiple nanobubbles. For that purpose we develop a finite difference scheme for the diffusion equation with the appropriate boundary conditions and with the immersed boundary method used to represent the growing or shrinking bubbles. After validation of the scheme against the exact results of Epstein and Plesset for a bulk bubble [J. Chem. Phys., 1950, 18, 1505] and of Lohse and Zhang for a surface bubble, the framework of these simulations is used to describe the coarsening process of competitively growing nanobubbles. The coarsening process for such diffusively interacting nanobubbles slows down with advancing time and increasing bubble distance. The present results for surface nanobubbles are also applicable for immersed surface nanodroplets, for which better controlled experimental results of the coarsening process exist.

Effects of mitral chordae tendineae on the flow in the left heart ventricle

In this paper a computational model for the ventricular flow with a mitral valve and modeled chordae tendineae is presented. The results are compared with an analogous case in which the chordae are not included and their presence is replaced by kinematic boundary conditions. The problem is studied using direct numerical simulation of the Navier-Stokes equations, two-way coupled with a structural solver for the ventricle and mitral valve dynamics. An experimental validation of the model is performed by a comparison of the results with a companion dedicated experiment. It is found that the inclusion of the chordae tendineae makes the model self-consistent thus avoiding the use of ad hoc kinematic constraints to mimic their effect. In this way it is possible to simulate the correct system dynamics without user-defined parameters. More in detail, the results have shown that the mitral valve dynamics can be described also without chordae with the help of ad hoc kinematic constrains, whereas the changes produced in the intra-ventricular flow need the explicit consideration of the chordae in the model. On the other hand, the computational load increases owing to the presence of additional structures that, being thin filaments, are also demanding for the spatial resolution requirements. Since the presence of the chordae tendineae produces only specific differences in the overall flow structure, we conclude that their explicit modeling should be limited to those cases in which their presence is unavoidable.

Controlling Heat Transport and Flow Structures in Thermal Turbulence Using Ratchet Surfaces

In this combined experimental and numerical study on thermally driven turbulence in a rectangular cell, the global heat transport and the coherent flow structures are controlled with an asymmetric ratchetlike roughness on the top and bottom plates. We show that, by means of symmetry breaking due to the presence of the ratchet structures on the conducting plates, the orientation of the large scale circulation roll (LSCR) can be locked to a preferred direction even when the cell is perfectly leveled out. By introducing a small tilt to the system, we show that the LSCR orientation can be tuned and controlled. The two different orientations of LSCR give two quite different heat transport efficiencies, indicating that heat transport is sensitive to the LSCR direction over the asymmetric roughness structure. Through a quantitative analysis of the dynamics of thermal plume emissions and the orientation of the LSCR over the asymmetric structure, we provide a physical explanation for these findings. The current work has important implications for passive and active flow control in engineering, biofluid dynamics, and geophysical flows.

Roughness-Facilitated Local 1/2 Scaling Does Not Imply the Onset of the Ultimate Regime of Thermal Convection

In thermal convection, roughness is often used as a means to enhance heat transport, expressed in Nusselt number. Yet there is no consensus on whether the Nusselt vs Rayleigh number scaling exponent (Nu∼Ra^{β}) increases or remains unchanged. Here we numerically investigate turbulent Rayleigh-Bénard convection over rough plates in two dimensions, up to Ra≈10^{12}. Varying the height and wavelength of the roughness elements with over 200 combinations, we reveal the existence of two universal regimes. In the first regime, the local effective scaling exponent can reach up to 1/2. However, this cannot be explained as the attainment of the so-called ultimate regime as suggested in previous studies, because a further increase in Ra leads to the second regime, in which the scaling saturates back to a value close to the smooth wall case. Counterintuitively, the transition from the first to the second regime corresponds to the competition between bulk and boundary layer flow: from the bulk-dominated regime back to the classical boundary-layer-controlled regime. Our study demonstrates that the local 1/2 scaling does not necessarily signal the onset of ultimate turbulence.

Confined Rayleigh-Bénard, Rotating Rayleigh-Bénard, and Double Diffusive Convection: A Unifying View on Turbulent Transport Enhancement through Coherent Structure Manipulation

Many natural and engineering systems are simultaneously subjected to a driving force and a stabilizing force. The interplay between the two forces, especially for highly nonlinear systems such as fluid flow, often results in surprising features. Here we reveal such features in three different types of Rayleigh-Bénard (RB) convection, i.e., buoyancy-driven flow with the fluid density being affected by a scalar field. In the three cases different stabilizing forces are considered, namely (i) horizontal confinement, (ii) rotation around a vertical axis, and (iii) a second stabilizing scalar field. Despite the very different nature of the stabilizing forces and the corresponding equations of motion, at moderate strength we counterintuitively but consistently observe an enhancement in the flux, even though the flow motion is weaker than the original RB flow. The flux enhancement occurs in an intermediate regime in which the stabilizing force is strong enough to alter the flow structures in the bulk to a more organized morphology, yet not too strong to severely suppress the flow motions. Near the optimal transport enhancements all three systems exhibit a transition from a state in which the thermal boundary layer (BL) is nested inside the momentum BL to the one with the thermal BL being thicker than the momentum BL. The observed optimal transport enhancement is explained through an optimal coupling between the suction of hot or fresh fluid and the corresponding scalar fluctuations.

Vertically Bounded Double Diffusive Convection in the Finger Regime: Comparing No-Slip versus Free-Slip Boundary Conditions

Vertically bounded fingering double diffusive convection is numerically investigated, focusing on the influences of different velocity boundary conditions, i.e., the no-slip condition, which is inevitable in the lab-scale experimental researches, and the free-slip condition, which is an approximation for the interfaces in many natural environments, such as the oceans. For both boundary conditions the flow is dominated by fingers and the global responses follow the same scaling laws, with enhanced prefactors for the free-slip cases. Therefore, the laboratory experiments with the no-slip boundaries serve as a good model for the finger layers in the ocean. Moreover, in the free-slip case, although the tangential shear stress is eliminated at the boundaries, the local dissipation rate in the near-wall region may exceed the value found in the no-slip cases, which is caused by the stronger vertical motions of horizontally focused fingers and sheet structures near the free-slip boundaries. This counterintuitive result might be relevant for properly estimating and modeling the mixing and entrainment phenomena at free surfaces and interfaces widespread in oceans and geophysical flows.

Logarithmic Mean Temperature Profiles and Their Connection to Plume Emissions in Turbulent Rayleigh-Bénard Convection

Two-dimensional simulations of Rayleigh-Bénard convection at Ra=5×10^{10} show that vertical logarithmic mean temperature profiles can be observed in regions of the boundary layer where thermal plumes are emitted. The profile is logarithmic only in these regions and not in the rest of the boundary layer where it is sheared by the large-scale wind and impacted by plumes. In addition, the logarithmic behavior is not visible in the horizontal average. The findings reveal that the temperature profiles are strongly connected to thermal plume emission, and they support a perception that parts of the boundary layer can be turbulent while others are not. The transition to the ultimate regime, in which the boundary layers are considered to be fully turbulent, can therefore be understood as a gradual increase in the fraction of the plume-emitting ("turbulent") regions of the boundary layer.

Effect of velocity boundary conditions on the heat transfer and flow topology in two-dimensional Rayleigh-Bénard convection

The effect of various velocity boundary condition is studied in two-dimensional Rayleigh-Bénard convection. Combinations of no-slip, stress-free, and periodic boundary conditions are used on both the sidewalls and the horizontal plates. For the studied Rayleigh numbers Ra between 10(8) and 10(11) the heat transport is lower for Γ=0.33 than for Γ=1 in case of no-slip sidewalls. This is, surprisingly, the opposite for stress-free sidewalls, where the heat transport increases for a lower aspect ratio. In wider cells the aspect-ratio dependence is observed to disappear for Ra ≥ 10(10). Two distinct flow types with very different dynamics can be seen, mostly dependent on the plate velocity boundary condition, namely roll-like flow and zonal flow, which have a substantial effect on the dynamics and heat transport in the system. The predominantly horizontal zonal flow suppresses heat flux and is observed for stress-free and asymmetric plates. Low aspect-ratio periodic sidewall simulations with a no-slip boundary condition on the plates also exhibit zonal flow. In all the other cases, the flow is roll like. In two-dimensional Rayleigh-Bénard convection, the velocity boundary conditions thus have large implications on both roll-like and zonal flow that have to be taken into consideration before the boundary conditions are imposed.

Annular dilatation and loss of sino-tubular junction in aneurysmatic aorta: implications on leaflet quality at the time of surgery. A finite element study

OBJECTIVES

In the belief that stress is the main determinant of leaflet quality deterioration, we sought to evaluate the effect of annular and/or sino-tubular junction dilatation on leaflet stress. A finite element computer-assisted stress analysis was used to model four different anatomic conditions and analyse the consequent stress pattern on the aortic valve.

METHODS

Theoretical models of four aortic root configurations (normal, with dilated annulus, with loss of sino-tubular junction and with both dilatation simultaneously) were created with computer-aided design technique. The pattern of stress and strain was then analysed by means of finite elements analysis, when a uniform pressure of 100 mmHg was applied to the model. Analysis produced von Mises charts (colour-coded, computational, three-dimensional stress-pattern graphics) and bidimensional plots of compared stress on arc-linear line, which allowed direct comparison of stress in the four different conditions.

RESULTS

Stresses both on the free margin and on the 'belly' of the leaflet rose from 0.28 MPa (normal conditions) to 0.32 MPa (+14%) in case of isolated dilatation of the sino-tubular junction, while increased to 0.42 MPa (+67%) in case of isolated annular dilatation, with no substantial difference whether sino-tubular junction dilatation was present or not.

CONCLUSIONS

Annular dilatation is the key element determining an increased stress on aortic leaflets independently from an associated sino-tubular junction dilatation. The presence of annular dilatation associated with root aneurysm greatly decreases the chance of performing a valve sparing procedure without the need for additional manoeuvres on leaflet tissue. This information may lead to a refinement in the optimal surgical strategy.

Spatial distribution of heat flux and fluctuations in turbulent Rayleigh-Bénard convection

We numerically investigate the radial dependence of the velocity and temperature fluctuations and of the time-averaged heat flux j ¯(r) in a cylindrical Rayleigh-Bénard cell with aspect ratio Γ=1 for Rayleigh numbers Ra between 2×10^{6} and 2×10^{9} at a fixed Prandtl number Pr=5.2. The numerical results reveal that the heat flux close to the sidewall is larger than in the center and that, just as the global heat transport, it has an effective power law dependence on the Rayleigh number, j ¯(r)∝Ra{γ{j}(r)}. The scaling exponent γ{j}(r) decreases monotonically from 0.43 near the axis (r≈0) to 0.29 close to the sidewalls (r≈D/2). The effective exponents near the axis and the sidewall agree well with the measurements of Shang et al. [Phys. Rev. Lett. 100, 244503 (2008)] and the predictions of Grossmann and Lohse [Phys. Fluids 16, 1070 (2004)]. Extrapolating our results to large Rayleigh number would imply a crossover at Ra≈10^{15}, where the heat flux near the axis would begin to dominate. In addition, we find that the local heat flux is more than twice as high at the location where warm or cold plumes go up or down than in plume depleted regions.

Logarithmic temperature profiles in turbulent Rayleigh-Bénard convection

We report results for the temperature profiles of turbulent Rayleigh-Bénard convection (RBC) in the interior of a cylindrical sample of aspect ratio Γ≡D/L=0.50 (D and L are the diameter and height, respectively). Both in the classical and in the ultimate state of RBC we find that the temperature varies as A×ln(z/L)+B, where z is the distance from the bottom or top plate. In the classical state, the coefficient A decreases in the radial direction as the distance from the side wall increases. For the ultimate state, the radial dependence of A has not yet been determined. These findings are based on experimental measurements over the Rayleigh-number range 4×10(12)≲Ra≲10(15) for a Prandtl number Pr≃0.8 and on direct numerical simulation at Ra=2×10(12), 2×10(11), and 2×10(10), all for Pr=0.7.

Thermal boundary layer profiles in turbulent Rayleigh-Bénard convection in a cylindrical sample

We numerically investigate the structures of the near-plate temperature profiles close to the bottom and top plates of turbulent Rayleigh-Bénard flow in a cylindrical sample at Rayleigh numbers Ra = 10(8) to Ra = 2 × 10(12) and Prandtl numbers Pr = 6.4 and Pr = 0.7 with the dynamical frame method [Zhou and Xia, Phys. Rev. Lett. 104, 104301 (2010)], thus extending previous results for quasi-two-dimensional systems to three-dimensional systems. The dynamical frame method shows that the measured temperature profiles in the spatially and temporally local frame are much closer to the temperature profile of a laminar, zero-pressure gradient boundary layer (BL) according to Pohlhausen than in the fixed reference frame. The deviation between the measured profiles in the dynamical reference frame and the laminar profiles increases with decreasing Pr, where the thermal BL is more exposed to the bulk fluctuations due to the thinner kinetic BL, and increasing Ra, where more plumes are passing the measurement location.

Effect of vapor bubbles on velocity fluctuations and dissipation rates in bubbly Rayleigh-Bénard convection

Numerical results for kinetic and thermal energy dissipation rates in bubbly Rayleigh-Bénard convection are reported. Bubbles have a twofold effect on the flow: on the one hand, they absorb or release heat to the surrounding liquid phase, thus tending to decrease the temperature differences responsible for the convective motion; but on the other hand, the absorbed heat causes the bubbles to grow, thus increasing their buoyancy and enhancing turbulence (or, more properly, pseudoturbulence) by generating velocity fluctuations. This enhancement depends on the ratio of the sensible heat to the latent heat of the phase change, given by the Jakob number, which determines the dynamics of the bubble growth.

Evaluation of prosthetic-valved devices by means of numerical simulations

The in vivo evaluation of prosthetic device performance is often difficult, if not impossible. In particular, in order to deal with potential problems such as thrombosis, haemolysis, etc., which could arise when a patient undergoes heart valve replacement, a thorough understanding of the blood flow dynamics inside the devices interacting with natural or composite tissues is required. Numerical simulation, combining both computational fluid and structure dynamics, could provide detailed information on such complex problems. In this work, a numerical investigation of the mechanics of two composite aortic prostheses during a cardiac cycle is presented. The numerical tool presented is able to reproduce accurately the flow and structure dynamics of the prostheses. The analysis shows that the vortical structures forming inside the two different grafts do not influence the kinematics of a bileaflet valve or the main coronary flow, whereas major differences are present for the stress status near the suture line of the coronaries to the prostheses. The results are in agreement with in vitro and in vivo observations found in literature.

Heat transfer mechanisms in bubbly Rayleigh-Bénard convection

The heat transfer mechanism in Rayleigh-Bénard convection in a liquid with a mean temperature close to its boiling point is studied through numerical simulations with pointlike vapor bubbles, which are allowed to grow or shrink through evaporation and condensation and which act back on the flow both thermally and mechanically. It is shown that the effect of the bubbles is strongly dependent on the ratio of the sensible heat to the latent heat as embodied in the Jakob number Ja. For very small Ja the bubbles stabilize the flow by absorbing heat in the warmer regions and releasing it in the colder regions. With an increase in Ja, the added buoyancy due to the bubble growth destabilizes the flow with respect to single-phase convection and considerably increases the Nusselt number.

Prandtl-, Rayleigh-, and Rossby-number dependence of heat transport in turbulent rotating Rayleigh-Bénard convection

Experimental and numerical data for the heat transfer as a function of the Rayleigh, Prandtl, and Rossby numbers in turbulent rotating Rayleigh-Bénard convection are presented. For relatively small Ra approximately 10(8) and large Pr modest rotation can enhance the heat transfer by up to 30%. At larger Ra there is less heat-transfer enhancement, and at small Pr less, similar 0.7 there is no heat-transfer enhancement at all. We suggest that the small-Pr behavior is due to the breakdown of the heat-transfer-enhancing Ekman pumping because of larger thermal diffusion.

Numerical and experimental investigation of structure-function scaling in turbulent Rayleigh-Bénard convection

Direct numerical simulation and stereoscopic particle image velocimetry of turbulent convection are used to gather spatial data for the calculation of structure functions. We wish to add to the ongoing discussion in the literature whether temperature acts as an active or passive scalar in turbulent convection, with consequences for structure-function scaling. The simulation results show direct confirmation of the scalings derived by Bolgiano and Obukhov for turbulence with an active scalar for both velocity and temperature statistics. The active-scalar range shifts to larger scales when the forcing parameter (Rayleigh number) is increased. Furthermore, a close inspection of local turbulent length scales (Kolmogorov and Bolgiano lengths) confirms conjectures from earlier studies that the oft-used global averages are not suited for the interpretation of structure functions. In the experiment, a characterization of the domain-filling large-scale circulation of confined convection is carried out for comparison with other studies. The measured velocity fields are also used to calculate velocity structure functions, further confirming the Bolgiano-Obukhov scalings when interpreted with the local turbulent length scales found in the simulations. An extended self-similarity analysis shows that the relative scalings are different for the Kolmogorov and Bolgiano-Obukhov regimes.