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

Richard J. A. M. Stevens

Richard J. A. M. Stevens 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.ao-phphysics.comp-ph

Catalog footprint

What is connected

47works
3topics
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

47 published item(s)

preprint2022arXiv

A three-dimensional dynamic mode decomposition analysis of wind farm flow aerodynamics

High-fidelity large-eddy simulations are suitable to obtain insight into the complex flow dynamics in extended wind farms. In order to better understand these flow dynamics, we use dynamic mode decomposition (DMD) to analyze and reconstruct the flow field in large-scale numerically simulated wind farms by large-eddy simulations (LES). Different wind farm layouts are considered, and we find that a combination of horizontal and vertical staggering leads to improved wind farm performance compared to traditional horizontal staggering. We analyze the wind farm flows using the amplitude selection (AP) and sparsity-promoting (SP method) DMD approach. We find that the AP method tends to select modes with a small length scale and a high frequency, while the SP method selects large coherent structures with low frequency. The latter are somewhat reminiscent of modes obtained using proper orthogonal decomposition (POD). We find that a relatively limited number of SP-DMD modes is sufficient to accurately reconstruct the flow field in the entire wind farm, whereas the AP-DMD method requires more modes to achieve an accurate reconstruction. Thus, the SP-DMD method has a smaller performance loss compared to the AP-DMD method in terms of the reconstruction of the flow field.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Evaluating the accuracy of the actuator line model against blade element momentum theory in uniform inflow

We evaluate the accuracy of the actuator line model (ALM) approach by performing simulations for the NREL~5~MW wind turbine in uniform inflow using three large eddy simulation codes. The power and thrust coefficients obtained using the three codes agrees within $1\%$ when the grid spacing $Δ_{\rm grid} \le 5.25$~m, and are cross-validated against blade element momentum (BEM) theory. We find that the results of ALM converge towards BEM theory without the need for tip correction when the numerical resolution is increased. For $Δ_{\rm grid}=0.98$~m the difference between the power and thrust coefficient obtained using ALM and BEM is $4.5\%$ and $2.1\%$, respectively, although we note that no absolute convergence between ALM and BEM can be obtained as both models use different assumptions, such as the use of a force projection method in the ALM. The difference in the local axial and tangential forces along the blades obtained from ALM simulations using $Δ_{\rm grid} = 1.97$~m and $Δ_{\rm grid} = 0.98$~m can be as large as $10\%$. The effect of the number of actuator points on the obtained turbine power and thrust coefficients is limited as the results converge when the spacing between the actuator points is about three times the grid spacing. This insight on the required number of blade points can be used to improve the efficiency of actuator line simulations.

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
preprint2021arXiv

Universal Wind Profile for Conventionally Neutral Atmospheric Boundary Layers

Conventionally neutral atmospheric boundary layers (CNBLs), which are characterized with zero surface potential temperature flux and capped by an inversion of potential temperature, are frequently encountered in nature. Therefore, predicting the wind speed profiles of CNBLs is relevant for weather forecasting, climate modeling, and wind energy applications. However, previous attempts to predict the velocity profiles in CNBLs have had limited success due to the complicated interplay between buoyancy, shear, and Coriolis effects. Here, we utilize ideas from the classical Monin-Obukhov similarity theory in combination with a local scaling hypothesis to derive an analytic expression for the stability correction function $ψ= -c_ψ(z/L)^{1/2}$, where $c_ψ= 4.2$ is an empirical constant, $z$ is the height above ground, and $L$ is the local Obukhov length based on potential temperature flux at that height, for CNBLs. An analytic expression for this flux is also derived using dimensional analysis and a perturbation method approach. We find that the derived profile agrees excellently with the velocity profile in the entire boundary layer obtained from high-fidelity large eddy simulations of typical CNBLs.

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

Interaction between low-level jets and wind farms in a stable atmospheric boundary layer

