Researcher profile

Joerg Schumacher

Joerg Schumacher contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2021arXiv

Connecting boundary layer dynamics with extreme bulk dissipation events in Rayleigh-Bénard flow

We study the connection between extreme events of thermal and kinetic energy dissipation rates in the bulk of three-dimensional Rayleigh-Bénard convection and the wall shear stress patterns at the top and the bottom planes that enclose the layer. Zero points of this two-dimensional vector field stand for detachments of strong thermal plumes. If their position at the opposite planes and a given time is close then they can be considered as precursors for high-amplitude bulk dissipation events triggered by plume collisions or close passings. This scenario requires a breaking of the synchronicity of the boundary layer dynamics at both plates which is found to be in line with a transition of the bulk derivative statistics from Gaussian to intermittent. Our studies are based on three-dimensional high-resolution direct numerical simulations for moderate Rayleigh numbers between $Ra=10^4$ and $5\times 10^5$.

preprint2020arXiv

Rayleigh-Bénard Convection in Strong Vertical Magnetic Field: Flow Structure and Verification of Numerical Method

Direct numerical simulations are performed to study turbulent Rayleigh-Bénard convection in a vertical cylindrical cavity with uniform axial magnetic field. Flows at high Hartmann and Rayleigh numbers are considered. The calculations reveal that, similarly to the behavior observed in Rayleigh-Bénard convection with strong rotation, flows at strong magnetic field develop a central vortex, while the heat transfer is suppressed.

preprint2020arXiv

Turbulent Rayleigh-Bénard convection in a strong vertical magnetic field

Direct numerical simulations are carried out to study flow structure and transport properties in turbulent Rayleigh-Bénard convection in a cylindrical cell of aspect ratio one with an imposed axial magnetic field. Flows at the Prandtl number 0.025 and the Rayleigh and Hartmann numbers up to $10^9$ and 1400 are considered. The results are consistent with those of earlier experimental and numerical data. As anticipated, the heat transfer rate and kinetic energy are suppressed by strong magnetic field. At the same time, their growth with the Rayleigh number is found to be faster in flows at high Hartmann numbers. This behaviour is attributed to the newly discovered flow regime characterized by prominent quasi two-dimensional structures reminiscent of vortex sheets observed earlier in simulations of magnetohydrodynamic turbulence. Rotating wall modes similar to those in the Rayleigh-Bénard convection with rotation are found in flows near the Chandrasekhar linear stability limit. Detailed analysis of the spatial structure of the flows and its effect of global transport properties is reported.

preprint2013arXiv

Resolving the fine-scale structure in turbulent Rayleigh-Benard convection

We present high-resolution direct numerical simulation studies of turbulent Rayleigh-Benard convection in a closed cylindrical cell with an aspect ratio of one. The focus of our analysis is on the finest scales of convective turbulence, in particular the statistics of the kinetic energy and thermal dissipation rates in the bulk and the whole cell. The fluctuations of the energy dissipation field can directly be translated into a fluctuating local dissipation scale which is found to develop ever finer fluctuations with increasing Rayleigh number. The range of these scales as well as the probability of high-amplitude dissipation events decreases with increasing Prandtl number. In addition, we examine the joint statistics of the two dissipation fields and the consequences of high-amplitude events. We also have investigated the convergence properties of our spectral element method and have found that both dissipation fields are very sensitive to insufficient resolution. We demonstrate that global transport properties, such as the Nusselt number, and the energy balances are partly insensitive to insufficient resolution and yield correct results even when the dissipation fields are under-resolved. Our present numerical framework is also compared with high-resolution simulations which use a finite difference method. For most of the compared quantities the agreement is found to be satisfactory.

preprint2012arXiv

Boundary layer structure in turbulent Rayleigh-Benard convection

The structure of the boundary layers in turbulent Rayleigh-Benard convection is studied by means of three-dimensional direct numerical simulations. We consider convection in a cylindrical cell at an aspect ratio one for Rayleigh numbers of Ra=3e+9 and 3e+10 at fixed Prandtl number Pr=0.7. Similar to the experimental results in the same setup and for the same Prandtl number, the structure of the laminar boundary layers of the velocity and temperature fields is found to deviate from the prediction of the Prandtl-Blasius-Pohlhausen theory. Deviations decrease when a dynamical rescaling of the data with an instantaneously defined boundary layer thickness is performed and the analysis plane is aligned with the instantaneous direction of the large-scale circulation in the closed cell. Our numerical results demonstrate that important assumptions which enter existing classical laminar boundary layer theories for forced and natural convection are violated, such as the strict two-dimensionality of the dynamics or the steadiness of the fluid motion. The boundary layer dynamics consists of two essential local dynamical building blocks, a plume detachment and a post-plume phase. The former is associated with larger variations of the instantaneous thickness of velocity and temperature boundary layer and a fully three-dimensional local flow. The post-plume dynamics is connected with the large-scale circulation in the cell that penetrates the boundary region from above. The mean turbulence profiles taken in localized sections of the boundary layer for both dynamical phases are also compared with solutions of perturbation expansions of the boundary layer equations of forced or natural convection towards mixed convection. Our analysis of both boundary layers shows that the near-wall dynamics combines elements of forced Blasius-type and natural convection.

