四章节大洋环流模式设计初步学习资料.pptVIP

第四讲大洋环流模式设计初步; 引言: 模式设计过程 基本近似, 方程组和边界条件 海表高度预报方程, 表面重力波 C-F-L条件，计算稳定性 “刚盖”近似, B-C-S模式 自由面模式, 模态分解算法 ‘耗散’和‘频散’，典型差分方案;海表高度预报方程 和 表面重力波;年平均海表高度 (LICOM1.0, cm);z0(x,y,t);;;正压流方程的推导;u-动量方程右端项的分析;;压力的 external / internal 分量;典型重力波过程;表面重力波—正压速度和海表压力的相互调整;表面波 (External gravity waves);? 缓慢变化的海表高度 ? 快速传播的表面波？ ? 表面波的计算问题？;对正压模态方程组的分析;垂直积分流 U 满足的方程;? 在Boussinesq近似下, 海表高度的变化由体积守恒决定. ? 海表高度的大尺度分布取决于风应力和斜压压力梯度. ? 表面重力波和‘地转调整’.;练习题. 证明恒等式 ：;垂直积分流严格满足的方程;;Courant-Friedrichs-Lewy 条件 差分格式的计算稳定性;模型问题: 一维行波方程;一般的波动方程;一维波动方程解的d’Alembert 公式;;;;Courant-Friedrichs-Lewy条件;‘迎风’格式的稳定性分析;;对‘迎风’格式而言, CFL条件;;Rigid-lid 近似 Bryan-Cox-Semtner 模式;Type;计算稳定性条件;Bryan, K., 1969: A numerical method for the study of the circulation of the World Ocean. J. Comp. Phys., 4, 347-376. (‘rigid-lid’模式);海表高度 z0 的预报方程;Vertical displacements of the ocean surface are not allowed. This constraint is called the ‘rigid-lid’ approximation. (Bryan,1969);垂直积分流速U满足的方程;正压无辐散流可以用流函数表示;正压流函数的预报方程;风应力强迫下的正压涡度方程;为什么rigid-lid模式较早获得成功？ ? 滤去了表面波，模式只包含慢过程，可用大时间步长积分，计算量较小 ; ? 虽然要求解关于正压流函数倾向的椭圆型方程，但早在1960年代椭圆型方程的迭代解法已经比较成熟. ;;“刚盖”近似模式的局限性 ? 不能直接预报海表高度, 因而不适用于海表高度资料的同化; ? 滤去了表面重力波，不适用于潮汐研究. —有必要发展自由面模式.;自由面模式 和 模态分解算法;? 将表面波模态(external mode)从模式中分离出来, 用小步长积分; ? 对描写垂直结构的模态(internal mode)用大步长积分; ? 正确描写两种模态的相互作用.;流速的斜压和正压分量;垂直积分流 U 满足的方程;对正压模态方程组的分析;刚盖近似模式;Blumburg, A.F., and G.L. Mellor, 1987: A description of a three-dimensional coastal ocean circulation model (POM). 刘海龙等：LICOM1.0用户手册;;差分格式的耗散和频散 (dissipation / dispersion);迎风格式的‘单调’性和强耗散性; ‘单调’性的充要条件;高精度格式的频散效应;Lax-Wendroff 格式(二阶精度, 低耗散 );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）;克服高精度格式频散效应的其他方法举例: Leap-frog格式+时间滤波;;第四讲的回顾; 大洋环流模式是建立在什么基础之上的？ 波动过程, 特别是和海表高度变化有关的表面重力波.