自适应网格技术中权函数的问题地研究.pdfVIP

自适应网格技术中权函数问题的研究+ 康红文王鹏云徐祥德 (中国气象科学研究院) 提要 利用变分原理的自适应网格技术被应用劲两个有解析解的问题上。权函数的选取耐网格的合理调度至关重要。对 运动激渡算例·权函数考虑二阶导数项比考虑梯度值更能有效地减小误差；对气旋锋生算例，仅以梯度作为权函数并 不能提高精度，考虑速度场、锋生函数分布的权函数能更台理地安排网格；通过与解析解对比发现，自适应网格技术 在提高精度、节省内存方面显示了优良的性能。 关键词：自适应网格，权函数，数值模拟． 1 引言 中尺度数值模式在预报具有强对流的天气系统时，不可避免地要遇到多尺度结构的问题．比如： 强风暴、锋面降水等，一方面要考虑其中的中小尺度结构特征：另一方面，还要考虑影响它们的天气 背景的大尺度结构流场及其反馈作用。目前．解决这一问题广泛采用的是套网格技术。然而，二十多 年来-套网格技术中粗细网格之间插值所产生的误差及内边界所存在的伪反射问题，一直令气象学家 们头疼。 本文所介绍的自适应网格技术正是针对包含多尺度结构流场进行数值模拟提出的。自适应网格技 术首先从航空航天领域发展起来的，它首先根据各网格间数值解情况取一个权函数W，并使得 W·Ax=常数。这样，某网格间上权函数越大，相应的△工就越小。为达到这一目的．先将其转化为 一个变分问题，通过求解其欧拉方程来得到新网格点位置…。然后通过坐标变换，把物理域上不规整 网格变换到计算域上规整的网格，使得阿格的移动写进控制方程中。并在每一步的求解中得以实现， 所提出。1980年代在计算流体力学领域已有较广泛的应用，井在1990年代应用到气象学上，比如我国 学者刘卓在1991年曾对这一技术在气象学上应用有过专门地研究pJ。美国也将这一技术列入下一代模 式的研制中。然而到目前为止，这一技术在气象业务上还没得到使用．原因主要有两个：一是在求解 网格方程时．需用迭代法。由于天气形势的变化，每预报一步．均需重新计算一次新网格点位置．迭 代法的使用，将使计算代价大大增加：二是权函数的确定。由于选择不同的权函数，将有不同的网格 分布。能否选选择一个好的权函数，从而充分发挥自适应网格技术的效率，至关重要。本文将对此进 行探讨。 2两个算例 本文选择两个具有解析解的算例来进行试验。第一个算倒是运动激波问题。其控制方程 汐比 Ou 1护2“ 万蜘万2面万了 初始条件 当r：o，“：J。 7 1， 扣o 【o，0工≤l 边界条件 f1，X=0 肛托x：l 其精确解为14 ’资助课题；田家重点基础研究发展规划项目G1998040907。IO国家攀登专项95．专4)3 哇H蓑矿-ed蓑)} ㈤ 的带动下的对流扩散问题。 第二个算铡是气旋锋生算倒。为研究大气环流中常见的气旋和锋的生消过程，bosewe¨(1984) 分析了一个由于涡旋流场的作用而产生的气旋锋生成过程【51。 锋区由一个任意的标量场口(如位温)的密集带代表。考虑一个二维情况，其初始场在X方向无变 化，在Y方向有一梯度，即 Q=一tanhy (3) 假定Q为守恒量，(当0为位温时，即假定为绝热过程)．则 塑+“塑+v塑：o 西 西 西 (4) 水平风场给定为一个定常无辐散的气旋性涡旋。即风向沿围绕中心点的圆的切线逆时针方向，风 速沿径向方向变化，径向风速和切向风速分别为： U=0 吁=＆曲2rtanhr