大洋环流模式设计初步.ppt

Lax-Wendroff 格式不是‘单调’格式 ? 满足Courant条件的‘迎风’格式是稳定和正定的,但耗散太强； ? Lax-Wendroff 格式具有二阶精度，耗散很小,但频散太强； ? 推测:上述两种格式的适当组合可以构成高精度的保形格式. 小 结 ? 先用Lax-Wendroff 格式积分一步得到 ‘试探’解, 用于检测出‘单调’性有可能被破坏的格点； 正式积分时, 在那些被检测出来的格点上改用‘迎风’格式. 两种格式混合使用, 扬长避短. 两步‘保形’方案的要点： 宇如聪, 1994: A Two-Step Shape-Preserving Advection Scheme Adv. Atmos. Sci., 11, 479-490. 对两步保形平流格式的检验（Yu, R.-C., 1994） 初值 “迎风差分” Lax-Wendroff TSPAS after 10 cycles 克服高精度格式频散效应的其他方法举例: Leap-frog格式+时间滤波 第四讲的回顾 大洋环流模式是建立在什么基础之上的？ 波动过程, 特别是和海表高度变化有关的表面重力波. ‘迎风’格式和波动过程差分格式的基本问题. 我们了解了那些模式原理？ ‘刚盖’近似模式和自由面模式—大洋环流模式的历史. 差分格式设计追求的主要目标之一—低耗散、弱频散. 对常用的海洋模式 (如MOM) 的原理有所了解. ? Equations of the model* ? Arrangement of variables ? General method of finite differencing Bryan (1969) 的设计要点 *其中包括‘rigid-lid’近似. 差分方法:用有限的格点系统逼近连续系统 0? 180? 360? 赤道 南极 北极 B-Grid Kirk Bryan GFDL 模式的主要研究者 K. Bryan M. Cox A. Semtner R. Pacanowski S. Griffies Bryan-Cox-Semtner model MOM 2004年12月26日印度洋海啸 /wiki/Tsunami * 5 xi tn 收敛的‘迎风’格式 和 CFL 条件 j j-1 j+1 n+1 n 差分方程解的依赖区域 Courant-Friedrichs-Lewy条件 保证了‘迎风’格式解的依赖区间包含微分方程解的依赖区间, 这是近似解收敛的必要条件. ‘迎风’格式的稳定性分析 α 1 0 1 f 对‘迎风’格式而言, CFL条件 也是计算稳定性条件, 它表明:时间步长要受到波速和空间网格距的限制. Rigid-lid 近似 Bryan-Cox-Semtner 模式 Type Atmos. Ocean Gravity waves External First internal mode 300 100 200 3 Currents Jets Interior 150 --- 1.5 0.2 Velocities of physical phenomena which may limit the time step of a numerical integration in an atmospheric or ocean model in units of ms-1. From Bryan (1984) 计算稳定性条件 C=200m/s (表面波) ? t 4 min C=3m/s (内波) ? t 4.6 hrs ? x = 50 km Bryan, K., 1969: A numerical method for the study of the circulation of the World Ocean. J. Comp. Phys., 4, 347-376. (‘rigid-lid’模式) Crowley, W.P., 1968: A Global numerical ocean model: Part I, J. Comp. Phys., 3, 111-147. (free-surface模式) 两篇海洋模式设计的论文 海表高度 z0 的预报方程 ‘刚盖’近似~正压流无辐散. = 0 U, V : 垂直积分流速(正压流) Vertical displacements of the ocean surface are not allowed. This constraint is called the ‘rigid-lid’ approximation. (Bryan,1969) 垂直积分流速U满足的方程 其中X只依赖于平