Researcher profile

Roberto Verzicco

Roberto Verzicco contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
29works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2026arXiv

Transition from classical to ultimate melting

Melting is omnipresent in nature and technology, with applications ranging from metallurgy, biology, food science, and latent thermal energy storage to oceanography, geophysics, and climate science, and occurring on all scales from sub-millimeter to global scales. The key objective is to understand the rate at which an object melts as a function of its size and of the ambient conditions. To achieve this it is important to be able to extrapolate from small scale experiments and observations to large or even global scales. This is done by scaling laws. However, these are only meaningful if there is no transition from one scaling relation to another one. Here we show, however, that for both fixed and freely-advected melting objects immersed in a turbulent flow a melting transition does exist, namely from slow melting at the small scales to fast melting at the large scales. We do so by controlled melting experiments and corresponding direct numerical simulations, covering four orders of magnitude in scale. The transition corresponds to the transition from a laminar-type boundary layer around the melting object to a turbulent-type boundary layer, i.e., from so-called classical turbulence to ultimate turbulence, with its enhanced transport properties. Our results thus provide a quantitative understanding of the flow physics of the melting process and thereby enable a better extrapolation and prediction of melt rates on large scales such as relevant in geophysics, oceanography, and climate science.

preprint2025arXiv

Collective effects of neighbouring melting ice objects

We present a study on the melting dynamics of neighbouring ice bodies by means of idealised simulations, focusing on collective effects, with the goal of obtaining fundamental insight into how collective interactions influence the melting of ice. Two neighbouring (vertically or horizontally aligned), square-shaped, and equally sized ice objects (size on the order of centimetres) are immersed in quiescent fresh water at a temperature of 20°C. By performing two-dimensional direct numerical simulations, and using the phase-field method to model the phase change, the collective melting of these objects is studied. When the objects are horizontally aligned, no significant influence of the neighbouring object on the melting time is observed. On the other hand, when vertically aligned, though the melting of the upper object is mostly unaffected, the melting time and the morphology of the lower ice body strongly depends on the initial inter-object distance. We report that the melting of the bottom object can be enhanced by more than 10%, or delayed more than 20%, displaying a non-monotonic dependence on the initial object size. We show that this behaviour results from a non-trivial competition between layering of cold fluid, which lowers the heat transfer, and convective flows, which favour mixing and heat transfer. For this melting in mixed convection, we were able to collapse our data onto a single curve.

preprint2024arXiv

Rising and settling 2D cylinders with centre-of-mass offset

Rotational effects are commonly neglected when considering the dynamics of freely rising or settling isotropic particles. Here, we demonstrate that particle rotations play an important role for rising as well as for settling cylinders in situations when mass eccentricity, and thereby a new pendulum timescale, is introduced to the system. We employ two-dimensional simulations to study the motion of a single cylinder in a quiescent unbounded incompressible Newtonian fluid. This allows us to vary the Galileo number, density ratio, relative moment of inertia, and Centre-Of-Mass offset (COM) systematically and beyond what is feasible experimentally. For certain buoyant density ratios, the particle dynamics exhibit a resonance mode, during which the coupling via the Magnus lift force causes a positive feedback between translational and rotational motions. This mode results in vastly different trajectories with significantly larger rotational and translational amplitudes and an increase of the drag coefficient easily exceeding a factor two. We propose a simple model that captures how the occurrence of the COM offset induced resonance regime varies, depending on the other input parameters, specifically the density ratio, the Galileo number, and the relative moment of inertia. Remarkably, depending on the input parameters, resonance can be observed for centre-of-mass offsets as small as a few percent of the particle diameter, showing that the particle dynamics can be highly sensitive to this parameter.

preprint2023arXiv

Double-diffusive transport in multicomponent vertical convection

Motivated by the ablation of vertical ice faces in salt water, we use three-dimensional direct numerical simulations to investigate the heat and salt fluxes in two-scalar vertical convection. For parameters relevant to ice-ocean interfaces in the convection-dominated regime, we observe that the salinity field drives the convection and that heat is essentially transported as a passive scalar. By varying the diffusivity ratio of heat and salt (i.e., the Lewis number $Le$), we identify how the different molecular diffusivities affect the scalar fluxes through the system. Away from the walls, we find that the heat transport is determined by a turbulent Prandtl number of $Pr_t\approx 1$ and that double-diffusive effects are practically negligible. However, the difference in molecular diffusivities plays an important role close to the boundaries. In the (unrealistic) case where salt diffused faster than heat, the ratio of salt-to-heat fluxes would scale as $Le^{1/3}$, consistent with classical nested scalar boundary layers. However, in the realistic case of faster heat diffusion (relative to salt), we observe a transition towards a $Le^{1/2}$ scaling of the ratio of the fluxes. This coincides with the thermal boundary layer width growing beyond the thickness of the viscous boundary layer. We find that this transition is not determined by a critical Lewis number, but rather by a critical Prandtl number $Pr\approx 10$, slightly below that for cold seawater where $Pr=14$. We compare our results to similar studies of sheared and double-diffusive flow under ice shelves, and discuss the implications for fluxes in large-scale ice-ocean models. By coupling our results to ice-ocean interface thermodynamics, we describe how the flux ratio impacts the interfacial salinity, and hence the strength of solutal convection and the ablation rate.

