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

Erwin P. van der Poel

Erwin P. van der Poel 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-ph

Catalog footprint

What is connected

11works
2topics
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

11 published item(s)

preprint2016arXiv

Penetrative internally heated convection in two and three dimensions

Convection of an internally heated fluid, confined between top and bottom plates of equal temperature, is studied by direct numerical simulation in two and three dimensions. The unstably stratified upper region drives convection that penetrates into the stably stratified lower region. The fraction of produced heat escaping across the bottom plate, which is one half without convection, initially decreases as convection strengthens. Entering the turbulent regime, this decrease reverses in two dimensions but continues monotonically in three dimensions. The mean fluid temperature, which grows proportionally to the heating rate ($H$) without convection, grows proportionally to $H^{4/5}$ when convection is strong in both two and three dimensions. The ratio of the heating rate to the fluid temperature is likened to the Nusselt number of Rayleigh-Bénard convection. Simulations are reported for Prandtl numbers between 0.1 and 10 and for Rayleigh numbers (defined in terms of the heating rate) up to $5\times10^{10}$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Translational and rotational dynamics of a large buoyant sphere in turbulence

We report experimental measurements of the translational and rotational dynamics of a large buoyant sphere in isotropic turbulence. We introduce an efficient method to simultaneously determine the position and (absolute) orientation of a spherical body from visual observation. The method employs a minimization algorithm to obtain the orientation from the 2D projection of a specific pattern drawn onto the surface of the sphere. This has the advantages that it does not require a database of reference images, is easily scalable using parallel processing, and enables accurate absolute orientation reference. Analysis of the sphere's translational dynamics reveals clear differences between the streamwise and transverse directions. The translational auto-correlations and PDFs provide evidence for periodicity in the particle's dynamics even under turbulent conditions. The angular autocorrelations show weak periodicity. The angular accelerations exhibit wide tails, however without a directional dependence.

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

Comparison between two and three dimensional Rayleigh-Bénard convection

Two dimensional (2D) and three dimensional (3D) Rayleigh-Bénard convection is compared using results from direct numerical simulations and prior experiments. The explored phase diagrams for both cases are reviewed. The differences and similarities between 2D and 3D are studied using Nu(Ra) for Pr = 4.38 and Pr = 0.7 and Nu(Pr) for Ra up to $10^8$. In the Nu(Ra) scaling at higher Pr, 2D and 3D are very similar; differing only by a constant factor up to $\text{Ra} = 10^{10}$. In contrast, the difference is large at lower Pr, due to the strong roll state dependence of Nu in 2D. The behaviour of Nu(Pr) is similar in 2D and 3D at large Pr. However, it differs significantly around Pr = 1. The Reynolds number values are consistently higher in 2D and additionally converge at large Pr. Finally, the thermal boundary layer profiles are compared in 2D and 3D.

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
preprint2013arXiv

The unifying theory of scaling in thermal convection: The updated prefactors

The unifying theory of scaling in thermal convection (Grossmann & Lohse (2000)) (henceforth the GL theory) suggests that there are no pure power laws for the Nusselt and Reynolds numbers as function of the Rayleigh and Prandtl numbers in the experimentally accessible parameter regime. In Grossmann & Lohse (2001) the dimensionless parameters of the theory were fitted to 155 experimental data points by Ahlers & Xu (2001) in the regime $3\times 10^7 \le Ra \le 3 \times 10^{9}$ and $4\le Pr \le 34$ and Grossmann & Lohse (2002) used the experimental data point from Qiu & Tong (2001) and the fact that Nu(Ra,Pr) is independent of the parameter a, which relates the dimensionless kinetic boundary thickness with the square root of the wind Reynolds number, to fix the Reynolds number dependence. Meanwhile the theory is on one hand well confirmed through various new experiments and numerical simulations. On the other hand these new data points provide the basis for an updated fit in a much larger parameter space. Here we pick four well established (and sufficiently distant) Nu(Ra,Pr) data points and show that the resulting Nu(Ra,Pr) function is in agreement with almost all established experimental and numerical data up to the ultimate regime of thermal convection, whose onset also follows from the theory. One extra Re(Ra,Pr) data point is used to fix Re(Ra,Pr). As Re can depend on the definition and the aspect ratio the transformation properties of the GL equations are discussed in order to show how the GL coefficients can easily be adapted to new Reynolds number data while keeping Nu(Ra,Pr) unchanged.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Flow states in two-dimensional Rayleigh-Bénard convection as a function of aspect-ratio and Rayleigh number

In this numerical study on two-dimensional Rayleigh-Bénard convection we consider $10^7 \leq Ra \leq 10^{12}$ in aspect ratio $0.23 \leq Γ\leq 13$ samples. We focus on several cases. First we consider small aspect ratio cells, where at high Ra number we find a sharp transition from a low Ra number branch towards a high Ra number branch, due to changes in the flow structure. Subsequently, we show that the influence of the aspect ratio on the heat transport decreases with increasing aspect ratio, although even at very large aspect ratio of $Γ\approx10$ variations up to 2.5% in the heat transport as a function of Γare observed. Finally, we observe long-lived transients up to at least $Ra=10^9$, as in certain aspect ratio cells we observe different flow states that are stable for thousands of turnover times.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Connecting flow structures and heat flux in turbulent Rayleigh-Bénard convection

The aspect ratio (Γ) dependence of the heat transfer (Nusselt number Nu in dimensionless form) in turbulent (two-dimensional) Rayleigh-Bénard convection is numerically studied in the regime $0.4 \le Γ\le 1.25$ for Rayleigh numbers $10^7 \le Ra \le Ra^{9}$ and Prandtl numbers Pr =0.7 (gas) and 4.3 (water). Nu (Γ) shows a very rich structure with sudden jumps and sharp transitions. We connect these structures to the way the flow organizes itself in the sample and explain why the aspect ratio dependence of Nu is more pronounced for small Pr. Even for fixed Γdifferent turbulent states (with different resulting Nu) can exist, between which the flow can or cannot switch. In the latter case the heat transfer thus depends on the initial conditions.

0 citations0 reviewsNo code linkedOpen access