第四讲大洋环流模式设计初步.pptVIP

第四讲大洋环流模式设计初步 张学洪 (zxh@lasg.iap.ac.cn) 引言: 模式设计过程 基本近似, 方程组和边界条件 海表高度预报方程, 表面重力波 C-F-L条件，计算稳定性 “刚盖”近似, B-C-S模式 自由面模式, 模态分解算法 ‘耗散’和‘频散’，典型差分方案 大洋环流模式是建立在什么基础之上的？ 波动过程, 特别是和海表高度变化有关的表面重力波. ‘迎风’格式和波动过程差分格式的基本问题. ‘刚盖’近似模式和自由面模式—大洋环流模式的历史. 差分格式设计追求的主要目标之一—低耗散、弱频散. 对常用的海洋模式 (如MOM) 的原理有所了解. 宇如聪, 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格式+时间滤波 第四讲的回顾 我们了解了那些模式原理？ α 1 0 1 f 对‘迎风’格式而言, CFL条件 也是计算稳定性条件, 它表明:时间步长要受到波速和空间网格距的限制. Rigid-lid 近似 Bryan-Cox-Semtner 模式 1.5 0.2 150 --- Currents Jets Interior 200 3 300 100 Gravity waves External First internal mode Ocean Atmos. Type 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) 计算稳定性条件 ? t 4.6 hrs C=3m/s (内波) ? t 4 min C=200m/s (表面波) ? 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只依赖于平流和水平湍流粘性. 正压无辐散流可以用流函数表示 正压流函数的预报方程 椭圆型算子 风应力强迫下的正压涡度方程 假定: 1. H=const 是无运动深度; 2. 平流项和水平湍流项可略 Sverdrup balance 为什么rigid-lid模式较早获得成功？ ? 滤去了表面波，模式只包含慢过程，可用大时间步长积分，计算量较小 ; ? 虽然要求解关于正压流函数倾向的椭圆型方程，但早在1960年代椭圆型方程的迭代解法已经比较成熟. L4, 4X5 (1989) L20, 4X5 (1994) L30,T63 (1999) LICOM (2004) LASG Semtner (1997) POP MOM “刚盖”近似模式的局限性 ? 不能直接预报海表高度, 因而不适用于海表高度资料的同化; ? 滤去了表面重力波，不适用于潮汐研究. —有必要发展自由面模式. 自由面模式 和 模态分解算法 ? 将表面波模态(external mode)从模式中分离出来, 用小步长积分; ? 对描写垂直结构的模态(internal mode)用大步长积分; ? 正确描写两种模态的相互作用. 如何减少自由面模式计算量？ 流速的斜压和正压分量 U , V — external mode 的速度场. 垂直积分流 U 满足的方