球坐标系下MASNUM海浪数值模式的建立及其应用.pdf

第 27 卷 　第 2 期 海 　　洋 　　学 　　报 Vol . 27 , No . 2 2005 年 3 月 A C TA OCEAN OL O GICA SIN ICA March 2005 球坐标系下 MASNUM 海浪数值模式的建立及其应用 杨永增1 ,2 ,乔方利1 ,2 ,赵 　伟1 ,2 ,滕 　涌1 ,2 ,袁业立1 ,2 ( 1. 国家海洋局 第一海洋研究所 ,山东 青岛 26606 1 ; 2 . 海洋环境科学与数值模拟国家海洋局重点实验室 , 山东 青 岛 26606 1) 摘要 :为开展海浪对海洋上混合层的搅拌混合作用及其对海气界面通量的影响等研究 ,在L A GFD WA M 区域海浪数值模式基础上建立了球坐标系下的全球海浪数值模式. 重点导出了球坐标系下 的海浪能量谱平衡方程及其复杂特征线方程 ,该组方程包含了背景流场对波动传播的调整 、波动沿 大圆传播的折射等. 数值积分则采用复杂特征线嵌入计算格式. 初步数值模拟结果表明 ,该海浪全 球数值模式能够较为精确地刻画海浪的动力过程. 关键词 :海浪数值模式 ;球坐标 ;特征线计算格式 中图分类号 : P73 122 　　　文献标识码 : A 　　　文章编号 : 02534 193 (2005) 02000 107 [ 5 ] 能量平衡方程积分计算. 杨春成等 开展了该模式 1 　引言 在深水大洋区域与 WA M 模式在近海与浅水区域 海浪过程是海洋动力环境演变的重要环节 ,对 的嵌套计算研究. 尹宝树等[ 6 ] 建立了 YWSW P 第 海气界面动量 、热量及物质通量具有重要的影响 ,它 二代海浪数值预报模式 ,利用文氏理论风浪谱进行 也是海洋上混合层与跃层形成的关键因素[ 1 ] . 海浪 非线性能量转移项的计算. 针对 WA M 第三代海浪 各统计要素是海洋工程设计和海上运输等经济活动 数值模式对大浪模拟偏低的问题 ,袁业立等[7 ,8 ] 建 的重要参数. 立了L A GFDWA M 第三代海浪数值模式 ,破碎耗 现今海浪数值计算方法大都是在能谱概念基础 散源函数利用 Yuan 等[9 ] 改进的参数化方案 ,非线 上建立的 ,且依据能谱传播 、增长与衰减的物理控制 性波 - 波相互作用源函数参照 Ha sselmann 等[ 10 ] 的 演化方程数值积分来进行模拟 ,不同模式之间所采 参数化方案 ,在数值模式中首次考虑了波流相互作 用的能量输入 、耗散源函数略有不同 ,数值计算格式 用源函数 ,数值计算格式上则提出了物理上更为合 也不尽一致. 国际上具有较广泛影响的是第三代海 理直观的复杂特征线嵌入计算格式 ,该模式广泛应 浪数值模式 WA M [2 ] ,破碎耗散源函数采用 Ko men 用到 中国近海 的数值模拟和海洋工程 中. Yang 等[3 ] 的参数化方案 ,数值计算采用差分方法 ,现在已 等[ 11 ,12 ] 在该模式基础上建立了海浪变分同化模式. 发展到 Cycle 4 版本. 国内关于海浪数值模式的