高级检索
李江浩, 彭新东. 阴阳网格上质量守恒计算性能分析[J]. 大气科学, 2013, 37(4): 852-862. DOI: 10.3878/j.issn.1006-9895.2012.12060
引用本文: 李江浩, 彭新东. 阴阳网格上质量守恒计算性能分析[J]. 大气科学, 2013, 37(4): 852-862. DOI: 10.3878/j.issn.1006-9895.2012.12060
LI Jianghao, PENG Xindong. Analysis of Computational Performance of Conservative Constraint on the Yin-Yang Grid[J]. Chinese Journal of Atmospheric Sciences, 2013, 37(4): 852-862. DOI: 10.3878/j.issn.1006-9895.2012.12060
Citation: LI Jianghao, PENG Xindong. Analysis of Computational Performance of Conservative Constraint on the Yin-Yang Grid[J]. Chinese Journal of Atmospheric Sciences, 2013, 37(4): 852-862. DOI: 10.3878/j.issn.1006-9895.2012.12060

阴阳网格上质量守恒计算性能分析

Analysis of Computational Performance of Conservative Constraint on the Yin-Yang Grid

  • 摘要: 质量守恒数值计算是球面准均匀阴阳网格构造全球大气环流模式的重要条件,也是提高阴阳网格应用质量的重要技术手段。本文针对通量形式平流方程,在球面坐标上采用多种理想数值试验对阴阳网格上的三种守恒计算方案和边界插值非守恒计算方案进行了比较检验。发现,质量守恒方案不仅对全球数值积分重要,还影响数值计算精度,满足局地守恒条件的全球强迫守恒方法可以获得较高的精度;网格内质量均匀分布的阴阳网格边界通量一致性守恒强迫计算方案,实现了在不增加计算误差条件下保证局地和全球守恒的目的,且具有很小的计算负担,可以作为阴阳网格上全球质量强迫守恒的有效计算方案;而网格质量的线性分布可以有效提高阴阳网格的数值积分计算精度,但在一定程度上会增加计算负担。

     

    Abstract: Mass conservation is an important factor in the construction of the general circulation model on the quasi-uniform overset spherical Yin-Yang grid, and is also a key issue when applying the Yin-Yang grid to numerical models. With the help of a flux-form advection equation, we carried out a series of idealized numerical tests on a sphere to assess three conservative constraints in addition to the non-conservative boundary interpolation scheme. The test cases included a standard solid-body advection, a benchmark deformational flow, and a sine-wave advection. The results showed that the conservative constraint of mass property is important for global numerical integration and has an effect on the computational accuracy. The global conservative constraint that ensures local mass conservation gave accurate results. The boundary flux identification based on cell-wise mass constant distribution on the Yin-Yang grid achieved local and global mass conservation without either decreasing the computational accuracy or increasing computing costs. This algorithm is computationally efficient and can be an effective numerical scheme for mass conservative constraints on the overset Yin-Yang grid. However, while the boundary flux identification based on cell-wise linear distribution can improve the computational accuracy, it requires much more calculation time, which will limit its application in a dynamical model.

     

/

返回文章
返回