简介:为提高数值求解大气方程的速度,研究了预处理JFNK(Jacobian—FreeNewton—Krylov)方法及其在大气方程中的应用。这是一种非线性外循环Newton迭代与线性内循环Krylobv迭代相结合的快速算法,其优点是进行外循环Newton迭代时不要求Jacobian矩阵的形成和存储;它的有效性取决于内循环中线性系统的预处理。首先介绍了JFNK算法,然后以浅水波方程为例,描述了非线性残值的形成、预处理矩阵的构造及其在JFNK算法中的应用。试验结果表明:对内循环线性系统进行适当的预处理,能大幅度提高JFNK算法的运算速度。因而,JFNK是一种值得在大气方程中推广应用的方法。
简介:利用北京大学校园地区PMIO质量浓度观测资料、中国科学院大气物理研究所325rn气象塔气象要素梯度和湍流观测资料,分析了北京地区2010年3月20-22日两次强沙尘暴过程微气象学要素和沙尘参量的时空演变以及湍流输送特征,为理解北京地区强沙尘暴天气沙尘输送规律和微气象学特征提供参考。结果表明:3月20-22日强沙尘暴过程前后不同高度温度先升后降,气压和相对湿度则相反。强沙尘暴来临时,高层风速先迅速增大,低层风速增加略有滞后,风切变明显加强,PMIO浓度最大值和风速极大值出现时间较吻合。强沙尘暴过境时,不同高度向下的湍流动量输送、向上的湍流热量输送和湍流动能明显加强。与3月21日非沙尘暴日相比,强沙尘暴过程湍流动量通量增加,有利于沙尘粒子的水平和垂直输送过程;由于冷锋过境,水平热通量增大;垂直热通量因白天温度垂直梯度减小而减小,夜间因逆温层被破坏而增加;水平湍流动能对湍流动能占主要贡献,垂直湍流动能仅占水平湍流动能的10%-25%。