高精度紧致平流算法的应用和区域天气精细化数值预报试验

来源 :兰州大学 | 被引量 : 0次 | 上传用户:baino1
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
平流是大气运动中最基本的动力过程,其数值计算对大气数值模式和数值天气预报起着至关重要的作用。平流计算方案的好坏不仅决定了模式大气运动计算的精度,还决定了数值模式积分的时间步长和计算效率。平流方案的特性对于天气系统的预报和模拟效果有显著的影响,特别对于一些强天气过程(如锋面、台风和强烈发展的温带气旋或冷涡等),物理量的不连续变化对数值模式的高精度平流计算和数值预报提出了挑战,因为通常高阶多项式构建的平流方案由于构建插值函数的模板较宽,往往对不连续物理量的平流计算缺乏正确的描述能力。另一方面,当采用低阶计算方案时会产生数值耗散、产生计算频散误差;二阶或更高阶的平流计算方案不可避免地产生数值振荡,对于模式中比湿等浓度计算出现负值振荡是致命的。随着数值预报模式分辨率的日益精细,模式变量不连续分布问题会因分辨率提高而更加突出,提高对锋面等不连续界面的计算和描述能力变得愈加重要。除大气数值模式中常用的有限差分平流计算方法外,还有许多高性能计算方法可以供气象模式使用,其中CIP(Constrained Interpolation Profiles)类的平流计算方法就是一个选择。这种方案采用同一网格点上配置多变量矩的Hermit类型插值方案构造半拉格朗日平流算法,所以可以在缩短插值函数模板的情况下构造高阶插值函数,从而在一个网格距内构建高阶连续函数,提高计算精度。本论文主要围绕CIP方案开展研究,包括一维和二维理想平流试验、调整植入MM5模式的CIP计算方案,并在我国一个暴雨过程个例数值模拟中对其影响进行分析。另外,结合模式计算需要,应用CIP方案开发了新型重力快波显式快速计算方法;论文的另一个主要工作就是结合实际数值预报的需求,利用中尺度非静力模式WRF探讨高分辨率精细化数值预报的初、边值应用方法。利用一维和二维简单理想流体的平流问题,对CIP法进行了计算和检验,在确信方法和计算精度要求的基础上,对Peng等(2003)引入MM5中的新型紧致CIP平流方案进行了适当的调整,并应用于我国南方2006年的一次梅雨暴雨过程的模拟研究。在对个例的模拟中,对比分析了CIP半拉格朗日方法与原有中央差欧拉平流方案的表现,证实CIP方案对这次暴雨过程具有较好的计算效果,更加真实地再现了梅雨锋降水的中尺度结构,改进了模式降水和锋面系统的模拟效果。对梅雨锋内对流运动和热力结构的模拟,显示了CIP对刻画大梯度分布和强对流运动计算的良好性能,平流的高精度计算对提高这类系统的模拟精度有明显效果。重力波是天气变化过程中有意义的动力过程,由于重力波激发一些重要天气过程,与地形性降水和对流天气过程中的能量传播有关,在数值模式中不可将这类快波滤掉,但由于重力波的传播速度较快,在数值模式积分过程中限制了模式积分时间步长。在大气数值模式中通常采用半隐式的算法处理重力波项,以扩展积分时间步长,但必须面对矩阵运算和变换,这对于高分辨模式而言是一个沉重的负担。由于CIP紧致平流方案属半拉格朗日型算法,不受积分线性稳定性条件限制,因此在提高计算精度的同时,还可以提高模式的计算效率。本研究在CIP法基础上开发的球面CIP特征线法计算重力波传播,以半拉格朗日算法代替半隐式算法,避免了矩阵运算,同时有效地克服了重力波对模式积分的线性稳定性限制,扩大时间步长。文中以球面浅水模式为例,直接使用半拉格朗日平流方案代替半隐式方案求解,并借助球面浅水模式的理想问题,给出了球面特征线法求解重力波的数值结果。采用CIP半拉格朗日方案,保证了模式数值积分对任意柯朗数的稳定性,节省大量CPU时间,运算在保证所需计算精度的同时,获得了极高的计算效率。作为示例,这里给出了在球面阴阳网格系统中理想化试验的数值解,以便评估。数值试验结果显示了CIP特征线方法在高分辨率模式中处理快重力波精确、经济和方便的特点,并且在理想试验中都获得了高分辨率大柯朗数(CFL>1.0)时稳定而精确的计算,这对开发和应用高分辨数值模式有重要价值。高分辨率、高精度中尺度数值预报和模拟是数值模式应用的目的。在目前计算能力不断提高和观测、分析资料不断丰富的条件下如何进行精细化的数值预报和数值模拟是我们面临的一个难题。中尺度数值预报受到初值、边值、模式动力和物理过程诸方面的影响,本研究的另一部分工作就是针对模式初、边值的影响来研究精细化中尺度数值预报的一个方面。运用WRF中尺度模式和高分辨率的T639全球预报结果,通过不同的初、边值形成方法来研究和寻找较好的预报结果。选择黄河下游的一次暴雨过程,对T639预报结果插值得到WRF初边值,分别采用连续积分和分段T639结果插值冷启动的方法进行降尺度处理,结果显示WRF模式连续积分的高分辨率嵌套数值预报结果较好,而T639预报结果分段插值冷启动预报较差。反映了具有丰富物理过程和动力过程的高分辨率WRF模式可以更好的获得天气过程发展的主要信息,尽管其边界仍然受到全球模式的影响,也说明T639静力模式对中尺度天气系统的描述较WRF中尺度模式粗糙,从而造成分段冷启动结果降水不足和偏离。这里呈现的仅为一个个例研究的结果,不具普遍性,有待更多实例模拟验证。
其他文献
椭圆问题是偏微分方程研究中最主要的问题之一,属于核心数学的研究范畴.早在1900年D.Hilbert提出的著名的23个问题中就有3个与椭圆问题有关.过去的一个多世纪椭圆问题得到了丰富的研究,获得了丰硕的成果,形成了庞大的理论体系,其在几何学、流体力学、热力学等自然科学和工程技术的各个领域都得到了广泛应用并有力地推动了这些学科的发展.本篇博士论文主要研究了一类从MEMS(micro electro
利用卫星观测资料、再分析资料,结合大气化学气候模式(WACCM),研究了热带上对流层下平流层(UTLS)区域一氧化碳(CO)时空变化和热带平流层CO准两年振荡(QBO)位相变化特征的形成机理,分析了亚洲夏季风反气旋环流造成的副热带地区向热带地区的CO水平输送对热带UTLS区域CO时空变化的影响,最后还分析了不同区域排放的大气示踪物向UTLS区域动力传输的特征。所得主要结论如下:(1)利用全球化学气
体心立方结构的铁及铁合金在核(原子能)工业领域是常见的结构材料。其中低活性铁素体/马氏体钢是聚变堆结构材料的候选材料。在14MeV中子的辐照下,由于(n,α)核反应,结构材料中将产生大量的He原子。一般结构材料是由多晶组成,包含大量的晶界,所以研究He泡在铁及铁合金晶界中的行为是发展耐辐照结构材料的重要研究课题之一。为了进一步理解多晶材料的辐照损伤效应,本文通过计算机来模拟缺陷在晶界中的演化。论文
细胞分裂素影响植物生长和发育的很多方面。在低温条件下植物会出现生长迟缓的现象.这里我们分离鉴定了一个在不同温度条件下生长速度有所差异的突变体cytokinin and temperature Sensitive Mutant (ckts)。突变体表现出在所有温度下,特别是在低温条件下主根生长缓慢,在低温生长环境中对外源生长素(IAA、2,4-D)敏感,并且使用细胞分裂素(6-BA、ZT)处理发现c
在本篇博士学位论文中,主要考虑在对数压力坐标下,关于三维大尺度大气运动原始简化方程的如下初边值问题:在区域上弱解的全局存在性,强解的全局存在唯一性以及在条件下强解全局吸引子的存在性.这里(x,y, z)∈Ω,其中:z=-Hs log(p/pm),Hs是定常高度尺度,p*是定常参考压力.并且记Ω的边界为:其中:对简化方程组(1),我们首先定义了水平速度v的平均v以及它的波动v并给出了它们的一些性质,
花粉发育及花粉管生长是有花植物有性生殖的关键步骤,受到很多蛋白和信号分子的调控,是一个非常复杂的生物学过程。膜联蛋白(Annexin)家族是一类进化上保守的且能结合钙和磷脂的多基因家族,该家族成员是否参与调控花粉发育和生长的研究还未见报道。在花粉萌发和花粉管生长过程中,微丝骨架在不同的时间和空间呈现出特定的结构并行使不同功能。ACTDIN DEPOLYMERIZING FACTOR (ADF)是一
线性系统的求解在数学、物理学、统计学、工程学甚至社会科学中的很多问题求解时都占有重要的地位.比如物理中的线性偏微分方程的近似求解就需要转化为线性系统的求解,所以线性系统数值求解历来都是数值线性代数研究的主要问题.随着大规模科学计算发展的需要,线性系统的求解越来越多地受到人们的关注.众所周知,求解线性系统的数值方法一般分为两类:直接法和迭代法.直接法的计算量小,但计算格式较复杂,因而适用于阶数不太高
图论的研究始于1736年,Euler用图的方法解决了哥尼斯堡(Konigsberg)七桥问题,并发表了第一篇关于图论的学术论文.从此,图论这门新的学科诞生了.20世纪60年代以来,图论在科学界异军突起,迅猛发展.同时,众多有趣的图论问题也应运而生,比如哈密尔顿问题,中国邮递员问题,染色问题,连通性问题,匹配问题等.并且,应用图论的方法解决化学,生物学,信息与计算科学等学科中的问题已显示出极大的优越
新疆大部分地区处于干旱区,导致人们对该地区的干旱问题习以为常,因而对其极端干旱问题重视不够。然而新疆是我国最大的内陆省份,占全国陆地总面积的1/6,其区域性极端干旱事件对当地的农牧业生产、水资源、民族团结和地区稳定以及丝绸之路经济带的繁荣都有重要的影响,对其极端干旱问题进行研究,既有重要的科学价值,又有很大的现实意义。本文首先针对干旱表征量SPI指数的适用性问题进行了研究,筛选出适合新疆降水量分布
从内质网转运囊泡到中央溶解液泡或者蛋白贮藏囊泡是植物中的一个主要的囊泡运输途径。囊泡运输需要SNARE复合体介导膜的融合。最近的研究表明,募集复合体对于内膜运输非常关键,推测它在运输囊泡和靶标膜的最初靠近过程中起作用。拟南芥MAG2蛋白是酵母Dsllp复合体TIP20亚基和哺乳动物syntaxin 18复合体RINT-1亚基的同源蛋白,参与在内质网和高尔基体之间贮藏蛋白前体的运输。文献报道囊泡运输