preprint2022arXiv

Enhancing Heat Transport in Multiphase Rayleigh-Bénard Turbulence by Changing the Plate-Liquid Contact Angles

This numerical study presents a simple but extremely effective way to considerably enhance heat transport in turbulent multiphase flows, namely by using oleophilic walls. As a model system, we pick the Rayleigh-Bénard setup, filled with an oil-water mixture. For oleophilic walls, e.g. using only $10\%$ volume fraction of oil in water, we observe a remarkable heat transport enhancement of more than $100\%$ as compared to the pure water case. In contrast, for oleophobic walls, the enhancement is then only about $20\%$ as compared to pure water. The physical explanation of the highly-efficient heat transport for oleophilic walls is that thermal plumes detach from the oil-rich boundary layer and are transported together with the oil phase. In the bulk, the oil-water interface prevents the plumes to mix with the turbulent water bulk. To confirm this physical picture, we show that the minimum amount of oil to achieve the maximum heat transport is set by the volume fraction of the thermal plumes. Our findings provide guidelines of how to optimize heat transport in thermal turbulence. Moreover, the physical insight of how coherent structures are coupled with one phase of a two-phase system has very general applicability for controlling transport properties in other turbulent multiphase flows.

preprint2022arXiv

Turbulent Rayleigh-Bénard convection with bubbles attached to the plate

We numerically investigate turbulent Rayleigh-Bénard convection with gas bubbles attached to the hot plate, mimicking a core feature in electrolysis, catalysis, or boiling. The existence of bubbles on the plate reduces the global heat transfer due to the much lower thermal conductivity of gases as compared to liquids and changes the structure of the boundary layers. The numerical simulations are performed in 3D at Prandtl number Pr=4.38 (water) and Rayleigh number $10^7\le Ra \le10^8$. For simplicity, we assume the bubbles to be equally-sized and having pinned contact lines. We vary the total gas-covered area fraction $0.18 \le S_0 \le 0.62$, the relative bubble height $0.02\le h/H \le0.05$ (where $H$ is the height of the Rayleigh-Bénard cell), the bubble number $40 \le n \le 144$, and their spatial distribution. In all cases, asymmetric temperature profiles are observed, which we quantitatively explain based on the heat flux conservation at each horizontal section. We further propose the idea of using an equivalent single-phase setup to mimic the system with attached bubbles. Based on this equivalence, we can calculate the heat transfer. Without introducing any free parameter, the predictions for the Nusselt number, the upper and lower thermal boundary layer thicknesses, and the mean centre temperature well agree with the numerical results. Finally, our predictions also work for the cases with much larger Pr (e.g. $400$), which indicates that our results can also be applied to predict the mass transfer in water electrolysis with bubbles attached to the electrode surface or in catalysis.

preprint2021arXiv

A fast computational model for the electrophysiology of the whole human heart

In this study we present a novel computational model for unprecedented simulations of the whole cardiac electrophysiology. According to the heterogeneous electrophysiologic properties of the heart, the whole cardiac geometry is decomposed into a set of coupled conductive media having different topology and electrical conductivities: (i) a network of slender bundles comprising a fast conduction atrial network, the AV-node and the ventricular bundles; (ii) the Purkinje network; and (iii) the atrial and ventricular myocardium. The propagation of the action potential in these conductive media is governed by the bidomain/monodomain equations, which are discretized in space using an in-house finite volume method and coupled to three different cellular models, the Courtemanche model [1] for the atrial myocytes, the Stewart model [2] for the Purkinje Network and the ten Tusscher-Panfilov model [3] for the ventricular myocytes. The developed numerical model correctly reproduces the cardiac electrophysiology of the whole human heart in healthy and pathologic conditions and it can be tailored to study and optimize resynchronization therapies or invasive surgical procedures. Importantly, the whole solver is GPU-accelerated using CUDA Fortran providing an unprecedented speedup, thus opening the way for systematic parametric studies and uncertainty quantification analyses.

