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.

Multiple states and transport properties of double-diffusive convection turbulence.

When fluid stratification is induced by the vertical gradients of two scalars with different diffusivities, double-diffusive convection (DDC) may occur and play a crucial role in mixing. Such a process exists in many natural and engineering environments. Especially in the ocean, DDC is omnipresent since the seawater density is affected by temperature and salinity. The most intriguing phenomenon caused by DDC is the thermohaline staircase, i.e., a stack of alternating well-mixed convection layers and sharp interfaces with very large gradients in both temperature and salinity. Here we investigate DDC and thermohaline staircases in the salt finger regime, which happens when warm saltier water lies above cold fresher water and is commonly observed in the (sub)tropic regions. By conducting direct numerical simulations over a large range of parameters, we reveal that multiple equilibrium states exist in fingering DDC and staircases even for the same control parameters. Different states can be established from different initial scalar distributions or different evolution histories of the flow parameters. Hysteresis appears during the transition from a staircase to a single salt finger interface. For the same local density ratio, salt finger interfaces in the single-layer state generate very different fluxes compared to those within staircases. However, the salinity flux for all salt finger interfaces follows the same dependence on the salinity Rayleigh number of the layer and can be described by an effective power law scaling. Our findings have direct applications to oceanic thermohaline staircases.

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.

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.

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 .

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.

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.

From convection rolls to finger convection in double-diffusive turbulence.

Double-diffusive convection (DDC), which is the buoyancy-driven flow with fluid density depending on two scalar components, is ubiquitous in many natural and engineering environments. Of great interests are scalars' transfer rate and flow structures. Here we systematically investigate DDC flow between two horizontal plates, driven by an unstable salinity gradient and stabilized by a temperature gradient. Counterintuitively, when increasing the stabilizing temperature gradient, the salinity flux first increases, even though the velocity monotonically decreases, before it finally breaks down to the purely diffusive value. The enhanced salinity transport is traced back to a transition in the overall flow pattern, namely from large-scale convection rolls to well-organized vertically oriented salt fingers. We also show and explain that the unifying theory of thermal convection originally developed by Grossmann and Lohse for Rayleigh-Bénard convection can be directly applied to DDC flow for a wide range of control parameters (Lewis number and density ratio), including those which cover the common values relevant for ocean flows.

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.

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.

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.

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.

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.

Heat transport in bubbling turbulent convection.

Boiling is an extremely effective way to promote heat transfer from a hot surface to a liquid due to numerous mechanisms, many of which are not understood in quantitative detail. An important component of the overall process is that the buoyancy of the bubble compounds with that of the liquid to give rise to a much-enhanced natural convection. In this article, we focus specifically on this enhancement and present a numerical study of the resulting two-phase Rayleigh-Bénard convection process in a cylindrical cell with a diameter equal to its height. We make no attempt to model other aspects of the boiling process such as bubble nucleation and detachment. The cell base and top are held at temperatures above and below the boiling point of the liquid, respectively. By keeping this difference constant, we study the effect of the liquid superheat in a Rayleigh number range that, in the absence of boiling, would be between 2 × 10(6) and 5 × 10(9). We find a considerable enhancement of the heat transfer and study its dependence on the number of bubbles, the degree of superheat of the hot cell bottom, and the Rayleigh number. The increased buoyancy provided by the bubbles leads to more energetic hot plumes detaching from the cell bottom, and the strength of the circulation in the cell is significantly increased. Our results are in general agreement with recent experiments on boiling Rayleigh-Bénard convection.

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.

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.