Researcher profile

Herman J. H. Clercx

Herman J. H. Clercx contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
1topics
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

16 published item(s)

preprint2022arXiv

Discontinuous Transitions Towards Vortex Condensates in Buoyancy-Driven Rotating Turbulence

Using direct numerical simulations of rotating Rayleigh-Bénard convection, we explore the transitions between turbulent states from a 3D flow state towards a quasi-2D condensate known as the large-scale vortex (LSV). We vary the Rayleigh number $Ra$ as control parameter and study the system response (strength of the LSV) in terms of order parameters assessing the energetic content in the flow and the upscale energy flux. By sensitively probing the boundaries of the region of existence of the LSV in parameter space, we find discontinuous transitions and we identify the presence of a hysteresis loop as well as memoryless abrupt growth dynamics. We show furthermore that the creation of the condensate state coincides with a discontinuous transition of the energy transport into the largest mode of the system.

preprint2022arXiv

Regime transitions in stratified shear flows: the link between horizontal and inclined ducts

We present the analytical solution for the two-dimensional velocity and density fields within an approximation for laminar stratified inclined duct (SID) flows where diffusion dominates over inertia in the along-channel momentum equation but it is negligible in the density transport equation. We refer to this approximation as the hydrostatic/gravitational/viscous in momentum and advective in density (HGV-A) approximation due to the leading balances in the governing equations. The analytical solution is valid for laminar flows in a two-layer configuration in the limit of long ducts. Under such conditions, the non-dimensional volume flux is given by the Froude number $Fr^* =Re_g/(A\,K)$ with $Re_g$ the gravitational Reynolds number, $A$ the aspect ratio of the duct, and $K$ a geometrical parameter that depends on the tilt of the duct and is obtained from the analytical solution. The analytical solution in the HGV-A approximation is validated against results from laboratory experiments, and allows us to gain new insight into the dynamics and properties of SID flows. Most importantly, constant values of $Fr^*$ describe, in both horizontal and inclined ducts, the transitions between increasingly turbulent flow regimes: from laminar flow, to interfacial waves, to intermittent turbulence and sustained turbulence.

preprint2020arXiv

Turbulent rotating convection confined in a slender cylinder: the sidewall circulation

Recent studies of rotating Rayleigh-Bénard convection at high rotation rates and strong thermal forcing have shown a significant discrepancy in total heat transport between experiments on a confined cylindrical domain on the one hand and simulations on a laterally unconfined periodic domain on the other. This paper addresses this discrepancy using direct numerical simulations on a cylindrical domain. An analysis of the flow field reveals a region of enhanced convection near the wall, the sidewall circulation. The sidewall circulation rotates slowly within the cylinder in anticyclonic direction. It has a convoluted structure, illustrated by mean flow fields in horizontal cross-sections of the flow where instantaneous snapshots are compensated for the orientation of the sidewall circulation before averaging. Through separate analysis of the sidewall region and the inner bulk flow, we find that for higher values of the thermal forcing the heat transport in the inner part of the cylindrical domain, outside the sidewall circulation region, coincides with the heat transport on the unconfined periodic domain. Thus the sidewall circulation accounts for the differences in heat transfer between the two considered domains, while in the bulk the turbulent heat flux is the same as that of a laterally unbounded periodic domain. Therefore, experiments, with their inherent confinement, can still provide turbulence akin to the unbounded domains of simulations, and at more extreme values of the governing parameters for thermal forcing and rotation. We also provide experimental evidence for the existence of the sidewall circulation that is in close agreement with the simulation results.

preprint2015arXiv

Clustering of vertically constrained passive particles in homogeneous, isotropic turbulence

We analyze the dynamics of small particles vertically confined, by means of a linear restoring force, to move within a horizontal fluid slab in a three-dimensional (3D) homogeneous isotropic turbulent velocity field. The model that we introduce and study is possibly the simplest description for the dynamics of small aquatic organisms that, due to swimming, active regulation of their buoyancy, or any other mechanism, maintain themselves in a shallow horizontal layer below the free surface of oceans or lakes. By varying the strength of the restoring force, we are able to control the thickness of the fluid slab in which the particles can move. This allows us to analyze the statistical features of the system over a wide range of conditions going from a fully 3D incompressible flow (corresponding to the case of no confinement) to the extremely confined case corresponding to a two-dimensional slice. The background 3D turbulent velocity field is evolved by means of fully resolved direct numerical simulations. Whenever some level of vertical confinement is present, the particle trajectories deviate from that of fluid tracers and the particles experience an effectively compressible velocity field. Here, we have quantified the compressibility, the preferential concentration of the particles, and the correlation dimension by changing the strength of the restoring force. The main result is that there exists a particular value of the force constant, corresponding to a mean slab depth approximately equal to a few times the Kolmogorov length scale, that maximizes the clustering of the particles.

