skip to main content


This content will become publicly available on June 10, 2024

Title: The Matsuno–Gill model on the sphere

We extend the Matsuno–Gill model, originally developed on the equatorial$\beta$-plane, to the surface of the sphere. While on the$\beta$-plane the non-dimensional model contains a single parameter, the damping rate$\gamma$, on a sphere the model contains a second parameter, the rotation rate$\epsilon ^{1/2}$(Lamb number). By considering the different combinations of damping and rotation, we are able to characterize the solutions over the$(\gamma, \epsilon ^{1/2})$plane. We find that the$\beta$-plane approximation is accurate only for fast rotation rates, where gravity waves traverse a fraction of the sphere's diameter in one rotation period. The particular solutions studied by Matsuno and Gill are accurate only for fast rotation and moderate damping rates, where the relaxation time is comparable to the time on which gravity waves traverse the sphere's diameter. Other regions of the parameter space can be described by different approximations, including radiative relaxation, geostrophic, weak temperature gradient and non-rotating approximations. The effect of the additional parameter introduced by the sphere is to alter the eigenmodes of the free system. Thus, unlike the solutions obtained by Matsuno and Gill, where the long-term response to a symmetric forcing consists solely of Kelvin and Rossby waves, the response on the sphere includes other waves as well, depending on the combination of$\gamma$and$\epsilon ^{1/2}$. The particular solutions studied by Matsuno and Gill apply to Earth's oceans, while the more general$\beta$-plane solutions are only somewhat relevant to Earth's troposphere. In Earth's stratosphere, Venus and Titan, only the spherical solutions apply.

 
more » « less
Award ID(s):
2004572
NSF-PAR ID:
10468677
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Cambridge University Press
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
964
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

    This paper will study almost everywhere behaviors of functions on partition spaces of cardinals possessing suitable partition properties. Almost everywhere continuity and monotonicity properties for functions on partition spaces will be established. These results will be applied to distinguish the cardinality of certain subsets of the power set of partition cardinals.

    The following summarizes the main results proved under suitable partition hypotheses.

    If$\kappa $is a cardinal,$\epsilon < \kappa $,${\mathrm {cof}}(\epsilon ) = \omega $,$\kappa \rightarrow _* (\kappa )^{\epsilon \cdot \epsilon }_2$and$\Phi : [\kappa ]^\epsilon _* \rightarrow \mathrm {ON}$, then$\Phi $satisfies the almost everywhere short length continuity property: There is a club$C \subseteq \kappa $and a$\delta < \epsilon $so that for all$f,g \in [C]^\epsilon _*$, if$f \upharpoonright \delta = g \upharpoonright \delta $and$\sup (f) = \sup (g)$, then$\Phi (f) = \Phi (g)$.

    If$\kappa $is a cardinal,$\epsilon $is countable,$\kappa \rightarrow _* (\kappa )^{\epsilon \cdot \epsilon }_2$holds and$\Phi : [\kappa ]^\epsilon _* \rightarrow \mathrm {ON}$, then$\Phi $satisfies the strong almost everywhere short length continuity property: There is a club$C \subseteq \kappa $and finitely many ordinals$\delta _0, ..., \delta _k \leq \epsilon $so that for all$f,g \in [C]^\epsilon _*$, if for all$0 \leq i \leq k$,$\sup (f \upharpoonright \delta _i) = \sup (g \upharpoonright \delta _i)$, then$\Phi (f) = \Phi (g)$.

    If$\kappa $satisfies$\kappa \rightarrow _* (\kappa )^\kappa _2$,$\epsilon \leq \kappa $and$\Phi : [\kappa ]^\epsilon _* \rightarrow \mathrm {ON}$, then$\Phi $satisfies the almost everywhere monotonicity property: There is a club$C \subseteq \kappa $so that for all$f,g \in [C]^\epsilon _*$, if for all$\alpha < \epsilon $,$f(\alpha ) \leq g(\alpha )$, then$\Phi (f) \leq \Phi (g)$.

    Suppose dependent choice ($\mathsf {DC}$),${\omega _1} \rightarrow _* ({\omega _1})^{\omega _1}_2$and the almost everywhere short length club uniformization principle for${\omega _1}$hold. Then every function$\Phi : [{\omega _1}]^{\omega _1}_* \rightarrow {\omega _1}$satisfies a finite continuity property with respect to closure points: Let$\mathfrak {C}_f$be the club of$\alpha < {\omega _1}$so that$\sup (f \upharpoonright \alpha ) = \alpha $. There is a club$C \subseteq {\omega _1}$and finitely many functions$\Upsilon _0, ..., \Upsilon _{n - 1} : [C]^{\omega _1}_* \rightarrow {\omega _1}$so that for all$f \in [C]^{\omega _1}_*$, for all$g \in [C]^{\omega _1}_*$, if$\mathfrak {C}_g = \mathfrak {C}_f$and for all$i < n$,$\sup (g \upharpoonright \Upsilon _i(f)) = \sup (f \upharpoonright \Upsilon _i(f))$, then$\Phi (g) = \Phi (f)$.

    Suppose$\kappa $satisfies$\kappa \rightarrow _* (\kappa )^\epsilon _2$for all$\epsilon < \kappa $. For all$\chi < \kappa $,$[\kappa ]^{<\kappa }$does not inject into${}^\chi \mathrm {ON}$, the class of$\chi $-length sequences of ordinals, and therefore,$|[\kappa ]^\chi | < |[\kappa ]^{<\kappa }|$. As a consequence, under the axiom of determinacy$(\mathsf {AD})$, these two cardinality results hold when$\kappa $is one of the following weak or strong partition cardinals of determinacy:${\omega _1}$,$\omega _2$,$\boldsymbol {\delta }_n^1$(for all$1 \leq n < \omega $) and$\boldsymbol {\delta }^2_1$(assuming in addition$\mathsf {DC}_{\mathbb {R}}$).

     
    more » « less
  3. Except in the trivial case of spatially uniform flow, the advection–diffusion operator of a passive scalar tracer is linear and non-self-adjoint. In this study, we exploit the linearity of the governing equation and present an analytical eigenfunction approach for computing solutions to the advection–diffusion equation in two dimensions given arbitrary initial conditions, and when the advecting flow field at any given time is a plane parallel shear flow. Our analysis illuminates the specific role that the non-self-adjointness of the linear operator plays in the solution behaviour, and highlights the multiscale nature of the scalar mixing problem given the explicit dependence of the eigenvalue–eigenfunction pairs on a multiscale parameter$q=2{\rm i}k\,{\textit {Pe}}$, where$k$is the non-dimensional wavenumber of the tracer in the streamwise direction, and${\textit {Pe}}$is the Péclet number. We complement our theoretical discussion on the spectra of the operator by computing solutions and analysing the effect of shear flow width on the scale-dependent scalar decay of tracer variance, and characterize the distinct self-similar dispersive processes that arise from the shear flow dispersion of an arbitrarily compact tracer concentration. Finally, we discuss limitations of the present approach and future directions.

     
    more » « less
  4. Abstract

    We study a class of ordinary differential equations with a non-Lipschitz point singularity that admits non-unique solutions through this point. As a selection criterion, we introduce stochastic regularizations depending on a parameter$\nu $: the regularized dynamics is globally defined for each$\nu> 0$, and the original singular system is recovered in the limit of vanishing$\nu $. We prove that this limit yields aunique statistical solutionindependent of regularization when the deterministic system possesses a chaotic attractor having a physical measure with the convergence to equilibrium property. In this case, solutions become spontaneously stochastic after passing through the singularity: they are selected randomly with an intrinsic probability distribution.

     
    more » « less
  5. The effects of surfactants on a mechanically generated plunging breaker are studied experimentally in a laboratory wave tank. Waves are generated using a dispersively focused wave packet with a characteristic wavelength of$\lambda _0 = 1.18$m. Experiments are performed with two sets of surfactant solutions. In the first set, increasing amounts of the soluble surfactant Triton X-100 are mixed into the tank water, while in the second set filtered tap water is left undisturbed in the tank for wait times ranging from 15 min to 21 h. Increasing Triton X-100 concentrations and longer wait times lead to surfactant-induced changes in the dynamic properties of the free surface in the tank. It is found that low surface concentrations of surfactants can dramatically change the wave breaking process by changing the shape of the jet and breaking up the entrained air cavity at the time of jet impact. Direct numerical simulations (DNS) of plunging breakers with constant surface tension are used to show that there is significant compression of the free surface near the plunging jet tip and dilatation elsewhere. To explore the effect of this compression/dilatation, the surface tension isotherm is measured in all experimental cases. The effects of surfactants on the plunging jet are shown to be primarily controlled by the surface tension gradient ($\Delta \mathcal {E}$) while the ambient surface tension of the undisturbed wave tank ($\sigma _0$) plays a secondary role.

     
    more » « less