9 results
Exploiting self-organized criticality in strongly stratified turbulence
- Gregory P. Chini, Guillaume Michel, Keith Julien, Cesar B. Rocha, Colm-cille P. Caulfield
-
- Journal:
- Journal of Fluid Mechanics / Volume 933 / 25 February 2022
- Published online by Cambridge University Press:
- 23 December 2021, A22
-
- Article
- Export citation
-
A multiscale reduced description of turbulent free shear flows in the presence of strong stabilizing density stratification is derived via asymptotic analysis of the Boussinesq equations in the simultaneous limits of small Froude and large Reynolds numbers. The analysis explicitly recognizes the occurrence of dynamics on disparate spatiotemporal scales, yielding simplified partial differential equations governing the coupled evolution of slow large-scale hydrostatic flows and fast small-scale isotropic instabilities and internal waves. The dynamics captured by the coupled reduced equations is illustrated in the context of two-dimensional strongly stratified Kolmogorov flow. A noteworthy feature of the reduced model is that the fluctuations are constrained to satisfy quasilinear (QL) dynamics about the comparably slowly varying large-scale fields. Crucially, this QL reduction is not invoked as an ad hoc closure approximation, but rather is derived in a physically relevant and mathematically consistent distinguished limit. Further analysis of the resulting slow–fast QL system shows how the amplitude of the fast stratified-shear instabilities is slaved to the slowly evolving mean fields to ensure the marginal stability of the latter. Physically, this marginal stability condition appears to be compatible with recent evidence of self-organized criticality in both observations and simulations of stratified turbulence. Algorithmically, the slaving of the fluctuation fields enables numerical simulations to be time-evolved strictly on the slow time scale of the hydrostatic flow. The reduced equations thus provide a solid mathematical foundation for future studies of three-dimensional strongly stratified turbulence in extreme parameter regimes of geophysical relevance and suggest avenues for new sub-grid-scale parametrizations.
Steady Rayleigh–Bénard convection between stress-free boundaries
- Baole Wen, David Goluskin, Matthew LeDuc, Gregory P. Chini, Charles R. Doering
-
- Journal:
- Journal of Fluid Mechanics / Volume 905 / 25 December 2020
- Published online by Cambridge University Press:
- 04 November 2020, R4
-
- Article
- Export citation
-
Steady two-dimensional Rayleigh–Bénard convection between stress-free isothermal boundaries is studied via numerical computations. We explore properties of steady convective rolls with aspect ratios ${\rm \pi} /5\leqslant \varGamma \leqslant 4{\rm \pi}$, where $\varGamma$ is the width-to-height ratio for a pair of counter-rotating rolls, over eight orders of magnitude in the Rayleigh number, $10^3\leqslant Ra\leqslant 10^{11}$, and four orders of magnitude in the Prandtl number, $10^{-2}\leqslant Pr\leqslant 10^2$. At large $Ra$ where steady rolls are dynamically unstable, the computed rolls display $Ra \rightarrow \infty$ asymptotic scaling. In this regime, the Nusselt number $Nu$ that measures heat transport scales as $Ra^{1/3}$ uniformly in $Pr$. The prefactor of this scaling depends on $\varGamma$ and is largest at $\varGamma \approx 1.9$. The Reynolds number $Re$ for large-$Ra$ rolls scales as $Pr^{-1} Ra^{2/3}$ with a prefactor that is largest at $\varGamma \approx 4.5$. All of these large-$Ra$ features agree quantitatively with the semi-analytical asymptotic solutions constructed by Chini & Cox (Phys. Fluids, vol. 21, 2009, 083603). Convergence of $Nu$ and $Re$ to their asymptotic scalings occurs more slowly when $Pr$ is larger and when $\varGamma$ is smaller.
A self-sustaining process theory for uniform momentum zones and internal shear layers in high Reynolds number shear flows
- Brandon Montemuro, Christopher M. White, Joseph C. Klewicki, Gregory P. Chini
-
- Journal:
- Journal of Fluid Mechanics / Volume 901 / 25 October 2020
- Published online by Cambridge University Press:
- 02 September 2020, A28
-
- Article
- Export citation
-
Many exact coherent states (ECS) arising in wall-bounded shear flows have an asymptotic structure at extreme Reynolds number $Re$ in which the effective Reynolds number governing the streak and roll dynamics is $\mathit {O}(1)$. Consequently, these viscous ECS are not suitable candidates for quasi-coherent structures away from the wall that necessarily are inviscid in the mean. Specifically, viscous ECS cannot account for the singular nature of the inertial domain, where the flow self-organizes into uniform momentum zones (UMZs) separated by internal shear layers and the instantaneous streamwise velocity develops a staircase-like profile. In this investigation, a large-$Re$ asymptotic analysis is performed to explore the potential for a three-dimensional, short streamwise- and spanwise-wavelength instability of the embedded shear layers to sustain a spatially distributed array of much larger-scale, effectively inviscid streamwise roll motions. In contrast to other self-sustaining process theories, the rolls are sufficiently strong to differentially homogenize the background shear flow, thereby providing a mechanistic explanation for the formation and maintenance of UMZs and interlaced shear layers that respects the leading-order balance structure of the mean dynamics.
Strong wave–mean-flow coupling in baroclinic acoustic streaming
- Guillaume Michel, Gregory P. Chini
-
- Journal:
- Journal of Fluid Mechanics / Volume 858 / 10 January 2019
- Published online by Cambridge University Press:
- 06 November 2018, pp. 536-564
-
- Article
- Export citation
-
The interaction of an acoustic wave with a stratified fluid can drive strong streaming flows owing to the baroclinic production of fluctuating vorticity, as recently demonstrated by Chini et al. (J. Fluid Mech., 744, 2014, pp. 329–351). In the present investigation, a set of wave/mean-flow interaction equations is derived that governs the coupled dynamics of a standing acoustic-wave mode of characteristic (small) amplitude $\unicode[STIX]{x1D716}$ and the streaming flow it drives in a thin channel with walls maintained at differing temperatures. Unlike classical Rayleigh streaming, the resulting mean flow arises at $O(\unicode[STIX]{x1D716})$ rather than at $O(\unicode[STIX]{x1D716}^{2})$. Consequently, fully two-way coupling between the waves and the mean flow is possible: the streaming is sufficiently strong to induce $O(1)$ rearrangements of the imposed background temperature and density fields, which modifies the spatial structure and frequency of the acoustic mode on the streaming time scale. A novel Wentzel–Kramers–Brillouin–Jeffreys analysis is developed to average over the fast wave dynamics, enabling the coupled system to be integrated strictly on the slow time scale of the streaming flow. Analytical solutions of the reduced system are derived for weak wave forcing and are shown to reproduce results from prior direct numerical simulations (DNS) of the compressible Navier–Stokes and heat equations with remarkable accuracy. Moreover, numerical simulations of the reduced system are performed in the regime of strong wave/mean-flow coupling for a fraction of the computational cost of the corresponding DNS. These simulations shed light on the potential for baroclinic acoustic streaming to be used as an effective means to enhance heat transfer.
A uniform momentum zone–vortical fissure model of the turbulent boundary layer
- Juan Carlos Cuevas Bautista, Alireza Ebadi, Christopher M. White, Gregory P. Chini, Joseph C. Klewicki
-
- Journal:
- Journal of Fluid Mechanics / Volume 858 / 10 January 2019
- Published online by Cambridge University Press:
- 06 November 2018, pp. 609-633
-
- Article
- Export citation
-
Recent studies reveal that at large friction Reynolds number $\unicode[STIX]{x1D6FF}^{+}$ the inertially dominated region of the turbulent boundary layer is composed of large-scale zones of nearly uniform momentum segregated by narrow fissures of concentrated vorticity. Experiments show that, when scaled by the boundary-layer thickness, the fissure thickness is $\mathit{O}(1/\sqrt{\unicode[STIX]{x1D6FF}^{+}})$, while the dimensional jump in streamwise velocity across each fissure scales in proportion to the friction velocity $u_{\unicode[STIX]{x1D70F}}$. A simple model that exploits these essential elements of the turbulent boundary-layer structure at large $\unicode[STIX]{x1D6FF}^{+}$ is developed. First, a master wall-normal profile of streamwise velocity is constructed by placing a discrete number of fissures across the boundary layer. The number of fissures and their wall-normal locations follow scalings informed by analysis of the mean momentum equation. The fissures are then randomly displaced in the wall-normal direction, exchanging momentum as they move, to create an instantaneous velocity profile. This process is repeated to generate ensembles of streamwise velocity profiles from which statistical moments are computed. The modelled statistical profiles are shown to agree remarkably well with those acquired from direct numerical simulations of turbulent channel flow at large $\unicode[STIX]{x1D6FF}^{+}$. In particular, the model robustly reproduces the empirically observed sub-Gaussian behaviour for the skewness and kurtosis profiles over a large range of input parameters.
The influence of a poroelastic till on rapid subglacial flooding and cavity formation
- Duncan R. Hewitt, Gregory P. Chini, Jerome A. Neufeld
-
- Journal:
- Journal of Fluid Mechanics / Volume 855 / 25 November 2018
- Published online by Cambridge University Press:
- 24 September 2018, pp. 1170-1207
-
- Article
- Export citation
-
We develop a model of the rapid propagation of water at the contact between elastic glacial ice and a poroelastic subglacial till, motivated by observations of the rapid drainage of supraglacial lakes in Greenland. By treating the ice as an elastic bending beam, the fluid dynamics of contact with the subglacial hydrological network, which is modelled as a saturated poroelastic till, can be examined in detail. The model describes the formation and dynamics of an axisymmetric subglacial cavity, and the spread of pore pressure, in response to injection of fluid. A combination of numerical simulation and asymptotic analysis is used to describe these dynamics for both a rigid and a deformable porous till, and for both laminar and turbulent fluid flow. For constant injection rates and laminar flow, the cavity is isostatic and its spread is controlled by bending of the ice and suction of pore water in the vicinity of the ice–till contact. For a deformable till, this control can be modified: generically, a flexural wave that is initially trapped in advance of the contact point relaxes over time by diffusion of pore pressure ahead of the cavity. While the dynamics are found to be relatively insensitive to the properties of the subglacial till during injection with a constant flux, significant dependence on the till properties is manifest during the subsequent spread of a constant volume. A simple hybrid turbulent–laminar model is presented to account for fast injection rates of water: in this case, self-similar turbulent propagation can initially control the spread of the cavity, but there is a transition to laminar control in the vicinity of the ice–till contact point as the flow slows. Finally, the model results are compared with recent geophysical observations of the rapid drainage of supraglacial lakes in Greenland; the comparison provides qualitative agreement and raises suggestions for future quantitative comparison.
Inclined porous medium convection at large Rayleigh number
- Baole Wen, Gregory P. Chini
-
- Journal:
- Journal of Fluid Mechanics / Volume 837 / 25 February 2018
- Published online by Cambridge University Press:
- 05 January 2018, pp. 670-702
-
- Article
- Export citation
-
High-Rayleigh-number ($Ra$) convection in an inclined two-dimensional porous layer is investigated using direct numerical simulations (DNS) and stability and variational upper-bound analyses. When the inclination angle $\unicode[STIX]{x1D719}$ of the layer satisfies $0^{\circ }<\unicode[STIX]{x1D719}\lesssim 25^{\circ }$, DNS confirm that the flow exhibits a three-region wall-normal asymptotic structure in accord with the strictly horizontal ($\unicode[STIX]{x1D719}=0^{\circ }$) case, except that as $\unicode[STIX]{x1D719}$ is increased the time-mean spacing between neighbouring interior plumes also increases substantially. Both DNS and upper-bound analysis indicate that the heat transport enhancement factor (i.e. the Nusselt number) $Nu\sim CRa$ with a $\unicode[STIX]{x1D719}$-dependent prefactor $C$. When $\unicode[STIX]{x1D719}>\unicode[STIX]{x1D719}_{t}$, however, where $30^{\circ }<\unicode[STIX]{x1D719}_{t}<32^{\circ }$ independently of $Ra$, the columnar flow structure is completely broken down: the flow transitions to a large-scale travelling-wave convective roll state, and the heat transport is significantly reduced. To better understand the physics of inclined porous medium convection at large $Ra$ and modest inclination angles, a spatial Floquet analysis is performed, yielding predictions of the linear stability of numerically computed, fully nonlinear steady convective states. The results show that there exist two types of instability when $\unicode[STIX]{x1D719}\neq 0^{\circ }$: a bulk-mode instability and a wall-mode instability, consistent with previous findings for $\unicode[STIX]{x1D719}=0^{\circ }$ (Wen et al., J. Fluid Mech., vol. 772, 2015, pp. 197–224). The background flow induced by the inclination of the layer intensifies the bulk-mode instability during its subsequent nonlinear evolution, thereby favouring increased spacing between the interior plumes relative to that observed in convection in a horizontal porous layer.
Structure and stability of steady porous medium convection at large Rayleigh number
- Baole Wen, Lindsey T. Corson, Gregory P. Chini
-
- Journal:
- Journal of Fluid Mechanics / Volume 772 / 10 June 2015
- Published online by Cambridge University Press:
- 05 May 2015, pp. 197-224
-
- Article
- Export citation
-
A systematic investigation of unstable steady-state solutions of the Darcy–Oberbeck–Boussinesq equations at large values of the Rayleigh number $\mathit{Ra}$ is performed to gain insight into two-dimensional porous medium convection in domains of varying aspect ratio $L$. The steady convective states are shown to transport less heat than the statistically steady ‘turbulent’ flow realised at the same parameter values: the Nusselt number $\mathit{Nu}\sim \mathit{Ra}$ for turbulent porous medium convection, while $\mathit{Nu}\sim \mathit{Ra}^{0.6}$ for the maximum heat-transporting steady solutions. A key finding is that the lateral scale of the heat-flux-maximising solutions shrinks roughly as $L\sim \mathit{Ra}^{-0.5}$, reminiscent of the decrease of the mean inter-plume spacing observed in turbulent porous medium convection as the thermal forcing is increased. A spatial Floquet analysis is performed to investigate the linear stability of the fully nonlinear steady convective states, extending a recent study by Hewitt et al. (J. Fluid Mech., vol. 737, 2013, pp. 205–231) by treating a base convective state, and secondary stability modes, that satisfy appropriate boundary conditions along plane parallel walls. As in that study, a bulk instability mode is found for sufficiently small-aspect-ratio base states. However, the growth rate of this bulk mode is shown to be significantly reduced by the presence of the walls. Beyond a certain critical $\mathit{Ra}$-dependent aspect ratio, the base state is most strongly unstable to a secondary mode that is localised near the heated and cooled walls. Direct numerical simulations, strategically initialised to investigate the fully nonlinear evolution of the most dangerous secondary instability modes, suggest that the (long time) mean inter-plume spacing in statistically steady porous medium convection results from a balance between the competing effects of these two types of instability.
Wall to wall optimal transport
- Pedram Hassanzadeh, Gregory P. Chini, Charles R. Doering
-
- Journal:
- Journal of Fluid Mechanics / Volume 751 / 25 July 2014
- Published online by Cambridge University Press:
- 24 June 2014, pp. 627-662
-
- Article
- Export citation
-
The calculus of variations is employed to find steady divergence-free velocity fields that maximize transport of a tracer between two parallel walls held at fixed concentration for one of two constraints on flow strength: a fixed value of the kinetic energy (mean square velocity) or a fixed value of the enstrophy (mean square vorticity). The optimizing flows consist of an array of (convection) cells of a particular aspect ratio $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\varGamma $. We solve the nonlinear Euler–Lagrange equations analytically for weak flows and numerically – as well as via matched asymptotic analysis in the fixed energy case – for strong flows. We report the results in terms of the Nusselt number ${\mathit{Nu}}$, a dimensionless measure of the tracer transport, as a function of the Péclet number ${\mathit{Pe}}$, a dimensionless measure of the strength of the flow. For both constraints, the maximum transport ${\mathit{Nu}}_{\mathit{MAX}}({\mathit{Pe}})$ is realized in cells of decreasing aspect ratio $\varGamma _{\mathit{opt}}({\mathit{Pe}})$ as ${\mathit{Pe}}$ increases. For the fixed energy problem, ${\mathit{Nu}}_{\mathit{MAX}} \sim {\mathit{Pe}}$ and $\varGamma _{\mathit{opt}} \sim {\mathit{Pe}}^{-1/2}$, while for the fixed enstrophy scenario, ${\mathit{Nu}}_{\mathit{MAX}} \sim {\mathit{Pe}}^{10/17}$ and $\varGamma _{\mathit{opt}} \sim {\mathit{Pe}}^{-0.36}$. We interpret our results in the context of buoyancy-driven Rayleigh–Bénard convection problems that satisfy the flow intensity constraints, enabling us to investigate how the transport scalings compare with upper bounds on ${\mathit{Nu}}$ expressed as a function of the Rayleigh number ${\mathit{Ra}}$. For steady convection in porous media, corresponding to the fixed energy problem, we find ${\mathit{Nu}}_{\mathit{MAX}} \sim {\mathit{Ra}}$ and $\varGamma _{\mathit{opt}} \sim {\mathit{Ra}}^{-1/2}$, while for steady convection in a pure fluid layer between stress-free isothermal walls, corresponding to fixed enstrophy transport, ${\mathit{Nu}}_{\mathit{MAX}} \sim {\mathit{Ra}}^{5/12}$ and $\varGamma _{\mathit{opt}} \sim {\mathit{Ra}}^{-1/4}$.