非线性大气化学动力学方程组数值解法比较.pdfVIP

第30卷 第4期 气 象 科 学 Vo1．30，No．4 2010年 8月 SCIENTIA METE0R0L0GICA SINI‘CA Aug．，2010 张欣 ，王体健 ，沈凡卉，等．非线性大气化学动力学方程组数值解法的比较．气象科学，2010，30(4)：427-437． ZhangXin， WangTijian，ShenFanhui，eta1．Thecomparisonofnumericalschemesfornonlinearatmosphericchemicalkineticequations．Scien— tiaMeteorologicaSinica，2010，30(4)：427’437． 非线性大气化学动力学方程组数值解法的比较 张 欣 王体健 沈凡卉 赵雨舟 (1南京大学 大气科学学院，南京210093)(2南京大学数学系，南京210093) 摘 要 大气化学模式中描述化学反应动力学的是一组高度耦合、刚性、非线性的常微分方程 组，高效精确的计算方法在空气质量模拟预测、气候化学相互作用研究中尤为重要。本文介绍了求 解此类非线性常微分方程组 (ODEs)的七种常见的数值计算方法，包括 QSSA、MQSSA、HYBRID、 MHYBRID、LSODE、GongCho以及TWOSTEP方法，并基于同一反应机理进行24h的数值模拟， 通过对模拟结果的讨论分析，比较了各种方法的优劣。研究结果表明，就计算效率而言，QSSA方 法最为省时。就计算精度而言，LSODE方法最优。若综合考虑计算精度和效率，HYBRID、MHY— BRID是比较适 中的数值解法。 关键词 大气化学动力学方程组 非线性 刚性 数值解法 分类号 P421．2 文献标识码 A Thecomparisonofnumericalschemesfornonlinear atmosphericchem icalkineticequations ZhangXin WangTijian ShenFanhui ZhaoYuzhou (1SchoolofAtmosphericSciences，NanjingUniversity，Nanjing210093，China) (2DepartmentofMathematics，Na ngUniversity，Na ng210093，China) Abstract Thechemicalreactionsinatmosphericchemistrymodelsaredescribedbyagroupofhigh— lycoupled，stiffandnon—linearordinarydifferentialequations(ODEs)，towhichahigh—accuracyand fastnumericalmethodisparticularlyimportant．Inthispaper，wepresentandcomparetheresultsofsev— eralmethods，includingQSSA，MQSSA，HYBRID，MHYBRID，LSODE，GongChoandTWOSTEP bya24hnumericalsimulationusingthesamechemicalmechanism．Amongallthosemethods，QSSA， TWOSTEP，MHYBRID，HYBRIDandMQSSAarerelativelyfasterthanLSODEandGong＆Cho，while LSODEisthebestintermsofaccuracyandQSSAcoststheleastamountofcup—time．Inbalance，we recommendHYBRIDandMHYBRID tobethebetterchoiceswithsatisfactory performanceinbothpreci-