An error evaluation on the vertical velocity algorithm in POM
-
摘要: 时间分裂技术是若干自由面海洋模式的常用技术。内外模态方程的不同的截断误差要求进行必要的数值调整,以保证算法能够正确地满足连续性方程,且使得示踪量守恒。POM模式应用了一种简单的将内模态流速的垂向平均调整为与垂向积分的外模态流速一致。但是,由于控制数值不稳定的Asselin时间平滑算法的引入,POM的速度调整方法不再能保证连续性方程严格成立,即使它采用了特殊的处理根据外模态水位来定义内模态的水位。误差被证明是Asselin平滑算子的二阶项。该误差的在数值模式中的一个影响为海底的运动学边界条件不再满足。通过一个区域模拟和一个准全球模拟实验,对该误差的量级进行了估计,并进行了若干敏感性实验。文章分析了该误差的特征,并提出了两个减少误差的算法。
-
关键词:
- 模态分离,计算误差,POM模式
Abstract: A time splitting technique is common to many free surface ocean models. The different truncation errors in the equations of the internal and external modes require a numerical adjustment to make sure that algorithms correctly satisfy continuity equations and conserve tracers quantities. The princeton ocean model (POM) has applied a simple method of adjusting the vertical mean of internal velocities to external velocities at each internal time step. However, due to the Asselin time filter method adopted to prevent the numerical instability, the method of velocity adjustment used in POM can no longer guarantee the satisfaction of the continuity equation in the internal mode, though a special treatment is used to relate the surface elevation of the internal mode with that of the external mode. The error is proved to be a second-order term of the coefficient in the Asselin filter. One influence of this error in the numerical model is the failure of the kinetic boundary condition at the sea floor. By a regional experiment and a quasi-global experiment, the magnitudes of this error are evaluated, and several sensitivity tests of this error are performed. The characteristic of this error is analyzed and two alternative algorithms are suggested to reduce the error.-
Key words:
- mode splitting
-
Courant R, Friedrichs K, Lewy H. 1967. On the partial difference equations of mathematical physics. IBM Journal, 3: 215-234 Da Silva A, Young C, Levitus S. 1994. Atlas of Surface Marine Data 1994, vol. 1, Algorithms and Procedures. Washington D C, US: Dep of Commer, 74 Ezer T, Arango H, Shchepetkin A F. 2002. Developments in terrainfollowing ocean models: intercomparison of numerical aspects. Ocean Modell, 4: 249-267 Higdon R L, Bennett A F. 1996. Modeling, stability analysis of operator for large-scale ocean modeling. J Comput Phys, 123: 311-329 Higdon R L, de Szoeke R A. 1997. Barotropic-baroclinic time splitting for ocean circulation modeling. J Comput Phys, 135: 31-53 Kantha L H, Clayson C A. 2000. Numerical Models of Oceans and Oceanic Processes. San Diego: Academic Press, 750 Mellor G L. 2003. Users Guide for a Three-dimensional, Primitive Equation, Numerical Ocean Model. June 2003 version. Princeton: Princeton University, 56 Munk W. 1966. Abyssal recipes. Deep-Sea Res, 13: 707-730 Shchepetkin A F, McWilliams J C. 2005. The regional ocean modeling system: A split-explicit, free-surface, topography following coordinates ocean model. Ocean Modell, 9: 347-404 Wunsch C, Ferrari R. 2004. Vertical mixing, energy, and the general circulation of the oceans. Annu Rev Fluid Mech, 36: 281-314 Xia Changshui, Qiao Fangli, Yang Yongzeng, et al. 2006. Threedimensional structure of the summertime circulation in the Yellow Sea from a wave-tide-circulation coupled model. J Geophys Res, 111: C11S-C13S Xia Changshui, Qiao Fangli, Zhang Qinghua, et al. 2004. Numerical modeling of the quasi-global ocean circulation based on POM. J Hydrodyn: Ser B, 16: 537-543 Zhou Weidong. 2002. A proper time integration with split stepping for the explicit free-surface modeling. Adv Atmos Sci, 19: 255-265
点击查看大图
计量
- 文章访问数: 1440
- HTML全文浏览量: 57
- PDF下载量: 1488
- 被引次数: 0