Low-level jets (LLJs) are the wind maxima in the lower regions of the atmosphere with a high wind energy potential. Here we use large-eddy simulations to study the effect of LLJ height on the flow dynamics in a wind farm with $10\times4$ turbines. We change the LLJ height and atmospheric thermal stratification by varying the surface cooling rate. We find that the first row power production is higher in the presence of a LLJ compared to a neutral reference case without LLJ. Besides, we show that the first row power production increases with decreasing LLJ height. Due to the higher turbulence intensity, the wind turbine wakes recover faster in a neutral boundary layer than in a stably stratified one. However, for strong thermal stratification with a low-height LLJ, the wake recovery can be faster than for the neutral reference case as energy can be entrained from the LLJ. Flow visualizations reveal that under stable stratification the growth of wind farm's internal boundary layer is restricted and the wind flows around the wind farm. Wind farms extract energy from LLJs through wake meandering and turbulent entrainment depending on the LLJ height. Both effects are advantageous for wake recovery, which is beneficial for the performance of downwind turbines. This finding is confirmed by an energy budget analysis, which reveals a significant increase in the kinetic energy flux in the presence of a LLJ. The jet strength reduces as it passes through consecutive turbine rows. For strong stratification, the combined effect of buoyancy destruction and turbulence dissipation is larger than the turbulent entrainment. Therefore, the power production of turbines in the back of the wind farm is relatively low for strong atmospheric stratifications.

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

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

Coherence of temperature and velocity superstructures in turbulent Rayleigh-Bénard flow

We investigate the interplay between large-scale patterns, so-called superstructures, in the fluctuation fields of temperature $θ$ and vertical velocity $w$ in turbulent Rayleigh-Bénard convection at large aspect ratios. Earlier studies suggested that velocity superstructures were smaller than their thermal counterparts in the center of the domain. However, a scale-by-scale analysis of the correlation between the two fields employing the linear coherence spectrum reveals that superstructures of the same size exist in both fields, which are almost perfectly correlated. The issue is further clarified by the observation that in contrast to the temperature, and unlike assumed previously, superstructures in the vertical velocity field do not result in a peak in the power spectrum of $w$. The origin of this difference is traced back to the production terms of the $θ$- and $w$-variance. These results are confirmed for a range of Rayleigh numbers $Ra = 10^5$--$10^9$, the superstructure size is seen to increase monotonically with $Ra$. Furthermore, the scale distribution of particularly the temperature fluctuations is pronouncedly bimodal. In addition to the large-scale peak caused by the superstructures, there exists a strong small-scale peak. This `inner peak' is most intense at a distance of $δ_θ$ from the wall and associated with structures of size $\approx 10 δ_θ$, where $δ_θ$ is the thermal boundary layer thickness. Finally, based on the vertical coherence relative to a reference height of $δ_θ$, a self-similar structure is identified in the velocity field (vertical and horizontal components) but not in the temperature.

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

Wind Turbine Large-Eddy Simulations on Very Coarse Grid Resolutions using an Actuator Line Model

In this work the accuracy of the Actuator Line Model (ALM) in Large Eddy Simulations of wind turbine flow is studied under the specific conditions of very coarse spatial resolutions. For finely-resolved conditions, it is known that ALM provides better accuracy compared to the standard Actuator Disk Model (ADM) without rotation. However, we show here that on very coarse resolutions, flow induction occurring at rotor scales can affect the predicted inflow angle and can adversely affect the ALM predictions. We first provide an illustration of coarse LES to reproduce wind tunnel measurements. The resulting flow predictions are good, but the challenges in predicting power outputs from the detailed ALM motivate more detailed analysis on a case with uniform inflow. We present a theoretical framework to compare the filtered quantities that enter the Large-Eddy Simulation equations as body forces with a scaling relation between the filtered and unfiltered quantities. The study aims to apply the theoretical derivation to the simulation framework and improve the current results for an ALM, especially in the near wake where the largest differences are observed.

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

Coupled wake boundary layer model of wind-farms

We present and test the coupled wake boundary layer (CWBL) model that describes the distribution of the power output in a wind-farm. The model couples the traditional, industry-standard wake model approach with a "top-down" model for the overall wind-farm boundary layer structure. This wake model captures the effect of turbine positioning, while the "top-down" portion of the model adds the interactions between the wind-turbine wakes and the atmospheric boundary layer. Each portion of the model requires specification of a parameter that is not known a-priori. For the wake model, the wake expansion coefficient is required, while the "top-down" model requires an effective spanwise turbine spacing within which the model's momentum balance is relevant. The wake expansion coefficient is obtained by matching the predicted mean velocity at the turbine from both approaches, while the effective spanwise turbine spacing depends on turbine positioning and thus can be determined from the wake model. Coupling of the constitutive components of the CWBL model is achieved by iterating these parameters until convergence is reached. We illustrate the performance of the model by applying it to both developing wind-farms including entrance effects and to fully developed (deep-array) conditions. Comparisons of the CWBL model predictions with results from a suite of large eddy simulations (LES) shows that the model closely represents the results obtained in these high-fidelity numerical simulations. A comparison with measured power degradation at the Horns Rev and Nysted wind-farms shows that the model can also be successfully applied to real wind-farms.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spatio-temporal spectra in the logarithmic layer of wall turbulence: large-eddy simulations and simple models

Motivated by the need to characterize the spatio-temporal structure of turbulence in wall-bounded flows, we study wavenumber-frequency spectra of the streamwise velocity component based on large-eddy simulation (LES) data. The LES data are used to measure spectra as a function of the two wall-parallel wavenumbers and the frequency in the equilibrium (logarithmic) layer. We then reformulate one of the simplest models that is able to reproduce the observations: the random sweeping model with a Gaussian large-scale fluctuating velocity and with additional mean flow. Comparison with LES data shows that the model captures the observed temporal decorrelation, which is related to the Doppler broadening of frequencies. We furthermore introduce a parameterization for the entire wavenumber-frequency spectrum $E_{11}(k_1,k_2,ω;z)$, where $k_1$, $k_2$ are the streamwise and spanwise wavenumbers, $ω$ is the frequency and $z$ is the distance to the wall. The results are found to be in good agreement with LES data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Turbulent Rayleigh-Bénard convection described by projected dynamics in phase space

Rayleigh-Bénard convection, i.e. the flow of a fluid between two parallel plates that is driven by a temperature gradient, is an idealised setup to study thermal convection. Of special interest are the statistics of the turbulent temperature field, which we are investigating and comparing for three different geometries, namely convection with periodic horizontal boundary conditions in three and two dimensions as well as convection in a cylindrical vessel, in order to work out similarities and differences. To this end, we derive an exact evolution equation for the temperature probability density function (PDF). Unclosed terms are expressed as conditional averages of velocities and heat diffusion, which are estimated from direct numerical simulations. This framework lets us identify the average behaviour of a fluid particle by revealing the mean evolution of fluid of different temperatures in different parts of the convection cell. We connect the statistics to the dynamics of Rayleigh-Bénard convection, giving deeper insights into the temperature statistics and transport mechanisms. We find that the average behaviour is described by closed cycles in phase space that reconstruct the typical Rayleigh-Bénard cycle of fluid heating up at the bottom, rising up to the top plate, cooling down and falling down again. The detailed behaviour shows subtle differences between the three cases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A concurrent precursor inflow method for Large Eddy Simulations and applications to finite length wind farms

In order to enable simulations of developing wind turbine array boundary layers with highly realistic inflow conditions a concurrent precursor method for Large Eddy Simulations is proposed. In this method we consider two domains simultaneously, i.e. in one domain a turbulent Atmospheric Boundary Layer (ABL) without wind turbines is simulated in order to generate the turbulent inflow conditions for a second domain in which the wind turbines are placed. The benefit of this approach is that a) it avoids the need for large databases in which the turbulent inflow conditions are stored and the correspondingly slow I/O operations and b) we are sure that the simulations are not negatively affected by statically swept fixed inflow fields or synthetic fields lacking the proper ABL coherent structures. Sample applications are presented, in which, in agreement with field data a strong decrease of the power output of downstream wind-turbines with respect to the first row of wind-turbines is observed for perfectly aligned inflow.

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

Large Eddy Simulation studies of the effects of alignment and wind farm length

Large eddy simulations of wind farms are performed to study the effects of wind turbine row alignment with respect to the incoming flow direction. Various wind farms with fixed stream-wise spacing (7.85 rotor diameters) and varying lateral displacements and span-wise turbine spacings are considered, for a fixed inflow direction. Simulations show that, contrary to common belief, a perfectly staggered (checker-board) configuration does not necessarily give the highest average power output. Instead, the highest mean wind farm power output is found to depend on several factors, the most important one being the alignment that leads to minimization of wake effects from turbines in several upstream rows. This alignment typically occurs at significantly smaller angles than those corresponding to perfect staggering. The observed trends have implications for wind farm designs, especially in sites with a well-defined prevailing wind direction.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Large-eddy simulation study of the logarithmic law for second and higher-order moments in turbulent wall-bounded flow

The logarithmic law for the mean velocity in turbulent boundary layers has long provided a valuable and robust reference for comparison with theories, models, and large-eddy simulations (LES) of wall-bounded turbulence. More recently, analysis of high-Reynolds number experimental boundary layer data has shown that also the variance and higher-order moments of the streamwise velocity fluctuations $u'^{+}$ display logarithmic laws. Such experimental observations motivate the question whether LES can accurately reproduce the variance and the higher-order moments, in particular their logarithmic dependency on distance to the wall. In this study we perform LES of very high Reynolds number wall-modeled channel flow and focus on profiles of variance and higher-order moments of the streamwise velocity fluctuations. In agreement with the experimental data, we observe an approximately logarithmic law for the variance in the LES, with a `Townsend-Perry' constant of $A_1\approx 1.25$. The LES also yields approximate logarithmic laws for the higher-order moments of the streamwise velocity. Good agreement is found between $A_p$, the generalized `Townsend-Perry' constants for moments of order $2p$, from experiments and simulations. Both are indicative of sub-Gaussian behavior of the streamwise velocity fluctuations. The near-wall behavior of the variance, the ranges of validity of the logarithmic law and in particular possible dependencies on characteristic length scales such as the roughness scale $z_0$, the LES grid scale $Δ$, and sub-grid scale (SGS) mixing length $C_sΔ$ are examined. We also present LES results on moments of spanwise and wall-normal fluctuations of velocity.

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
preprint2014arXiv

