Numerical Dispersion of Gravity Waves

Schroeder, Guido; Schlünzen, K. Heinke
December 2009
Monthly Weather Review;Dec2009, Vol. 137 Issue 12, p4344
Academic Journal
When atmospheric gravity waves are simulated in numerical models, they are not only dispersive for physical but also for numerical reasons. Their wave properties (e.g., damping or propagation speed and direction) can depend on grid spacing as well as on the numerical schemes. In this work numerical dispersion relations for atmospheric gravity waves are theoretically derived as well as experimentally measured using the anelastic Mesoscale Transport and Stream model (METRAS). Both the theoretical solution and the numerical model show a retardation of gravity waves with decreasing grid resolution. Furthermore, the influence of a Shapiro seven-point filter is analyzed. The Shapiro seven-point filter causes damping of the shorter waves. Therefore, shorter waves can better be simulated without the seven-point filter. The influence of different advection schemes is analyzed by prescribing a background wind. A first-order upstream scheme and second- and third-order flux integrated essentially nonoscillatory (FIENO) schemes are used. As expected, the damping is the smaller the higher the order of the scheme. The numerical dispersion has severe consequences, when nonuniform grid spacing is used. Waves moving from the fine grid to the coarse are reflected because of numerical dispersion if they are only poorly resolved on the coarse grid. In tests with different refinement factors and wave lengths the reflection is found to be the larger the greater the refinement factor. The results show that refinement factors larger than 3 should not be used with nonuniform grid spacing or two-way nested grids.


Related Articles

  • Numerical study of double diffusive convection in a cylindrical annular geometry Validation and mesh sensibility. Cheddadi, Bouchra; Aboulaich, Rajae; Cheddadi, Abdelkhalek // International Journal of Computer Science Issues (IJCSI);Jan2014, Vol. 11 Issue 1, p96 

    In this paper, we consider the problem of double diffusion in an annular cylindrical space with an aspect ratio R=2. The ADI method has been used to describe the two-dimensional symmetric solutions in terms of isotherms, streamlines and iso-concentrations. This numerical study has been validated...

  • Cracks, stress risers no problem for superconverging FEA code. Mraz, Stephen J. // Machine Design;9/26/96, Vol. 68 Issue 17, p32 

    Reports on Superconvergent Adaptive General Elements, a finite-element-analysis package. Capabilities of the new package; Handling of singularities of high-order elements and advanced algorithms.

  • Efficient List Ranking on the Reconfigurable Mesh, with Applications. Hayashi, T.; Nakano, K.; Olariu, S. // Theory of Computing Systems;1998, Vol. 31 Issue 5, p593 

    Finding a vast array of applications, the list-ranking problem has emerged as one of the fundamental techniques in parallel algorithm design. Surprisingly, the best previously known algorithm to rank a list of n items on a reconfigurable mesh of size n × n was running in O(log n) time. It...

  • An adaptive solution of the 3-D Euler equations on an unstructured grid. Zhu, Z.Q.; Wang, P.; Tuo, S.F. // Acta Mechanica;2002, Vol. 155 Issue 3/4, p215 

    Presents an adaptive algorithm of the three-dimensional surface unstructured grid. Spatial discretization of the cell-centered finite volume formulation on the unstructured grid; Application of a novel cell reconstruction procedure; Euler solution for an adaptive tetrahedral grid.

  • Design of an Atmospheric Model Based on a Generalized Vertical Coordinate. Konor, Celal S.; Arakawa, Akio // Monthly Weather Review;Jul97, Vol. 125 Issue 7, p1649 

    Although there are important advantages in the use of an isentropic vertical coordinate in atmospheric models, it requires overcoming computational difficulties associated with intersections of coordinate surfaces with the earth's surface. In this paper, the authors present a model based on...

  • Quadrilateral Mesh. Ming, Pingbing; Shi, Zhongci // Chinese Annals of Mathematics;Apr2002, Vol. 23 Issue 2, p235 

    Several quadrilateral shape regular mesh conditions commonly used in the finite element method are proven to be equivalent. Their influence on the finite element interpolation error and the consistency error committed by nonconforming finite elements are investigated. The effect of the...

  • Grid Computing. Foster, Ian // AIP Conference Proceedings;2001, Vol. 583 Issue 1, p51 

    The term “Grid Computing” refers to the use, for computational purposes, of emerging distributed Grid infrastructures: that is, network and middleware services designed to provide on-demand and high-performance access to all important computational resources within an organization or...

  • Worldwide computing working group summary. Gaines, Irwin // AIP Conference Proceedings;2001, Vol. 583 Issue 1, p357 

    The worldwide computing working group brought together a number of experts in this rapidly emerging field. A combination of presented talks and discussion identified common interests and issues to be pursued in the future. © 2001 American Institute of Physics.

  • THE ENTIRE DOMAIN AFTER CREATION OF ITS MESH.  // Diagram;May/Jun2012, Vol. 12 Issue 3, p3 

    An image of a mesh from the 1991 publication of "Automatic Mesh Generation" by P. L. George is presented.


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics