资源预览内容
第1页 / 共40页
第2页 / 共40页
第3页 / 共40页
第4页 / 共40页
第5页 / 共40页
第6页 / 共40页
第7页 / 共40页
第8页 / 共40页
第9页 / 共40页
第10页 / 共40页
亲,该文档总共40页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述
博士开题报告 有限区域三维流场分解方法 及其在中尺度系统中的 应用研究 指导老师:高守亭 研究员 许 秦 研究员 学 生: 曹 洁(2005级硕博连读) 2007年6月 1 选题的目的和意义 2 相关研究回顾 3 拟解决的关键问题和 4 已完成的工作 5 待完成的工作 6 本研究的主要创新点 研究手段 全风场看不清楚确切的中尺度信息 资料同化处理 数值模拟准确性 完全分解的旋转辐散和垂直运动信息 然而 : 中尺度运动学和动力学研究中: 风场信息在中尺度系统中有着重要的地位 雷达、卫星等高时空分辨率探测仪 中尺度数值模式的发展 问题得到一定缓解 常规的观测风场: 时空精度不够 选题的目的和意义(续 1) 选题的出发点 如何简单、有效地从全风场信息中 提取中尺度流场? 意义 1)提高灾害性天气监测水平; 2)深入理解中尺度系统结构机理等研究问题; 3)提高中尺度系统业务预报水平。 基本思想 水平风场 旋转风辐散风 三维流场 地转风强迫的非地转风 完全的非地转风 2-1)流函数、速度势的一般计算方法 2-2)调和正(余)弦方法和有限区域 水平风场分解方法 2-3)准地转近似下的三维流场分解理论 2-4)半平衡近似下的三维流场分解理论 相关研究回顾(续 1) 2-1)流函数、速度势的一般计算方法 对全球大气, 即周期条件下Poisson方程的求解: ok 对有限区域问题,即耦合边条件下Poisson方程的求解: 历史上有代表性的几个方法(图示)的基本思想: 在保证求得的满足D 分布,满足分布的条件下 ,人为地简化耦合边条件的形式(纯数学处理) 难以处理 1, 边条件 2, 解不唯一 相关研究回顾(续 2) Shukla and Saha(1974) 两种迭代方案 历史上有代表性的几个方法(图示) Phillips(1958) Brown and Neilon (1961) Sangster(1960 ) 相关研究回顾(续 3) 历史方法常导致: 计算不稳定、原始风场无法还原、边界上的系统缺 失、迭代次数多、结果很不准确 2-2)H-S(C)方法和有限区域风场分解理论 Chen等(1992a,b)提出的调和-正(余)弦展开法: 函数fi (x, y)满足: 完全由内部变量决定, 故称为内部部分 函数fh (x, y)满足: 完全由边界条件决定, 故称为外部部分 相关研究回顾(续 4) 有限区域风场分解和重建流程图 相关研究回顾(续 5) 2-3)准地转近似下的三维流场分解理论 基本思想(Xu,1992,1993): 从准地转近似下基本运动方程出发, 三维流场 准地转分量 由准地转分量激发的非地转A分量 包括重力波等信息的完全的非地转N分量 2-4)半平衡近似下的三维流场分解理论 对于大尺度锋面系统,准地转假设近似成立; 然而对于高度弯曲的锋面和旋涡, Allen(1990,1991), Xu(1994)等提出了更高阶近似下的三维流场分解理论。 资料:NCEP再分析资料,雷达资料,模式输出资料等。 研究手段: 算法改进, 诊断分析, 动力推导, 数值模拟, 综合集成化等 拟解决的关键问题: 1.改进传统天气图分析方法 中尺度系统的运动学诊断新思路 2.建立一套包含地形、热力影响的三维流场控制方程 分离天气尺度及中尺度系统的动力信息 3.运用能量学、动力学等途径找寻中尺度系统信息 为中尺度系统相关动力研究和业务应用提供新途径 本论文具体研究问题与常用的张驰法在实际个例中的效果比较 地形坐标系下辐散风和旋转风动能收支方程的研究 方程的理论推导和物理解释;一次次天气尺度涡旋生消过程的能量学研究 水平风场分解在天气分析中的应用 对一次华北暴雨过程的分析;从风场分解角度比较两次台风过程 准地转近似下的三维流场分解理论的流程化 地形坐标下准地转近似三维Psi方程的理论求解方法及其在 梅雨锋系统动力分析中的应用 方程的理论推导;二维风场分解方法对梅雨锋系统水平结构的运动学分析; 包含地形影响的准地转近似下的三维风场分解理论对梅雨锋系统三维结构、 发生发展机理、与周围环境场相互作用的动力分析 包含地形、热力效应的更高价近似下三维流场分解 理论及其在台风系统动力分析中的应用 方程的理论推导和算法构建;在ARPS模拟的一次台风个例中的应用 有限区域三维流场分解方法的业务化框架 4.1 改进的H-C算法及其在实际个例中的检验 从内部分量中分离出由平均涡度散度和齐次边值决定的 和 ; 后面求解内部分量时,以原始涡度与平均涡度之差作为方程右项 ,散度类似,得到了内部分量的另一部分: 和 ; 改进算法的基本思路 求解外部分量方法相同,得到 和 。 对范围任意小的区域可能就会出现问题了 H-C方法要求满足以下的必要性条件: 对稍大区域 一般都满足 改进的必要 性 已完成的工作(续 1)1979年3月30日00UTC, 500hPa, (41N,130E), 120km, NCEP资料 B CD A 风场的流线图 已完成的工作(续 2) 风场的 分量 U BA CD 已完成的工作(续 3) 风场的 分量 V AB CD 已完成的工作(续 4) 更多试验: 不同时次,不同层次,不同空间分辨率,不同 环流背景,不同资料来源, 结论: 改进的H-C方法能以很高的计算效率实现更高 的计算精度,完成风场分解和重建。 计算效率 及精度: 里查逊法:2000次迭代,9.209933 加速利布曼法:1000次迭代,9.202168 改进的H-C方法:3次迭代,0.4114308 已完成的工作(续5) 4.2 水平风场分解在天气分析中的应用 基本思路: 利用改进的H-C算法分解水平全风场,得到旋转风和 辐散风两个分量场信息; 在此基础上,应用传统天气图分析,试图得到更多的 流场信息。 改进的传统天气图分析方法 实验设计: 一、对2004年8月12-13日一次华北暴雨的诊断分析 ; 二、对2006年台风“桑美”和“碧利斯”的比较研 究。 着眼点:水平风场分解 已完成的工作(续6) 850hPa 全风场 旋转风 辐散风 12日08UTC 华北暴雨 已完成的工作(续7) 全风场 旋转风 辐散风 700hPa 12日08UTC 华北暴雨 已完成的工作(续 13) “桑美” 7日18UTC 全风场 旋转风 “碧利斯” 11日00UTC 已完成的工作(续 14) “桑美” 9日12UTC 全风场 辐散风 地面 高空 已完成的工作(续 15) “碧利斯” 11日06UTC 全风场 辐散风 地面 高空 已完成的工作(续 11) a)全风场,b)旋转风场,c)辐散风场,卫星云图,雷达图像 “桑美” 10日12UTC 已完成的工作(续 12) a)全风场,b)旋转风场,c)辐散风场,卫星云图,雷达图像 “碧利斯” 14日06UTC 已完成的工作(续 18) “桑美”:8月11日18UTC“碧利斯”:7月14日12UTC 低层水汽输送旋转风场 已完成的工作(续 14) 辐散风和旋转风能显示原始全风场看不 到的天气尺度系统和中尺度系统信息 优势 不足 1. 单纯的水平风场分析不足以提供中尺度系统中占 重要地位的垂直运动; 2. 对旋转风辐散风相互作用的表现不很明显。 从风场分解角度改进传统天气图分析方法: 能量学研究 三维流场分解 已完成的工作(续 15) 4.3 地形坐标系下辐散风和旋转风动能收支 方程的研究 过去的能量学研究表明: 辐散风起着 “催化剂”的作用;但对辐散风和旋转风的计算 欠精确,且出发方程不包括地形效应。 地形坐标系下的动力理论: 在模式构建中很有效;但直接应用等面上的位势高度得不 到与天气系统有关的信息 。 目的: 为中尺度研究提供新思路和新途径; 地形坐标系的新应用。 出发点: 辐散风和旋转风相互作用在传统的涡度方程和散度方程中 很难体现。 已完成的工作(续 16)有限区域旋转风和辐散风动能收支方程 已完成的工作(续 17) 坐标系下有限区域系统动能收支示意图 已完成的工作(续 18)一次次天气尺度涡旋生消过程的能量学研究 2004年8月12日在0.825上的旋转风场 06UTC09UTC 10UTC 14UTC 已完成的工作(续 19) 涡生消过程中区域- 垂直平均的四种动能 时间变率随时间变化 a) 涡发展阶段b) 涡衰亡阶段 已完成的工作(续 20)一些有趣的现象和结论: 两种风分量能反映天气尺度和中尺度系统的运动学差异 辐散风有“催化剂”般作用 考虑地形效应的必要性 风场分解准确性的重要性 动能转换项C(KD, KR)在涡旋发展到了顶峰之后发生反向; 科氏力项在低层起主要作用,垂直相关项随高度变化有很大差异。 GR和HKR是DKR稳定的能量源,HKD是DKD的汇,INVKD是主要源项; 其他一些大值项都与垂直运动有关,垂直梯度大,积分后贡献小。 INTR和GD在涡旋由发展向衰亡转变时变号,二者与辐散流有关。 CR和CD量值很小,与地面气压倾向有关; 地面气压倾向项有利于动能的增加,在涡衰亡阶段更加明显。 已完成的工作(续 21) 4.4 准地转三维流场分解理论的流程化 求解齐次边条件下的绝热、无摩擦、f平面的准地转 方程,产生的问题是: 得到的地转流强迫的部分非地转解不完整 为得到完整的有限区域三维非地转A分量环流的完整 解,Xu和Keyser(1993a,b)提出了两种不同的解决方法。 但是,存在几点不足: 1.没有算法化 2.缺少地形、热力影响 已完成的工作(续 22) 直接求解非齐次边条件下 的三维Psi方程 针对不同类型的边条件,求解得 到有限区域矢量Psi流函数 由 得到 分别求解非地转流的斜压 和正压分量 把地转风激发的非地转风分量 分成 齐次侧边条件决定的 非齐次边条件决定的 分别求解,由 得到 从 分解出位势高度的地转部分 ; 利用地转关系 ,得到地转风 ; 用不同的方法得到地转风激发的非地转风分量 的边条件: 和 由已知的观测风场 ,计算得到的地转风 及其激发的 非地转风分量 之差,得到完全的非地转风分量 。 1 地形坐标下准地转近似的三维Psi方程的 理论推导 正在进行; 难点:针对不同的模式有不同的出发方程 2 包含地形影响的准地转三维流场分解理 论在实际梅雨锋个例中的诊断和动力分析 即将开始个例的选取和数值模拟; 难点:模拟效果 待完成的工作(续 1) 3 探讨包含了地形、热力效应的更高价近似下 三维流场分解理论的可行性和求解方法 难点:理论推导和物理解释 4 更高价近似下的三维流场分解方法在实际个 例中的应用 难点:中尺度模式ARPS调试; 动力知识积累 5 完善本研究的理论和算法,建立一套可供业 务单位应用的分析系统 本研究期望达到的最终目标 发展普适的有限区域水平风场分解 算法和应用方法 为传统天气图分析提供新途径 运动学 理论化并算法化高阶近似下的三维 流场分解方法 为中尺度动力学研究提供新思路 动力学 系统化并业务应用中尺度系统动力 分析新方法 业务应用 Gao S T, and Cao J: Physical basis of generalized potential temperature and its application to cyclone tracks in non-unif
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号