Abstract Heavy particles sink straight in water, while buoyant bubbles and spheres may zigzag or spiral as they rise. The precise conditions that trigger such path-instabilities are still not completely understood. For a buoyant rising sphere, two parameters are believed to govern the development of unsteady dynamics: the particle’s density relative to the fluid, and its Galileo number. Consequently, with these parameters specified, the opportunities for variation in particle dynamics appear limited. In contrast to this picture, here we demonstrate that vigorous path-oscillations can be triggered by modulating a spherical particle’s moment of inertia (MoI). For a buoyant sphere rising in a turbulent flow, MoI reduction triggers a tumble–flutter transition, while in quiescent liquid, it induces a modification of the sphere wake resulting in large-amplitude path-oscillations. The present finding opens the door for control of particle path- and wake-instabilities, with potential for enhanced mixing and heat transfer in particle-laden and dispersed multiphase environments.
The growth and evolution of microbial populations is often subjected to advection by fluid flows in spatially extended environments, with immediate consequences for questions of spatial population genetics in marine ecology, planktonic diversity and origin of life scenarios. Here, we review recent progress made in understanding this rich problem in the simplified setting of two competing genetic microbial strains subjected to fluid flows. As a pedagogical example we focus on antagonsim, i.e., two killer microorganism strains, each secreting toxins that impede the growth of their competitors (competitive exclusion), in the presence of stationary fluid flows. By solving two coupled reaction-diffusion equations that include advection by simple steady cellular flows composed of characteristic flow motifs in two dimensions (2d), we show how local flow shear and compressibility effects can interact with selective advantage to have a dramatic influence on genetic competition and fixation in spatially distributed populations. We analyze several 1d and 2d flow geometries including sources, sinks, vortices and saddles, and show how simple analytical models of the dynamics of the genetic interface can be used to shed light on the nucleation, coexistence and flow-driven instabilities of genetic drops. By exploiting an analogy with phase separation with nonconserved order parameters, we uncover how these genetic drops harness fluid flows for novel evolutionary strategies, even in the presence of number fluctuations, as confirmed by agent-based simulations as well.
A numerical study of decaying wall-bounded turbulence is made by stopping the forcing in a Taylor-Couette system. The decay progressively exhausts all modes through which perturbations grow, from most unstable to least unstable until a final stage of self-similar viscous decay is reached.
In this combined experimental and numerical study on thermally driven turbulence in a rectangular cell, the global heat transport and the coherent flow structures are controlled with an asymmetric ratchetlike roughness on the top and bottom plates. We show that, by means of symmetry breaking due to the presence of the ratchet structures on the conducting plates, the orientation of the large scale circulation roll (LSCR) can be locked to a preferred direction even when the cell is perfectly leveled out. By introducing a small tilt to the system, we show that the LSCR orientation can be tuned and controlled. The two different orientations of LSCR give two quite different heat transport efficiencies, indicating that heat transport is sensitive to the LSCR direction over the asymmetric roughness structure. Through a quantitative analysis of the dynamics of thermal plume emissions and the orientation of the LSCR over the asymmetric structure, we provide a physical explanation for these findings. The current work has important implications for passive and active flow control in engineering, biofluid dynamics, and geophysical flows.
In the generalized quasilinear approximation (GQLA) (Marston et al. , Phys. Rev. Lett. , vol. 116, 2016, 214501), a threshold wavenumber ( $k_0$ ) in the direction of translational symmetry segregates the total into large- ( $l$ ) and small-scale ( $h$ ) fields. While the governing equation for the large-scale field is fully nonlinear, that for the small scales is linearized with respect to the large-scale field. In addition, some nonlinear triad interactions are omitted in the GQLA. Herein, the GQLA is applied to two-dimensional planar Rayleigh–Bénard convection (RBC). A scale separation between the large-scale convection rolls and small-scale turbulent fluctuations is typical in RBC. The present work explores the efficacy of GQLA in capturing the scale-by-scale energy transfer processes in RBC. The initial condition for the GQLA simulations was either the statistically stationary state obtained in direct numerical simulation (DNS) or random fluctuations superimposed on the linear conductive temperature profile and $\boldsymbol {u}=0$ . The GQLA simulations can capture the convection rolls for $k_0$ larger than or equal to the dominant wavenumber for thermal driving of the flow ( $k_0 \ge k_{\hat {Q}}$ ). Additionally, the GQLA emulates the fully nonlinear dynamics for $k_0$ larger than or equal to the first harmonic of the convection-roll wavenumber ( $k_0 \ge 2k_{roll}$ ). In the intermediate regime with $2k_{roll} > k_0 \ge k_{\hat {Q}}$ , the dynamics captured in GQLA simulations is different from the DNS. In DNS, two primary energy transfer processes dominate: (i) the energy transfer to/from the convection rolls and (ii) the scale-by-scale inverse kinetic energy and forward thermal energy cascades mediated by the convection rolls. The fully nonlinear dynamics is emulated by GQLA when these energy transfer processes are faithfully reproduced. Utilizing the framework of altering the triad interactions in GQLA, an additional intrusive calculation, including target triad interactions, is performed here to study their influence. This intrusive calculation shows that the convection rolls are not captured in GQLA for $k_0 < k_{\hat {Q}}$ because of the exclusion of the $h \to l \to l$ and $l \to h \to h$ triad interactions in GQLA. The inclusion of these triad interactions in the intrusive calculation yields the convection rolls, and the reproduced dynamics is similar to that of the intermediate GQLA regime with $2k_{roll} > k_0 \ge k_{\hat {Q}}$ .
We report on a three-dimensional direct numerical simulation study of flow structure and heat transport in the annular centrifugal Rayleigh–Bénard convection (ACRBC) system, with cold inner and hot outer cylinders corotating axially, for the Rayleigh number range $Ra \in [{10^6},{10^8}]$ and radius ratio range $\eta = {R_i}/{R_o} \in [0.3,0.9]$ ( $R_i$ and $R_o$ are the radius of the inner and outer cylinders, respectively). This study focuses on the dependence of flow dynamics, heat transport and asymmetric mean temperature fields on the radius ratio $\eta$ . For the inverse Rossby number $Ro^{-1} = 1$ , as the Coriolis force balances inertial force, the flow is in the inertial regime. The mechanisms of zonal flow revolving in the prograde direction in this regime are attributed to the asymmetric movements of plumes and the different curvatures of the cylinders. The number of roll pairs is smaller than the circular roll hypothesis as the convection rolls are probably elongated by zonal flow. The physical mechanism of zonal flow is verified by the dependence of the drift frequency of the large-scale circulation (LSC) rolls and the space- and time-averaged azimuthal velocity on $\eta$ . The larger $\eta$ is, the weaker the zonal flow becomes. We show that the heat transport efficiency increases with $\eta$ . It is also found that the bulk temperature deviates from the arithmetic mean temperature and the deviation increases as $\eta$ decreases. This effect can be explained by a simple model that accounts for the curvature effects and the radially dependent centrifugal force in ACRBC.
Abstract Solar flares can release coronal magnetic energy explosively and may impact the safety of near-Earth space environments. Their structures and properties on the macroscale have been interpreted successfully by the generally accepted 2D standard model, invoking magnetic reconnection theory as the key energy conversion mechanism. Nevertheless, some momentous dynamical features as discovered by recent high-resolution observations remain elusive. Here, we report a self-consistent high-resolution 3D magnetohydrodynamical simulation of turbulent magnetic reconnection within a flare current sheet. It is found that fragmented current patches of different scales are spontaneously generated with a well-developed turbulence spectrum at the current sheet, as well as at the flare loop-top region. The close coupling of tearing mode and Kelvin–Helmholtz instabilities plays a critical role in developing turbulent reconnection and in forming dynamical structures with synthetic observables in good agreement with realistic observations. The sophisticated modeling makes a paradigm shift from the traditional to a 3D turbulent reconnection model unifying flare dynamical structures of different scales.
High-Rayleigh number convective turbulence is ubiquitous in many natural phenomena and in industries, such as atmospheric circulations, oceanic flows, flows in the fluid core of planets, and energy generations. In this work, we present a novel approach to boost the Rayleigh number in thermal convection by exploiting centrifugal acceleration and rapidly rotating a cylindrical annulus to reach an effective gravity of 60 times Earth's gravity. We show that in the regime where the Coriolis effect is strong, the scaling exponent of Nusselt number versus Rayleigh number exceeds one-third once the Rayleigh number is large enough. The convective rolls revolve in prograde direction, signifying the emergence of zonal flow. The present findings open a new avenue on the exploration of high-Rayleigh number turbulent thermal convection and will improve the understanding of the flow dynamics and heat transfer processes in geophysical and astrophysical flows and other strongly rotating systems.
In turbulent Rayleigh–Bénard (RB) convection, a transition to the so-called ultimate regime, in which the boundary layers (BL) are of turbulent type, has been postulated. Indeed, at very large Rayleigh number $Ra\approx 10^{13}{-}10^{14}$ a transition in the scaling of the global Nusselt number $Nu$ (the dimensionless heat transfer) and the Reynolds number with $Ra$ has been observed in experiments and very recently in direct numerical simulations (DNS) of two-dimensional (2D) RB convection. In this paper, we analyse the local scaling properties of the lateral temperature structure functions in the BLs of this simulation of 2D RB convection, employing extended self-similarity (ESS) (i.e., plotting the structure functions against each other, rather than only against the scale) in the spirit of the attached-eddy hypothesis, as we have recently introduced for velocity structure functions in wall turbulence (Krug et al. , J. Fluid Mech. , vol. 830, 2017, pp. 797–819). We find no ESS scaling at $Ra$ below the transition and in the near-wall region. However, beyond the transition and for large enough wall distance $z^{+}>100$ , we find clear ESS behaviour, as expected for a scalar in a turbulent boundary layer. In striking correspondence to the $Nu$ scaling, the ESS scaling region is negligible at $Ra=10^{11}$ and well developed at $Ra=10^{14}$ , thus providing strong evidence that the observed transition in the global Nusselt number at $Ra\approx 10^{13}$ indeed is the transition from a laminar type BL to a turbulent type BL. Our results further show that the relative slopes for scalar structure functions in the ESS scaling regime are the same as for their velocity counterparts, extending their previously established universality. The findings are confirmed by comparing to scalar structure functions in three-dimensional turbulent channel flow.
'%3e%3ccircle r='20' fill='%23008'/%3e%3ccircle r='17.5' fill='%23fff'/%3e%3ccircle r='3.5' fill='%23008'/%3e%3cg id='d'%3e%3cg id='c'%3e%3cg id='b'%3e%3cg id='a' fill='%23008'%3e%3ccircle r='.875' transform='rotate(7.5 -8.75 133.5)'/%3e%3cpath d='M0 17.5.6 7 0 2l-.6 5L0 17.5z'/%3e%3c/g%3e%3cuse xlink:href='%23a' transform='rotate(15)'/%3e%3c/g%3e%3cuse xlink:href='%23b' transform='rotate(30)'/%3e%3c/g%3e%3cuse xlink:href='%23c' transform='rotate(60)'/%3e%3c/g%3e%3cuse xlink:href='%23d' transform='rotate(120)'/%3e%3cuse xlink:href='%23d' transform='rotate(-120)'/%3e%3c/g%3e%3c/svg%3e)
'/%3e%3cuse xlink:href='%23a' transform='rotate(23.036 .093 25.536)'/%3e%3cuse xlink:href='%23a' transform='rotate(45.87 1.273 16.18)'/%3e%3cuse xlink:href='%23a' transform='rotate(69.945 .996 12.078)'/%3e%3cuse xlink:href='%23a' transform='rotate(20.66 -19.689 31.932)'/%3e%3c/svg%3e)
