BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Roberto Verzicco

Roberto Verzicco appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynphysics.comp-phcond-mat.softnlin.CDphysics.ao-phphysics.geo-phphysics.med-phphysics.app-phphysics.hist-phphysics.soc-ph

Catalog footprint

What is connected

52works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

52 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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)$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 `$+$' 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Deformation and orientation statistics of neutrally buoyant sub-Kolmogorov ellipsoidal droplets in turbulent Taylor-Couette flow

The influence of the underlying flow topology on the shape and size of sub-Kolmogorov droplets dispersed in a turbulent flow is of considerable interest in many industrial and scientific applications. In this work we study the deformation and orientation statistics of sub-Kolmogorov droplets dispersed into a turbulent Taylor-Couette flow. Along with Direct Numerical Simulations (DNS) of the carrier phase and Lagrangian tracking of the dispersed droplets, we solve a phenomenological equation proposed by Maffettone and Minale (\emph{J. Fluid Mech.} 78, 227-241 (1998)) to track the shape evolution and orientation of approximately $10^5$ ellipsoidal droplets. By varying the capillary number $Ca$ and viscosity ratio $\hat μ$ of the droplets we find that the droplets deform more with increasing capillary number $Ca$ and this effect is more pronounced in the boundary layer regions. This indicates that along with a capillary number effect there is also a strong correlation between spatial position and degree of deformation of the droplet. Regardless of the capillary number $Ca$, the major-axis of the ellipsoids tends to align with the stream-wise direction and the extensional strain rate eigen direction in the boundary layer region while the distribution is highly isotropic in the bulk. When the viscosity ratio between the droplet and the carrier fluid is increased we find that there is no preferential stretched axis which is due to the increased influence of rotation over stretching and relaxation. Droplets in high viscosity ratio systems are thus less deformed and oblate (disk-like) as compared to highly deformed prolate (cigar-like) droplets in low viscosity ratio systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Direct numerical simulation of Taylor-Couette flow with grooved walls: torque scaling and flow structure

We present direct numerical simulations of Taylor-Couette flow with grooved walls at a fixed radius ratio $η=r_i/r_o=0.714$ with inner cylinder Reynolds number up to $Re_i=3.76\times10^4$, corresponding to Taylor number up to $Ta=2.15\times10^9$. The grooves are axisymmetric V-shaped obstacles attached to the wall with a tip angle of $90^\circ$. Results are compared to the smooth wall case in order to investigate the effects of grooves on Taylor-Couette flow. We focus on the effective scaling laws for the torque, flow structures, and boundary layers. It is found that, when the groove height is smaller than the boundary layer thickness, the torque is the same as that of the smooth wall cases. With increasing $Ta$, the boundary layer thickness becomes smaller than the groove height. Plumes are ejected from the tips of the grooves and secondary circulations between the latter are formed. This is associated to a sharp increase of the torque and thus the effective scaling law for the torque vs. $Ta$ becomes much steeper. Further increasing $Ta$ does not result in an additional slope increase. Instead, the effective scaling law saturates to the "ultimate" regime effective exponents seen for smooth walls. It is found that even though after saturation the slope is the same as for the smooth wall case, the absolute value of torque is increased, and the more the larger size of the grooves.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Disentangling the origins of torque enhancement through wall roughness in Taylor-Couette turbulence

Direct numerical simulations (DNSs) are performed to analyze the global transport properties of turbulent Taylor-Couette flow with inner rough wall up to Taylor number $Ta=10^{10}$. The dimensionless torque $Nu_ω$ shows an effective scaling of $Nu_ω\propto Ta^{0.42\pm0.01}$, which is steeper than the ultimate regime effective scaling $Nu_ω\propto Ta^{0.38}$ seen for smooth inner and outer walls. It is found that at the inner rough wall, the dominant contribution to the torque comes from the pressure forces on the radial faces of the rough elements; while viscous shear stresses on the rough surfaces contribute little to $Nu_ω$. Thus, the log layer close to the rough wall depends on the roughness length scale, rather than on the viscous length scale. We then separate the torque contributed from the smooth inner wall and the rough outer wall. It is found that the smooth wall torque scaling follows $Nu_s \propto Ta_s^{0.38\pm0.01}$, in excellent agreement with the case where both walls are smooth. In contrast, the rough wall torque scaling follows $Nu_r \propto Ta_r^{0.47\pm0.03}$, very close to the pure ultimate regime scaling $Nu_ω\propto Ta^{1/2}$. The energy dissipation rate at the wall of inner rough cylinder decreases significantly as a consequence of the wall shear stress reduction caused by the flow separation at the rough elements. On the other hand, the latter shed vortices in the bulk that are transported towards the outer cylinder and dissipated. Compared to the purely smooth case, the inner wall roughness renders the system more bulk dominated and thus increases the effective scaling exponent.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Scaling laws and flow structures of double diffusive convection in the finger regime

Direct numerical simulations are conducted for double diffusive convection (DDC) bounded by two parallel plates, with fluid properties similar to the values of seawater. The DDC flow is driven by an unstable salinity difference and stabilized at the same time by a temperature difference. For these conditions the flow can be in the finger regime. We develop scaling laws for three key response parameters of the system: The non-dimensional salinity flux $Nu_S$ mainly depends on the salinity Rayleigh number $Ra_S$, which measures the strength of the salinity difference, and exhibits a very weak dependence on the density ratio $Λ$, which is the ratio of the buoyancy forces induced by two scalar differences. The non-dimensional flow velocity $Re$ and the non-dimensional heat flux $Nu_T$ are dependent on both $Ra_S$ and $Λ$. However, the rescaled Reynolds number $ReΛ^{α^{\rm eff}_u}$ and the rescaled convective heat flux $(Nu_T-1)Λ^{α^{\rm eff}_T}$ depend only on $Ra_S$. The two exponents are dependent on the fluid properties and are determined from the numerical results. Moreover, the behaviors of $Nu_S$ and $ReΛ^{α^{\rm eff}_u}$ agree with the predictions of the Grossmann-Lohse theory which was originally developed for the Rayleigh-Bénard flow. The non-dimensional salt-finger width and the thickness of the velocity boundary layers, after being rescaled by $Λ^{α^{\rm eff}_u/2}$, collapse and obey a similar power-law scaling relation with $Ra_S$. When $Ra_S$ is large enough, salt fingers do not extend from one plate to the other and horizontal zonal flows emerge in the bulk region. We then show that the current scaling strategy can be successfully applied to the experimental results of a heat-copper-ion system~(Hage and Tilgner, Phys. Fluids, 22, 076603, 2010).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Vertically bounded double diffusive convection in the fingering regime: comparing no-slip vs free-slip boundary conditions

Vertically bounded fingering double diffusive convection (DDC) 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 research, 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 fingers and sheet structures near the free-slip boundaries. This counter intuitive result might be relevant for properly estimating and modelling the mixing and entrainment phenomena at free-surfaces and interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A multiple--resolution strategy for Direct Numerical Simulation of scalar turbulence

In this paper a numerical procedure to simulate low diffusivity scalar turbulence is presented. The method consists of using a grid for the advected scalar with a higher spatial resolutions than that of the momentum. The latter usually requires a less refined mesh and integrating both fields on a single grid tailored to the most demanding variable, produces an unnecessary computational overhead. A multiple resolution approach is used also in the time integration in order to maintain the stability of the scalars on the finer grid. The method is the more advantageous the less diffusive the scalar is with respect to momentum, therefore it is particularly well suited for large Prandtl or Schmidt number flows. However, even in the case of equal diffusivities the present procedure gives CPU time and memory occupation savings. The reason is that the absence of the pressure term in the scalar equation leads to much steeper gradients in the scalar field as compared to the velocity field.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A pencil distributed finite difference code for strongly turbulent wall-bounded flows

We present a numerical scheme geared for high performance computation of wall-bounded turbulent flows. The number of all-to-all communications is decreased to only six instances by using a two-dimensional (pencil) domain decomposition and utilizing the favourable scaling of the CFL time-step constraint as compared to the diffusive time-step constraint. As the CFL condition is more restrictive at high driving, implicit time integration of the viscous terms in the wall-parallel directions is no longer required. This avoids the communication of non-local information to a process for the computation of implicit derivatives in these directions. We explain in detail the numerical scheme used for the integration of the equations, and the underlying parallelization. The code is shown to have very good strong and weak scaling to at least 64K cores.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Deformation statistics of sub-Kolmogorov-scale ellipsoidal neutrally buoyant drops in isotropic turbulence

Small droplets in turbulent flows can undergo highly variable deformations and orientational dynamics. For neutrally buoyant droplets smaller than the Kolmogorov scale, the dominant effects from the surrounding turbulent flow arise through Lagrangian time histories of the velocity gradient tensor. Here we study the evolution of representative droplets using a model that includes rotation and stretching effects from the surrounding fluid, and restoration effects from surface tension including a constant droplet volume constraint, while assuming that the droplets maintain an ellipsoidal shape. The model is combined with Lagrangian time histories of the velocity gradient tensor extracted from DNS of turbulence to obtain simulated droplet evolutions. These are used to characterize the size, shape and orientation statistics of small droplets in turbulence. A critical capillary number, $Ca_c$ is identified associated with unbounded growth of one or two of the droplet's semi-axes. Exploiting analogies with dynamics of polymers in turbulence, the $Ca_c$ number can be predicted based on the large deviation theory for the largest Finite Time Lyapunov exponent. Also, for sub-critical $Ca$ the theory enables predictions of the slope of the power-law tails of droplet size distributions in turbulence. For cases when the viscosities of droplet and outer fluid differ in a way that enables vorticity to decorrelate the shape from the straining directions, the large deviation formalism based on the stretching properties of the velocity gradient tensor loses validity and its predictions fail. Even considering the limitations of the assumed ellipsoidal droplet shape, the results highlight the complex coupling between droplet deformation, orientation and the local fluid velocity gradient tensor to be expected when small viscous drops interact with turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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 enviroments. Of great interests are scalers 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-organised 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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\times10^{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 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 increases in fraction of the plume-emitting ('turbulent') regions of the boundary layer.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Salinity transfer in bounded double diffusive convection

The double diffusive convection between two parallel plates is numerically studied for a series of parameters. The flow is driven by the salinity difference and stabilized by the thermal field. Our simulations are directly compared to experiments by Hage and Tilgner (\emph{Phys. Fluids} 22, 076603 (2010)) for several sets of parameters and reasonable agreement is found. This in particular holds for the salinity flux and its dependence on the salinity Rayleigh number. Salt fingers are present in all simulations and extend through the entire height. The thermal Rayleigh number seems to have minor influence on salinity flux but affects the Reynolds number and the morphology of the flow. Next to the numerical calculation, we apply the Grossmann-Lohse theory for Rayleigh-Bénard flow to the current problem without introducing any new coefficients. The theory successfully predicts the salinity flux both with respect to the scaling and even with respect to the absolute value for the numerical and experimental results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The 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 opposite for stress-free sidewalls, where the heat transport increases for lower aspect-ratio. In wider cells the aspect-ratio dependence is observed to disappear for $\text{Ra} \ge 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 horizontal 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Heat transport in boiling turbulent Rayleigh-Bénard convection

Boiling is an extremely effective way to promote heat transfer from a hot surface to a liquid due to several mechanisms many of which are not understood in quantitative detail. An important component of the overall process is that the buoyancy of the bubbles compounds with that of the liquid to give rise to a much enhanced natural convection. In this paper we focus specifically on this enhancement and present a numerical study of the resulting two-phase Rayleigh-Bénard convection process. We make no attempt to model other aspects of the boiling process such as bubble nucleation and detachment. We consider a cylindrical cell with a diameter equal to its height. The cell base and top are held at temperatures above and below the boiling point of the liquid, respectively. By keeping the temperature difference constant and changing the liquid pressure we study the effect of the liquid superheat in a Rayleigh number range that, in the absence of boiling, would be between $2\times10^6$ and $5\times10^9$. We find a considerable enhancement of the heat transfer and study its dependence on the bubble number, 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 hot cell bottom and, as a consequence, the strength of the circulation in the cell is significantly increased. Our results are in general agreement with recent experimental results of Zhong et al., Phys. Rev. Lett. {\bf 102}, 124501 (2009) for boiling Rayleigh-Bénard convection.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sidewall effects in Rayleigh-Bénard convection

We investigate the influence of the temperature boundary conditions at the sidewall on the heat transport in Rayleigh-Bénard (RB) convection using direct numerical simulations. For relatively low Rayleigh numbers Ra the heat transport is higher when the sidewall is isothermal, kept at a temperature $T_c+Δ/2$ (where $Δ$ is the temperature difference between the horizontal plates and $T_c$ the temperature of the cold plate), than when the sidewall is adiabatic. The reason is that in the former case part of the heat current avoids the thermal resistance of the fluid layer by escaping through the sidewall that acts as a short-circuit. For higher Ra the bulk becomes more isothermal and this reduces the heat current through the sidewall. Therefore the heat flux in a cell with an isothermal sidewall converges to the value obtained with an adiabatic sidewall for high enough Ra ($\simeq 10^{10}$). However, when the sidewall temperature deviates from $T_c+Δ/2$ the heat transport at the bottom and top plates is different from the value obtained using an adiabatic sidewall. In this case the difference does not decrease with increasing Ra thus indicating that the ambient temperature of the experimental apparatus can influence the heat transfer. A similar behavior is observed when only a very small sidewall region close to the horizontal plates is kept isothermal, while the rest of the sidewall is adiabatic. The reason is that in the region closest to the horizontal plates the temperature difference between the fluid and the sidewall is highest. This suggests that one should be careful with the placement of thermal shields outside the fluid sample to minimize spurious heat currents.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mixing in thermally stratified nonlinear spin-up with sources and sinks

Stratified spin-up experiments in enclosed cylinders have reported the presence of small pockets of well-mixed fluids but quantitative measurements of the mixedness of the fluid has been lacking. Previous numerical simulations have not addressed these measurements. Here we present numerical simulations that address how the combined effect of spin-up and thermal boundary conditions enhances or hinders mixing of a fluid in a cylinder. Measurements of efficiency of mixing are based on the variance of temperature and explained in terms of the potential energy available. The numerical simulations of the Navier--Stokes equations for the problem with different sets of thermal boundary conditions at the horizontal walls helped shed some light on the physical mechanisms of mixing, for which a clear explanation was lacking.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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 $Γ\equiv D/L = 0.50$ ($D$ and $L$ are the diameter and height respectively). Results from experiment over the Rayleigh number range $4\times 10^{12} \alt Ra \alt 10^{15}$ for a Prandtl number $\Pra \simeq 0.8$ and from direct numerical simulation (DNS) at $Ra = 2 \times 10^{12}$ for $\Pra = 0.7$ are presented. We find that the temperature varies as $A*ln(z/L) + B$ where $z$ is the distance from the bottom or top plate. This is the case in the classical as well as in the ultimate state of RBC. From DNS we find that $A$ in the classical state decreases in the radial direction as the distance from the side wall increases and becomes small near the sample center.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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 $\bar{j}(r)$ in a cylindrical Rayleigh-Bénard cell with aspect ratio Γ=1 for Rayleigh numbers Ra between $2 \times 10^6$ and $2\times 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, $\bar{j}(r)\propto Ra^{γ_j(r)}$. The scaling exponent $γ_j(r)$ decreases monotonically from 0.43 near the axis ($r \approx 0$) to 0.29 close to the side walls ($r \approx D/2$). The effective exponents near the axis and the side wall agree well with the measurements of Shang et al. (Phys.\ Rev.\ Lett.\ \textbf{100}, 244503, 2008) and the predictions of Grossmann and Lohse (Phys.\ Fluids \textbf{16}, 1070, 2004). Extrapolating our results to large Rayleigh number would imply a crossover at $Ra\approx 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 the plume depleted regions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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\times10^{12} and Prandtl numbers Pr=6.4 and Pr=0.7 with the dynamical frame method [Q. Zhou and K.-Q. Xia, Phys. Rev. Lett. 104, 104301 (2010)] thus extending previous results for quasi-2-dimensional systems to 3D systems for the first time. 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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Axially-homogeneous Rayleigh-Benard convection in a cylindrical cell

Previous numerical studies have shown that the "ultimate regime of thermal convection" can be attained in a Rayleigh-Benard cell when the kinetic and thermal boundary layers are eliminated by replacing the walls with periodic boundary conditions (homogeneous Rayleigh-Benard convection). Then, the heat transfer scales like Nu ~ Ra^{1/2} and turbulence intensity as Re ~ Ra^{1/2}, where the Rayleigh number Ra indicates the strength of the driving force. However, experiments never operate in unbounded domains and it is important to understand how confinement might alter the approach to this ultimate regime. Here we consider homogeneous Rayleigh-Benard convection in a laterally confined geometry - a small aspect-ratio vertical cylindrical cell - and show evidence of the ultimate regime as Ra is increased: In spite of the confinement and the resulting kinetic boundary layers, we still find Nu ~ Re ~ Ra^{1/2}. The system supports exact solutions composed of modes of exponentially growing vertical velocity and temperature fields, with Ra as the critical parameter determining the properties of these modes. Counterintuitively, in the low Ra regime, or for very narrow cylinders, the numerical simulations are susceptible to these solutions which can dominate the dynamics and lead to very high and unsteady heat transfer. As Ra is increased, interaction between modes stabilizes the system, evidenced by the increasing homogeneity and reduced fluctuations in the r.m.s. velocity and temperature fields. We also test that physical results become independent of the periodicity length of the cylinder, a purely numerical parameter, as the aspect ratio is increased.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Prandtl and Rayleigh number dependence of heat transport in high Rayleigh number thermal convection

Results from direct numerical simulation for three-dimensional Rayleigh-Bénard convection in samples of aspect ratio $Γ=0.23$ and $Γ=0.5$ up to Rayleigh number $Ra=2\times10^{12}$ are presented. The broad range of Prandtl numbers $0.5<Pr<10$ is considered. In contrast to some experiments, we do not see any increase in $Nu/Ra^{1/3}$, neither due to $Pr$ number effects, nor due to a constant heat flux boundary condition at the bottom plate instead of constant temperature boundary conditions. Even at these very high $Ra$, both the thermal and kinetic boundary layer thicknesses obey Prandtl-Blasius scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Radial boundary layer structure and Nusselt number in Rayleigh-Benard convection

Results from direct numerical simulations for three dimensional Rayleigh-Benard convection in a cylindrical cell of aspect ratio 1/2 and Pr=0.7 are presented. They span five decades of Ra from $2\times 10^6$ to $2 \times10^{11}$. Good numerical resolution with grid spacing $\sim$ Kolmogorov scale turns out to be crucial to accurately calculate the Nusselt number, which is in good agreement with the experimental data by Niemela et al., Nature, 404, 837 (2000). In underresolved simulations the hot (cold) plumes travel further from the bottom (top) plate than in the fully resolved case, because the thermal dissipation close to the sidewall (where the grid cells are largest) is insufficient. We compared the fully resolved thermal boundary layer profile with the Prandtl-Blasius profile. We find that the boundary layer profile is closer to the Prandtl Blasius profile at the cylinder axis than close to the sidewall, due to rising plumes in that region.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Heat transfer mechanisms in bubbly Rayleigh-Benard convection

The heat transfer mechanism in Rayleigh-Benard convection in a liquid with a mean temperature close to its boiling point is studied through numerical simulations with point-like 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 Jacob 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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Prandtl-, Rayleigh-, and Rossby-number dependence of heat transport in turbulent rotating Rayleigh-Benard convection

Experimental and numerical data for the heat transfer as a function of the Rayleigh-, Prandtl-, and Rossby numbers in turbulent rotating Rayleigh-Benard convection are presented. For relatively small $Ra ~ 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 = 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.

0 citations0 reviewsNo code linkedOpen access