Researcher profile

Olga Shishkina

Olga Shishkina contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
8works
0followers
1topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

8 published item(s)

preprint2022arXiv

Data-driven identification of the spatio-temporal structure of turbulent flows by streaming Dynamic Mode Decomposition

Streaming Dynamic Mode Decomposition (sDMD) (Hemati et al., Phys. Fluids 26(2014)) is a low-storage version of Dynamic Mode Decomposition (DMD) (Schmid, J. Fluid Mech. 656 (2010)), a data-driven method to extract spatio-temporal flow patterns. Streaming DMD avoids storing the entire data sequence in memory by approximating the dynamic modes through incremental updates with new available data. In this paper, we use sDMD to identify and extract dominant spatio-temporal structures of different turbulent flows, requiring the analysis of large datasets. First, the efficiency and accuracy of sDMD are compared to the classical DMD, using a publicly available test dataset that consists of velocity field snapshots obtained by direct numerical simulation of a wake flow behind a cylinder. Streaming DMD not only reliably reproduces the most important dynamical features of the flow; our calculations also highlight its advantage in terms of the required computational resources. We subsequently use sDMD to analyse three different turbulent flows that all show some degree of large-scale coherence: rapidly rotating Rayleigh--Bénard convection, horizontal convection and the asymptotic suction boundary layer. Structures of different frequencies and spatial extent can be clearly separated, and the prominent features of the dynamics are captured with just a few dynamic modes. In summary, we demonstrate that sDMD is a powerful tool for the identification of spatio-temporal structures in a wide range of turbulent flows.

preprint2020arXiv

Elliptical Instability and Multi-Roll Flow Modes of the Large-scale Circulation in Confined Turbulent Rayleigh-Bénard Convection

Turbulent Rayleigh-Bénard convection in slender cylindrical cells exhibits rich dynamics of the large-scale circulation (LSC), with several rolls stacked on top of each other. We propose that the elliptical instability is the mechanism which causes the twisting and breaking of the LSC into multiple rolls and that the volume-averaged heat and momentum transport, represented by the Nusselt number and Reynolds number, is generally weaker for larger number $n$ of the LSC rolls. This is supported by direct numerical simulations for $Ra=5\times10^5$, $Pr=0.1$, $H=5D$ and $1\leq~n\leq4$.

preprint2020arXiv

Multiple states in turbulent large-aspect ratio thermal convection: What determines the number of convection rolls?

Recent findings suggest that wall-bounded turbulent flow can take different statistically stationary turbulent states, with different transport properties, even for the very same values of the control parameters. What state the system takes depends on the initial conditions. Here we analyze the multiple states in large-aspect ratio ($Γ$) two-dimensional turbulent Rayleigh--Bénard flow with no-slip plates and horizontally periodic boundary conditions as model system. We determine the number $n$ of convection rolls, their mean aspect ratios $Γ_r = Γ/n$, and the corresponding transport properties of the flow (i.e., the Nusselt number $Nu$), as function of the control parameters Rayleigh ($Ra$) and Prandtl number. The effective scaling exponent $β$ in $Nu \sim Ra^β$ is found to depend on the realized state and thus $Γ_r$, with a larger value for the smaller $Γ_r$. By making use of a generalized Friedrichs inequality, we show that the elliptical instability and viscous damping determine the $Γ_r$-window for the realizable turbulent states. The theoretical results are in excellent agreement with our numerical finding $2/3 \le Γ_r \le 4/3$, where the lower threshold is approached for the larger $Ra$. Finally, we show that the theoretical approach to frame $Γ_r$ also works for free-slip boundary conditions.

preprint2020arXiv

Scaling in internally heated convection: a unifying theory

We offer a unifying theory for turbulent purely internally heated convection, generalizing the unifying theories of Grossmann and Lohse (2000, 2001) for Rayleigh--Bénard turbulence and of Shishkina, Grossmann and Lohse (2016) for turbulent horizontal convection, which are both based on the splitting of the kinetic and thermal dissipation rates in respective boundary and bulk contributions. We obtain the mean temperature of the system and the Reynolds number (which are the response parameters) as function of the control parameters, namely the internal thermal driving strength (called, when nondimensionalized, the Rayleigh--Roberts number) and the Prandtl number. The results of the theory are consistent with our direct numerical simulations.

preprint2019arXiv

Boundary Zonal Flow in Rotating Turbulent Rayleigh-Bénard Convection

For rapidly rotating turbulent Rayleigh--Bénard convection in a slender cylindrical cell, experiments and direct numerical simulations reveal a boundary zonal flow (BZF) that replaces the classical large-scale circulation. The BZF is located near the vertical side wall and enables enhanced heat transport there. Although the azimuthal velocity of the BZF is cyclonic (in the rotating frame), the temperature is an anticyclonic traveling wave of mode one whose signature is a bimodal temperature distribution near the radial boundary. The BZF width is found to scale like $Ra^{1/4}Ek^{2/3}$ where the Ekman number $Ek$ decreases with increasing rotation rate.

preprint2015arXiv

Thermal Boundary Layer Equation for Turbulent Rayleigh-Benard Convection

We report a new thermal boundary layer equation for turbulent Rayleigh-Benard convection for Prandtl number Pr>1 that takes into account the effect of turbulent fluctuations. These fluctuations are neglected in existing equations, which are based on steady-state and laminar assumptions. Using this new equation, we derive analytically the mean temperature profiles in two limits: (a) Pr>1, Pr~1 and (b) Pr>>1. These two theoretical predictions are in excellent agreement with the results of our direct numerical simulations for Pr=4.38 (water) and Pr=2547.9 (glycerol) respectively.

preprint2014arXiv

Toroidal and poloidal energy in rotating Rayleigh-Bénard convection

We consider rotating Rayleigh-Bénard convection of a fluid with a Prandtl number of $Pr = 0.8$ in a cylindrical cell with an aspect ratio $Γ= 1/2$. Direct numerical simulations were performed for the Rayleigh number range $10^5 \leq Ra \leq 10^9$ and the inverse Rossby number range $0 \leq 1/Ro \leq 20$. We propose a method to capture regime transitions based on the decomposition of the velocity field into toroidal and poloidal parts. We identify four different regimes. First, a buoyancy dominated regime occurring as long as the toroidal energy $e_{tor}$ is not affected by rotation and remains equal to that in the non-rotating case, $e^0_{tor}$. Second, a rotation influenced regime, starting at rotation rates where $e_{tor} > e^0_{tor}$ and ending at a critical inverse Rossby number $1/Ro_{cr}$ that is determined by the balance of the toroidal and poloidal energy, $e_{tor} = e_{pol}$. Third, a rotation dominated regime, where the toroidal energy $e_{tor}$ is larger than both, $e_{pol}$ and $e^0_{tor}$. Fourth, a geostrophic turbulence regime for high rotation rates where the toroidal energy drops below the value of non-rotating convection.

preprint2011arXiv

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

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