preprint2012arXiv

Cloud microphysical effects of turbulent mixing and entrainment

Turbulent mixing and entrainment at the boundary of a cloud is studied by means of direct numerical simulations that couple the Eulerian description of the turbulent velocity and water vapor fields with a Lagrangian ensemble of cloud water droplets that can grow and shrink by condensation and evaporation, respectively. The focus is on detailed analysis of the relaxation process of the droplet ensemble during the entrainment of subsaturated air, in particular the dependence on turbulence time scales, droplet number density, initial droplet radius and particle inertia. We find that the droplet evolution during the entrainment process is captured best by a phase relaxation time that is based on the droplet number density with respect to the entire simulation domain and the initial droplet radius. Even under conditions favoring homogeneous mixing, the probability density function of supersaturation at droplet locations exhibits initially strong negative skewness, consistent with droplets near the cloud boundary being suddenly mixed into clear air, but rapidly approaches a narrower, symmetric shape. The droplet size distribution, which is initialized as perfectly monodisperse, broadens and also becomes somewhat negatively skewed. Particle inertia and gravitational settling lead to a more rapid initial evaporation, but ultimately only to slight depletion of both tails of the droplet size distribution. The Reynolds number dependence of the mixing process remained weak over the parameter range studied, most probably due to the fact that the inhomogeneous mixing regime could not be fully accessed when phase relaxation times based on global number density are considered.

preprint2012arXiv

Moist turbulent Rayleigh-Benard convection with Neumann and Dirichlet boundary conditions

Turbulent Rayleigh-Benard convection with phase changes in an extended layer between two parallel impermeable planes is studied by means of three-dimensional direct numerical simulations for Rayleigh numbers between 10^4 and 1.5\times 10^7 and for Prandtl number Pr=0.7. Two different sets of boundary conditions of temperature and total water content are compared: imposed constant amplitudes which translate into Dirichlet boundary conditions for the scalar field fluctuations about the quiescent diffusive equilibrium and constant imposed flux boundary conditions that result in Neumann boundary conditions. Moist turbulent convection is in the conditionally unstable regime throughout this study for which unsaturated air parcels are stably and saturated air parcels unstably stratified. A direct comparison of both sets of boundary conditions with the same parameters requires to start the turbulence simulations out of differently saturated equilibrium states. Similar to dry Rayleigh-Benard convection the differences in the turbulent velocity fluctuations, the cloud cover and the convective buoyancy flux decrease across the layer with increasing Rayleigh number. At the highest Rayleigh numbers the system is found in a two-layer regime, a dry cloudless and stably stratified layer with low turbulence level below a fully saturated and cloudy turbulent one which equals classical Rayleigh-Benard convection layer. Both are separated by a strong inversion that gets increasingly narrower for growing Rayleigh number.

preprint2011arXiv

Low-dimensional model of turbulent Rayleigh-Benard convection in a Cartesian cell with square domain

A low-dimensional model (LDM) for turbulent Rayleigh-Benard convection in a Cartesian cell with square domain, based on the Galerkin projection of the Boussinesq equations onto a finite set of empirical eigenfunctions, is presented. The empirical eigenfunctions are obtained from a joint Proper Orthogonal Decomposition (POD) of the velocity and temperature fields using the Snapshot Method on the basis of a direct numerical simulation (DNS). The resulting LDM is a quadratic inhomogeneous system of coupled ordinary differential equations which we use to describe the long-time temporal evolution of the large-scale mode amplitudes for a Rayleigh number of 1e5 and a Prandtl number of 0.7. The truncation to a finite number of degrees of freedom, that does not exceed a number of 310 for the present, requires the additional implementation of an eddy viscosity-diffusivity to capture the missing dissipation of the small-scale modes. The magnitude of this additional dissipation mechanism is determined by requiring statistical stationarity and a total dissipation that corresponds with the original DNS data. We compare the performance of two models, a constant so-called Heisenberg viscosity--diffusivity and a mode-dependent or modal one. The latter viscosity--diffusivity model turns out to reproduce the large-scale properties of the turbulent convection qualitatively well, even for a model with only a few hundred POD modes.

preprint2011arXiv

Rayleigh-Benard convection with phase changes in a Galerkin model

