Filtering inertia-gravity waves from the initial conditions of the linear shallow water equations

2007 ◽  
Vol 19 (3-4) ◽  
pp. 204-218 ◽  
Author(s):  
Alexander Barth ◽  
Jean-Marie Beckers ◽  
Aida Alvera-Azcárate ◽  
Robert H. Weisberg
2020 ◽  
Vol 35 (6) ◽  
pp. 355-366
Author(s):  
Vladimir V. Shashkin ◽  
Gordey S. Goyman

AbstractThis paper proposes the combination of matrix exponential method with the semi-Lagrangian approach for the time integration of shallow water equations on the sphere. The second order accuracy of the developed scheme is shown. Exponential semi-Lagrangian scheme in the combination with spatial approximation on the cubed-sphere grid is verified using the standard test problems for shallow water models. The developed scheme is as good as the conventional semi-implicit semi-Lagrangian scheme in accuracy of slowly varying flow component reproduction and significantly better in the reproduction of the fast inertia-gravity waves. The accuracy of inertia-gravity waves reproduction is close to that of the explicit time-integration scheme. The computational efficiency of the proposed exponential semi-Lagrangian scheme is somewhat lower than the efficiency of semi-implicit semi-Lagrangian scheme, but significantly higher than the efficiency of explicit, semi-implicit, and exponential Eulerian schemes.


2013 ◽  
Vol 20 (1) ◽  
pp. 25-34 ◽  
Author(s):  
A. Wirth

Abstract. Using a fine resolution numerical model (40002 × 2 grid-points) of the two-layer shallow-water equations of the mid-latitude β-plane dynamics, it is shown that there is no sudden breakdown of balance in the turbulent enstrophy cascade but a faint and continuous emission of inertia–gravity waves. The wave energy accumulates in the equator-ward region of the domain due to the Coriolis parameter depending on the latitude and dispersion relation of inertia–gravity waves.


2015 ◽  
Vol 137 (11) ◽  
Author(s):  
Abdellah Kharicha ◽  
Jan Bohacek ◽  
Andreas Ludwig ◽  
Menghuai Wu

A numerical model based on the shallow water equations (SWE) was proposed to simulate the two-dimensional (2D) average flow dynamics of the liquid metal spreading inside a horizontally rotating mold. The SWE were modified to account for the forces, such as the centrifugal force, Coriolis force, shear force with the mold wall, and gravity. In addition, inherent vibrations caused by a poor roundness of the mold and the mold deformation due to temperature gradients were applied explicitly by perturbing the gravity and the axis bending, respectively. Several cases were studied with the following initial conditions: a constant average height of the liquid metal (5, 10, 20, 30, and 40 mm) patched as a flat or a perturbed surface. The angular frequency Ω of the mold (∅1150–3200) was 71.2 (or 30) rad/s. Results showed various wave patterns propagating on the free surface. In early stages, a single longitudinal wave moved around the circumference. As the time proceeded, it slowly diminished and waves traveled mainly in the axial direction. It was found that the mean amplitude of the oscillations grows with the increasing liquid height.


2012 ◽  
Vol 69 (1) ◽  
pp. 303-322 ◽  
Author(s):  
Mark D. Fruman ◽  
Ulrich Achatz

Abstract The three-dimensionalization of turbulence in the breaking of nearly vertically propagating inertia–gravity waves is investigated numerically using singular vector analysis applied to the Boussinesq equations linearized about three two-dimensional time-dependent basic states obtained from nonlinear simulations of breaking waves: a statically unstable wave perturbed by its leading transverse normal mode, the same wave perturbed by its leading parallel normal mode, and a statically stable wave perturbed by a leading transverse singular vector. The secondary instabilities grow through interaction with the buoyancy gradient and velocity shear in the basic state. Which growth mechanism predominates depends on the time-dependent structure of the basic state and the wavelength of the secondary perturbation. The singular vectors are compared to integrations of the linear model using random initial conditions, and the leading few singular vectors are found to be representative of the structures that emerge in the randomly initialized integrations. A main result is that the length scales of the leading secondary instabilities are an order of magnitude smaller than the wavelength of the initial wave, suggesting that the essential dynamics of the breaking might be captured by tractable nonlinear three-dimensional simulations in a relatively small triply periodic domain.


2007 ◽  
Vol 64 (12) ◽  
pp. 4417-4431 ◽  
Author(s):  
Chris Snyder ◽  
David J. Muraki ◽  
Riwal Plougonven ◽  
Fuqing Zhang