preprint2021arXiv

An efficient phase-field method for turbulent multiphase flows

With the aim of efficiently simulating three-dimensional multiphase turbulent flows with a phase-field method, we propose a new discretization scheme for the biharmonic term (the 4th-order derivative term) of the Cahn-Hilliard equation. This novel scheme can significantly reduce the computational cost while retaining the same accuracy as the original procedure. Our phase-field method is built on top of a direct numerical simulation solver, named AFiD (www.afid.eu) and open-sourced by our research group. It relies on a pencil distributed parallel strategy and a FFT-based Poisson solver. To deal with large density ratios between the two phases, a pressure split method [1] has been applied to the Poisson solver. To further reduce computational costs, we implement a multiple-resolution algorithm which decouples the discretizations for the Navier-Stokes equations and the scalar equation: while a stretched wall-resolving grid is used for the Navier-Stokes equations, for the Cahn-Hilliard equation we use a fine uniform mesh. The present method shows excellent computational performance for large-scale computation: on meshes up to 8 billion nodes and 3072 CPU cores, a multiphase flow needs only slightly less than 1.5 times the CPU time of the single-phase flow solver on the same grid. The present method is validated by comparing the results to previous studies for the cases of drop deformation in shear flow, including the convergence test with mesh refinement, and breakup of a rising buoyant bubble with density ratio up to 1000. Finally, we simulate the breakup of a big drop and the coalescence of O(10^3) drops in turbulent Rayleigh-Bénard convection at a Rayleigh number of $10^8$, observing good agreement with theoretical results.

preprint2021arXiv

Extended lifetime of respiratory droplets in a turbulent vapour 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 10um) 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 William F. 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 two meters advection range of the respiratory droplets within one second. 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.

preprint2021arXiv

FSEI-GPU: GPU accelerated simulations of the fluid-structure-electrophysiology interaction in the left heart

The reliability of cardiovascular computational models depends on the accurate solution of the hemodynamics, the realistic characterization of the hyperelastic and electric properties of the tissues along with the correct description of their interaction. The resulting fluid-structure-electrophysiology interaction (FSEI) thus requires an immense computational power, usually available in large supercomputing centers, and requires long time to obtain results even if multi-CPU processors are used (MPI acceleration). In recent years, graphics processing units (GPUs) have emerged as a convenient platform for high performance computing, as they allow for considerable reductions of the time-to-solution. This approach is particularly appealing if the tool has to support medical decisions that require solutions within reduced times and possibly obtained by local computational resources. Accordingly, our multi-physics solver has been ported to GPU architectures using CUDA Fortran to tackle fast and accurate hemodynamics simulations of the human heart without resorting to large-scale supercomputers. This work describes the use of CUDA to accelerate the FSEI on heterogeneous clusters, where both the CPUs and GPUs are used in synergistically with minor modifications of the original source code. The resulting GPU accelerated code solves a single heartbeat within a few hours (from three to ten depending on the grid resolution) running on premises computing facility made of few GPU cards, which can be easily installed in a medical laboratory or in a hospital, thus opening towards a systematic computational fluid dynamics (CFD) aided diagnostic.

preprint2021arXiv

Heat transfer in turbulent Rayleigh-Bénard convection within two immiscible fluid layers

We numerically investigate turbulent Rayleigh-Bénard convection within two immiscible fluid layers, aiming to understand how the layer thickness and fluid properties affect the heat transfer (characterized by the Nusselt number $Nu$) in two-layer systems. Both two- and three-dimensional simulations are performed at fixed global Rayleigh number $Ra=10^8$, Prandtl number $Pr=4.38$, and Weber number $We=5$. We vary the relative thickness of the upper layer between $0.01 \le α\le 0.99$ and the thermal conductivity coefficient ratio of the two liquids between $0.1 \le λ_k \le 10$. Two flow regimes are observed: In the first regime at $0.04\leα\le0.96$, convective flows appear in both layers and $Nu$ is not sensitive to $α$. In the second regime at $α\le0.02$ or $α\ge0.98$, convective flow only exists in the thicker layer, while the thinner one is dominated by pure conduction. In this regime, $Nu$ is sensitive to $α$. To predict $Nu$ in the system in which the two layers are separated by a unique interface, we apply the Grossmann-Lohse theory for both individual layers and impose heat flux conservation at the interface. Without introducing any free parameter, the predictions for $Nu$ and for the temperature at the interface well agree with our numerical results and previous experimental data.

preprint2021arXiv

