Multilayer shallow water equations with complete Coriolis force. Part 2. Linear plane waves

2011 ◽  
Vol 690 ◽  
pp. 16-50 ◽  
Author(s):  
Andrew L. Stewart ◽  
Paul J. Dellar

AbstractWe investigate the behaviour of linear plane waves in multilayer shallow water equations that include a complete treatment of the Coriolis force. These equations improve upon the conventional shallow water equations, based on the traditional approximation, that include only the part of the Coriolis force due to the locally vertical component of the rotation vector. Including the complete Coriolis force leads to dramatic changes in the structure of long linear plane waves. It allows subinertial waves to exist with frequencies below the inertial frequency, the minimum frequency for which waves exist under the traditional approximation. These subinertial waves are characterized by a distinguished limit in which the horizontal pressure gradient becomes comparable to the upwellings and downwellings driven by the non-traditional Coriolis term in the vertical momentum equation. The subinertial waves connect wave modes that remain separate in the conventional multilayer shallow water equations, such as the surface and internal waves in a two-layer system. Eastward-propagating surface waves in a two-layer system connect with westward-propagating internal waves, and vice versa, via the long subinertial waves. The long subinertial waves cannot be classified as either surface or internal waves, due to the phase difference between the disturbances to the interfaces in these waves.

2013 ◽  
Vol 723 ◽  
pp. 289-317 ◽  
Author(s):  
Andrew L. Stewart ◽  
Paul J. Dellar

AbstractWe analyse the hyperbolicity of our multilayer shallow water equations that include the complete Coriolis force due to the Earth’s rotation. Shallow water theory represents flows in which the vertical shear is concentrated into vortex sheets between layers of uniform velocity. Such configurations are subject to Kelvin–Helmholtz instabilities, with arbitrarily large growth rates for sufficiently short-wavelength disturbances. These instabilities manifest themselves through a loss of hyperbolicity in the shallow water equations, rendering them ill-posed for the solution of initial value problems. We show that, in the limit of vanishingly small density difference between the two layers, our two-layer shallow water equations remain hyperbolic when the velocity difference remains below the same threshold that also ensures the hyperbolicity of the standard shallow water equations. Direct calculation of the domain of hyperbolicity becomes much less tractable for three or more layers, so we demonstrate numerically that the threshold for the velocity differences, below which the three-layer equations remain hyperbolic, is also unchanged by the inclusion of the complete Coriolis force. In all cases, the shape of the domain of hyperbolicity, which extends outside the threshold, changes considerably. The standard shallow water equations only lose hyperbolicity due to shear parallel to the direction of wave propagation, but the complete Coriolis force introduces another mechanism for loss of hyperbolicity due to shear in the perpendicular direction. We demonstrate that this additional mechanism corresponds to the onset of a transverse shear instability driven by the non-traditional components of the Coriolis force in a three-dimensional continuously stratified fluid.


Author(s):  
Vladimir Zeitlin

The derivation of the rotating shallow-water model by vertical averaging is carried on in the tangent plane approximation without neglecting the vertical component of the Coriolis force, and contributions of the vertical component of velocity in its horizontal component (‘non-traditional’ terms), leading to one- and two-layer ‘non-traditional’ rotating shallow-water models. A similar approach on the whole sphere encounters difficulties with conservation of angular momentum. Consistent ‘non-traditional’ rotating shallow-water equations in this case are obtained from the variational principle, which is first formulated for full primitive equations. It is shown that columnar motion hypothesis should be replaced by solid-angle motion one on the sphere. Two-layer non-traditional rotating shallow-water equations are used to analyse inertial instability of jets and compare the results with Chapter 10. It is shown that non-traditional terms can increase the growth rates up to 30% in some configurations and can also change the structure of the unstable modes.


2014 ◽  
Vol 748 ◽  
pp. 789-821 ◽  
Author(s):  
Marine Tort ◽  
Thomas Dubos ◽  
François Bouchut ◽  
Vladimir Zeitlin

AbstractConsistent shallow-water equations are derived on the rotating sphere with topography retaining the Coriolis force due to the horizontal component of the planetary angular velocity. Unlike the traditional approximation, this ‘non-traditional’ approximation captures the increase with height of the solid-body velocity due to planetary rotation. The conservation of energy, angular momentum and potential vorticity are ensured in the system. The caveats in extending the standard shallow-water wisdom to the case of the rotating sphere are exposed. Different derivations of the model are possible, being based, respectively, on (i) Hamilton’s principle for primitive equations with a complete Coriolis force, under the hypothesis of columnar motion, (ii) straightforward vertical averaging of the ‘non-traditional’ primitive equations, and (iii) a time-dependent change of independent variables in the primitive equations written in the curl (‘vector-invariant’) form, with subsequent application of the columnar motion hypothesis. An intrinsic, coordinate-independent form of the non-traditional equations on the sphere is then given, and used to derive hyperbolicity criteria and Rankine–Hugoniot conditions for weak solutions. The relevance of the model for the Earth’s atmosphere and oceans, as well as other planets, is discussed.


Author(s):  
Xiao-Hua Zhu ◽  
Xiao-Hua Zhu ◽  
Ze-Nan Zhu ◽  
Ze-Nan Zhu ◽  
Xinyu Guo ◽  
...  

A coastal acoustic tomography (CAT) experiment for mapping the tidal currents in the Zhitouyang Bay was successfully carried out with seven acoustic stations during July 12 to 13, 2009. The horizontal distributions of tidal current in the tomography domain are calculated by the inverse analysis in which the travel time differences for sound traveling reciprocally are used as data. Spatial mean amplitude ratios M2 : M4 : M6 are 1.00 : 0.15 : 0.11. The shallow-water equations are used to analyze the generation mechanisms of M4 and M6. In the deep area, velocity amplitudes of M4 measured by CAT agree well with those of M4 predicted by the advection terms in the shallow water equations, indicating that M4 in the deep area where water depths are larger than 60 m is predominantly generated by the advection terms. M6 measured by CAT and M6 predicted by the nonlinear quadratic bottom friction terms agree well in the area where water depths are less than 20 m, indicating that friction mechanisms are predominant for generating M6 in the shallow area. Dynamic analysis of the residual currents using the tidally averaged momentum equation shows that spatial mean values of the horizontal pressure gradient due to residual sea level and of the advection of residual currents together contribute about 75% of the spatial mean values of the advection by the tidal currents, indicating that residual currents in this bay are induced mainly by the nonlinear effects of tidal currents.


2021 ◽  
Vol 70 ◽  
pp. 31-44
Author(s):  
E. Audusse ◽  
V. Dubos ◽  
A. Duran ◽  
N. Gaveau ◽  
Y. Nasseri ◽  
...  

We investigate in this work a class of numerical schemes dedicated to the non-linear Shallow Water equations with topography and Coriolis force. The proposed algorithms rely on Finite Volume approximations formulated on collocated and staggered meshes, involving appropriate diffusion terms in the numerical fluxes, expressed as discrete versions of the linear geostrophic balance. It follows that, contrary to standard Finite-Volume approaches, the linear versions of the proposed schemes provide a relevant approximation of the geostrophic equilibrium. We also show that the resulting methods ensure semi-discrete energy estimates. Numerical experiments exhibit the efficiency of the approach in the presence of Coriolis force close to the geostrophic balance, especially at low Froude number regimes.


Sign in / Sign up

Export Citation Format

Share Document