藉不可逆热力学算子技术对飓风Andrew模拟精度的改进.pdfVIP

藉不可逆热力学算子技术对咫风Andrew模拟精度的改进 刘 英 柳崇健 康红文 (中国气象科学研究院，北京 100081) 1 引 言 大气动力学控制方程组中引人水平扩散项虽 动力享的一个组成部分，但最初是基于 计算技巧上的考虑，旨在克服非线性混淆或计算不稳定，同时光滑物理场。但从目前通用的二 阶/四阶扩散算子的数学表示以及实际计算/模拟的效果看，无论线性或非线性的偶阶扩散方案 均不能恰当地反映扩散的物理实质，从而影响到预报和模拟的精 文从扩散过程的物理耗 散特性出发，藉刻划多体系统不可逆过程演变方向的热力学第二 进了M5V3的水平扩 散方案，以一个具有解析解的著名的Burgers方程和给美国带来严重社会经济后果的咫风An drew(1992)例“具体说明本文建议的水平扩散方案的一些改进细节。〕 2 趋于定常流动的一维Burgers方程 一维的Burgers方程可写为 au十。au_李a2u (1) dt dx R,a, 式中u二 “(x,t)为流速，R。系雷诺数，本文取R。二100.0。这种无强迫的耗散流最终将趋 于定常流。例如，若取初始条件呈线性分布: u :二 一 x (t=0， 1x 1) (2) 则在边界条件 ! 1 1 X = 一 1 U - . J 1 1 (3) . ~ 1 X 之下，存在一个特解(参见忻孝康等1980): _ R u=一th(才x) (4) 显然，这是一个不含时间的定常解。若采用通常的(中央)差分方案，则以初始条件(2)出发积分 到一定时间(在2s之内)即会以一定的精度逼近此解。 然而，按通常的中央差分方案，最后求得的数值解在速度梯度突增的 “激波”附近出现了 Gibbs振荡(参见图1)。只有将(四阶)方案修正到符合热力学第二定律约束的物理扩散形式， 才获得了与解析解一致的数值解(图1)0 3 在美国佛罗里达州登陆的一个咫风实例 本文以开始迅速加强的1992年8月21日12时UTC的分析场作为初始条件用MM5V3 (参见Dudhia1993及Grelletal.1995)积分到登陆后的8月24日12时UTC，共进行了72h的 对比模拟试验。 188 一 . ，”“辛件苦行 门 」 0oa 一 . . a.oo 0.50 1- 图I趋于定常流动的一维Burgere方程采用新的光滑技术前后数值与精确解的比较(详见正文)。 下面拟着重分析对咫风描述而言极具重要性的地面低气压和地面最大风速在引人新的扩