Layering and vertical transport in sheared double diffusive convection in the diffusive regime

A sequence of two and three-dimensional simulations is conducted for the double diffusive convection (DDC) flows in the diffusive regime subjected to an imposed shear. The flow is confined between two horizontal plates which are maintained at different constant temperature, salinity, and different velocity, thus setting up a shear across the flow. The lower plate is fixed at higher temperature and salinity, while the overall (unperturbed) density gradient is statically stable. For a wide range of control parameters, and for sufficiently strong perturbation of the conductive initial state, we find that staircase-like structures spontaneously develop, with relatively well-mixed layers separated by sharp interfaces of enhanced scalar gradient. Such staircases appear to be robust even in the presence of strong shear over very long times, although we typically observe early time coarsening of the number of observed layers. For the same set of control parameters, different asymptotic layered states, with markedly different vertical scalar fluxes, can arise for different initial perturbation structures. The imposed shear does significantly spatio-temporally modify the vertical transport of the various scalars. The flux ratio (i.e., the ratio between the density fluxes due to the total (convective and diffusive) salt flux and the total heat flux) is found, at steady state, to be essentially equal to the square root of the ratio of the salt diffusivity to the thermal diffusivity, consistently with the physical model originally proposed by Linden and Shirtcliffe (1978) and the variational arguments presented by Stern (1982) for unsheared double diffusive convection.

preprint2021arXiv

Off-center gravity induces large-scale flow patterns in spherical Rayleigh-Bénard

Inspired by the hemispherical asymmetry observed in the Earth's inner core, we perform direct numerical simulations to study the effect of the gravity center offset in spherical Rayleigh-Bénard convection. We find that even a minimal shift of the gravity center has a pronounced influence on the flow structures. When the gravity center is shifted towards the South, the co-latitudinal buoyancy component creates an energetic jet on the Northern side of the inner sphere that is directed towards the outer sphere. As a result, a large-scale meridional circulation is formed. However, surprisingly, the global heat flux is not affected by the changes in the large-scale flow organization induced by the gravity center offset. Our results suggest that the hemispherical core asymmetry is key to model the flow phenomena in the Earth's outer core and mantle.

preprint2021arXiv

Strong alignment of prolate ellipsoids in Taylor-Couette flow

We report on the mobility and orientation of finite-size, neutrally buoyant prolate ellipsoids (of aspect ratio $Λ=4$) in Taylor-Couette flow, using interface resolved numerical simulations. The setup consists of a particle-laden flow in between a rotating inner and a stationary outer cylinder. We simulate two particle sizes $\ell/d=0.1$ and $\ell/d=0.2$, $\ell$ denoting the particle major axis and $d$ the gap-width between the cylinders. The volume fractions are $0.01\%$ and $0.07\%$, respectively. The particles, which are initially randomly positioned, ultimately display characteristic spatial distributions which can be categorised into four modes. Modes $(i)$ to $(iii)$ are observed in the Taylor vortex flow regime, while mode ($iv$) encompasses both the wavy vortex, and turbulent Taylor vortex flow regimes. Mode $(i)$ corresponds to stable orbits away from the vortex cores. Remarkably, in a narrow $\textit{Ta}$ range, particles get trapped in the Taylor vortex cores (mode ($ii$)). Mode $(iii)$ is the transition when both modes $(i)$ and $(ii)$ are observed. For mode $(iv)$, particles distribute throughout the domain due to flow instabilities. All four modes show characteristic orientational statistics. We find the particle clustering for mode ($ii$) to be size-dependent, with two main observations. Firstly, particle agglomeration at the core is much higher for $\ell/d=0.2$ compared to $\ell/d=0.1$. Secondly, the $\textit{Ta}$ range for which clustering is observed depends on the particle size. For this mode $(ii)$ we observe particles to align strongly with the local cylinder tangent. The most pronounced particle alignment is observed for $\ell/d=0.2$ around $\textit{Ta}=4.2\times10^5$. This observation is found to closely correspond to a minimum of axial vorticity at the Taylor vortex core ($\textit{Ta}=6\times10^5$) and we explain why.

preprint2020arXiv

Diffusive interaction of multiple surface nanobubbles and nanodroplets: 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 {\it single} surface nanobubble has been theoretically explained, i.e. contact line pinning and gas oversaturation are required to stabilize them against diffusive dissolution [Lohse and Zhang, Phys.\ Rev.\ E 91, 031003 (R) (2015)], here we focus on the {\it collective} diffusive interaction of {\it 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. 18, 1505 (1950)] 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 thus 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.

preprint2020arXiv