Abstract Vortex dipoles provide a simple representation of localized atmospheric jets. Numerical simulations of a synoptic-scale dipole in surface potential temperature are considered in a rotating, stratified fluid with approximately uniform potential vorticity. Following an initial period of adjustment, the dipole propagates along a slightly curved trajectory at a nearly steady rate and with a nearly fixed structure for more than 50 days. Downstream from the jet maximum, the flow also contains smaller-scale, upward-propagating inertia–gravity waves that are embedded within and stationary relative to the dipole. The waves form elongated bows along the leading edge of the dipole. Consistent with propagation in horizontal deformation and vertical shear, the waves’ horizontal scale shrinks and the vertical slope varies as they approach the leading stagnation point in the dipole’s flow. Because the waves persist for tens of days despite explicit dissipation in the numerical model that would otherwise damp the waves on a time scale of a few hours, they must be inherent features of the dipole itself, rather than remnants of imbalances in the initial conditions. The wave amplitude varies with the strength of the dipole, with waves becoming obvious once the maximum vertical vorticity in the dipole is roughly half the Coriolis parameter. Possible mechanisms for the wave generation are spontaneous wave emission and the instability of the underlying balanced dipole.


2012 ◽  
Vol 706 ◽  
pp. 71-107 ◽  
Author(s):  
Noé Lahaye ◽  
Vladimir Zeitlin

AbstractWe study formation and properties of new coherent structures: ageostrophic modons in the two-layer rotating shallow water model. The ageostrophic modons are obtained by ‘ageostrophic adjustment’ of the exact modon solutions of the two-layer quasi-geostrophic equations with the free surface, which are used to initialize the full two-layer shallow water model. Numerical simulations are performed using a well-balanced high-resolution finite volume numerical scheme. For large enough Rossby numbers, the initial configurations undergo ageostrophic adjustment towards asymmetric ageostrophic quasi-stationary coherent dipoles. This process is accompanied by substantial emission of inertia–gravity waves. The resulting dipole is shown to be robust and survives frontal collisions. It contains captured inertia–gravity waves and, for higher Rossby numbers and weak stratification, carries a (baroclinic) hydraulic jump at its axis. For stronger stratifications and high enough Rossby numbers, ‘rider’ coherent structures appear as a result of adjustment, with a monopole in one layer and a dipole in another. Other ageostrophic coherent structures, such as two-layer tripoles and two-layer modons with nonlinear scatter plot, result from the collisions of ageostrophic modons. They are shown to be long-living and robust, and to capture waves.


Author(s):  
Rick Salmon

This second chapter offers a brief introduction to geophysical fluid dynamics—the dynamics of rotating, stratified flows. We start with the shallow water equations, which govern columnar motion in a thin layer of homogeneous fluid. Roughly speaking, the solutions of the shallow-water equations comprise two types of motion: ageostrophic motions, including inertia-gravity waves, on the one hand, and nearly geostrophic motions on the other. In rapidly rotating flow, these two types of motion may, in some sense, decouple. We seek simpler equations that describe only the nearly geostrophic motion. The simplest such equations are the quasigeostrophic equations. In the quasigcostrophic equations, potential vorticity plays the key role: The potential vorticity completely determines the velocity field that transports it, thereby controlling the whole dynamics. We begin by generalizing our previously derived fluid equations to a rotating coordinate frame.


2008 ◽  
Vol 65 (5) ◽  
pp. 1622-1637 ◽  
Author(s):  
J. Vanneste

Abstract This paper discusses some of the mechanisms whereby fast inertia–gravity waves can be generated spontaneously by slow, balanced atmospheric and oceanic flows. In the small Rossby number regime relevant to midlatitude dynamics, high-accuracy balanced models, which filter out inertia–gravity waves completely, can in principle describe the evolution of suitably initialized flows up to terms that are exponentially small in the Rossby number ɛ, that is, of the form exp(−α/ɛ) for some α > 0. This suggests that the mechanisms of inertia–gravity wave generation, which are not captured by these balanced models, are also exponentially weak. This has been confirmed by explicit analytical results obtained for a few highly simplified models. These results are reviewed, and some of the exponential-asymptotic techniques that have been used in their derivation are presented. Two types of mechanisms are examined: spontaneous-generation mechanisms, which generate exponentially small waves from perfectly balanced initial conditions, and unbalanced instability mechanisms, which amplify unbalanced initial perturbations of steady flows. The relevance of the results to realistic flows is discussed.


Sign in / Sign up

Export Citation Format

Share Document