导读:本文包含了气溶胶热力学平衡论文开题报告文献综述及选题提纲参考文献,主要关键词:气溶胶预测,气溶胶动力学,气溶胶热力学平衡,HDMR方法
气溶胶热力学平衡论文文献综述
程余[1](2010)在《气溶胶动力学与热力学平衡预测的有效算法》一文中研究指出气溶胶是液态或固态微粒在空气中的悬浮体系。这些空气中悬浮颗粒的相关直径可以从几纳米到几十乃至几千微米。由于大气气溶胶在全球气候变化以及人类健康的问题中有着突出的作用,气溶胶模型也在大气环境预测中扮演了越来越重要的角色。气溶胶能够散射以及吸收外来的太阳日照和辐射,影响云层的热效应,进而影响到全球的气候。与此同时,由于气溶胶是非常微小的颗粒物质,那些直径在微米级别的细小颗粒容易被人体吸入,并进入人体内部,从而损害到人类的健康。除了这些不利的影响,另一方面,气溶胶微粒物质被广泛应用于工业生产中。这些微粒物质主要指的是细颗粒,它们在生产染料、黑炭、光纤、硅制品以及陶瓷粉等工业产品中都发挥了非常积极的作用。而这些消极和积极的作用都跟气溶胶颗粒的大小以及成分有着密不可分的关系。考虑气溶胶动力学以及热力学平衡问题时,模拟以及预测气溶胶的矢径分布和气溶胶的成分分布是研究气溶胶的两个关键性的工具(Moya et al.,2002).首先,在第一个研究课题中,越来越多的目光集中到了研究在气溶胶化学反应和动力过程中如何预测其矢径分布。从气溶胶渐渐受到人们重视发展至今,学者们提出了非常多的数值方法以及化学方法来求解气溶胶动力学方程,其中比较经典的方法有区域方法(sectional method) (Gelbard et al.,1980)、动量法(moment method) (Brock et al.,1987;Seo et al.,1990)、模式法(modal method)(Ackermann et al.,1998; Whitby et al.,1997)、随机法(stochastic approach)(Debry et al.,2003)等等。在这些方法中,传统的区域方法对扩散项有着较大的限制并且数值精度比较低;模式方法虽然有着较高的数值精度但是缺乏物理意义并且不能有效解决多气溶胶成分时的矢径分布问题;动量方法其本质是一种化学分析方法,它是基于单种气溶胶成分的物理以及化学属性所得到的微粒大小分布,但是它也同样不适用于多成分的气溶胶矢径分布预测问题;而随机方法的最大限制就是最终得到数值结果的误差精度得不到一个满意的保证。在近期的研究成果中,Sandu和Borden(Sandu et al.,2003)构建出关于气溶胶动力学方程有限元方法的框架。有限元方法在解决此类偏微分方程无论在时间方向还是在颗粒大小方向上都有着较高的精度。进一步地,Sandu(Sandu et al.,2006)成功的将Runge-Kutta方法和分片线性多项式相结合提出时间方向更高阶的数值算法。Liang(Liang etal.,2008)在其发表的论文中提出了分片小波函数应用于气溶胶动力方程的时间以及颗粒大小这两个方向,得到了很好的精度和效果。但是上述方法都忽略了气溶胶动力方程的重要特点:第一,气溶胶的增长过程产生了方程中的对流项;第二,气溶胶的凝并过程在方程中是一个非线性项。这两项的结合构成了气溶胶动力学方程最大的特征。除此之外,根据气溶胶的观测结果,气溶胶颗粒大小的分布在欧拉坐标下变化非常剧烈而在对数坐标下呈正态分布。众所周知,对具有陡峭峰值的对流占优问题进行数值逼近是一件较为困难的事情。很多求解对流占优问题的数值算法虽然达到了相应的精度,但是是条件收敛的,不能计算所有的增长因子和凝并因子,或者对时间步长有所限制;而另一些算法则是精度得不到很好的保证。因此,高阶有效的数值方法去处理气溶胶动力学方程对流项以及非线性凝并项成为学者们热衷研究的问题。由于气溶胶动力学方程是非线性的积分-微分方程,并且由于对流项和非线性凝并项的作用,在欧拉坐标下分布变化剧烈,呈对数正态分布,因此我们考虑特征线方法。由于偏微分方程中的对流项具有双曲特性,Douglas提出的修正特征线方法(MMOC)(Douglas et al.,1982)能够很好的求解对流扩散问题。特征线方法主要用于解决对流占优问题,它的主要思路是沿特征线方向由前一层网格的结果推出后一层网格的数值结果,这样就可以很好的解决让人难以接受的非物理震荡以及过多的数值弥散,从而磨平了尖锐的移动前沿。这种修正特征有限元方法被成功并且广泛的应用于多孔介质流问题的计算中如油藏数值模拟、海水侵入数值模拟等(Ewing et al.,1983; Russell et al.,1985).但是Douglas所提出的特征方法(Douglas et al.,1982)在时间方向只有一阶精度。为了使时间方向精度有所提高,Bermudez在他的论文中(Bermudez et al.,2006)关于线性的对流-扩散-反应问题提出了特征有限元方法,该方法在时间方向达到了二阶精度。因此提出新的改进的有效高阶特征算法,成为模拟求解气溶胶动力方程非常有意义的研究课题。其次,在第二个研究课题中,过去很多的文章研究并发展了气溶胶热力学平衡模型的建立,比如MARS (Saxena et al.,1986), SEQUILIB (Pilinis and Seinfeld,1987), SCAPE (Kim et al.,1993),以及ISORROPIA (Nenes et al.,1998,1999).以上所有的气溶胶热力学平衡模型都是基于热力学平衡方程,因此所有的计算这些平衡模型都是用了迭代算法。与其它模型所不同的,由于考虑了共同潮解度(mutual deliquescencehumidity) ISORROPIA方法(Nenes et al.,1998,1999;Makar et al.,2003;Metzger et al.,2002)被认为是使用最广泛并且预测气溶胶热力学平衡较为准确的一个模型方法。而该模型方法的不足之处与其它方法类似,在模型计算中使用了迭代技巧去解决非线性气/液气溶胶平衡方程。其中每一次迭代求解气溶胶热力学平衡方程都依赖于这次迭代时气溶胶的成分以及外部环境(如温度,相对湿度),而这些迭代计算将耗费预测气溶胶热力学平衡模型时的大量CPU时间。由于ISORROPIA模型的计算占用大量CPU时间,因此提出关于多态多相气溶胶的热平衡学输入输出算法在预测空气质量AQ (air quality)问题中显得尤为重要。高维数值模拟方法(high dimensional model representation (HDMR))是处理、评估、分析大量数据模型的一个有效方法,它的主要思想是捕捉高维输入以及输出变量之间的关系,文献(Rabitz et al.,1999; Li et al.,2003a,b,2004; Rabitz and Alis, 1999; Alis and Rabitz,2001)都对该方法进行了研究和扩展。HDMR方法在最近的研究成果中,被广泛应用于输入输出模型(input-output (10) system)用以减轻CPU负担提高计算效率减少计算时间中去。事实上,高维数值模拟HDMR方法的技巧类似于黑盒效应,通过实验数据或者观测数据建立高维数学模型中的输入与输出量之间的关系。而该技巧所得到的优势是,能够很好的避免严重的计算迭代并且能够大大的降低计算时间。HDMR方法主要包括随机抽样HDMR方法(random sampling HDMR (RS-HDMR)) (Rabitz et al.,1999; Li et al.,2003a,b)以及切割点HDMR方法(cut point HDMR (cut-HDMR)) (Rabitz and Alis,1999;Li et al.,2004).根据近几年的文献,HDMR的计算技巧已经被广泛应用于各种化学以及物理模拟实验计算中,成为一种解决多维问题的行之有效的办法。因此,利用改进后的HDMR方法模拟气溶胶热力学平衡模型成为一种新的尝试。在导师王文洽教授和梁栋教授的悉心指导下,本文作者关于气溶胶动力学和热力学平衡模型做了部分研究工作。在气溶胶动力学的方面,我们提出了有效的二阶特征有限元方法解决气溶胶动力方程。高阶的特征线方法被用于处理气溶胶的增长过程,同时二阶沿特征线方向的外推格式被用于解决气溶胶的非线性凝并过程。与标准的特征有限元方法比较而言,数值算例很好地证明了该方法在时间方向达到了二阶精度;在气溶胶体积方向(空间方向)达到了最优阶。在另一方面,我们提出了一种有效的高维方法预测气溶胶热力学平衡。结合移动切割点的技巧,我们利用HDMR方法来解决多态多相的高维气溶胶输入输出平衡模型。该方法能够在大范围的气溶胶浓度环境下(10-10(mol/m3)-10-6(mol/m3)),很好地模拟大气气溶胶平衡模型。数值算例显示这两种方法在不同的气溶胶预测领域中都有非常好的表现。全文共分叁章。在第一章中,我们主要介绍了气溶胶相关背景知识,包括气溶胶动力学以及气溶胶热力学平衡模型。首先给出气溶胶的简单定义以及在大气中发挥的主要作用。在1.2节中,我们给出了气溶胶动力方程的具体表达形式,其中每一项的气溶胶过程我们都作出了详细解释。在1.3节中,我们给出了气溶胶热力学平衡方程。同时,我们列出了气溶胶热力学平衡模型所涉及到的气溶胶成分输入量以及输出量,其中考虑了内陆地区和海洋地区两种情况。在第二章中,根据在模拟气溶胶矢径分布时所遇到的实际问题和需要,我们关于非线性积分-微分气溶胶动力学方程提出了有效的二阶特征有限元算法。我们所提出的二阶特征有限元方法分为两步,首先,将时间方向导数与对流项相结合,把它们转换为沿特征线方向的方向导数,然后沿特征线方向做差分,这样完成气溶胶动力方程半离散有限元格式。最后,在处理方程右边的非线性凝并项上,我们提出了沿特征线方向的二阶外推格式,其中当前层的数值解由沿特征线方向前两层的数值解所得到。结合上述的两种数值技巧以及有限元方法,我们提出了关于时间方向高阶的特征有限元方法。这种方法能够很好的适应气溶胶动力方程的特点,在数值算例中也显示出其特定的优势。由于这种方法的最大特点是将传统的特征有限元格式在时间方向的精度进行了提高,因此在模拟气溶胶矢径分布时可以在时间上采用大步长并且同时保证数值逼近的精度。同时该方法继承了特征线法的优点,最大程度的减少了过多的数值弥散以及很好的克服了在尖锐解处的非物理震荡。第二章给出的数值算例包括存在精确解的数值问题以及实际气溶胶矢径分布问题,与此同时我们也在对数坐标下给出了相应的数值算例。所有的数值算例都表明,我们所提出的二阶特征有限元方法能够很好的求解对数正态分布的气溶胶动力学方程。第二章中的内容已发表在高水平SCI杂志International Journal for Numerical Methods in Engineering中(2009年影响因子2.229)。在第叁章中,对于多态多相高维的气溶胶热力学平衡模型预测,我们提出了结合移动切割点技术的有效的高维数值模拟HDMR方法。在经典的cut-HDMR方法中,有且仅有一个切割点被应用。但是,当输入值的数据偏离这个切割点较远的时候,那么预测结果可能会出现比较大的误差。为了解决计算误差的问题,相应地,多切割点高维数值模拟方法(multicut-HDMR method (Li et al.,2004))在最近的研究成果中被提出来。但是,在一个非常庞大的高维输入定义系统中,如何去选取多切割点(multi-cut points)使得误差精度得到较大的提高成为一项困难的研究课题。在这一章中,我们结合了移动切割点,根据气溶胶热力学平衡物理以及化学特性,提出了移动切割点HDMR方法去模拟和预测气溶胶热力学平衡结果。这种方法的主要思路是结合多移动切割点解决高维的输入输出系统。与经典的cut-HDMR方法比较而言,我们所提出的方法能够很大程度上的提高数值模拟的误差精度。数值实验表明我们的方法不仅在计算误差上能够达到较为满意的效果,并且在CPU计算时间上也能得到非常大的改善。在与经典的气溶胶热力学平衡模型算法ISORROPIA的比较中,我们所用的计算时间要远远低于ISORROPIA所用的时间。对于实际的气溶胶热力学平衡预测,我们的方法也取得了比较好的效果.首先,我们在大气溶胶摩尔浓度范围内10-10到10-6 mol/m3,计算结果能够得到较好的保证;而经典的HDMR方法一般输入量最多涉及到两个量级的变化。进一步地,我们利用移动切割点HDMR方法去做了关于PM (particulate matter)浓度的预测,在与其它气溶胶热力学平衡模型的比较中比如MARS、SEQUILIB、EQUISOLV以及ISORROPIA,也取得了很好的效果。北京地区一天之中的气溶胶热力学平衡模型成分浓度预测也将在出现在我们的模拟实验中。第叁章中的内容已发表在国际顶尖环境SCI期刊Atmospheric Environment中(2009年影响因子2.89)。本文研究课题为国家基础研究973项目的子课题,课题编号2006CB403703。文中所涉及到部分气溶胶实际预测的数据由973项目子课题组提供。(本文来源于《山东大学》期刊2010-03-10)
李树,高丽洁,王体健,张美根[2](2005)在《利用耦合的气相化学模式和热力学平衡模式模拟硝酸盐气溶胶》一文中研究指出将热力学平衡模式ISORROPIA和气相化学模式相结合,通过敏感性试验,探讨了对流层硝酸盐形成的理化机制和影响因子。研究表明,硝酸盐浓度分别随太阳高度角增加、温度上升、湿度增加、SO2浓度、 NO2浓度和NH3浓度上升都有不同程度的增加,但关系是非线性的,它们是影响硝酸盐气溶胶浓度的主要因子。(本文来源于《大气气溶胶科学技术研究进展——第八届全国气溶胶会议暨第二届海峡两岸气溶胶科技研讨会文集》期刊2005-11-01)
气溶胶热力学平衡论文开题报告
(1)论文研究背景及目的
此处内容要求:
首先简单简介论文所研究问题的基本概念和背景,再而简单明了地指出论文所要研究解决的具体问题,并提出你的论文准备的观点或解决方法。
写法范例:
将热力学平衡模式ISORROPIA和气相化学模式相结合,通过敏感性试验,探讨了对流层硝酸盐形成的理化机制和影响因子。研究表明,硝酸盐浓度分别随太阳高度角增加、温度上升、湿度增加、SO2浓度、 NO2浓度和NH3浓度上升都有不同程度的增加,但关系是非线性的,它们是影响硝酸盐气溶胶浓度的主要因子。
(2)本文研究方法
调查法:该方法是有目的、有系统的搜集有关研究对象的具体信息。
观察法:用自己的感官和辅助工具直接观察研究对象从而得到有关信息。
实验法:通过主支变革、控制研究对象来发现与确认事物间的因果关系。
文献研究法:通过调查文献来获得资料,从而全面的、正确的了解掌握研究方法。
实证研究法:依据现有的科学理论和实践的需要提出设计。
定性分析法:对研究对象进行“质”的方面的研究,这个方法需要计算的数据较少。
定量分析法:通过具体的数字,使人们对研究对象的认识进一步精确化。
跨学科研究法:运用多学科的理论、方法和成果从整体上对某一课题进行研究。
功能分析法:这是社会科学用来分析社会现象的一种方法,从某一功能出发研究多个方面的影响。
模拟法:通过创设一个与原型相似的模型来间接研究原型某种特性的一种形容方法。
气溶胶热力学平衡论文参考文献
[1].程余.气溶胶动力学与热力学平衡预测的有效算法[D].山东大学.2010
[2].李树,高丽洁,王体健,张美根.利用耦合的气相化学模式和热力学平衡模式模拟硝酸盐气溶胶[C].大气气溶胶科学技术研究进展——第八届全国气溶胶会议暨第二届海峡两岸气溶胶科技研讨会文集.2005