Source author record

Siegfried Grossmann

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

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

18 published item(s)

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.

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.

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.

preprint2013arXiv

Velocity profiles in strongly turbulent Taylor-Couette flow

We derive the velocity profiles in strongly turbulent Taylor-Couette flow for the general case of independently rotating cylinders. The theory is based on the Navier-Stokes equations in the appropriate (cylinder) geometry. In particular, we derive the axial and the angular velocity profiles as functions of distance from the cylinder walls and find that both follow a logarithmic profile, with downwards-bending curvature corrections, which are more pronounced for the angular velocity profile as compared to the axial velocity profile, and which strongly increase with decreasing ratio $η$ between inner and outer cylinder radius. In contrast, the azimuthal velocity does not follow a log-law. We then compare the angular and azimuthal velocity profiles with the recently measured profiles in the ultimate state of (very) large Taylor numbers. Though the {\em qualitative} trends are the same -- down-bending for large wall distances and (properly shifted and non-dimensionalized) angular velocity profile $ω^+(r)$ being closer to a log-law than (properly shifted and non-dimensionalized) azimuthal velocity profile $u^+_φ(r)$ -- {\em quantitative} deviations are found for large wall distances. We attribute these differences to the Taylor rolls and the height dependence of the profiles, neither of which are considered in the theoretical approach.

preprint2012arXiv

Logarithmic temperature profiles in the ultimate regime of thermal convection

We report on the theory of logarithmic temperature profiles in very strongly developed thermal convection in the geometry of a Rayleigh-Benard cell with aspect ratio one and discuss the degree of agreement with the recently measured profiles in the ultimate state of very large Rayleigh number flow. The parameters of the log-profile are calculated and compared with the measure ones. Their physical interpretation as well as their dependence on the radial position are discussed.

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.

preprint2012arXiv

Optimal Taylor-Couette turbulence

Strongly turbulent Taylor-Couette flow with independently rotating inner and outer cylinders with a radius ratio of η= 0.716 is experimentally studied. From global torque measurements, we analyse the dimensionless angular velocity flux Nu_ω(Ta, a) as a function of the Taylor number Ta and the angular velocity ratio a = -ω_o/ω_i in the large-Taylor-number regime 10^{11} \lesssim Ta \lesssim 10^{13}. We analyse the data with the common power-law ansatz for the dimensionless angular velocity transport flux Nu_ω(Ta, a) = f(a)Ta^γ, with an amplitude f(a) and an exponent γ. The data are consistent with one effective exponent γ= 0.39\pm0.03 for all a. The amplitude of the angular velocity flux f(a) = Nu_ω(Ta, a)/Ta^0.39 is measured to be maximal at slight counter-rotation, namely at an angular velocity ratio of a_opt = 0.33\pm0.04. This value is theoretically interpreted as the result of a competition between the destabilizing inner cylinder rotation and the stabilizing but shear-enhancing outer cylinder counter-rotation. With the help of laser Doppler anemometry, we provide angular velocity profiles and identify the radial position r_n of the neutral line. While for moderate counter-rotation -0.40 ω_i \lesssim ω_o < 0, the neutral line still remains close to the outer cylinder and the probability distribution function (p.d.f.) of the bulk angular velocity is observed to be monomodal. For stronger counter-rotation the neutral line is pushed inwards towards the inner cylinder; in this regime the p.d.f. of the bulk angular velocity becomes bimodal, reflecting intermittent bursts of turbulent structures beyond the neutral line into the outer flow domain, which otherwise is stabilized by the counter-rotating outer cylinder. Finally, a hypothesis is offered allowing a unifying view for all these various results.

preprint2012arXiv

Spatial distribution of heat flux and fluctuations in turbulent Rayleigh-Bénard convection

We numerically investigate the radial dependence of the velocity and temperature fluctuations and of the time-averaged heat flux $\bar{j}(r)$ in a cylindrical Rayleigh-Bénard cell with aspect ratio Γ=1 for Rayleigh numbers Ra between $2 \times 10^6$ and $2\times 10^{9}$ at a fixed Prandtl number Pr = 5.2. The numerical results reveal that the heat flux close to the sidewall is larger than in the center and that, just as the global heat transport, it has an effective power law dependence on the Rayleigh number, $\bar{j}(r)\propto Ra^{γ_j(r)}$. The scaling exponent $γ_j(r)$ decreases monotonically from 0.43 near the axis ($r \approx 0$) to 0.29 close to the side walls ($r \approx D/2$). The effective exponents near the axis and the side wall agree well with the measurements of Shang et al. (Phys.\ Rev.\ Lett.\ \textbf{100}, 244503, 2008) and the predictions of Grossmann and Lohse (Phys.\ Fluids \textbf{16}, 1070, 2004). Extrapolating our results to large Rayleigh number would imply a crossover at $Ra\approx 10^{15}$, where the heat flux near the axis would begin to dominate. In addition, we find that the local heat flux is more than twice as high at the location where warm or cold plumes go up or down, than in the plume depleted regions.

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.

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.

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.

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.