preprint2013arXiv

Geometrical statistics of the vorticity vector and the strain rate tensor in rotating turbulence

We report results on the geometrical statistics of the vorticity vector obtained from experiments in electromagnetically forced rotating turbulence. A range of rotation rates $Ω$ is considered, from non-rotating to rapidly rotating turbulence with a maximum background rotation rate of $Ω=5$ rad/s (with Rossby number much smaller than unity). Typically, in our experiments ${\rm{Re}}_λ\approx 100$. The measurement volume is located in the centre of the fluid container above the bottom boundary layer, where the turbulent flow can be considered locally statistically isotropic and horizontally homogeneous for the non-rotating case, see van Bokhoven et al., Phys. Fluids 21, 096601 (2009). Based on the full set of velocity derivatives, measured in a Lagrangian way by 3D Particle Tracking Velocimetry, we have been able to quantify statistically the effect of system rotation on several flow properties. The experimental results show how the turbulence evolves from almost isotropic 3D turbulence ($Ω\lesssim 0.2$ rad/s) to quasi-2D turbulence ($Ω\approx 5.0$ rad/s) and how this is reflected by several statistical quantities. In particular, we have studied the orientation of the vorticity vector with respect to the three eigenvectors of the local strain rate tensor and with respect to the vortex stretching vector. Additionally, we have quantified the role of system rotation on the self-amplification terms of the enstrophy and strain rate equations and the direct contribution of the background rotation on these evolution equations. The main effect is the strong reduction of extreme events and related (strong) reduction of the skewness of PDFs of several quantities such as, for example, the intermediate eigenvalue of the strain rate tensor and the enstrophy self-amplification term.

preprint2013arXiv

Lagrangian acceleration of passive tracers in statistically-steady rotating turbulence

The statistical properties of the Lagrangian acceleration vector of passive tracers in statistically-steady rotating turbulence is studied by Particle Tracking Velocimetry. Direct effects of the background rotation are the suppression of high-acceleration events parallel to the (vertical) rotation axis, the enhancement of high-acceleration events for the horizontal acceleration, and the strong amplification of the auto-correlation of the acceleration component perpendicular to both the rotation vector Ω and local velocity vector u. The auto-correlation of the acceleration component in the plane set up by Ω and u is only mildly enhanced.

preprint2013arXiv

Lagrangian velocity auto-correlations in statistically-steady rotating turbulence

Lagrangian statistics of passive tracers in rotating turbulence is investigated by Particle Tracking Velocimetry. A confined and steadily-forced turbulent flow is subjected to five different rotation rates. The PDFs of the velocity components clearly reveal the anisotropy induced by background rotation. Although the statistical properties of the horizontal turbulent flow field are approximately isotropic, in agreement with previously reported results by van Bokhoven and coworkers [Phys. Fluids 21, 096601 (2009)], the velocity component parallel to the (vertical) rotation axis gets strongly reduced (compared to the horizontal ones) while the rotation is increased. The auto-correlation coefficients of all three components are progressively enhanced for increasing rotation rates, although the vertical one shows a tendency to decrease for slow rotation rates. The decorrelation is approximately exponential. Lagrangian data compare favourably with previously reported Eulerian data for horizontal velocity components, but show different behaviour for the vertical velocity components at higher rotation rates.

preprint2012arXiv

Breakdown of the large-scale wind in Γ=1/2 rotating Rayleigh-Bénard flow

