Showing 1–50 of 62 results for author: Grossmann, S

OpenLB User Guide: Associated with Release 1.6 of the Code

Authors: Adrian Kummerländer, Samuel J. Avis, Halim Kusumaatmaja, Fedor Bukreev, Michael Crocoll, Davide Dapelo, Simon Großmann, Nicolas Hafen, Shota Ito, Julius Jeßberger, Eliane Kummer, Jan E. Marquardt, Johanna Mödl, Tim Pertzel, František Prinz, Florian Raichle, Martin Sadric, Maximilian Schecher, Dennis Teutscher, Stephan Simonis, Mathias J. Krause

Abstract: OpenLB is an object-oriented implementation of LBM. It is the first implementation of a generic platform for LBM programming, which is shared with the open source community (GPLv2). Since the first release in 2007, the code has been continuously improved and extended which is documented by thirteen releases as well as the corresponding release notes which are available on the OpenLB website (https… ▽ More OpenLB is an object-oriented implementation of LBM. It is the first implementation of a generic platform for LBM programming, which is shared with the open source community (GPLv2). Since the first release in 2007, the code has been continuously improved and extended which is documented by thirteen releases as well as the corresponding release notes which are available on the OpenLB website (https://www.openlb.net). The OpenLB code is written in C++ and is used by application programmers as well as developers, with the ability to implement custom models OpenLB supports complex data structures that allow simulations in complex geometries and parallel execution using MPI, OpenMP and CUDA on high-performance computers. The source code uses the concepts of interfaces and templates, so that efficient, direct and intuitive implementations of the LBM become possible. The efficiency and scalability has been checked and proved by code reviews. This user manual and a source code documentation by DoxyGen are available on the OpenLB project website. △ Less

Submitted 7 August, 2024; v1 submitted 17 May, 2023; originally announced July 2023.

MSC Class: 97N80; 74S30; 76M25; 80M25; 82C80; ACM Class: G.4

Dielectric Modeling of Oil-paper Insulation Systems at High DC Voltage Stress Using a Charge-carrier-based Approach

Authors: Tobias Gabler, Karsten Backhaus, Steffen Großmann, Ronny Fritsche

Abstract: It is state-of-the-art to describe the dielectric behavior of an insulation material by its permittivity and its specific electric conductivity in order to estimate the dielectric stress of an insulation system. Thus, the electric field at DC voltage stress is determined according to the stationary electrical conduction field with the electric conductivity of the insulation materials. However, at… ▽ More It is state-of-the-art to describe the dielectric behavior of an insulation material by its permittivity and its specific electric conductivity in order to estimate the dielectric stress of an insulation system. Thus, the electric field at DC voltage stress is determined according to the stationary electrical conduction field with the electric conductivity of the insulation materials. However, at oil-insulated arrangements a higher field strength in front of bare metal electrodes at high DC voltage stress occurs, which is not explainable with this model. Therefore, a charge carrier-based approach is presented to describe the dielectric behavior of the oil-paper insulation. It describes the movement of charge carriers and their effect on the electric field strength. Their drift leads to an accumulation of charge carriers in front of electrodes which results in a higher field strength in these areas, which can be calculated numerically using the Poisson-Nernst-Planck equation system. Compared to the conductivity-based model fundamental differences can be shown. Breakdown experiments qualitatively confirm the expectations according to the charge carrier-based approach. The results show that the charge carrier-based field distribution has to be considered for modelling the electrical field strength distribution at high DC voltage stress. They also show, that the dielectric behavior of these arrangements cannot be explained according to the state-of-the-art model. △ Less

Submitted 21 November, 2019; originally announced November 2019.

Comments: This accepted version of the contribution was published in IEEE Transactions on Dielectrics and Electrical Insulation Vol. 26, No. 5; October 2019 with DOI 10.1109/TDEI.2019.008175, see https://doi.org/10.1109/TDEI.2019.008175

Journal ref: IEEE Trans. Dielectr. Electr. Insul. 26-5 (2019) 1549-1557

Nonclassical optical properties of mesoscopic gold

Authors: Swen Großmann, Daniel Friedrich, Michael Karolak, René Kullock, Enno Krauss, Monika Emmerling, Giorgio Sangiovanni, Bert Hecht

Abstract: Gold nanostructures have important applications in nanoelectronics, nano-optics as well as in precision metrology due to their intriguing opto-electronic properties. These properties are governed by the bulk band structure but to some extend are tunable via geometrical resonances. Here we show that the band structure of gold itself exhibits significant size-dependent changes already for mesoscopic… ▽ More Gold nanostructures have important applications in nanoelectronics, nano-optics as well as in precision metrology due to their intriguing opto-electronic properties. These properties are governed by the bulk band structure but to some extend are tunable via geometrical resonances. Here we show that the band structure of gold itself exhibits significant size-dependent changes already for mesoscopic critical dimensions below 30 nm. To suppress the effects of geometrical resonances and grain boundaries, we prepared atomically flat ultrathin films of various thicknesses by utilizing large chemically grown single-crystalline gold platelets. We experimentally probe thickness-dependent changes of the band structure by means of two-photon photoluminescence and observe a surprising 100-fold increase of the nonlinear signal when the gold film thickness is reduced below 30 nm allowing us to optically resolve single-unit-cell steps. The effect is well explained by density functional calculations of the thickness-dependent 2D band structure of gold. △ Less

Submitted 23 May, 2019; originally announced May 2019.

High Reynolds number Taylor-Couette turbulence

Authors: Siegfried Grossmann, Detlef Lohse, Chao Sun

Abstract: Taylor-Couette flow -- the flow between two coaxial co- or counter-rotating cylinders -- is one of the paradigmatic systems in physics of fluids. The (dimensionless) control parameters are the Reynolds numbers of the inner and outer cylinder, the ratio of the cylinder radii, and the aspect ratio. The response of the system is the torque required to retain constant angular velocities, which can be… ▽ More Taylor-Couette flow -- the flow between two coaxial co- or counter-rotating cylinders -- is one of the paradigmatic systems in physics of fluids. The (dimensionless) control parameters are the Reynolds numbers of the inner and outer cylinder, the ratio of the cylinder radii, and the aspect ratio. The response of the system is the torque required to retain constant angular velocities, which can be connected to the angular velocity transport through the gap. While the low Reynolds number regime has been very well explored in the '80s and '90s of the last century, in the fully turbulent regime major research activity only developed in the last decade. In this paper we review this recent progress in our understanding of fully developed Taylor-Couette turbulence, from the experimental, numerical, and theoretical point of view. We will focus on the parameter dependence of the global torque and on the local flow organisation, including velocity profiles and boundary layers. Next, we will discuss transitions between different (turbulent) flow states. We will also elaborate on the relevance of this system for astrophysical disks (Keplerian flows). The review ends with a list of challenges for future research on turbulent Taylor-Couette flow. The published version in ARFM: Annu. Rev. Fluid Mech. 48:53-80 (2016). △ Less

Submitted 30 March, 2019; originally announced April 2019.

A Spin-Optical Nano Device

Authors: Enno Krauss, Gary Razinskas, Dominik Köck, Swen Grossmann, Bert Hecht

Abstract: The photon spin is an important resource for quantum information processing as is the electron spin in spintronics. However, for subwavelength confined optical excitations, polarization as a global property of a mode cannot be defined. Here, we show that any polarization state of a plane-wave photon can reversibly be mapped to a pseudo-spin embodied by the two fundamental modes of a subwavelength… ▽ More The photon spin is an important resource for quantum information processing as is the electron spin in spintronics. However, for subwavelength confined optical excitations, polarization as a global property of a mode cannot be defined. Here, we show that any polarization state of a plane-wave photon can reversibly be mapped to a pseudo-spin embodied by the two fundamental modes of a subwavelength plasmonic two-wire transmission line. We design a device in which this pseudo-spin evolves in a well-defined fashion throughout the device reminiscent of the evolution of photon polarization in a birefringent medium and the behaviour of electron spins in the channel of a spin field-effect transistor. The significance of this pseudo-spin is enriched by the fact that it is subject to spin-orbit locking. Combined with optically active materials to exert external control over the pseudo-spin precession, our findings could enable spin-optical transistors, i.e. the routing and processing of quantum information with light on a subwavelength scale. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: 21 pages, 10 figures

MSC Class: 78-05

Evidence for cascaded third harmonic generation in non-centrosymmetric gold nanoantennas

Authors: Michele Celebrano, Andrea Locatelli, Lavinia Ghirardini, Giovanni Pellegrini, Paolo Biagioni, Attilio Zilli, Xiaofei Wu, Swen Grossmann, Luca Carletti, Costantino De Angelis, Lamberto Duò, Bert Hecht, Marco Finazzi

Abstract: The optimization of nonlinear optical processes at the nanoscale is a crucial step for the development of nanoscale photon sources for quantum-optical networks. The development of innovative plasmonic nanoantenna designs and hybrid nanostructures to enhance optical nonlinearities in very small volumes represents one of the most promising routes. In such systems, the upconversion of photons can be… ▽ More The optimization of nonlinear optical processes at the nanoscale is a crucial step for the development of nanoscale photon sources for quantum-optical networks. The development of innovative plasmonic nanoantenna designs and hybrid nanostructures to enhance optical nonlinearities in very small volumes represents one of the most promising routes. In such systems, the upconversion of photons can be achieved with high efficiencies via third-order processes, such as third harmonic generation (THG), thanks to the resonantly-enhanced volume currents. Conversely, second-order processes, such as second harmonic generation (SHG), are often inhibited by the symmetry of metal lattices and of common nanoantenna geometries. SHG and THG processes in plasmonic nanostructures are generally treated independently, since they both represent a small perturbation in the light-matter interaction mechanisms. In this work, we demonstrate that this paradigm does not hold in general, by providing evidence of a cascaded process in THG, which is fueled by SHG and sizably contributes to the overall yield. We address this mechanism by unveiling an anomalous fingerprint in the polarization state of the nonlinear emission from non-centrosymmetric gold nanoantennas and point out that such cascaded processes may also appear for structures that exhibit only moderate SHG yields - signifying its general relevance in plasmon-enhanced nonlinear optics. The presence of this peculiar mechanism in THG from plasmonic nanoantennas at telecommunication wavelengths allows gaining further insight on the physics of plasmon-enhanced nonlinear optical processes. This could be crucial in the realization of nanoscale elements for photon conversion and manipulation operating at room-temperature. △ Less

Submitted 2 July, 2019; v1 submitted 9 March, 2018; originally announced March 2018.

Comments: 25 pages, 4 figures

Scaling relations in large-Prandtl-number natural thermal convection

Authors: Olga Shishkina, Mohammad S. Emran, Siegfried Grossmann, Detlef Lohse

Abstract: In this study we follow Grossmann and Lohse, Phys. Rev. Lett. 86 (2001), who derived various scalings regimes for the dependence of the Nusselt number $Nu$ and the Reynolds number $Re$ on the Rayleigh number $Ra$ and the Prandtl number $Pr$. We focus on theoretical arguments as well as on numerical simulations for the case of large-$Pr$ natural thermal convection. Based on an analysis of self-simi… ▽ More In this study we follow Grossmann and Lohse, Phys. Rev. Lett. 86 (2001), who derived various scalings regimes for the dependence of the Nusselt number $Nu$ and the Reynolds number $Re$ on the Rayleigh number $Ra$ and the Prandtl number $Pr$. We focus on theoretical arguments as well as on numerical simulations for the case of large-$Pr$ natural thermal convection. Based on an analysis of self-similarity of the boundary layer equations, we derive that in this case the limiting large-$Pr$ boundary-layer dominated regime is I$_\infty^<$, introduced and defined in [1], with the scaling relations $Nu\sim Pr^0\,Ra^{1/3}$ and $Re\sim Pr^{-1}\,Ra^{2/3}$. Our direct numerical simulations for $Ra$ from $10^4$ to $10^9$ and $Pr$ from 0.1 to 200 show that the regime I$_\infty^<$ is almost indistinguishable from the regime III$_\infty$, where the kinetic dissipation is bulk-dominated. With increasing $Ra$, the scaling relations undergo a transition to those in IV$_u$ of reference [1], where the thermal dissipation is determined by its bulk contribution. △ Less

Submitted 27 September, 2017; v1 submitted 22 July, 2017; originally announced July 2017.

Journal ref: Phys. Rev. Fluids 2, 103502 (2017)

Electromechanically Tunable Suspended Optical Nano-antenna

Authors: Kai Chen, Gary Razinskas, Thorsten Feichtner, Swen Grossmann, Silke Christiansen, Bert Hecht

Abstract: Coupling mechanical degrees of freedom with plasmonic resonances has potential applications in optomechanics, sensing, and active plasmonics. Here we demonstrate a suspended two-wire plasmonic nano-antenna acting like a nano-electrometer. The antenna wires are supported and electrically connected via thin leads without disturbing the antenna resonance. As a voltage is applied, equal charges are in… ▽ More Coupling mechanical degrees of freedom with plasmonic resonances has potential applications in optomechanics, sensing, and active plasmonics. Here we demonstrate a suspended two-wire plasmonic nano-antenna acting like a nano-electrometer. The antenna wires are supported and electrically connected via thin leads without disturbing the antenna resonance. As a voltage is applied, equal charges are induced on both antenna wires. The resulting equilibrium between the repulsive Coulomb force and the restoring elastic bending force enables us to precisely control the gap size. As a result the resonance wavelength and the field enhancement of the suspended optical nano-antenna (SONA) can be reversibly tuned. Our experiments highlight the potential to realize large bandwidth optical nanoelectromechanical systems (NEMS). △ Less

Submitted 30 January, 2016; originally announced February 2016.

Salinity transfer in bounded double diffusive convection

Authors: Yantao Yang, Erwin P. van der Poel, Rodolfo Ostilla-Mónico, Chao Sun, Roberto Verzicco, Siegfried Grossmann, Detlef Lohse

Abstract: 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 fo… ▽ More 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. △ Less

Submitted 6 October, 2015; originally announced October 2015.

Journal ref: Journal of fluid mechanics 768 (2015): 476-491

Logarithmic mean temperature profiles and their connection to plume emissions in turbulent Rayleigh-Bénard convection

Authors: Erwin P. van der Poel, Rodolfo Ostilla-Mónico, Roberto Verzicco, Siegfried Grossmann, Detlef Lohse

Abstract: 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 l… ▽ More 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. △ Less

Submitted 6 October, 2015; originally announced October 2015.

Comments: 6 pages

Journal ref: Phys. Rev. Lett. 115, 154501 (2015)

The near-wall region of highly turbulent Taylor-Couette flow

Authors: Rodolfo Ostilla-Mónico, Roberto Verzicco, Siegfried Grossmann, Detlef Lohse

Abstract: Direct numerical simulations of the Taylor-Couette (TC) problem, the flow between two coaxial and independently rotating cylinders, have been performed. The study focuses on TC flow with mild curvature (small gap) with a radius ratio of $η=r_i/r_o=0.909$, an aspect ratio of $Γ=L/d=2π/3$, and a stationary outer cylinder. Three inner cylinder Reynolds of $1\cdot10^5$, $2\cdot10^5$ and $3\cdot 10^5$… ▽ More Direct numerical simulations of the Taylor-Couette (TC) problem, the flow between two coaxial and independently rotating cylinders, have been performed. The study focuses on TC flow with mild curvature (small gap) with a radius ratio of $η=r_i/r_o=0.909$, an aspect ratio of $Γ=L/d=2π/3$, and a stationary outer cylinder. Three inner cylinder Reynolds of $1\cdot10^5$, $2\cdot10^5$ and $3\cdot 10^5$ were simulated, corresponding to frictional Reynolds numbers between $Re_τ\approx 1400$ and $Re_τ\approx 4000$. An additional case with a large gap, $η=0.5$ and driving of $Re=2\cdot10^5$ was also performed. Small-gap TC was found to be dominated by spatially-fixed large-scale structures, known as Taylor rolls (TRs). TRs are attached to the boundary layer, and are active, i.e. they transport angular velocity through Reynolds stresses. An additional simulation with inner cylinder Reynolds number of $Re=1\cdot10^5$ and fixed outer cylinder with an externally imposed axial flow of comparable strength as the wind of the TRs was also conducted. The axial flow was found to convect the TRs without any weakening effect. For small-gap TC, evidence for the existence of logarithmic velocity fluctuations, and of an overlap layer, in which the velocity fluctuations collapse in outer units, was found. Profiles consistent with a logarithmic dependence were also found for the angular velocity in large-gap TC, albeit in a very reduced range of scales. Finally, the behaviour of both small- and large-gap TC was compared to other canonical flows. Small-gap TC has similar behaviour in the near-wall region to other canonical flows, while large-gap TC displays very different behaviour. △ Less

Submitted 23 April, 2017; v1 submitted 13 January, 2015; originally announced January 2015.

Journal ref: Journal of Fluid Mechanics, (788), p. 95-117, 2016

Mode-matching in multiresonant plasmonic nanoantennas for enhanced second harmonic generation

Authors: Michele Celebrano, Xiaofei Wu, Milena Baselli, Swen Großmann, Paolo Biagioni, Andrea Locatelli, Costantino De Angelis, Giulio Cerullo, Roberto Osellame, Bert Hecht, Lamberto Duò, Franco Ciccacci, Marco Finazzi

Abstract: Boosting nonlinear frequency conversion in extremely confined volumes remains a key challenge in nano-optics, nanomedicine, photocatalysis, and background-free biosensing. To this aim, field enhancements in plasmonic nanostructures are often exploited to effectively compensate for the lack of phase-matching at the nanoscale. Second harmonic generation (SHG) is, however, strongly quenched by the hi… ▽ More Boosting nonlinear frequency conversion in extremely confined volumes remains a key challenge in nano-optics, nanomedicine, photocatalysis, and background-free biosensing. To this aim, field enhancements in plasmonic nanostructures are often exploited to effectively compensate for the lack of phase-matching at the nanoscale. Second harmonic generation (SHG) is, however, strongly quenched by the high degree of symmetry in plasmonic materials at the atomic scale and in nanoantenna designs. Here, we devise a plasmonic nanoantenna lacking axial symmetry, which exhibits spatial and frequency mode overlap at both the excitation and the SHG wavelengths. The effective combination of these features in a single device allows obtaining unprecedented SHG conversion efficiency. Our results shed new light on the optimization of SHG at the nanoscale, paving the way to new classes of nanoscale coherent light sources and molecular sensing devices based on nonlinear plasmonic platforms. △ Less

Submitted 10 December, 2014; v1 submitted 1 December, 2014; originally announced December 2014.

Comments: 14 pages, 4 figures

Exploring the phase diagram of fully turbulent Taylor-Couette flow

Authors: Rodolfo Ostilla Mónico, Erwin P. van der Poel, Roberto Verzicco, Siegfried Grossmann, Detlef Lohse

Abstract: Direct numerical simulations of Taylor-Couette flow (TC). Shear Reynolds numbers of up to $3\cdot10^5$, corresponding to Taylor numbers of $Ta=4.6\cdot10^{10}$, were reached. Effective scaling laws for the torque are presented. The transition to the ultimate regime, in which asymptotic scaling laws (with logarithmic corrections) for the torque are expected to hold up to arbitrarily high driving, i… ▽ More Direct numerical simulations of Taylor-Couette flow (TC). Shear Reynolds numbers of up to $3\cdot10^5$, corresponding to Taylor numbers of $Ta=4.6\cdot10^{10}$, were reached. Effective scaling laws for the torque are presented. The transition to the ultimate regime, in which asymptotic scaling laws (with logarithmic corrections) for the torque are expected to hold up to arbitrarily high driving, is analysed for different radius ratios $η$, different aspect ratios $Γ$ and different rotation ratios $Ro$. It is shown that the transition is approximately independent of $Ro$ and $Γ$, but depends significantly on $η$. We furthermore calculate the local angular velocity profiles and visualize different flow regimes that depend both on the shearing of the flow, and the Coriolis force originating from the outer cylinder rotation. Two main regimes are distinguished, based on the magnitude of the Coriolis force, namely the co-rotating and weakly counter-rotating regime dominated by Rayleigh-unstable regions, and the strongly counter-rotating regime where a mixture of stable and unstable regions exist. Furthermore, an analogy between $η$ and outer-cylinder rotation is revealed, namely that smaller gaps behave like a wider gap with co-rotating cylinders, and that wider gaps behave like smaller gaps with weakly counter-rotating cylinders. Finally, the effect of $Γ$ on the effective torque versus $Ta$ scaling is analysed and it is shown that different branches of the torque-versus-$Ta$ relationships associated to different aspect ratios are found to cross within $15%$ of the $Re$ associated to the transition to the ultimate regime. The paper culminates in phase diagram in the inner vs outer $Re$ number parameter space and in the $Ta$ vs $Ro$ parameter space, which can be seen as the extension of the Andereck \emph{et al.} phase diagram towards the ultimate regime. △ Less

Submitted 15 October, 2014; v1 submitted 1 May, 2014; originally announced May 2014.

Comments: 23 pages, submitted to JFM

Turbulence decay towards the linearly-stable regime of Taylor-Couette flow

Authors: Rodolfo Ostilla Mónico, Roberto Verzicco, Siegfried Grossmann, Detlef Lohse

Abstract: Taylor-Couette (TC) flow is used to probe the hydrodynamical stability of astrophysical accretion disks. Experimental data on the subcritical stability of TC are in conflict about the existence of turbulence (cf. Ji et al. Nature, 444, 343-346 (2006) and Paoletti et al., A$\&$A, 547, A64 (2012)), with discrepancies attributed to end-plate effects. In this paper we numerically simulate TC flow with… ▽ More Taylor-Couette (TC) flow is used to probe the hydrodynamical stability of astrophysical accretion disks. Experimental data on the subcritical stability of TC are in conflict about the existence of turbulence (cf. Ji et al. Nature, 444, 343-346 (2006) and Paoletti et al., A$\&$A, 547, A64 (2012)), with discrepancies attributed to end-plate effects. In this paper we numerically simulate TC flow with axially periodic boundary conditions to explore the existence of sub-critical transitions to turbulence when no end-plates are present. We start the simulations with a fully turbulent state in the unstable regime and enter the linearly stable regime by suddenly starting a (stabilizing) outer cylinder rotation. The shear Reynolds number of the turbulent initial state is up to $Re_s \sim10^5$ and the radius ratio is $η=0.714$. The stabilization causes the system to behave as a damped oscillator and correspondingly the turbulence decays. The evolution of the torque and turbulent kinetic energy is analysed and the periodicity and damping of the oscillations are quantified and explained as a function of shear Reynolds number. Though the initially turbulent flow state decays, surprisingly, the system is found to absorb energy during this decay. △ Less

Submitted 11 November, 2013; originally announced November 2013.

Comments: Preprint submitted to PRL, 12 pages, 5 figures

Velocity profiles in strongly turbulent Taylor-Couette flow

Authors: Siegfried Grossmann, Detlef Lohse, Chao Sun

Abstract: 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 downw… ▽ More 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^+_{\varphi}(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. △ Less

Submitted 23 October, 2013; originally announced October 2013.

Boundary layer dynamics at the transition between the classical and the ultimate regime of Taylor-Couette flow

Authors: Rodolfo Ostilla Monico, Erwin P. van der Poel, Roberto Verzicco, Siegfried Grossmann, Detlef Lohse

Abstract: Direct numerical simulations of turbulent Taylor-Couette flow are performed up to inner cylinder Reynolds numbers of {Re_i=10^5} for a radius ratio of {η=r_i/r_o=0.714} between the inner and outer cylinder. With increasing {Re_i}, the flow undergoes transitions between three different regimes: (i) a flow dominated by large coherent structures, (ii) an intermediate transitional regime, and (iii) a… ▽ More Direct numerical simulations of turbulent Taylor-Couette flow are performed up to inner cylinder Reynolds numbers of {Re_i=10^5} for a radius ratio of {η=r_i/r_o=0.714} between the inner and outer cylinder. With increasing {Re_i}, the flow undergoes transitions between three different regimes: (i) a flow dominated by large coherent structures, (ii) an intermediate transitional regime, and (iii) a flow with developed turbulence. In the first regime the large--scale rolls completely drive the meridional flow while in the second one the coherent structures recover only on average. The presence of a mean flow allows for the coexistence of laminar and turbulent boundary layer dynamics. In the third regime the mean flow effects fade away and the flow becomes dominated by plumes. The effect of the local driving on the azimuthal and angular velocity profiles is quantified, in particular we show when and where those profiles develop. △ Less

Submitted 29 October, 2013; v1 submitted 12 June, 2013; originally announced June 2013.

Comments: 22 pages, submitted to PoF

Journal ref: Phys. Fluids 26, 015114 (2014)

Optimal Taylor-Couette flow: Radius ratio dependence

Authors: Rodolfo Ostilla Monico, Sander G. Huisman, Tim J. G. Jannink, Dennis P. M. van Gils, Roberto Verzicco, Siegfried Grossmann, Chao Sun, Detlef Lohse

Abstract: Taylor-Couette flow with independently rotating inner (i) and outer (o) cylinders is explored numerically and experimentally to determine the effects of the radius ratio η on the system response. Numerical simulations reach Reynolds numbers of up to Re_i=9.5 x 10^3 and Re_o=5x10^3, corresponding to Taylor numbers of up to Ta=10^8 for four different radius ratios η=r_i/r_o between 0.5 and 0.909. Th… ▽ More Taylor-Couette flow with independently rotating inner (i) and outer (o) cylinders is explored numerically and experimentally to determine the effects of the radius ratio η on the system response. Numerical simulations reach Reynolds numbers of up to Re_i=9.5 x 10^3 and Re_o=5x10^3, corresponding to Taylor numbers of up to Ta=10^8 for four different radius ratios η=r_i/r_o between 0.5 and 0.909. The experiments, performed in the Twente Turbulent Taylor-Couette (T^3C) setup, reach Reynolds numbers of up to Re_i=2x10^6$ and Re_o=1.5x10^6, corresponding to Ta=5x10^{12} for η=0.714-0.909. Effective scaling laws for the torque J^ω(Ta) are found, which for sufficiently large driving Ta are independent of the radius ratio η. As previously reported for η=0.714, optimum transport at a non-zero Rossby number Ro=r_i|ω_i-ω_o|/[2(r_o-r_i)ω_o] is found in both experiments and numerics. Ro_opt is found to depend on the radius ratio and the driving of the system. At a driving in the range between {Ta\sim3\cdot10^8} and {Ta\sim10^{10}}, Ro_opt saturates to an asymptotic η-dependent value. Theoretical predictions for the asymptotic value of Ro_{opt} are compared to the experimental results, and found to differ notably. Furthermore, the local angular velocity profiles from experiments and numerics are compared, and a link between a flat bulk profile and optimum transport for all radius ratios is reported. △ Less

Submitted 29 October, 2013; v1 submitted 23 April, 2013; originally announced April 2013.

Comments: Submitted to JFM, 28 pages, 17 figures

Journal ref: Journal of Fluid Mechanics / Volume 747 / May 2014, pp 1- 29

Superfluorescent emission in electrically pumped semiconductor laser

Authors: D. L. Boiko, X. Zeng, T. Stadelmann, S. Grossmann, A. Hoogerwerf, T. Weig, U. T. Schwarz, L. Sulmoni, J. -M. Lamy, N. Grandjean

Abstract: We report superfluorescent (SF) emission in electrically pumped InGaN/InGaN QW lasers with saturable absorber. In particular, we observe a superlinear growth of the peak power of SF pulses with increasing amplitude of injected current pulses and attribute it to cooperative pairing of electron-hole (e-h) radiative recombinations. The phase transitions from amplified spontaneous emission to superflu… ▽ More We report superfluorescent (SF) emission in electrically pumped InGaN/InGaN QW lasers with saturable absorber. In particular, we observe a superlinear growth of the peak power of SF pulses with increasing amplitude of injected current pulses and attribute it to cooperative pairing of electron-hole (e-h) radiative recombinations. The phase transitions from amplified spontaneous emission to superfluorescence and then to lasing regime is confirmed by observing (i) abrupt peak power growth accompanied by spectral broadening, (ii) spectral shape with hyperbolic secant envelope and (iii) red shift of central wavelength of SF emission pulse. The observed red shift of SF emission is shown to be caused by the pairing of e-h pairs in an indirect cooperative X-transition. △ Less

Submitted 11 February, 2013; v1 submitted 1 February, 2013; originally announced February 2013.

Comments: 18 pages, 7 figures

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

Authors: Richard J. A. M. Stevens, Erwin P. van der Poel, Siegfried Grossmann, Detlef Lohse

Abstract: 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… ▽ More 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. △ Less

Submitted 20 June, 2013; v1 submitted 29 January, 2013; originally announced January 2013.

Comments: 14 pages, 8 figures

Journal ref: J. Fluid Mech. 730 , 295-308 (2013)

Logarithmic temperature profiles in the ultimate regime of thermal convection

Authors: Siegfried Grossmann, Detlef Lohse

Abstract: 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… ▽ More 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. △ Less

Submitted 13 August, 2012; originally announced August 2012.

Comments: 14 pages, no figure

Optimal Taylor-Couette flow: direct numerical simulations

Authors: Rodolfo Ostilla Mónico, Richard J. A. M. Stevens, Siegfried Grossmann, Roberto Verzicco, Detlef Lohse

Abstract: We numerically simulate turbulent Taylor-Couette flow for independently rotating inner and outer cylinders, focusing on the analogy with turbulent Rayleigh-Bénard flow. Reynolds numbers of $Re_i=8\cdot10^3$ and $Re_o=\pm4\cdot10^3$ of the inner and outer cylinders, respectively, are reached, corresponding to Taylor numbers Ta up to $10^8$. Effective scaling laws for the torque and other system res… ▽ More We numerically simulate turbulent Taylor-Couette flow for independently rotating inner and outer cylinders, focusing on the analogy with turbulent Rayleigh-Bénard flow. Reynolds numbers of $Re_i=8\cdot10^3$ and $Re_o=\pm4\cdot10^3$ of the inner and outer cylinders, respectively, are reached, corresponding to Taylor numbers Ta up to $10^8$. Effective scaling laws for the torque and other system responses are found. Recent experiments with the Twente turbulent Taylor-Couette ($T^3C$) setup and with a similar facility in Maryland at very high Reynolds numbers have revealed an optimum transport at a certain non-zero rotation rate ratio $a = -ω_o / ω_i$ of about $a_{opt}=0.33-0.35$. For large enough $Ta$ in the numerically accessible range we also find such an optimum transport at non-zero counter-rotation. The position of this maximum is found to shift with the driving, reaching a maximum of $a_{opt}=0.15$ for $Ta=2.5\cdot10^7$. An explanation for this shift is elucidated, consistent with the experimental result that $a_{opt}$ becomes approximately independent of the driving strength for large enough Reynolds numbers. We furthermore numerically calculate the angular velocity profiles and visualize the different flow structures for the various regimes. By writing the equations in a frame co-rotating with the outer cylinder a link is found between the local angular velocity profiles and the global transport quantities. △ Less

Submitted 31 October, 2012; v1 submitted 10 July, 2012; originally announced July 2012.

Comments: Under consideration for publication in JFM, 31 pages, 25 figures

Journal ref: J. Fluid Mech. 719, 14-46 (2013)

Logarithmic temperature profiles in turbulent Rayleigh-Bénard convection

Authors: Guenter Ahlers, Eberhard Bodenschatz, Denis Funfschilling, Siegfried Grossmann, Xiaozhou He, Detlef Lohse, Richard J. A. M. Stevens, Roberto Verzicco

Abstract: 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… ▽ More 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. △ Less

Submitted 29 April, 2012; originally announced April 2012.

Journal ref: Phys. Rev. Lett. 109, 114501 (2012)

Atomic-scale confinement of optical fields

Authors: Johannes Kern, Swen Grossmann, Nadezda V. Tarakina, Tim Häckel, Monika Emmerling, Martin Kamp, Jer-Shing Huang, Paolo Biagioni, Jord C. Prangsma, Bert Hecht

Abstract: In the presence of matter there is no fundamental limit preventing confinement of visible light even down to atomic scales. Achieving such confinement and the corresponding intensity enhancement inevitably requires simultaneous control over atomic-scale details of material structures and over the optical modes that such structures support. By means of self-assembly we have obtained side-by-side al… ▽ More In the presence of matter there is no fundamental limit preventing confinement of visible light even down to atomic scales. Achieving such confinement and the corresponding intensity enhancement inevitably requires simultaneous control over atomic-scale details of material structures and over the optical modes that such structures support. By means of self-assembly we have obtained side-by-side aligned gold nanorod dimers with robust atomically-defined gaps reaching below 0.5 nm. The existence of atomically-confined light fields in these gaps is demonstrated by observing extreme Coulomb splitting of corresponding symmetric and anti-symmetric dimer eigenmodes of more than 800 meV in white-light scattering experiments. Our results open new perspectives for atomically-resolved spectroscopic imaging, deeply nonlinear optics, ultra-sensing, cavity optomechanics as well as for the realization of novel quantum-optical devices. △ Less

Submitted 24 May, 2012; v1 submitted 21 December, 2011; originally announced December 2011.

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

Authors: Rajaram Lakkaraju, Richard. J. A. M. Stevens, Roberto Verzicco, Siegfried Grossmann, Andrea Prosperetti, Chao Sun, Detlef Lohse

Abstract: 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 c… ▽ More 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. △ Less

Submitted 16 April, 2012; v1 submitted 2 December, 2011; originally announced December 2011.

Comments: 10 pages, 14 figures, 2 tables

Optimal Taylor-Couette turbulence

Authors: Dennis P. M. van Gils, Sander G. Huisman, Siegfried Grossmann, Chao Sun, Detlef Lohse

Abstract: 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}… ▽ More 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. △ Less

Submitted 2 October, 2012; v1 submitted 27 November, 2011; originally announced November 2011.

Comments: 30 pages, 22 figures

Journal ref: J. Fluid Mech. 706, 118 (2012)

Ultimate Turbulent Taylor-Couette Flow

Authors: Sander G. Huisman, Dennis P. M. van Gils, Siegfried Grossmann, Chao Sun, Detlef Lohse

Abstract: 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 tur… ▽ More 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. △ Less

Submitted 31 October, 2011; originally announced November 2011.

Journal ref: Phys. Rev. Lett. 108, 024501 (2012)

Thermal boundary layer profiles in turbulent Rayleigh-Bénard convection in a cylindrical sample

Authors: Richard J. A. M. Stevens, Quan Zhou, Siegfried Grossmann, Roberto Verzicco, Ke-Qing Xia, Detlef Lohse

Abstract: 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-di… ▽ More 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. △ Less

Submitted 8 February, 2012; v1 submitted 28 October, 2011; originally announced October 2011.

Comments: 5 pages, 2 figures

Journal ref: Phys. Rev. E 85, 027301 (2012)

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

Authors: Quan Zhou, Kazuyasu Sugiyama, Richard J. A. M. Stevens, Siegfried Grossmann, Detlef Lohse, Ke-Qing Xia

Abstract: 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 mo… ▽ More 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. △ Less

Submitted 30 September, 2011; originally announced September 2011.

Comments: 16 pages, 16 figures

Journal ref: Phys. Fluids 23, 125104 (2011)

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

Authors: Olga Shishkina, Richard J. A. M. Stevens, Siegfried Grossmann, Detlef Lohse

Abstract: 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… ▽ More 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. △ Less

Submitted 30 September, 2011; originally announced September 2011.

Journal ref: New J. Phys. 12 075022 (2010)

Flow reversals in thermally driven turbulence

Authors: Kazuyasu Sugiyama, Rui Ni, Richard J. A. M. Stevens, Tak Shing Chan, Sheng-Qi Zhou, Heng-Dong Xi, Chao Sun, Siegfried Grossmann, Ke-Qing Xia, Detlef Lohse

Abstract: 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 t… ▽ More 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. △ Less

Submitted 30 September, 2011; originally announced September 2011.

Comments: 4 pages, 4 figures

Journal ref: Phys. Rev. Lett. 105, 034503 (2010)

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

Authors: Quan Zhou, Richard J. A. M. Stevens, Kazuyasu Sugiyama, Siegfried Grossmann, Detlef Lohse, Ke-Qing Xia

Abstract: 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 Prandt… ▽ More 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. △ Less

Submitted 5 February, 2010; originally announced February 2010.

Comments: 10 pages, 6 figures

Journal ref: J. Fluid Mech. 664, 297-312 (2010)

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

Authors: Kazuyasu Sugiyama, Enrico Calzavarini, Siegfried Grossmann, Detlef Lohse

Abstract: 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; the… ▽ More 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. △ Less

Submitted 20 December, 2008; originally announced December 2008.

Comments: 31 pages, 22 figures

Journal ref: J. Fluid Mech. 637, pp. 105-135, (2009)

Large scale dynamics in turbulent Rayleigh-Benard convection

Authors: Guenter Ahlers, Siegfried Grossmann, Detlef Lohse

Abstract: 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 l… ▽ More 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. △ Less

Submitted 4 November, 2008; originally announced November 2008.

Comments: Review article, 34 pages, 13 figures, Rev. Mod. Phys. 81, in press (2009)

Report number: Rev.Mod.Phys.81-2009

Non-Oberbeck-Boussinesq effects in turbulent thermal convection in ethane close to the critical point

Authors: Guenter Ahlers, Enrico Calzavarini, Francisco Fontenele Araujo, Denis Funfschilling, Siegfried Grossmann, Detlef Lohse, Kazuyasu Sugiyama

Abstract: As shown in earlier work (Ahlers et al., J. Fluid Mech. 569, p.409 (2006)), non-Oberbeck Boussinesq (NOB) corrections to the center temperature in turbulent Rayleigh-Benard convection in water and also in glycerol are governed by the temperature dependences of the kinematic viscosity and the thermal diffusion coefficient. If the working fluid is ethane close to the critical point the origin of n… ▽ More As shown in earlier work (Ahlers et al., J. Fluid Mech. 569, p.409 (2006)), non-Oberbeck Boussinesq (NOB) corrections to the center temperature in turbulent Rayleigh-Benard convection in water and also in glycerol are governed by the temperature dependences of the kinematic viscosity and the thermal diffusion coefficient. If the working fluid is ethane close to the critical point the origin of non-Oberbeck-Boussinesq corrections is very different, as will be shown in the present paper. Namely, the main origin of NOB corrections then lies in the strong temperature dependence of the isobaric thermal expansion coefficient β(T). More precisely, it is the nonlinear T-dependence of the density ρ(T) in the buoyancy force which causes another type of NOB effect. We demonstrate that through a combination of experimental, numerical, and theoretical work, the latter in the framework of the extended Prandtl-Blasius boundary layer theory developed in Ahlers et al., J. Fluid Mech. 569, p.409 (2006). The latter comes to its limits, if the temperature dependence of the thermal expension coefficient β(T) is significant. △ Less

Submitted 22 February, 2008; v1 submitted 2 December, 2007; originally announced December 2007.

Comments: 18 pages, 15 figures, 3 tables

Journal ref: Phys. Rev. E 77, 046302 (2008)

Non-Oberbeck-Boussinesq effects in two-dimensional Rayleigh-Benard convection in glycerol

Authors: Kazuyasu Sugiyama, Enrico Calzavarini, Siegfried Grossmann, Detlef Lohse

Abstract: We numerically analyze Non-Oberbeck-Boussinesq (NOB) effects in two-dimensional Rayleigh-Benard flow in glycerol, which shows a dramatic change in the viscosity with temperature. The results are presented both as functions of the Rayleigh number (Ra) up to $10^8$ (for fixed temperature difference between the top and bottom plates) and as functions of "non-Oberbeck-Boussinesqness'' or "NOBness''… ▽ More We numerically analyze Non-Oberbeck-Boussinesq (NOB) effects in two-dimensional Rayleigh-Benard flow in glycerol, which shows a dramatic change in the viscosity with temperature. The results are presented both as functions of the Rayleigh number (Ra) up to $10^8$ (for fixed temperature difference between the top and bottom plates) and as functions of "non-Oberbeck-Boussinesqness'' or "NOBness'' ($Δ$) up to 50 K (for fixed Ra). For this large NOBness the center temperature $T_c$ is more than 5 K larger than the arithmetic mean temperature $T_m$ between top and bottom plate and only weakly depends on Ra. To physically account for the NOB deviations of the Nusselt numbers from its Oberbeck-Boussinesq values, we apply the decomposition of $Nu_{NOB}/Nu_{OB}$ into the product of two effects, namely first the change in the sum of the top and bottom thermal BL thicknesses, and second the shift of the center temperature $T_c$ as compared to $T_m$. While for water the origin of the $Nu$ deviation is totally dominated by the second effect (cf. Ahlers et al., J. Fluid Mech. 569, pp. 409 (2006)) for glycerol the first effect is dominating, in spite of the large increase of $T_c$ as compared to $T_m$. △ Less

Submitted 2 October, 2007; v1 submitted 1 July, 2007; originally announced July 2007.

Comments: 6 pages, 7 figures

Journal ref: Europhys. Lett. 80, 34002 (2007)

Some generic properties of level spacing distributions of 2D real random matrices

Authors: S. Grossmann, M. Robnik

Abstract: We study the level spacing distribution $P(S)$ of 2D real random matrices both symmetric as well as general, non-symmetric. In the general case we restrict ourselves to Gaussian distributed matrix elements, but different widths of the various matrix elements are admitted. The following results are obtained: An explicit exact formula for $P(S)$ is derived and its behaviour close to S=0 is studied… ▽ More We study the level spacing distribution $P(S)$ of 2D real random matrices both symmetric as well as general, non-symmetric. In the general case we restrict ourselves to Gaussian distributed matrix elements, but different widths of the various matrix elements are admitted. The following results are obtained: An explicit exact formula for $P(S)$ is derived and its behaviour close to S=0 is studied analytically, showing that there is linear level repulsion, unless there are additional constraints for the probability distribution of the matrix elements. The constraint of having only positive or only negative but otherwise arbitrary non-diagonal elements leads to quadratic level repulsion with logarithmic corrections. These findings detail and extend our previous results already published in a preceding paper. For the {\em symmetric} real 2D matrices also other, non-Gaussian statistical distributions are considered. In this case we show for arbitrary statistical distribution of the diagonal and non-diagonal elements that the level repulsion exponent $ρ$ is always $ρ= 1$, provided the distribution function of the matrix elements is regular at zero value. If the distribution function of the matrix elements is a singular (but still integrable) power law near zero value of $S$, the level spacing distribution $P(S)$ is a fractional exponent pawer law at small $S$. The tail of $P(S)$ depends on further details of the matrix element statistics. We explicitly work out four cases: the constant (box) distribution, the Cauchy-Lorentz distribution, the exponential distribution and, as an example for a singular distribution, the power law distribution for $P(S)$ near zero value times an exponential tail. △ Less

Submitted 27 June, 2007; originally announced June 2007.

Comments: 21 pages, no figures, submitted to Zeitschrift fuer Naturforschung A

Report number: CAMTP Preprint 5 April 2007

Extended phase diagram of the Lorenz model

Authors: Holger R. Dullin, Sven Schmidt, Peter H. Richter, Siegfried K. Grossmann

Abstract: The parameter dependence of the various attractive solutions of the three variable nonlinear Lorenz model equations for thermal convection in Rayleigh-Bénard flow is studied. Its bifurcation structure has commonly been investigated as a function of r, the normalized Rayleigh number, at fixed Prandtl number σ. The present work extends the analysis to the entire (r,σ) parameter plane. An onion lik… ▽ More The parameter dependence of the various attractive solutions of the three variable nonlinear Lorenz model equations for thermal convection in Rayleigh-Bénard flow is studied. Its bifurcation structure has commonly been investigated as a function of r, the normalized Rayleigh number, at fixed Prandtl number σ. The present work extends the analysis to the entire (r,σ) parameter plane. An onion like periodic pattern is found which is due to the alternating stability of symmetric and non-symmetric periodic orbits. This periodic pattern is explained by considering non-trivial limits of large r and σ. In addition to the limit which was previously analyzed by Sparrow, we identify two more distinct asymptotic regimes in which either σ/r or σ^2/r is constant. In both limits the dynamics is approximately described by Airy functions whence the periodicity in parameter space can be calculated analytically. Furthermore, some observations about sequences of bifurcations and coexistence of attractors, periodic as well as chaotic, are reported. △ Less

Submitted 12 April, 2005; originally announced April 2005.

Comments: 36 pages, 20 figures

Journal ref: Internat. J. Bifur. Chaos Appl. Sci. Engrg. 17, 3013-3033 (2007)

Wind reversals in turbulent Rayleigh-Benard convection

Authors: Francisco Fontenele Araujo, S. Grossmann, D. Lohse

Abstract: The phenomenon of irregular cessation and subsequent reversal of the large-scale circulation in turbulent Rayleigh-Bénard convection is theoretically analysed. The force and thermal balance on a single plume detached from the thermal boundary layer yields a set of coupled nonlinear equations, whose dynamics is related to the Lorenz equations. For Prandtl and Rayleigh numbers in the range… ▽ More The phenomenon of irregular cessation and subsequent reversal of the large-scale circulation in turbulent Rayleigh-Bénard convection is theoretically analysed. The force and thermal balance on a single plume detached from the thermal boundary layer yields a set of coupled nonlinear equations, whose dynamics is related to the Lorenz equations. For Prandtl and Rayleigh numbers in the range $10^{-2} \leq \Pr \leq 10^{3}$ and $10^{7} \leq \Ra \leq 10^{12}$, the model has the following features: (i) chaotic reversals may be exhibited at Ra $\geq 10^{7}$; (ii) the Reynolds number based on the root mean square velocity scales as $\Re_{rms} \sim \Ra^{[0.41 ... 0.47]}$ (depending on Pr), and as $\Re_{rms} \sim \Pr^{-[0.66 ... 0.76]}$ (depending on Ra); and (iii) the mean reversal frequency follows an effective scaling law $ω/ (νL^{-2}) \sim \Pr^{-(0.64 \pm 0.01)} \Ra^{0.44 \pm 0.01}$. The phase diagram of the model is sketched, and the observed transitions are discussed. △ Less

Submitted 29 August, 2005; v1 submitted 13 July, 2004; originally announced July 2004.

Comments: 4 pages, 5 figures

Journal ref: Phys. Rev. Lett. 95, 084502 (2005)

Response maxima in modulated turbulence Part II: Numerical simulations

Authors: A. von der Heydt, S. Grossmann, D. Lohse

Abstract: Numerical simulations of fully developed turbulence driven by a modulated energy input rate or driving force are performed within two dynamical cascade models, the GOY shell model and a reduced wave vector set approximation of the Navier-Stokes equation (REWA). The frequency behavior of the system response is studied and compared with predictions from a variable range mean-field theory, which ex… ▽ More Numerical simulations of fully developed turbulence driven by a modulated energy input rate or driving force are performed within two dynamical cascade models, the GOY shell model and a reduced wave vector set approximation of the Navier-Stokes equation (REWA). The frequency behavior of the system response is studied and compared with predictions from a variable range mean-field theory, which excludes turbulent fluctuations. In agreement with the mean-field approach we find a constant response amplitude for low driving frequencies and a 1/w-decay of the amplitude for high frequencies. In the mean-field theory, the finite cascade time scale had lead to an oscillating behavior of the response amplitude as a function of the driving frequency. In the simulations of both models we observe the main maximum. The higher maxima and minima are completely washed out by fluctuations. △ Less

Submitted 23 May, 2003; originally announced May 2003.

Comments: 12 pages, 10 figures, submitted to Phys.Rev. E

Journal ref: Phys. Rev. E 68, 066302 (2003)

Response maxima in modulated turbulence

Authors: A. von der Heydt, S. Grossmann, D. Lohse

Abstract: Isotropic and homogeneous turbulence driven by an energy input modulated in time is studied within a variable range mean-field theory. The response of the system, observed in the second order moment of the large-scale velocity difference D(L,t)=<<(\u(x+L)-\u(x))^2>>~Re(t)^2$, is calculated for varying modulation frequencies w and weak modulation amplitudes. For low frequencies the system follows… ▽ More Isotropic and homogeneous turbulence driven by an energy input modulated in time is studied within a variable range mean-field theory. The response of the system, observed in the second order moment of the large-scale velocity difference D(L,t)=<<(\u(x+L)-\u(x))^2>>~Re(t)^2$, is calculated for varying modulation frequencies w and weak modulation amplitudes. For low frequencies the system follows the modulation of the driving with almost constant amplitude, whereas for higher driving frequencies the amplitude of the response decreases on average 1/w. In addition, at certain frequencies the amplitude of the response either almost vanishes or is strongly enhanced. These frequencies are connected with the frequency scale of the energy cascade and multiples thereof. △ Less

Submitted 23 May, 2003; originally announced May 2003.

Comments: 11 pages, 6 figures

Journal ref: Phys. Rev. E 67, 046308 (2003)

On geometry effects in Rayleigh-Benard convection

Authors: Siegfried Grossmann, Detlef Lohse

Abstract: Various recent experiments hint at a geometry dependence of scaling relations in Rayleigh-Bénard convection. Aspect ratio and shape dependences have been found. In this paper a mechanism is offered which can account for such dependences. It is based on Prandtl's theory for laminar boundary layers and on the conservation of volume flux of the large scale wind. The mechanism implies the possibilit… ▽ More Various recent experiments hint at a geometry dependence of scaling relations in Rayleigh-Bénard convection. Aspect ratio and shape dependences have been found. In this paper a mechanism is offered which can account for such dependences. It is based on Prandtl's theory for laminar boundary layers and on the conservation of volume flux of the large scale wind. The mechanism implies the possibility of different thicknesses of the kinetic boundary layers at the sidewalls and the top/bottom plates, just as experimentally found by Qiu and Xia (Phys. Rev. E58, 486 (1998)), and also different $Ra$-scaling of the wind measured over the plates and at the sidewalls. In the second part of the paper a scaling argument for the velocity and temperature fluctuations in the bulk is developeVarious recent experiments hint at a geometry dependence of scaling relations in Rayleigh-Benard convection. Aspect ratio and shape dependences have been found. In this paper a mechanism is offered which can account for such dependences. It is based on Prandtl's theory for laminar boundary layers and on the conservation of volume flux of the large scale wind. The mechanism implies the possibility of different thicknesses of the kinetic boundary layers at the sidewalls and the top/bottom plates, just as experimentally found by Qiu and Xia (Phys. Rev. E58, 486 (1998)), and also different $Ra$-scaling of the wind measured over the plates and at the sidewalls. In the second part of the paper a scaling argument for the velocity and temperature fluctuations in the bulk is develope △ Less

Submitted 5 December, 2002; originally announced December 2002.

Comments: 4 pages, 1 figure

Thermal convection for large Prandtl numbers

Authors: Siegfried Grossmann, Detlef Lohse

Abstract: The Rayleigh-Benard theory by Grossmann and Lohse [J. Fluid Mech. 407, 27 (2000)] is extended towards very large Prandtl numbers Pr. The Nusselt number Nu is found here to be independent of Pr. However, for fixed Rayleigh numbers Ra> 10^{10} a maximum around Pr= 2 in the Nu(Pr)-dependence is predicted which is absent for lower Ra. We moreover offer the full functional dependences of Nu(Ra,Pr) an… ▽ More The Rayleigh-Benard theory by Grossmann and Lohse [J. Fluid Mech. 407, 27 (2000)] is extended towards very large Prandtl numbers Pr. The Nusselt number Nu is found here to be independent of Pr. However, for fixed Rayleigh numbers Ra> 10^{10} a maximum around Pr= 2 in the Nu(Pr)-dependence is predicted which is absent for lower Ra. We moreover offer the full functional dependences of Nu(Ra,Pr) and Re(Ra,Pr) within this extended theory, rather than only giving the limiting power laws as done previously. This enables us to more realistically describe the transitions between the various scaling regimes, including their widths. △ Less

Submitted 8 November, 2000; originally announced November 2000.

Comments: 7 pages

Journal ref: Phys. Rev. Lett. 86, 3316-3319 (2001)

Scaling of global momentum transport in Taylor-Couette and pipe flow

Authors: Bruno Eckhardt, Siegfried Grossmann, Detlef Lohse

Abstract: We interpret measurements of the Reynolds number dependence of the torque in Taylor-Couette flow by Lewis and Swinney [Phys. Rev. E 59, 5457 (1999)] and of the pressure drop in pipe flow by Smits and Zagarola, [Phys. Fluids 10, 1045 (1998)] within the scaling theory of Grossmann and Lohse [J. Fluid Mech. 407, 27 (2000)], developed in the context of thermal convection. The main idea is to split t… ▽ More We interpret measurements of the Reynolds number dependence of the torque in Taylor-Couette flow by Lewis and Swinney [Phys. Rev. E 59, 5457 (1999)] and of the pressure drop in pipe flow by Smits and Zagarola, [Phys. Fluids 10, 1045 (1998)] within the scaling theory of Grossmann and Lohse [J. Fluid Mech. 407, 27 (2000)], developed in the context of thermal convection. The main idea is to split the energy dissipation into contributions from a boundary layer and the turbulent bulk. This ansatz can account for the observed scaling in both cases if it is assumed that the internal wind velocity $U_w$ introduced through the rotational or pressure forcing is related to the the external (imposed) velocity U, by $U_w/U \propto Re^ξ$ with xi = -0.051 and xi = -0.041 for the Taylor-Couette (U inner cylinder velocity) and pipe flow (U mean flow velocity) case, respectively. In contrast to the Rayleigh-Benard case the scaling exponents cannot (yet) be derived from the dynamical equations. △ Less

Submitted 25 September, 2000; originally announced September 2000.

Comments: 7 pages, 4 ps figures with 4 program files included in the source. European Physical Journal B, accepted

Journal ref: Eur. Phys. J. B 18, 541-544 (2000).

Scaling in thermal convection: A unifying theory

Authors: Siegfried Grossmann, Detlef Lohse

Abstract: A systematic theory for the scaling of the Nusselt number $Nu$ and of the Reynolds number $Re$ in strong Rayleigh-Benard convection is suggested and shown to be compatible with recent experiments. It assumes a coherent large scale convection roll (``wind of turbulence'') and is based on the dynamical equations both in the bulk and in the boundary layers. Several regimes are identified in the Ray… ▽ More A systematic theory for the scaling of the Nusselt number $Nu$ and of the Reynolds number $Re$ in strong Rayleigh-Benard convection is suggested and shown to be compatible with recent experiments. It assumes a coherent large scale convection roll (``wind of turbulence'') and is based on the dynamical equations both in the bulk and in the boundary layers. Several regimes are identified in the Rayleigh number versus Prandtl number phase space, defined by whether the boundary layer or the bulk dominates the global kinetic and thermal dissipation, respectively. The crossover between the regimes is calculated. In the regime which has most frequently been studied in experiment (Ra smaller than 10^{11}) the leading terms are $Nu\sim Ra^{1/4}Pr^{1/8}$, $Re \sim Ra^{1/2} Pr^{-3/4}$ for $Pr < 1$ and $Nu\sim Ra^{1/4}Pr^{-1/12}$, $Re \sim Ra^{1/2} Pr^{-5/6}$ for $Pr > 1$. In most measurements these laws are modified by additive corrections from the neighboring regimes so that the impression of a slightly larger (effective) Nu vs Ra scaling exponent can arise. -- The presented theory is best summarized in the phase diagram figure 1. △ Less

Submitted 21 September, 1999; originally announced September 1999.

Comments: 30 pages, latex, 7 figures, under review at Journal of Fluid Mech

Report number: MR-TW-99-3

Continued Fraction Representation of Temporal Multi Scaling in Turbulence

Authors: David Daems, Siegfried Grossmann, Victor S. L'vov, Itamar Procaccia

Abstract: 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… ▽ More 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. △ Less

Submitted 25 November, 1998; originally announced November 1998.

Comments: 7 pages, no figure

Analysis of Rayleigh-Plesset dynamics for sonoluminescing bubbles

Authors: Sascha Hilgenfeldt, Michael P. Brenner, Siegfried Grossmann, Detlef Lohse

Abstract: Recent work on single bubble sonoluminescence (SBSL) has shown that many features of this phenomenon, especially the dependence of SBSL intensity and stability on experimental parameters, can be explained within a hydrodynamic approach. More specifically, many important properties can already be derived from an analysis of bubble wall dynamics. This dynamics is conveniently described by the Rayl… ▽ More Recent work on single bubble sonoluminescence (SBSL) has shown that many features of this phenomenon, especially the dependence of SBSL intensity and stability on experimental parameters, can be explained within a hydrodynamic approach. More specifically, many important properties can already be derived from an analysis of bubble wall dynamics. This dynamics is conveniently described by the Rayleigh-Plesset (RP) equation. In this work we derive analytical approximations for RP dynamics and subsequent analytical laws for parameter dependences. These results include (i) an expression for the onset threshold of SL, (ii) an analytical explanation of the transition from diffusively unstable to stable equilibria for the bubble ambient radius (unstable and stable sonoluminescence), and (iii) a detailed understanding of the resonance structure of the RP equation. It is found that the threshold for SL emission is shifted to larger bubble radii and larger driving pressures if surface tension is enlarged, whereas even a considerable change in liquid viscosity leaves this threshold virtually unaltered. As an enhanced viscosity stabilizes the bubbles against surface oscillations, we conclude that the ideal liquid for violently collapsing, surface stable SL bubbles should have small surface tension and large viscosity, although too large viscosity (>40 times the viscosity of water) will again preclude collapses. △ Less

Submitted 30 May, 1998; v1 submitted 18 May, 1998; originally announced May 1998.

Comments: 41 pages, 21 eps and ps figures; LaTeX stylefiles replaced because the PostScript file produced at the archive had misplaced and misscaled figures

Journal ref: J. Fluid Mech. 365, 171-204 (1998)

Variational bound on energy dissipation in turbulent shear flow

Authors: Rolf Nicodemus, Siegfried Grossmann, Martin Holthaus

Abstract: We present numerical solutions to the extended Doering-Constantin variational principle for upper bounds on the energy dissipation rate in plane Couette flow, bridging the entire range from low to asymptotically high Reynolds numbers. Our variational bound exhibits structure, namely a pronounced minimum at intermediate Reynolds numbers, and recovers the Busse bound in the asymptotic regime. The… ▽ More We present numerical solutions to the extended Doering-Constantin variational principle for upper bounds on the energy dissipation rate in plane Couette flow, bridging the entire range from low to asymptotically high Reynolds numbers. Our variational bound exhibits structure, namely a pronounced minimum at intermediate Reynolds numbers, and recovers the Busse bound in the asymptotic regime. The most notable feature is a bifurcation of the minimizing wavenumbers, giving rise to simple scaling of the optimized variational parameters, and of the upper bound, with the Reynolds number. △ Less

Submitted 24 April, 1998; originally announced April 1998.

Comments: 4 pages, RevTeX, 5 postscript figures are available as one .tar.gz file from holthaus@stat.physik.uni-marburg.de

Report number: TUR-97-05

Journal ref: Phys. Rev. Lett. 79, 4170 (1997)

Variational bound on energy dissipation in plane Couette flow

Authors: Rolf Nicodemus, Siegfried Grossmann, Martin Holthaus

Abstract: We present numerical solutions to the extended Doering-Constantin variational principle for upper bounds on the energy dissipation rate in turbulent plane Couette flow. Using the compound matrix technique in order to reformulate this principle's spectral constraint, we derive a system of equations that is amenable to numerical treatment in the entire range from low to asymptotically high Reynold… ▽ More We present numerical solutions to the extended Doering-Constantin variational principle for upper bounds on the energy dissipation rate in turbulent plane Couette flow. Using the compound matrix technique in order to reformulate this principle's spectral constraint, we derive a system of equations that is amenable to numerical treatment in the entire range from low to asymptotically high Reynolds numbers. Our variational bound exhibits a minimum at intermediate Reynolds numbers, and reproduces the Busse bound in the asymptotic regime. As a consequence of a bifurcation of the minimizing wavenumbers, there exist two length scales that determine the optimal upper bound: the effective width of the variational profile's boundary segments, and the extension of their flat interior part. △ Less

Submitted 24 April, 1998; originally announced April 1998.

Comments: 22 pages, RevTeX, 11 postscript figures are available as one uuencoded .tar.gz file from holthaus@stat.physik.uni-marburg.de

Report number: TUR-97-04

Journal ref: Phys. Rev. E 56, 6774 (1997)

Improved variational principle for bounds on energy dissipation in turbulent shear flow

Authors: Rolf Nicodemus, Siegfried Grossmann, Martin Holthaus

Abstract: We extend the Doering-Constantin approach to upper bounds on energy dissipation in turbulent flows by introducing a balance parameter into the variational principle. This parameter governs the relative weight of different contributions to the dissipation rate. Its optimization leads to improved bounds without entailing additional technical difficulties. For plane Couette flow, the high-Re-bounds… ▽ More We extend the Doering-Constantin approach to upper bounds on energy dissipation in turbulent flows by introducing a balance parameter into the variational principle. This parameter governs the relative weight of different contributions to the dissipation rate. Its optimization leads to improved bounds without entailing additional technical difficulties. For plane Couette flow, the high-Re-bounds obtainable with one-dimensional background flows are methodically lowered by a factor of 27/32. △ Less

Submitted 24 April, 1998; originally announced April 1998.

Comments: 15 pages, RevTeX, 3 postscript figures

Report number: TUR-97-01

Journal ref: Physica D 101, 178 (1997)

Nonperturbative renormalization group approach to turbulence

Authors: Alexander Esser, Siegfried Grossmann

Abstract: We suggest a new, renormalization group (RG) based, nonperturbative method for treating the intermittency problem of fully developed turbulence which also includes the effects of a finite boundary of the turbulent flow. The key idea is not to try to construct an elimination procedure based on some assumed statistical distribution, but to make an ansatz for possible RG transformations and to pose… ▽ More We suggest a new, renormalization group (RG) based, nonperturbative method for treating the intermittency problem of fully developed turbulence which also includes the effects of a finite boundary of the turbulent flow. The key idea is not to try to construct an elimination procedure based on some assumed statistical distribution, but to make an ansatz for possible RG transformations and to pose constraints upon those, which guarantee the invariance of the nonlinear term in the Navier-Stokes equation, the invariance of the energy dissipation, and other basic properties of the velocity field. The role of length scales is taken to be inverse to that in the theory of critical phenomena; thus possible intermittency corrections are connected with the outer length scale. Depending on the specific type of flow, we find different sets of admissible transformations with distinct scaling behavior: for the often considered infinite, isotropic, and homogeneous system K41 scaling is enforced, but for the more realistic plane Couette geometry no restrictions on intermittency exponents were obtained so far. △ Less

Submitted 28 December, 1997; originally announced January 1998.

Comments: 45 pages, 2 figures