Double maxima of angular momentum transport in $η=0.91$ TC turbulence

We use experiments and direct numerical simulations to probe the phase-space of low-curvature Taylor--Couette (TC) flow in the vicinity of the ultimate regime. The cylinder radius ratio is fixed at $η=r_i/r_o=0.91$. Non-dimensional shear drivings (Taylor numbers $\text{Ta}$) in the range $10^7\leq\text{Ta}\leq10^{11}$ are explored for both co- and counter-rotating configurations. In the $\text{Ta}$ range $10^8\leq\text{Ta}\leq10^{10}$, we observe two local maxima of the angular momentum transport as a function of the cylinder rotation ratio, which can be described as either as "co-" and "counter-rotating" due to their location or as "broad" and "narrow" due to their shape. We confirm that the broad peak is accompanied by the strengthening of the large-scale structures, and that the narrow peak appears once the driving (Ta) is strong enough. As first evidenced in numerical simulations by Brauckmann \emph{et al.}~(2016), the broad peak is produced by centrifugal instabilities and that the narrow peak is a consequence of shear instabilities. We describe how the peaks change with $\text{Ta}$ as the flow becomes more turbulent. Close to the transition to the ultimate regime when the boundary layers (BLs) become turbulent, the usual structure of counter-rotating Taylor vortex pairs breaks down and stable unpaired rolls appear locally. We attribute this state to changes in the underlying roll characteristics during the transition to the ultimate regime. Further changes in the flow structure around $\text{Ta}\approx10^{10}$ cause the broad peak to disappear completely and the narrow peak to move. This second transition is caused when the regions inside the BLs which are locally smooth regions disappear and the whole boundary layer becomes active.

preprint2020arXiv

Flow organization in laterally unconfined Rayleigh-Bénard turbulence

We investigate the large-scale circulation (LSC) of turbulent Rayleigh-Bénard convection in a large box of aspect ratio $Γ=32$ for Rayleigh numbers up to $Ra=10^9$ and at a fixed Prandtl number $Pr=1$. A conditional averaging technique allows us to extract statistics of the LSC even though the number and the orientation of the structures vary throughout the domain. We find that various properties of the LSC obtained here, such as the wall-shear stress distribution, the boundary layer thicknesses and the wind Reynolds number, do not differ significantly from results in confined domains ($Γ\approx 1$). This is remarkable given that the size of the structures (as measured by the width of a single convection roll) more than doubles at the highest $Ra$ as the confinement is removed. An extrapolation towards the critical shear Reynolds number of $Re_s^{\textrm{crit}} \approx 420$, at which the boundary layer (BL) typically becomes turbulent, predicts that the transition to the ultimate regime is expected at $Ra_{\textrm{crit}} \approx \mathcal{O}(10^{15})$ in unconfined geometries. This result is in line with the Göttingen experimental observations. Furthermore, we confirm that the local heat transport close to the wall is highest in the plume impacting region, where the thermal BL is thinnest, and lowest in the plume emitting region, where the thermal BL is thickest. This trend, however, weakens with increasing $Ra$.

preprint2020arXiv

From Rayleigh-Bénard convection to porous-media convection: how porosity affects heat transfer and flow structure

We perform a numerical study of the heat transfer and flow structure of Rayleigh-Bénard (RB) convection in (in most cases regular) porous media, which are comprised of circular, solid obstacles located on a square lattice. This study is focused on the role of porosity $ϕ$ in the flow properties during the transition process from the traditional RB convection with $ϕ=1$ (so no obstacles included) to Darcy-type porous-media convection with $ϕ$ approaching 0. Simulations are carried out in a cell with unity aspect ratio, for the Rayleigh number $Ra$ from $10^5$ to $10^{10}$ and varying porosities $ϕ$, at a fixed Prandtl number $Pr=4.3$, and we restrict ourselves to the two dimensional case. For fixed $Ra$, the Nusselt number $Nu$ is found to vary non-monotonously as a function of $ϕ$; namely, with decreasing $ϕ$, it first increases, before it decreases for $ϕ$ approaching 0. The non-monotonous behaviour of $Nu(ϕ)$ originates from two competing effects of the porous structure on the heat transfer. On the one hand, the flow coherence is enhanced in the porous media, which is beneficial for the heat transfer. On the other hand, the convection is slowed down by the enhanced resistance due to the porous structure, leading to heat transfer reduction. For fixed $ϕ$, depending on $Ra$, two different heat transfer regimes are identified, with different effective power-law behaviours of $Nu$ vs $Ra$, namely, a steep one for low $Ra$ when viscosity dominates, and the standard classical one for large $Ra$. The scaling crossover occurs when the thermal boundary layer thickness and the pore scale are comparable. The influences of the porous structure on the temperature and velocity fluctuations, convective heat flux, and energy dissipation rates are analysed, further demonstrating the competing effects of the porous structure to enhance or reduce the heat transfer.