The transition to turbulence in Rayleigh-Benard convection with phase changes and the resulting convective patterns are studied in a three-dimensional Galerkin model. Our study is focused to the conditionally unstable regime of moist convection in which the stratification is stable for unsaturated air parcels and unstable for saturated parcels. We perform a comprehensive statistical analysis of the transition to convection that samples the dependence of attractors (or fixed points) in the phase space of the model on the dimensionless parameters. Conditionally unstable convection can be initiated either from a fully unsaturated linearly stable equilibrium or a fully saturated linearly unstable equilibrium. Highly localized moist convection can be found in steady state, in oscillating recharge-discharge regime or turbulent in dependence on the aspect ratio and the degree of stable stratification of the unsaturated air. Our phase space analysis predicts parameter ranges for which self-sustained convective regimes in the case of subcritical conditional instability can be observed. The observed regime transitions for moist convection bear some similarities to transitions to turbulence in simple shear flows.

preprint2010arXiv

Aspect ratio dependence of heat transfer and large-scale flow in turbulent convection

The heat transport and corresponding changes in the large-scale circulation (LSC) in turbulent Rayleigh-Bénard convection are studied by means of three-dimensional direct numerical simulations as a function of the aspect ratio $Γ$ of a closed cylindrical cell and the Rayleigh number $Ra$. For small and moderate aspect ratios, the global heat transfer law $Nu=A\times Ra^β$ shows a power law dependence of both fit coefficients $A$ and $β$ on the aspect ratio. A minimum Nusselt number coincides with the point where the LSC undergoes a transition from a single-roll to a double-roll pattern. With increasing aspect ratio, we detect complex multi-roll LSC configurations. The aspect ratio dependence of the turbulent heat transfer for small and moderate $Γ$ is in line with a varying amount of energy contained in the LSC, as quantified by the Proper Orthogonal Decomposition analysis. For $Γ\gtrsim 8$ the heat transfer becomes independent of the aspect ratio.

preprint2010arXiv

Buoyancy statistics in moist turbulent Rayleigh-Benard convection

We study shallow moist Rayleigh-Benard convection in the Boussinesq approximation in three-dimensional direct numerical simulations. The thermodynamics of phase changes is approximated by a piecewise linear equation of state close to the phase boundary. The impact of phase changes on the turbulent fluctuations and the transfer of buoyancy through the layer is discussed as a function of the Rayleigh number and the ability to form liquid water. The enhanced buoyancy flux due to phase changes is compared with dry convection reference cases and related to the cloud cover in the convection layer. This study indicates that the moist Rayleigh-Benard problem offers a practical framework for the development and evaluation of parametrizations for atmospheric convection.

preprint2010arXiv

Lagrangian tracer dynamics in a closed cylindrical turbulent convection cell

Turbulent Rayleigh-Benard convection in a closed cylindrical cell is studied in the Lagrangian frame of reference with the help of three-dimensional direct numerical simulations. The aspect ratio of the cell Gamma is varied between 1 and 12, and the Rayleigh number Ra between 10^7 and 10^9. The Prandtl number Pr is fixed at 0.7. It is found that both the pair dispersion of the Lagrangian tracer particles and the statistics of the acceleration components measured along the particle trajectories depend on the aspect ratio for a fixed Rayleigh number for the parameter range covered in our studies. This suggests that large-scale circulations present in the convection cell affect the Lagrangian dynamics. Our findings are in qualitative agreement with existing Lagrangian laboratory experiments on turbulent convection.

preprint2009arXiv

Extreme vorticity growth in Navier-Stokes turbulence

According to statistical turbulence theory, the ensemble averaged squared vorticity rho_E is expected to grow not faster than drho_E/dt ~ rho_E^{3/2}. Solving a variational problem for maximal bulk enstrophy (E) growth, velocity fields were found for which the growth rate is as large as dE/dt ~ E^3. Using numerical simulations with well resolved small scales and a quasi-Lagrangian advection to track fluid subvolumes with rapidly growing vorticity, we study spatially resolved statistics of vorticity growth. We find that the volume ensemble averaged growth bound is satisfied locally to a remarkable degree of accuracy. Elements with dE/dt ~ E^3 can also be identified, but their growth tends to be replaced by the ensemble-averaged law when the intensities become too large.

preprint2009arXiv

Lagrangian views on turbulent mixing of passive scalars

The Lagrangian view of passive scalar turbulence has recently produced interesting results and interpretations. Innovations in theory, experiments, simulations and data analysis of Lagrangian turbulence are reviewed here in brief. Part of the review is closely related to the so-called Kraichnan model for the advection of the passive scalar in synthetic turbulence. Possible implications for a better understanding of the passive scalar mixing in Navier-Stokes turbulence are also discussed.