Temporal structure of aggregate power fluctuations in large-eddy simulations of extended wind-farms

Fluctuations represent a major challenge for the incorporation of electric power from large wind-farms into power grids. Wind farm power output fluctuates strongly in time, over various time scales. Understanding these fluctuations, especially their spatio-temporal characteristics, is particularly important for the design of backup power systems that must be readily available in conjunction with wind-farms. In this work we analyze the power fluctuations associated with the wind-input variability at scales between minutes to several hours, using large eddy simulations (LES) of extended wind-parks, interacting with the atmospheric boundary layer. LES studies enable careful control of parameters and availability of wind-velocities simultaneously across the entire wind-farm. The present study focuses on neutral atmospheric conditions and flat terrain, using actuator-disk representations of the individual wind-turbines. We consider power from various aggregates of wind-turbines such as the total average power signal, or signals from sub-averages within the wind-farm. Non-trivial correlations are observed due to the complex interactions between turbines placed downstream of each other, and they lead to noticeable spectral peaks at frequencies associated with the inter-turbine spacings when the wind-direction is completely fixed. In that case we observe that the frequency spectra of the total wind-farm output show a decay that follows approximately a $-5/3$ power-law scaling regime, qualitatively consistent with some observations made in field-scale operational wind-parks (Apt, 2007). We find that these features are still observed when the wind-speed varies in magnitude. However, significant changes in the wind-direction over time tend to smooth out the observed spectral peak and reduce the extent of the observed $-5/3$ power-law.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Heat transport and flow structure in rotating Rayleigh-Bénard convection

Here we summarize the results from our direct numerical simulations (DNS) and experimental measurements on rotating Rayleigh-Bénard (RB) convection. Our experiments and simulations are performed in cylindrical samples with an aspect ratio Γvarying from 1/2 to 2. Here Γ=D/L, where D and L are the diameter and height of the sample, respectively. When the rotation rate is increased, while a fixed temperature difference between the hot bottom and cold top plate is maintained, a sharp increase in the heat transfer is observed before the heat transfer drops drastically at stronger rotation rates. Here we focus on the question of how the heat transfer enhancement with respect to the non-rotating case depends on the Rayleigh number Ra, the Prandtl number Pr, and the rotation rate, indicated by the Rossby number Ro. Special attention will be given to the influence of the aspect ratio on the rotation rate that is required to get heat transport enhancement. In addition, we will discuss the relation between the heat transfer and the large scale flow structures that are formed in the different regimes of rotating RB convection and how the different regimes can be identified in experiments and simulations.

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

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.

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

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