preprint2011arXiv

Ultimate Turbulent Taylor-Couette Flow

The flow structure of strongly turbulent Taylor-Couette flow with Reynolds numbers up to Re_i = 2*10^6 of the inner cylinder is experimentally examined with high-speed particle image velocimetry (PIV). The wind Reynolds numbers Re_w of the turbulent Taylor-vortex flow is found to scale as Re_w ~ Ta^(1/2), exactly as predicted for the ultimate turbulence regime, in which the boundary layers are turbulent. The dimensionless angular velocity flux has an effective scaling of Nu_ω ~ Ta^0.38, also in correspondence with turbulence in the ultimate regime. The scaling of Nu_ω is confirmed by local angular velocity flux measurements extracted from high-speed PIV measurements: though the flux shows huge fluctuations, its spatial and temporal average nicely agrees with the result from the global torque measurements.

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.

preprint2008arXiv

Flow organization in non-Oberbeck-Boussinesq Rayleigh-Benard convection in water

Non-Oberbeck-Boussinesq (NOB) effects on the flow organization in two-dimensional Rayleigh-Benard turbulence are numerically analyzed. The working fluid is water. We focus on the temperature profiles, the center temperature, the Nusselt number, and on the analysis of the velocity field. Several velocity amplitudes (or Reynolds numbers) and several kinetic profiles are introduced and studied; these together describe the various features of the rather complex flow organization. The results are presented both as functions of the Rayleigh number Ra (with Ra up to 10^8) for fixed temperature difference (Delta) between top and bottom plates and as functions of Delta ("non-Oberbeck-Boussinesqness") for fixed Ra with Delta up to 60 K. All results are consistent with the available experimental NOB data for the center temperature Tc and the Nusselt number ratio Nu_{NOB}/Nu_{OB} (the label OB meaning that the Oberbeck-Boussinesq conditions are valid). Beyond Ra ~ 10^6 the flow consists of a large diagonal center convection roll and two smaller rolls in the upper and lower corners. In the NOB case the center convection roll is still characterized by only one velocity scale.

preprint2008arXiv

Large scale dynamics in turbulent Rayleigh-Benard convection

The progress in our understanding of several aspects of turbulent Rayleigh-Benard convection is reviewed. The focus is on the question of how the Nusselt number and the Reynolds number depend on the Rayleigh number Ra and the Prandtl number Pr, and on how the thicknesses of the thermal and the kinetic boundary layers scale with Ra and Pr. Non-Oberbeck-Boussinesq effects and the dynamics of the large-scale convection-roll are addressed as well. The review ends with a list of challenges for future research on the turbulent Rayleigh-Benard system.

preprint1998arXiv

Continued Fraction Representation of Temporal Multi Scaling in Turbulence

It was shown recently that the anomalous scaling of simultaneous correlation functions in turbulence is intimately related to the breaking of temporal scale invariance, which is equivalent to the appearance of infinitely many times scales in the time dependence of time-correlation functions. In this paper we derive a continued fraction representation of turbulent time correlation functions which is exact and in which the multiplicity of time scales is explicit. We demonstrate that this form yields precisely the same scaling laws for time derivatives and time integrals as the "multi-fractal" representation that was used before. Truncating the continued fraction representation yields the "best" estimates of time correlation functions if the given information is limited to the scaling exponents of the simultaneous correlation functions up to a certain, finite order. It is worth noting that the derivation of a continued fraction representation obtained here for an operator which is not Hermitian or anti-Hermitian may be of independent interest.

preprint1997arXiv

Maxwell's Demon at work: Two types of Bose condensate fluctuations in power-law traps

After discussing the key idea underlying the Maxwell's Demon ensemble, we employ this idea for calculating fluctuations of ideal Bose gas condensates in traps with power-law single-particle energy spectra. Two essentially different cases have to be distinguished. If the heat capacity remains continuous at the condensation point in the large-N-limit, the fluctuations of the number of condensate particles vanish linearly with temperature, independent of the trap characteristics. If the heat capacity becomes discontinuous, the fluctuations vanish algebraically with temperature, with an exponent determined by the trap. Our results are based on an integral representation that yields the solution to both the canonical and the microcanonical fluctuation problem in a singularly transparent manner.