简介:为提高数值求解大气方程的速度,研究了预处理JFNK(Jacobian—FreeNewton—Krylov)方法及其在大气方程中的应用。这是一种非线性外循环Newton迭代与线性内循环Krylobv迭代相结合的快速算法,其优点是进行外循环Newton迭代时不要求Jacobian矩阵的形成和存储;它的有效性取决于内循环中线性系统的预处理。首先介绍了JFNK算法,然后以浅水波方程为例,描述了非线性残值的形成、预处理矩阵的构造及其在JFNK算法中的应用。试验结果表明:对内循环线性系统进行适当的预处理,能大幅度提高JFNK算法的运算速度。因而,JFNK是一种值得在大气方程中推广应用的方法。
简介:利用2012年夏季山西省太原、大同、临汾3个温室气体观测站的近地面O_3及相关前体物NOx、NO2、NO、CO观测数据及同期气象观测资料,分析山西省夏季O_3体积分数与污染状况、O_3的时间变化特征及O_3与其前体物体积分数的相关性,同时分析气温、相对湿度、风速、降水和日照时数等气象因素对O_3体积分数的影响。结果表明:山西太原、大同、临汾3个城市的O_3小时体积分数与日最大8h平均体积分数都有超标情况发生,夏季O_3污染以临汾最严重、大同相对最轻;一日内O_3小时体积分数最大值出现在15:00左右,最小值出现在06:00左右,日变化呈单峰型分布,其中临汾O_3小时体积分数的昼夜变化振幅最大、大同最小;各城市NOx、NO2、NO、CO等前体物体积分数均呈现白天低、夜间高的日变化过程,与O_3日变化呈负相关;气温是影响夏季O_3体积分数的最重要因素,其次为相对湿度,风速贡献最小,夏季O_3体积分数高值多在高温低湿的午后,且太原、临汾O_3体积分数受局地气象因素影响比大同显著。