skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Characteristic scales during the onset of radiatively driven convection: linear analysis and simulations
This paper considers the initial stage of radiatively driven convection, when the perturbations from a quiescent but time-dependent background state are small. Radiation intensity is assumed to decay exponentially away from the surface, and we consider parameter regimes in which the depth of the water is greater than the decay scale of$$e$$of the radiation intensity. Both time-independent and time-periodic radiation are considered. In both cases, the background temperature profile of the water column is time-dependent. A linear analysis of the system is performed based on these time-dependent profiles. We find that the perturbations grow in time according to$$\exp [(\sigma (t) t)]$$, where$$\sigma (t)$$is a time-dependent growth rate. An appropriately defined Reynolds number is the primary dimensionless number characterising the system, determining the wavelength, vertical structure and growth rate of the perturbations. Simulations using a Boussinesq model (the Stratified Ocean Model with Adaptive Refinement) confirm the linear analysis.  more » « less
Award ID(s):
1829924
PAR ID:
10497613
Author(s) / Creator(s):
;
Publisher / Repository:
Camnbridge University Press
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
973
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Electrophoresis is the motion of a charged colloidal particle in an electrolyte under an applied electric field. The electrophoretic velocity of a spherical particle depends on the dimensionless electric field strength$$\beta =a^*e^*E_\infty ^*/k_B^*T^*$$, defined as the ratio of the product of the applied electric field magnitude$$E_\infty ^*$$and particle radius$$a^*$$, to the thermal voltage$$k_B^*T^*/e^*$$, where$$k_B^*$$is Boltzmann's constant,$$T^*$$is the absolute temperature, and$$e^*$$is the charge on a proton. In this paper, we develop a spectral element algorithm to compute the electrophoretic velocity of a spherical, rigid, dielectric particle, of fixed dimensionless surface charge density$$\sigma$$over a wide range of$$\beta$$. Here,$$\sigma =(e^*a^*/\epsilon ^*k_B^*T^*)\sigma ^*$$, where$$\sigma ^*$$is the dimensional surface charge density, and$$\epsilon ^*$$is the permittivity of the electrolyte. For moderately charged particles ($$\sigma ={O}(1)$$), the electrophoretic velocity is linear in$$\beta$$when$$\beta \ll 1$$, and its dependence on the ratio of the Debye length ($$1/\kappa ^*$$) to particle radius (denoted by$$\delta =1/(\kappa ^*a^*)$$) agrees with Henry's formula. As$$\beta$$increases, the nonlinear contribution to the electrophoretic velocity becomes prominent, and the onset of this behaviour is$$\delta$$-dependent. For$$\beta \gg 1$$, the electrophoretic velocity again becomes linear in field strength, approaching the Hückel limit of electrophoresis in a dielectric medium, for all$$\delta$$. For highly charged particles ($$\sigma \gg 1$$) in the thin-Debye-layer limit ($$\delta \ll 1$$), our computations are in good agreement with recent experimental and asymptotic results. 
    more » « less
  2. Abstract Cohesive powersof computable structures are effective analogs of ultrapowers, where cohesive sets play the role of ultrafilters. Let$$\omega $$,$$\zeta $$, and$$\eta $$denote the respective order-types of the natural numbers, the integers, and the rationals when thought of as linear orders. We investigate the cohesive powers of computable linear orders, with special emphasis on computable copies of$$\omega $$. If$$\mathcal {L}$$is a computable copy of$$\omega $$that is computably isomorphic to the usual presentation of$$\omega $$, then every cohesive power of$$\mathcal {L}$$has order-type$$\omega + \zeta \eta $$. However, there are computable copies of$$\omega $$, necessarily not computably isomorphic to the usual presentation, having cohesive powers not elementarily equivalent to$$\omega + \zeta \eta $$. For example, we show that there is a computable copy of$$\omega $$with a cohesive power of order-type$$\omega + \eta $$. Our most general result is that if$$X \subseteq \mathbb {N} \setminus \{0\}$$is a Boolean combination of$$\Sigma _2$$sets, thought of as a set of finite order-types, then there is a computable copy of$$\omega $$with a cohesive power of order-type$$\omega + \boldsymbol {\sigma }(X \cup \{\omega + \zeta \eta + \omega ^*\})$$, where$$\boldsymbol {\sigma }(X \cup \{\omega + \zeta \eta + \omega ^*\})$$denotes the shuffle of the order-types inXand the order-type$$\omega + \zeta \eta + \omega ^*$$. Furthermore, ifXis finite and non-empty, then there is a computable copy of$$\omega $$with a cohesive power of order-type$$\omega + \boldsymbol {\sigma }(X)$$. 
    more » « less
  3. We study the spaces of twisted conformal blocks attached to a$$\Gamma$$-curve$$\Sigma$$with marked$$\Gamma$$-orbits and an action of$$\Gamma$$on a simple Lie algebra$$\mathfrak {g}$$, where$$\Gamma$$is a finite group. We prove that if$$\Gamma$$stabilizes a Borel subalgebra of$$\mathfrak {g}$$, then the propagation theorem and factorization theorem hold. We endow a flat projective connection on the sheaf of twisted conformal blocks attached to a smooth family of pointed$$\Gamma$$-curves; in particular, it is locally free. We also prove that the sheaf of twisted conformal blocks on the stable compactification of Hurwitz stack is locally free. Let$$\mathscr {G}$$be the parahoric Bruhat–Tits group scheme on the quotient curve$$\Sigma /\Gamma$$obtained via the$$\Gamma$$-invariance of Weil restriction associated to$$\Sigma$$and the simply connected simple algebraic group$$G$$with Lie algebra$$\mathfrak {g}$$. We prove that the space of twisted conformal blocks can be identified with the space of generalized theta functions on the moduli stack of quasi-parabolic$$\mathscr {G}$$-torsors on$$\Sigma /\Gamma$$when the level$$c$$is divisible by$$|\Gamma |$$(establishing a conjecture due to Pappas and Rapoport). 
    more » « less
  4. This paper focuses on a 2D magnetohydrodynamic system with only horizontal dissipation in the domain$$\Omega = \mathbb {T}\times \mathbb {R}$$with$$\mathbb {T}=[0,\,1]$$being a periodic box. The goal here is to understand the stability problem on perturbations near the background magnetic field$$(1,\,0)$$. Due to the lack of vertical dissipation, this stability problem is difficult. This paper solves the desired stability problem by simultaneously exploiting two smoothing and stabilizing mechanisms: the enhanced dissipation due to the coupling between the velocity and the magnetic fields, and the strong Poincaré type inequalities for the oscillation part of the solution, namely the difference between the solution and its horizontal average. In addition, the oscillation part of the solution is shown to converge exponentially to zero in$$H^{1}$$as$$t\to \infty$$. As a consequence, the solution converges to its horizontal average asymptotically. 
    more » « less
  5. Large-eddy simulation was used to model turbulent atmospheric surface layer (ASL) flow over canopies composed of streamwise-aligned rows of synthetic trees of height,$$h$$, and systematically arranged to quantify the response to variable streamwise spacing,$$\delta _1$$, and spanwise spacing,$$\delta _2$$, between adjacent trees. The response to spanwise and streamwise heterogeneity has, indeed, been the topic of a sustained research effort: the former resulting in formation of Reynolds-averaged counter-rotating secondary cells, the latter associated with the$$k$$- and$$d$$-type response. No study has addressed the confluence of both, and results herein show secondary flow polarity reversal across ‘critical’ values of$$\delta _1$$and$$\delta _2$$. For$$\delta _2/\delta \lesssim 1$$and$$\gtrsim 2$$, where$$\delta$$is the flow depth, the counter-rotating secondary cells are aligned such that upwelling and downwelling, respectively, occurs above the elements. The streamwise spacing$$\delta _1$$regulates this transition, with secondary cell reversal occurring first for the largest$$k$$-type cases, as elevated turbulence production within the canopy necessitates entrainment of fluid from aloft. The results are interpreted through the lens of a benchmark prognostic closure for effective aerodynamic roughness,$$z_{0,{Eff.}} = \alpha \sigma _h$$, where$$\alpha$$is a proportionality constant and$$\sigma _h$$is height root mean square. We report$$\alpha \approx 10^{-1}$$, the value reported over many decades for a broad range of rough surfaces, for$$k$$-type cases at small$$\delta _2$$, whereas the transition to$$d$$-type arrangements necessitates larger$$\delta _2$$. Though preliminary, results highlight the non-trivial response to variation of streamwise and spanwise spacing. 
    more » « less