Experiments and simulations of rotating Rayleigh-Bénard convection in cylindrical samples have revealed an increase in heat transport with increasing rotation rate. This heat transport enhancement is intimately related to a transition in the turbulent flow structure from a regime dominated by a large-scale circulation (LSC), consisting of a single convection roll, at no or weak rotation to a regime dominated by vertically-aligned vortices at strong rotation. For a sample with an aspect ratio Γ= D/L = 1 (D is the sample diameter and L its height) the transition between the two regimes is indicated by a strong decrease in the LSC strength. In contrast, for Γ= 1/2 Weiss and Ahlers [J. Fluid Mech. {\bf{688}}, 461 (2011)] revealed the presence of a LSC-like sidewall temperature signature beyond the critical rotation rate. They suggested that this might be due to the formation of a two-vortex state, in which one vortex extends vertically from the bottom into the sample interior and brings up warm fluid, while another vortex brings down cold fluid from the top; this flow field would yield a sidewall temperature signature similar to that of the LSC. Here we show by direct numerical simulations for Γ= 1/2 and parameters that allow direct comparison with experiment that the spatial organization of the vertically-aligned vortical structures in the convection cell do indeed yield (for the time average) a sinusoidal variation of the temperature near the sidewall, as found in the experiment. This is also the essential and non-trivial difference with the Γ=1 sample, where the vertically-aligned vortices are distributed randomly.

preprint2011arXiv

Effect of aspect-ratio on vortex distribution and heat transfer in rotating Rayleigh-Bénard convection

Numerical and experimental data for the heat transfer as function of the Rossby number Ro in turbulent rotating Rayleigh-Bénard convection are presented for Prandtl number Pr=4.38 and Rayleigh number $Ra=2.91\times10^8$ up to $Ra=4.52\times10^9$. The aspect ratio Γ= D/L, where L is the height and D the diameter of the cylindrical sample, is varied between Γ=0.5 and Γ=2.0. Without rotation, where the aspect ratio influences the global large scale circulation, we see a small aspect-ratio dependence in the Nusselt number for $Ra=2.91\times10^8$. However, for stronger rotation, i.e. $1/Ro \gg 1/Ro_c$, the heat transport becomes independent of the aspect-ratio. We interpret this finding as follows: In the rotating regime the heat is mainly transported by vertically-aligned vortices. Since the vertically-aligned vortices are local, the aspect ratio has a negligible effect on the heat transport in the rotating regime. Indeed, a detailed analysis of vortex statistics shows that the fraction of the horizontal area that is covered by vortices is independent of the aspect ratio when $1/Ro \gg 1/Ro_c$. In agreement with the results of Weiss et al. Phys. Rev. Lett., vol 105, 224501 (2010) we find a vortex-depleted area close to the sidewall. Here, we in addition show that there is also an area with enhanced vortex concentration next to the vortex-depleted edge region and that the absolute widths of both regions are independent of the aspect ratio.

preprint2011arXiv

Effect of Plumes on Measuring the Large Scale Circulation in Turbulent Rayleigh-Bénard Convection

We studied the properties of the large-scale circulation (LSC) in turbulent Rayleigh-Bénard (RB) convection by using results from direct numerical simulations in which we placed a large number of numerical probes close to the sidewall. The LSC orientation is determined by either a cosine or a polynomial fit to the azimuthal temperature or azimuthal vertical velocity profile measured with the probes. We study the LSC in Γ=D/L=1/2 and Γ=1 samples, where D is the diameter and L the height. For Pr=6.4 in an aspect ratio Γ=1 sample at $Ra=1\times10^8$ and $5\times10^8$ the obtained LSC orientation is the same, irrespective of whether the data of only 8 or all 64 probes per horizontal plane are considered. In a Γ=1/2 sample with $Pr=0.7$ at $Ra=1\times10^8$ the influence of plumes on the azimuthal temperature and azimuthal vertical velocity profiles is stronger. Due to passing plumes and/or the corner flow the apparent LSC orientation obtained using a cosine fit can result in a misinterpretation of the character of the large-scale flow. We introduce the relative LSC strength, which we define as the ratio between the energy in the first Fourier mode and the energy in all modes that can be determined from the azimuthal temperature and azimuthal vertical velocity profiles, to further quantify the large-scale flow. For $Ra=1\times10^8$ we find that this relative LSC strength is significantly lower in a Γ=1/2 sample than in a Γ=1 sample, reflecting that the LSC is much more pronounced in a Γ=1 sample than in a Γ=1/2 sample. The determination of the relative LSC strength can be applied directly to available experimental data to study high Rayleigh number thermal convection and rotating RB convection.

preprint2011arXiv

Finite-size effects lead to supercritical bifurcations in turbulent rotating Rayleigh-Bénard convection

