天气动力学与诊断分析：任意形状有限区域调和-余弦谱分解方法.pptxVIP

任意形状有限区域调和-余弦谱分解方法 ;Method in Rectanglular Domain;Method in Rectanglular Domain;Method in Rectanglular Domain;Method in Rectanglular Domain;Method in Rectanglular Domain;规则区域;我国地形复杂： 复杂海岸线、山脉盆地 数值模式：非结构化网格;不规则区域;Capcitance Matrix Method Step 2 在矩形区域内求解Helmholtz方程 Step 3 求右端修正项 ;Capcitance Matrix Method Step 4 在矩形区域内求解修正Helmholtz方程 Dirichlet边界 将地形格点分为12类求线积分，首先求其中8类拐角边界格点的值 利用拐角格点值和两个拐角点之间线积分求其他边界格点的制 ;方法流程;Example 1;Example 1;Example 1;Example 1;Example 1;Example 2;Example 2;Example 2;Example 2;Example 2;优： 迭代次数少，计算精度高， 能有效还原风场。 劣： 谱展开，谱求导，一个点的资料有问题就会影响全区域的计算精度;为解决这个问题，美国XU Qin教授和大气所曹洁博士合作提出了一种新方法，分别求解内部、边界和资料洞对应的分量，然后求和。算法化流程是：从解析解出发，通过构造多种离散化方案提高精度， 设计自动识别内外边界方案，实现全程自动化。从右边两幅图可以看到：在这个有不规则形状资料洞的不规则有限区域里，重建流场完全还原了真实流场。真正实现任意形状有限区域风场分解的高效、高精度和自动化。 小结：提出并建立的用于任意形状有限区域（包括地形资料洞）的高效、高精度、且全程自动化的风场分解方法，可直接用于台风等中小尺度天气系统诊断分析中。;任意形状有限区域风场分解的差分计算方法及其诊断分析;涡度散度 外边界 内边界(地形洞);3°设计自动识别内外边界方案。;任意形状有限区域y, c计算; (c, y) = (ci, yi) + (ce, ye) + (cd, yd) 内场解(ci, yi) = (2p)-1 ∫Ddx’a(x’)ln|x-x’|; 无B.C. 资料空区解(cd, yd) = ?1(2p)-1ln|x-xk|∫Skds(vn-vni) 外场解(ce, ye)：B.C. 可解性条件： ∫ Skds?nce = ∫ Skds(vn - vni - vnd) =0 ∫ Skds?nye = ∫ Skds(vs - vsi - vsd) =0 解不唯一问题：设ce=0（或ye=0） 联合求解:复速度势we=je-iye=(i2p)-1∫Sdzwe(z)/(z’-z) 单一边条件求解：Dirichlet, Neumann, 切向导数 求解边界密度函数方程 → 线积分得到外场解; 离散化???案及数值检验 (ci, yi) ：ln|x-x’| -- 奇异点问题 交错格点、局部加密、分段连续函数积分 (cd, yd) ：边界及边界值确定；点源位置 两种边界自动搜索方案；空区的几何中心 (ce, ye)～we=je-iye：构造兼容的虚部?nye=?sje ye(s)=ye(s0)+ ∫s0sds’(-vne), je(s)=je(s0)+ ∫0sds’vse, 在S上； 分段连续函数积分方案： we(z) = (i2?)-1∑l ?l(z)?s ?l(z) = [?e(zl+1) - ?e(zl)]+??e(zl)+[?e(zl+1)-?e(zl)](z-zl)/(zl+1-zl)? ln[(zl+1 - z)/(zl - z)]; 内场解、外场解离散化方案示意图;(ce, ye)～边界密度函数形式解： Dirichlet B.C.：pm(xl) + ∑l′m(xl’)??ll’ = ye(xl)； ye(xij) = ∑l′m(xl′)??ijl′ 矩形区域角点的特殊处理：角点区域局部加密、 分段连续积分 m(xl’)??ll’ = ∫l′[m(xl′) + m(xl′+1)sl′/?s]dqll′ Neumann B.C.：pn(xl) - ∑l′n(xl’)??l’l = ?nye(xl)； ye(xij) = -∑l′n(xl′)(ln|xl′ - xij|)?sl 切向导数B.C.：-∑l′n(xl’)?rl′l/rl′l = ?sye(xl) 后两类B.C.问题都可转化到Dirichlet问题求解。; 数值试验结果 08UTC,