Breakdown of the slow manifold in the Shallow-Water equations

Abstract

Numerical solutions are obtained by implicit multigrid solvers for initial-value problems in the rotating Shallow-Water Equations (SWE) with spatially complex initial conditions. Companion solutions are also obtained with the Shallow-Water Balance Equations (SWBE), both to determine the initial conditions for the SWE and to provide a comparison solution that lies entirely on the slow, advective manifold. We make use of a control parameter (here the Rossby number, R) to regulate the degree of slowness and balance. While there are measurable discrepancies between the evolving SWE and SWBE solutions for all R, there is a distinct, spatially local breakdown both of the slow manifold in the SWE solution and in the closeness of correspondence between the SWE and SWBE solutions. This critical value for breakdown is only slightly smaller than the R values at which, first, the SWE evolution becomes singular (i.e., the fluid depth vanishes), or second, a consistent initial condition for the SWBE cannot be defined. This breakdown is most clearly evident in a sudden increase in vertical velocity near the center of a strong, cyclonic vortex; its behavior is primarily associated with an enhanced dissipation rather than an initiation of gravity-wave propagation. The numerical performance of the multigrid solvers is satisfactory even in the difficult circumstances near solution breakdown or singularity.     

An elliptic problem with integral constraints with application to large-scale geophysical flows

We present a weak formulation of a non-standard elliptic equation whose boundary values are determined in part by integral relations. Existence and uniqueness of its solution are proved, and a finite element discretization is described, analyzed, and implemented on a test problem. The equation is a generalization of one that is solved during integration of the three-dimensional Quasigeostrophic equations, which model large-scale rotating stratified flows, where the integral constraints represent conservation of physical properties.