Boundary layer structure in turbulent thermal convection and its consequences for the required numerical resolution

Results on the Prandtl-Blasius type kinetic and thermal boundary layer thicknesses in turbulent Rayleigh-Bénard convection in a broad range of Prandtl numbers are presented. By solving the laminar Prandtl-Blasius boundary layer equations, we calculate the ratio of the thermal and kinetic boundary layer thicknesses, which depends on the Prandtl number Pr only. It is approximated as $0.588Pr^{-1/2}$ for $Pr\ll Pr^*$ and as $0.982 Pr^{-1/3}$ for $Pr^*\ll\Pr$, with $Pr^*= 0.046$. Comparison of the Prandtl--Blasius velocity boundary layer thickness with that evaluated in the direct numerical simulations by Stevens, Verzicco, and Lohse (J. Fluid Mech. 643, 495 (2010)) gives very good agreement. Based on the Prandtl--Blasius type considerations, we derive a lower-bound estimate for the minimum number of the computational mesh nodes, required to conduct accurate numerical simulations of moderately high (boundary layer dominated) turbulent Rayleigh-Bénard convection, in the thermal and kinetic boundary layers close to bottom and top plates. It is shown that the number of required nodes within each boundary layer depends on Nu and Pr and grows with the Rayleigh number Ra not slower than $\sim\Ra^{0.15}$. This estimate agrees excellently with empirical results, which were based on the convergence of the Nusselt number in numerical simulations.

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Flow reversals in thermally driven turbulence

We analyze the reversals of the large scale flow in Rayleigh-Bénard convection both through particle image velocimetry flow visualization and direct numerical simulations (DNS) of the underlying Boussinesq equations in a (quasi) two-dimensional, rectangular geometry of aspect ratio 1. For medium Prandtl number there is a diagonal large scale convection roll and two smaller secondary rolls in the two remaining corners diagonally opposing each other. These corner flow rolls play a crucial role for the large scale wind reversal: They grow in kinetic energy and thus also in size thanks to plume detachments from the boundary layers up to the time that they take over the main, large scale diagonal flow, thus leading to reversal. Based on this mechanism we identify a typical time scale for the reversals. We map out the Rayleigh number vs Prandtl number phase space and find that the occurrence of reversals very sensitively depends on these parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Horizontal Structures of Velocity and Temperature Boundary Layers in 2D Numerical Turbulent Rayleigh-Bénard Convection

We investigate the structures of the near-plate velocity and temperature profiles at different horizontal positions along the conducting bottom (and top) plate of a Rayleigh-Bénard convection cell, using two-dimensional (2D) numerical data obtained at the Rayleigh number Ra=10^8 and the Prandtl number Pr=4.4 of an Oberbeck-Boussinesq flow with constant material parameters. The results show that most of the time, and for both velocity and temperature, the instantaneous profiles scaled by the dynamical frame method [Q. Zhou and K.-Q. Xia, Phys. Rev. Lett. 104, 104301 (2010) agree well with the classical Prandtl-Blasius laminar boundary layer (BL) profiles. Therefore, when averaging in the dynamical reference frames, which fluctuate with the respective instantaneous kinematic and thermal BL thicknesses, the obtained mean velocity and temperature profiles are also of Prandtl-Blasius type for nearly all horizontal positions. We further show that in certain situations the traditional definitions based on the time-averaged profiles can lead to unphysical BL thicknesses, while the dynamical method also in such cases can provide a well-defined BL thickness for both the kinematic and the thermal BLs.

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Prandtl-Blasius temperature and velocity boundary layer profiles in turbulent Rayleigh-Bénard convection

The shape of velocity and temperature profiles near the horizontal conducting plates in turbulent Rayleigh-Bénard convection are studied numerically and experimentally over the Rayleigh number range $10^8\lesssim Ra\lesssim3\times10^{11}$ and the Prandtl number range $0.7\lesssim Pr\lesssim5.4$. The results show that both the temperature and velocity profiles well agree with the classical Prandtl-Blasius laminar boundary-layer profiles, if they are re-sampled in the respective dynamical reference frames that fluctuate with the instantaneous thermal and velocity boundary-layer thicknesses.

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

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

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