preprint2020arXiv

From zonal flow to convection rolls in Rayleigh-Bénard convection with free-slip plates

Rayleigh-Bénard (RB) convection with free-slip plates and horizontally periodic boundary conditions is investigated using direct numerical simulations. Two configurations are considered, one is two-dimension (2D) RB convection and the other one three-dimension (3D) RB convection with a rotating axis parallel to the plate. We explore the parameter range of Rayleigh numbers Ra from $10^7 to $10^9$ and Prandtl numbers $Pr$ from $1$ to $100$. We show that zonal flow, which was observed, for example, by Goluskin \emph{et al}. \emph{J. Fluid. Mech.} 759, 360-385 (2014) for $Γ=2$, is only stable when $Γ$ is smaller than a critical value, which depends on $Ra$ and $Pr$. With increasing $Γ$, we find a second regime in which both zonal flow and different convection roll states can be statistically stable. For even larger $Γ$, in a third regime, only convection roll states are statistically stable and zonal flow is not sustained. For the 3D simulations, we fix $Ra=10^7$ and $Pr=0.71$, and compare the flow for $Γ=8$ and $Γ= 16$. We demonstrate that with increasing aspect ratio $Γ$, zonal flow, which was observed for small $Γ=2π$ by von Hardenberg \emph{et al}. \emph{Phys. Rev. Lett.} 15, 134501 (2015), completely disappears for $Γ=16$. For such large $Γ$ only convection roll states are statistically stable. In between, here for medium aspect ratio $Γ= 8$, the convection roll state and the zonal flow state are both statistically stable. What state is taken depends on the initial conditions, similarly as we found for the 2D case.

preprint2020arXiv

Multiple states in turbulent large-aspect ratio thermal convection: What determines the number of convection rolls?

Recent findings suggest that wall-bounded turbulent flow 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 \sim 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 instability 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 \le Γ_r \le 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.

preprint2020arXiv

Subcritical behaviour in double diffusive convection within the diffusive regime

We conduct two- and three-dimensional simulations for double diffusive convection in the diffusive regime, where the fluid flow is driven by a destabilizing temperature gradient and stabilized by a stably stratified salinity gradient. We study how the heat flux, Reynolds number, and flow structures change with the density ratio $Λ$, which is the ratio of the buoyancy force induced by the salinity gradient to that by the temperature gradient. When $Λ$ increases from zero, the flow first behaves similarly as in pure Rayleigh-Bénard (RB) convection, both with respect to flow structure and to heat transport. The linear stability analysis of Baines & Gill (J. Fluid Mech., vol. 37, 1969, pp. 289-306) had estimated the critical density ratio $Λ_c$, above which the flow becomes stable. However, here we show that by using a large-scale circulation as initial condition (rather than the linear profiles assumed in the linear stability analysis), DDC in the diffusive regime can exhibit subcritical behaviour when $Λ> Λ_c$, i.e., coexistence of states at the same control parameters. Even though the density ratio becomes thousands times that of the critical value $Λ_c$, there is still convection with strongly enhanced heat transfer properties compared to the pure conduction case. We reveal the corresponding flow structures and find an unstably-stratified region sandwiched between two stably-stratified layers. Our results demonstrate the importance of the initial condition for DDC in the diffusive regime, especially in the situation of a large density ratio, which occurs in high-latitude ocean regions.

preprint2020arXiv

The effect of Prandtl number on turbulent sheared thermal convection

In turbulent wall sheared thermal convection, there are three different flow regimes, depending on the relative relevance of thermal forcing and wall shear. In this paper we report the results of direct numerical simulations of such sheared Rayleigh-Bénard convection, at fixed Rayleigh number $Ra=10^6$, varying the wall Reynolds number in the range $0 \leq Re_w \leq 4000$ and Prandtl number $0.22 \leq Pr \leq 4.6$, extending our prior work by Blass et al. (2020), where $Pr$ was kept constant at unity and the thermal forcing ($Ra$) varied. We cover a wide span of bulk Richardson numbers $0.014 \leq Ri \leq 100$ and show that the Prandtl number strongly influences the morphology and dynamics of the flow structures. In particular, at fixed $Ra$ and $Re_w$, a high Prandtl number causes stronger momentum transport from the walls and therefore yields a greater impact of the wall shear on the flow structures, resulting in an increased effect of $Re_w$ on the Nusselt number. Furthermore, we analyse the thermal and kinetic boundary layer thicknesses and relate their behaviour to the resulting flow regimes. For the largest shear rates and $Pr$ numbers, we observe the emergence of a Prandtl-von Karman log-layer, signalling the onset of turbulent dynamics in the boundary layer. Finally, our results allow to extend the Grossmann-Lohse theory for heat transport in Rayleigh-Bénard convection to the sheared case, universally describing $Nu(Ra,Pr,Re_w)$.