In turbulent thermal convection in cylindrical samples of aspect ratio Γ= D/L (D is the diameter and L the height) the Nusselt number Nu is enhanced when the sample is rotated about its vertical axis, because of the formation of Ekman vortices that extract additional fluid out of thermal boundary layers at the top and bottom. We show from experiments and direct numerical simulations that the enhancement occurs only above a bifurcation point at a critical inverse Rossby number $1/\Ro_c$, with $1/\Ro_c \propto 1/Γ$. We present a Ginzburg-Landau like model that explains the existence of a bifurcation at finite $1/\Ro_c$ as a finite-size effect. The model yields the proportionality between $1/\Ro_c$ and $1/Γ$ and is consistent with several other measured or computed system properties.

preprint2011arXiv

The role of Stewartson and Ekman layers in turbulent rotating Rayleigh-Bénard convection

When the classical Rayleigh-Bénard (RB) system is rotated about its vertical axis roughly three regimes can be identified. In regime I (weak rotation) the large scale circulation (LSC) is the dominant feature of the flow. In regime II (moderate rotation) the LSC is replaced by vertically aligned vortices. Regime III (strong rotation) is characterized by suppression of the vertical velocity fluctuations. Using results from experiments and direct numerical simulations of RB convection for a cell with a diameter-to-height aspect ratio equal to one at $Ra \sim 10^8-10^9$ ($Pr=4-6$) and $0 \lesssim 1/Ro \lesssim 25$ we identified the characteristics of the azimuthal temperature profiles at the sidewall in the different regimes. In regime I the azimuthal wall temperature profile shows a cosine shape and a vertical temperature gradient due to plumes that travel with the LSC close to the sidewall. In regime II and III this cosine profile disappears, but the vertical wall temperature gradient is still observed. It turns out that the vertical wall temperature gradient in regimes II and III has a different origin than that observed in regime I. It is caused by boundary layer dynamics characteristic for rotating flows, which drives a secondary flow that transports hot fluid up the sidewall in the lower part of the container and cold fluid downwards along the sidewall in the top part.

preprint2010arXiv

Boundary layers in rotating weakly turbulent Rayleigh-Benard convection

The effect of rotation on the boundary layers (BLs) in a Rayleigh-Benard (RB) system at a relatively low Rayleigh number, i.e. $Ra = 4\times10^7$, is studied for different Pr by direct numerical simulations and the results are compared with laminar BL theory. In this regime we find a smooth onset of the heat transfer enhancement as function of increasing rotation rate. We study this regime in detail and introduce a model based on the Grossmann-Lohse theory to describe the heat transfer enhancement as function of the rotation rate for this relatively low Ra number regime and weak background rotation $Ro\gtrsim 1$. The smooth onset of heat transfer enhancement observed here is in contrast to the sharp onset observed at larger $Ra \gtrsim 10^8$ by Stevens {\it{et al.}} [Phys. Rev. Lett. {\bf{103}}, 024503, 2009], although only a small shift in the Ra-Ro-Pr phase space is involved.

preprint2009arXiv

Optimal Prandtl number for heat transfer in rotating Rayleigh-Benard convection

Numerical data for the heat transfer as a function of the Prandtl (Pr) and Rossby (Ro) numbers in turbulent rotating Rayleigh-Benard convection are presented for Rayleigh number Ra = 10^8. When Ro is fixed the heat transfer enhancement with respect to the non-rotating value shows a maximum as function of Pr. This maximum is due to the reduced efficiency of Ekman pumping when Pr becomes too small or too large. When Pr becomes small, i.e. for large thermal diffusivity, the heat that is carried by the vertical vortices spreads out in the middle of the cell, and Ekman pumping thus becomes less efficient. For higher Pr the thermal boundary layers (BLs) are thinner than the kinetic BLs and therefore the Ekman vortices do not reach the thermal BL. This means that the fluid that is sucked into the vertical vortices is colder than for lower Pr which limits the efficiency of the upwards heat transfer.

preprint2009arXiv

Transitions between turbulent states in rotating Rayleigh-Benard convection

Weakly-rotating turbulent Rayleigh-Benard convection was studied experimentally and numerically. With increasing rotation and large enough Rayleigh number an abrupt transition from a turbulent state with nearly rotation-independent heat transport to another turbulent state with enhanced heat transfer is observed at a critical inverse Rossby number $1/Ro_c \simeq 0.4$. Whereas for $1/Ro < 1/Ro_c$ the strength of the large-scale convection-roll is either enhanced or essentially unmodified depending on parameters, its strength is increasingly diminished beyond $1/Ro_c$ where it competes with Ekman vortices that cause vertical fluid transport and thus heat-transfer enhancement.

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.