An optimization of the icosahedral grid modified by spring dynamics.

*(English)*Zbl 1058.76591Summary: We have investigated an optimum form of the modified icosahedral grid that is generated by applying the spring dynamics to the standard icosahedral grid system. The spring dynamics can generate a more homogeneous grid system than the standard icosahedral grid system by tuning the natural spring length: as the natural spring length becomes longer, the ratio of maximum grid interval to minimum one becomes closer to unit. When the natural spring length is larger than a critical value, however, the spring dynamic system does not have a stable equilibrium. By setting the natural spring length to be the marginally critical value, we can obtain the most homogeneous grid system, which is most efficient in terms of the CFL condition. We have analyzed eigenmodes involved in the initial error of the geostrophic balance problem [test case 2 of D. L. Williamson et al., J. Comput. Phys. 102, 211–224 (1992; Zbl 0756.76060)]. Since the balance state in the discrete system differs slightly from the exact solution of the analytic system, the initial error field includes both the gravity wave mode and the Rossby wave mode. As the results of the analysis are based on Hough harmonics decompositions, we detected Rossby and gravity wave modes with zonal wavenumber 5, which are asymmetric against the equator. These errors are associated with icosahedral grid structure. The symmetric gravity wave mode with zonal wavenumber 0 also appears in the error field. To clarify the evolution of Rossby waves, we introduce divergence damping to reduce the gravity wave mode. From the simulated results of the geostrophic problem with various grid systems, we found that the spuriously generated Rossby wave mode is eliminated most effectively when the most homogeneously distributed grid system is used. It is therefore concluded that the most homogeneous grid system is the best choice from the viewpoint of numerical accuracy as well as computational efficiency.

##### MSC:

76M25 | Other numerical methods (fluid mechanics) (MSC2010) |

86A10 | Meteorology and atmospheric physics |

76B15 | Water waves, gravity waves; dispersion and scattering, nonlinear interaction |

76U05 | General theory of rotating fluids |

##### Software:

chammp
PDF
BibTeX
XML
Cite

\textit{H. Tomita} et al., J. Comput. Phys. 183, No. 1, 307--331 (2002; Zbl 1058.76591)

Full Text:
DOI

**OpenURL**

##### References:

[1] | Cheong, H., Double Fourier series on a sphere: applications to elliptic and vorticity equations, J. comput. phys., 157, 327, (2000) · Zbl 0961.76062 |

[2] | Cheong, H., Application of double Fourier series to the shallow-water equations on a sphere, J. comput. phys., 165, 261, (2000) · Zbl 0986.76064 |

[3] | Roache, P.J., computational fluid dynamics, (1976), Hermora Publishers |

[4] | Sadourny, R.; Arakawa, A.; Mintz, Y., Integration of the nondivergent barotropic vorticity equation with an icosahedral hexagonal grid for the sphere, Mon. wea. rev., 96, 351, (1968) |

[5] | Williamson, D.L., Integration of the barotropic vorticity equation on a spherical geodesic grid, Tellus, 20, 642, (1968) |

[6] | Ringler, T.D.; Heikes, R.H.; Randall, D.A., Modeling the atmospheric general circulation using a spherical geodesic grid: A new class of dynamical cores, Mon. wea. rev., 128, 2471, (2000) |

[7] | Majewski, D., The new global icosahedral-hexagonal gridpoint model GME of the deutscher wetterdienst, in proceeding of a seminar held at ECMWF on recent developments in numberical methods for atmospheric modelling, (Sept. 7-11, 1998), p. 172 |

[8] | Majewski, D.; Liermann, D.; Prohl, P.; Ritter, B.; Buchhold, M.; Hanisch, T.; Paul, G.; Wergen, W., The operational global icosahedral – hexagonal gridpoint model GME: description and high-resolution test, Mon. wea. rev., 130, 319, (2002) |

[9] | Stuhne, G.R.; Peltier, W.R., Vortex erosion and amalgamation in a new model of large scale flow on the shere, J. comput. phys., 128, 58, (1996) · Zbl 0861.76045 |

[10] | Masuda, Y.; Ohnishi, H., An integration scheme of the primitive equation model with an icosahedral-hexagonal grid system and its application to the shallow water equations, Short- and medium-range numerical weather prediction: collection of papers presented at the WMO/IUGG NWP symposium, Tokyo, 317, (Aug. 4-8, 1986) |

[11] | Cullen, M.J.P., Integration of the primitive equations on a sphere using the finite element method, Q. J. R. met. soc., 100, 555, (1974) · Zbl 0279.65090 |

[12] | Heikes, R.H.; Randall, D.A., Numerical integration of the shallow-water equations on a twisted icosahedral grid. I. basic design and results of tests, Mon. wea. rev., 123, 1862, (1995) |

[13] | Heikes, R.H.; Randall, D.A., Numerical integration of the shallow-water equations on a twisted icosahedral grid. II. A detailed description of the grid and analysis of numerical accuracy, Mon. wea. rev., 123, 1881, (1995) |

[14] | Stuhne, G.R.; Peltier, W.R., New icosahedral grid-point discretizations of the shallow water equations on the sphere, J. comput. phys., 148, 23, (1999) · Zbl 0930.76067 |

[15] | Thuburn, J., A PV-based-shallow-water model on a hexagonal-icosahedral grid, Mon. wea. rev., 125, 2328, (1997) |

[16] | Sadourny, R., Numerical integration of the primitive equations on a spherical grid with hexagonal cells, Proceedings of the WMO/IUGG symposium on numerical weather prediction in Tokyo, tech. rep. of JMA, VII45-VII52, (Nov. 26-Dec. 4, 1969) |

[17] | Williamson, D.L., Integration of the primitive barotropic model over a spherical geodesic grid, Mon. wea. rev., 98, 512, (1969) |

[18] | Cullen, M.J.P.; Hall, C.D., Forecasting and general circulation results from finite element models, Q. J. R. met. soc., 105, 571, (1979) |

[19] | Tomita, H.; Tsugawa, M.; Satoh, M.; Goto, K., Shallow water model on a modified icosahedral geodesic grid by using spring dynamics, J. comput. phys., 174, 579-613, (2001) · Zbl 1056.76058 |

[20] | Williamson, D.L.; Drake, J.B.; Hack, J.J.; Jacob, R.; Swarztrauber, P.N., A standard test set for numerical approximations to the shallow water equations in spherical geometry, J. comput. phys., 102, 211, (1992) · Zbl 0756.76060 |

[21] | Longuet-Higgins, M.S., The eigenfunctions of Laplace’s tidal equations over a sphere, Phil. trans. roy. soc. London, A262, 511, (1968) · Zbl 0207.27501 |

[22] | Swarztrauber, P.N.; Kasahara, A., The vector harmonic analysis of Laplace’s tidal equations, SIAM J. sci. stat. comput., 6, 464, (1985) · Zbl 0572.76019 |

[23] | Skamarock, W.C.; Klemp, J., The stability of time-split numerical method s for the hydrostatic and the nonhydrostatic elastic equations, Mon. wea. rev., 120, 2109, (1992) |

[24] | Jakob-Chien, R.; Hack, J.J.; Williamson, D.L., Spectral transform solutions to the shallow water test set, J. comput. phys., 119, 164, (1995) · Zbl 0878.76059 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.