preprint2020arXiv

Two-layer Thermally Driven Turbulence: Mechanisms for Interface Breakup

It is commonly accepted that the breakup criteria of drops or bubbles in turbulence is governed by surface tension and inertia. However, also {\it{buoyancy}} can play an important role at breakup. In order to better understand this role, here we numerically study Rayleigh-Bénard convection for two immiscible fluid layers, in order to identify the effects of buoyancy on interface breakup. We explore the parameter space spanned by the Weber number $5\leq We \leq 5000$ (the ratio of inertia to surface tension) and the density ratio between the two fluids $0.001 \leq Λ\leq 1$, at fixed Rayleigh number $Ra=10^8$ and Prandtl number $Pr=1$. At low $We$, the interface undulates due to plumes. When $We$ is larger than a critical value, the interface eventually breaks up. Depending on $Λ$, two breakup types are observed: The first type occurs at small $Λ\ll 1$ (e.g. air-water systems) when local filament thicknesses exceed the Hinze length scale. The second, strikingly different, type occurs at large $Λ$ with roughly $0.5 < Λ\le 1$ (e.g. oil-water systems): The layers undergo a periodic overturning caused by buoyancy overwhelming surface tension. For both types the breakup criteria can be derived from force balance arguments and show good agreement with the numerical results.

preprint2020arXiv

What rotation rate maximizes heat transport in rotating Rayleigh-Bénard convection with Prandtl number larger than one?

The heat transfer and flow structure in rotating Rayleigh-Bénard convection are strongly influenced by the Rayleigh ($Ra$), Prandtl ($Pr$), and Rossby ($Ro$) number. For $Pr\gtrsim 1$ and intermediate rotation rates, the heat transfer is increased compared to the non-rotating case. We find that the regime of increased heat transfer is subdivided into a low and a high $Ra$ number regime. For $Ra\lesssim 5\times10^8$ the heat transfer at a given $Ra$ and $Pr$ is highest at an optimal rotation rate, at which the thickness of the viscous and thermal boundary layer is about equal. From the scaling relations of the thermal and viscous boundary layer thicknesses, we derive that the optimal rotation rate scales as $1/Ro_\mathrm{opt} \approx 0.12 Pr^{1/2}Ra^{1/6}$. In the low $Ra$ regime the heat transfer is similar in a periodic domain and cylindrical cells with different aspect ratios, i.e.\ the ratio of diameter to height. This is consistent with the view that the vertically aligned vortices are the dominant flow structure. For $Ra\gtrsim 5\times10^8$ the above scaling for the optimal rotation rate does not hold anymore. It turns out that in the high $Ra$ regime, the flow structures at the optimal rotation rate are very different than for lower $Ra$. Surprisingly, the heat transfer in the high $Ra$ regime differs significantly for a periodic domain and cylindrical cells with different aspect ratios, which originates from the sidewall boundary layer dynamics and the corresponding secondary circulation.

preprint2019arXiv

Controlling secondary flow in Taylor-Couette turbulence through spanwise-varying roughness

Highly turbulent Taylor-Couette flow with spanwise-varying roughness is investigated experimentally and numerically (direct numerical simulations (DNS) with an immersed boundary method (IBM)) to determine the effects of the spacing and axial width $s$ of the spanwise varying roughness on the total drag and {on} the flow structures. We apply sandgrain roughness, in the form of alternating {rough and smooth} bands to the inner cylinder. Numerically, the Taylor number is $\mathcal{O}(10^9)$ and the roughness width is varied between $0.47\leq \tilde{s}=s/d \leq 1.23$, where $d$ is the gap width. Experimentally, we explore $\text{Ta}=\mathcal{O}(10^{12})$ and $0.61\leq \tilde s \leq 3.74$. For both approaches the radius ratio is fixed at $η=r_i/r_o = 0.716$, with $r_i$ and $r_o$ the radius of the inner and outer cylinder respectively. We present how the global transport properties and the local flow structures depend on the boundary conditions set by the roughness spacing $\tilde{s}$. Both numerically and experimentally, we find a maximum in the angular momentum transport as function of $\tilde s$. This can be atributed to the re-arrangement of the large-scale structures triggered by the presence of the rough stripes, leading to correspondingly large-scale turbulent vortices.

preprint2019arXiv

Convection-dominated dissolution for single and multiple immersed sessile droplets

We numerically investigate both single and multiple droplet dissolution with droplets consisting of lighter liquid dissolving in a denser host liquid. The significance of buoyancy is quantified by the Rayleigh number Ra which is the buoyancy force over the viscous damping force. In this study, Ra spans almost four decades from 0.1 to 400. We focus on how the mass flux, characterized by the Sherwood number Sh, and the flow morphologies depend on Ra. For single droplet dissolution, we first show the transition of the Sh(Ra) scaling from a constant value to $Sh\sim Ra^{1/4}$, which confirms the experimental results by Dietrich et al. (J. Fluid Mech., vol. 794, 2016, pp. 45--67). The two distinct regimes, namely the diffusively- and the convectively-dominated regime, exhibit different flow morphologies: when Ra>=10, a buoyant plume is clearly visible which contrasts sharply to the pure diffusion case at low Ra. For multiple droplet dissolution, the well-known shielding effect comes into play at low Ra so that the dissolution rate is slower as compared to the single droplet case. However, at high Ra, convection becomes more and more dominant so that a collective plume enhances the mass flux, and remarkably the multiple droplets dissolve faster than a single droplet. This has also been found in the experiments by Laghezza et al. (Soft Matter, vol. 12, 2016, pp. 5787--5796). We explain this enhancement by the formation of a single, larger plume rather than several individual plumes. Moreover, there is an optimal Ra at which the enhancement is maximized, because the single plume is narrower at larger Ra, which thus hinders the enhancement. Our findings demonstrate a new mechanism in collective droplet dissolution, which is the merging of the plumes, that leads to non-trivial phenomena, contrasting the shielding effect.

preprint2019arXiv

Flow organization and heat transfer in turbulent wall sheared thermal convection

We perform direct numerical simulations of wall sheared Rayleigh-Bénard (RB) convection for Rayleigh numbers up to $Ra=10^8$, Prandtl number unity, and wall shear Reynolds numbers up to $Re_w=10000$. Using the Monin-Obukhov length $L_{MO}$ we identify three different flow states, a buoyancy dominated regime ($L_{MO} \lesssim λ_θ$; with $λ_θ$ the thermal boundary layer thickness), a transitional regime ($0.5H \gtrsim L_{MO} \gtrsim λ_θ$; with $H$ the height of the domain), and a shear dominated regime ($L_{MO} \gtrsim 0.5H$). In the buoyancy dominated regime the flow dynamics are 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 $Ra$ 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. For $L_{MO} \gtrsim 0.5H$ the Nusselt number $Nu$ effectively scales as $Nu \sim Ra^α$, with $α\ll 1/3$ while we find $α\simeq 0.31$ in the buoyancy dominated regime. In the transitional regime the effective scaling exponent is $α> 1/3$, 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.

preprint2018arXiv

Direct numerical simulations of Taylor--Couette turbulence: the effect of sand grain roughness

Progress in roughness research, mapping any given roughness geometry to its fluid dynamic behaviour, has been hampered by the lack of accurate and direct measurements of skin-friction drag, especially in open systems. The Taylor--Couette (TC) system has the benefit of being a closed system, but its potential for characterizing irregular, realistic, 3-D roughness has not been previously considered in depth. Here, we present direct numerical simulations (DNSs) of TC turbulence with sand grain roughness mounted on the inner cylinder. The model proposed by Scotti (\textit{Phys. Fluids}, vol. 18, 031701, 2006) has been improved to simulate a random rough surface of monodisperse sand grains, which is characterized by the equivalent sand grain height $k_s$. Taylor numbers range from $Ta = 1.0\times 10^7$(corresponding to $Re_τ= 82$) to $Ta = 1.0\times 10^9$($Re_τ= 635$). We focus on the influence of the roughness height $k_s^+$ in the transitionally rough regime, through simulations of TC with rough surfaces, ranging from $k_s^+=5$ up to $k_s^+ = 92$, where the superscript `$+$&#39; indicates non-dimensionalization in viscous units. We find that the downwards shift of the logarithmic layer, due to transitionally rough sand grains exhibits remarkably similar behavior to that of the Nikuradse (\textit{VDI-Forschungsheft} 361, 1933) data of sand grain roughness in pipe flow, regardless of the Taylor number dependent constants of the logarithmic layer.