《煤粉燃烧数学模型PPT学习教案.pptx》由会员分享，可在线阅读，更多相关《煤粉燃烧数学模型PPT学习教案.pptx（121页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、会计学1煤粉燃烧煤粉燃烧(rnsho)数学模型数学模型第一页，共121页。燃烧过程模化的一般燃烧过程模化的一般(ybn)研究研究n n在过去的30多年的时间里，人们对燃烧过程的模化予以很大的重视，计算机的发展更使这一发展成为可能，目前，已有商业化应用程序的出现。现有的模型方法的发展已开始走向实用。n n 燃烧过程的模化不仅受到计算机储存能力和运算速度的限制，同时又缺乏(quf)估计评价这些计算模型的基本数据。因此近10年的研究中，人们在开发燃烧通用商业化程序的同时，将重点放在其涉及的煤的燃烧过程的机理研究方面，使各种机理模型更加接近实际的应用要求。第1页/共121页第二页，共121页。发展发展
2、(fzhn)煤燃烧模型的关键煤燃烧模型的关键在于在于1.能否用目前广泛通过实验研究的单颗煤粒的特性数据预示出整个颗粒群的特性。2.有限的、稳态的试验数据是否可用于预示非稳态或准稳态过程的特性。3.煤燃烧过程的主要控制因素是各个环节的速率，只要对这些速率控制过程进行描述就可对整个过程作出描述。4.对简单的计算系统作了估价的数值方法(fngf)可用于复杂系统，尽管精度尚需提高但仍能得到有益的结果。5.由于程序的广泛性，不可能对各个子模型给出充分、完整的评价，但其总体参数的评价仍是可靠的。第2页/共121页第三页，共121页。n n这些基础前提有可能是不准确的，而其正确性又得不到直接证明，因此在整个
3、模型的应用中，必须小心地结合试验数据，运用数值试验方法来估价模型的精确性与可靠性。n n从另一个角度看，虽然燃烧过程包含着复杂的微观过程，但其宏观(hnggun)特性却呈现出明显的规律性，包括宏观(hnggun)的温度场、速度场、浓度场、传热、传质、流动等特性，这表明用数学方法来描述这种过程是可能的。第3页/共121页第四页，共121页。获得试验方法无法获得试验方法无法(wf)得到的信息包得到的信息包括括1)确定炉膛，燃烧器的总的基本特性；2)解释和进行测量结果的分析，以及进行优化数值试验（CAT方法(fngf)）；3)确定敏感的变量；4)确定速率控制过程；5)确定需要深入研究的领域；6)控制
4、燃烧过程；7)进行设计和优化；8)帮助进行比例放大的工程设计；9)提出合适的优化运行方式；10)控制污染物排放的可行措施。第4页/共121页第五页，共121页。1)1)求出燃烧室中的温度分布和壁面热流分布，分析其热工况；求出燃烧室中的温度分布和壁面热流分布，分析其热工况；2)2)求求解解(qi(qi ji)ji)气气相相流流场场，分分析析流流动动工工况况（各各股股气气流流的的混混合合、流流速速、湍流和回流情况）；湍流和回流情况）；3)3)求出颗粒的反应经历、分析积灰、结渣和磨损过程；求出颗粒的反应经历、分析积灰、结渣和磨损过程；4)4)求出颗粒的反应经历，合理组织气流流动；求出颗粒的反应经历，
5、合理组织气流流动；5)5)求解求解(qi ji)(qi ji)气相组分分布，分析合理的反应混合情况；气相组分分布，分析合理的反应混合情况；6)6)辅助燃烧器的设计；辅助燃烧器的设计；7)7)对低对低NOxNOx控制的指导；控制的指导；8)8)了解燃料变化对锅炉总体运行的影响；了解燃料变化对锅炉总体运行的影响；9)9)进进行行计计算算机机数数值值试试验验（CATCAT），估估计计实实验验结结果果，帮帮助助进进行行按按比比例例放大的工程设计及优化运行。放大的工程设计及优化运行。第5页/共121页第六页，共121页。第6页/共121页第七页，共121页。煤粉燃烧主要煤粉燃烧主要(zhyo)物理、化物
6、理、化学现象及其模型学现象及其模型第7页/共121页第八页，共121页。综合的煤粉燃烧过程综合的煤粉燃烧过程(guchng)的系统模型的基础有三个方面的系统模型的基础有三个方面1)1)气气相相湍湍流流运运动动的的研研究究方方法法，这这一一方方面面的的研研究究由由于于基基于于湍湍流流模模型型的的发发展展，人人们们对对各各种种场场合合应应用用双双方方程程或或多多方方程程的的湍湍流流模模型型，获获得得了了可可靠靠的的结结果果。这这方方面面研研究究的的成成果果为为进进一一步步研研究究整整个个燃燃烧烧过过程程打下基础。打下基础。2)2)湍湍流流和和其其他他物物理理现现象象的的相相互互作作用用，最最重重要
7、要的的发发展展是是颗颗粒粒湍湍流流扩扩散散的的随随机机方方法法的的出出现现，使使得得液液滴滴，颗颗粒粒流流的的扩扩散散问问题题获获得得了了令令人人满满意意的的结结果果。反反映映湍湍流流-气气相相反反应应的的湍湍流流燃燃烧烧模模型型近近年年也也得得到到了了相当的发展。相当的发展。3)3)煤煤粉粉燃燃烧烧动动力力学学研研究究和和煤煤粉粉燃燃烧烧形形态态学学方方面面的的发发展展，为为煤煤粉粉颗颗粒粒本本身身的的模模化化打打下下了了基基础础。如如煤煤的的加加热热(ji(ji r)r)，热热解解（挥挥发发物物析析出出），焦炭反应过程中的形态，结团、破碎、膨胀、收缩等研究。，焦炭反应过程中的形态，结团、破
8、碎、膨胀、收缩等研究。第8页/共121页第九页，共121页。单颗煤粒经历单颗煤粒经历(jngl)模型模型n n煤粒的加热n n水份蒸发模型(mxng)n n挥发份析出模型(mxng)n n焦炭的非均相反应模型(mxng)n n煤粒在燃烧室中的其它经历模型(mxng)n n气体颗粒传质n n煤粉颗粒在反应中的物理变化第9页/共121页第十页，共121页。煤粒的加热煤粒的加热(ji r)第10页/共121页第十一页，共121页。水份蒸发水份蒸发(zhngf)模型模型第11页/共121页第十二页，共121页。挥发份析出挥发份析出(xch)模型模型第12页/共121页第十三页，共121页。焦炭焦炭(j
9、iotn)的非均相反应模的非均相反应模型型第13页/共121页第十四页，共121页。煤粒在燃烧室中的其它煤粒在燃烧室中的其它(qt)经历模型经历模型n n气体(qt)颗粒传质第14页/共121页第十五页，共121页。煤粉颗粒煤粉颗粒(kl)在反应在反应中的物理变化中的物理变化第15页/共121页第十六页，共121页。煤燃烧过程中流动煤燃烧过程中流动(lidng)、气相反应过程及其模型气相反应过程及其模型n n连续性方程是流体力学中质量守恒的表达式6。对于任何一个(y)化学组份K，其组份连续性方程第16页/共121页第十七页，共121页。动量动量(dngling)方程方程第17页/共121页第十
10、八页，共121页。能量能量(nngling)方程方程第18页/共121页第十九页，共121页。湍流湍流(tunli)模型模型n n时均湍流运动(yndng)方程的导出第19页/共121页第二十页，共121页。时均动量时均动量(dngling)方程方程第20页/共121页第二十一页，共121页。时均能量时均能量(nngling)方程方程第21页/共121页第二十二页，共121页。时均方程时均方程(fngchng)可用以可用以 为为标量参数的统一形式标量参数的统一形式第22页/共121页第二十三页，共121页。混合混合(hnh)长度模型长度模型n nBoussinesq建议(jiny)把湍流剪切应
11、力表示为第23页/共121页第二十四页，共121页。第24页/共121页第二十五页，共121页。单微分方程单微分方程(wi fn fn chn)模型模型第25页/共121页第二十六页，共121页。第26页/共121页第二十七页，共121页。k-双方双方(shungfng)程模型程模型第27页/共121页第二十八页，共121页。雷诺应力雷诺应力(yngl)模型模型第28页/共121页第二十九页，共121页。代数应力代数应力(yngl)模型模型第29页/共121页第三十页，共121页。气相燃烧气相燃烧(rnsho)第30页/共121页第三十一页，共121页。煤粉燃烧时炉内传热煤粉燃烧时炉内传热(c
12、hun r)的的模型及计算模型及计算第31页/共121页第三十二页，共121页。燃烧室的热辐射燃烧室的热辐射n n火焰(huyn)热辐射性质第32页/共121页第三十三页，共121页。辐射传热辐射传热(chun r)的模型的模型n n热流(rli)法(Heat flux)n n区域分析法(Zone Analysis)n n蒙特卡洛法(Monte-Carlo)第33页/共121页第三十四页，共121页。热流热流(rli)法法n n热流法的特点是将复杂的不均匀(jnyn)的多向的界面辐射热流和用均匀(jnyn)的界面辐射热流来代替，并取其平均值第34页/共121页第三十五页，共121页。区域区域(
13、qy)法法n n 最早由Hottel等提出的区域(qy)法，将整个燃烧室划分成若干个区域(qy)，把壁面划分成面积区，假定在区域(qy)内的温度和物性都是均匀的，按照各区域(qy)直接和周围进行空间辐射换热来计算。n n由于区域(qy)法考虑了一个区与其他所有区域(qy)的辐射换热关系，从原理上讲显然有其正确性，特别是当将区域(qy)划分足够小的时候。n n但是区域(qy)法带来的是一个巨大的工作量，因为需要求解各区域(qy)间的辐射交换面积，等于是直接求解多个重积分，占用了大量时间，所以对于只有可能分成少数区域(qy)的时候，才有较好的实用意义。第35页/共121页第三十六页，共121页。M
14、onte-Carlo概率概率(gil)模拟法模拟法n n热热流流法法忽忽略略了了空空间间运运算算而而获获得得方方程程的的简简化化，但但降降低低了了计计算算精精度，而区域法则直接计算积分度，而区域法则直接计算积分(jfn)(jfn)，工作量大。，工作量大。n n概概率率模模拟拟法法是是建建立立在在离离散散发发射射法法的的基基础础上上,所所谓谓的的离离散散发发射射法法即即是是把把每每一一个个微微元元体体向向周周围围发发射射的的辐辐射射能能量量按按空空间间角角分分为为若若干干等等分分(如如NiNi个个),),则则每每一一个个空空间间角角大大小小为为4 4/Ni/Ni，将将全全部部燃燃烧烧室室划划分分
15、为为MxMx个个体体积积微微元元和和MsMs个个面面积积微微元元(壁壁面面),),则则就就可可以以在在体体积积和和面面积积上上对对所所有有的的微微元元发发射射出出的的辐辐射射能能量量进进行行“吸吸收收”,”,从从而而获获得得每每个个微微元元的的能能量量方方程程辐辐射射换换热热项项，概概率率模模拟拟法法在在上上面面的的基基础础上上,在在所所有有方方向向上上“发发射射”并并进进行行跟跟踪踪，当当NiNi取取得得一一定定大大，则则模模拟拟可可以以相相当当的的正正确确,它它的的特特点点是是具具灵灵活活性性,通通用用性性,准准确确性性(可可以以随随NN增增加加而而与与区区域域法法相相匹匹美美),),而而
16、且且对对于于复复杂杂几几何何形形状状的的燃燃烧烧室室适适应应性性也也很很强强,因因此此越越来来越越受受到到重重视视,缺缺点点是是为为了了获获得得足足够够的的精精度度,仍仍需需要要大大量量的的计计算算次次数数,而而经经济性不如通量法。济性不如通量法。第36页/共121页第三十七页，共121页。实际煤火焰辐射传热实际煤火焰辐射传热(chun r)模模拟结果及分析拟结果及分析n n二维煤浆火炬辐射传热(chun r)计算第37页/共121页第三十八页，共121页。第38页/共121页第三十九页，共121页。第39页/共121页第四十页，共121页。第40页/共121页第四十一页，共121页。第41页
17、/共121页第四十二页，共121页。第42页/共121页第四十三页，共121页。第43页/共121页第四十四页，共121页。第44页/共121页第四十五页，共121页。三维煤粉火焰辐射传热计算三维煤粉火焰辐射传热计算(j sun)示例示例第45页/共121页第四十六页，共121页。第46页/共121页第四十七页，共121页。第47页/共121页第四十八页，共121页。第48页/共121页第四十九页，共121页。标准标准(biozhn)工况下工况下炉内温度场的数值模拟炉内温度场的数值模拟第49页/共121页第五十页，共121页。第50页/共121页第五十一页，共121页。第51页/共121页第五
18、十二页，共121页。第52页/共121页第五十三页，共121页。第53页/共121页第五十四页，共121页。第54页/共121页第五十五页，共121页。第55页/共121页第五十六页，共121页。第56页/共121页第五十七页，共121页。第57页/共121页第五十八页，共121页。二次风反切二次风反切(fnqi)时炉时炉内温度场的数值模拟内温度场的数值模拟第58页/共121页第五十九页，共121页。第59页/共121页第六十页，共121页。第60页/共121页第六十一页，共121页。第61页/共121页第六十二页，共121页。第62页/共121页第六十三页，共121页。一次风反切时炉内温度场
19、一次风反切时炉内温度场的数值的数值(shz)模拟模拟第63页/共121页第六十四页，共121页。第64页/共121页第六十五页，共121页。第65页/共121页第六十六页，共121页。第66页/共121页第六十七页，共121页。第67页/共121页第六十八页，共121页。煤粉颗粒煤粉颗粒(kl)扩散及两相流模扩散及两相流模型型n n描述两相流动(lidng)的基本方法n n一类是通过拉格朗日方法，一类是欧拉方法第68页/共121页第六十九页，共121页。稀相流动稀相流动(lidng)的基本分析的基本分析 容积比V颗粒/V汽1%；颗粒碰撞不频繁；单个颗粒的传热、传质、阻力系数不直接受邻近(ln
20、jn)颗粒影响。第69页/共121页第七十页，共121页。n n最早进行气-固两相问题的完整(wnzhng)尝试可以追溯到Migdal和Agoata8521967年给出的方法，拉格朗日方法提出的考虑两相耦合问题的方法主要是PSIC法，PSIC的主要思想是：n n首先解无颗粒存在的气相场，在此场中求解各组颗粒的轨迹，沿轨迹求颗粒、尺寸、速度、温度等经历变化，在气相单元边界上记录质量，动量和能量、组分、提供气相求解的源项，再解气相场，反复直到收敛。第70页/共121页第七十一页，共121页。第71页/共121页第七十二页，共121页。颗粒在气流颗粒在气流(qli)中的受力分中的受力分析析即是颗粒受
21、力的总和(zngh)，一般包括：稳态气动阻力，浮力，虚假质量效应，Basset力，Magnus效应，Saffman升力和热泳力等，在颗粒反应过程中由于表面质量流动不均匀引起的推力作用第72页/共121页第七十三页，共121页。颗粒的湍流颗粒的湍流(tunli)扩散扩散 离散相本身(bnshn)的湍流输运(湍流扩散)；由于离散相存在的输运特性对连续相湍流特性的修正；湍流脉动引起的相间传输率特性的修正。第73页/共121页第七十四页，共121页。两相模型两相模型(mxng)第74页/共121页第七十五页，共121页。数值求解数值求解(qi ji)方法方法各项相应称为时间导数(do sh)项、对流项
22、、扩散项和源项第75页/共121页第七十六页，共121页。n n.形式相同。这一点使我们有可能发展建立求解的通用程序。n n.非线性。即存在因变量或它们的导数的非一次项。它主要(zhyo)表现在对流项和化学反应源项。其它源项和扩散项也可能包含非线性，例如湍流脉动动能输运方程的源项，考虑交换系数为组分浓度和温度函数时的扩散项。在变密度的过程中，动量、能量和组分方程中的时间导数项也具有非线性。这表明基本方程一般来说无法用解析方法而必须用数值方法求解。n n.耦合。构成方程组的各方程不互相独立，因变量交错地存在于各个方程中。因此不能直接求解，而必须选代求解。第76页/共121页第七十七页，共121页
23、。离散离散(lsn)化方法化方法n n一个微分方程的数值解系一组的空间和时间分布的离散的数的集合。这与在实验室中进行的实验相类似。数值分析方法和实验方法两者所获得的数据都只能是一些有限数量的数值；就实际需要(xyo)而言，这个数量的数值是够实用就可以了。n n所谓的数值方法就是把计算域内有限数量位置(叫做网格结点)上的因变量当作为基本的未知量来处理。该方法的任务是提供一组关于这些未知量的代数方程并规定求解这组方程的算法。第77页/共121页第七十八页，共121页。离散离散(lsn)化的概念化的概念n n用离散的值取代了包含在微分方程精确解中的连续信息，就将注意力集中在网格结点处的值。这实际上已
24、经离散了的分布，这一类数值方法叫做离散化方法。n n这样，就从另一个(y)意义上遇到了离散化的概念。连续的计算域也被离散开了。这种对空间和因变量所作的系统的离散化使得有可能用比较容易求解的简单的代数方程取代前面提到过的控制微分方程。第78页/共121页第七十九页，共121页。离散化方程离散化方程(fngchng)的结构的结构n n一个离散化方程是连接一组网格结点处值的代数关系式。这样的关系式由的微分方程导出，并表示与该微分方程相同的物理信息。一定的离散化方程只与少数(shosh)的几个网格结点有关，这种情况是选取分段分布这一特性的结果。n n对于一个已知的微分方程，可能的离散化方程决不是唯一的
25、；尽管在网格结点数非常大的极限条件下，预计所有这些可能类型的离散化方程将会给出相同的解。离散化方程的不同形式起因于分布假设以及推导方法的不同。第79页/共121页第八十页，共121页。积分积分(jfn)区域的网区域的网格化格化n n积分区域即待求函数定义域的网格化是微分方程离散化的基础。网格化的方式影响微分方程离散化的难易，也关系到解的精确性、收敛性和求解的经济性。n n在一个积分区域中，网格的形式可以是单一的，也可以是几种形式网格的组合(zh)。例如，在计算内燃机气缸内的气体状态和活塞内的温度分布时，可以采用组合(zh)网格。在活塞不能进入的区域A取固定网格，在容积不断随活塞位置变化的区域B
26、取拉伸式网格，在活塞内部区域C取运动但相对位置固定的网格。第80页/共121页第八十一页，共121页。第81页/共121页第八十二页，共121页。第82页/共121页第八十三页，共121页。贴体坐标概念贴体坐标概念(ginin)初步初步n n对于非矩形域的计算问题，用矩形网格就无法(wf)获得正确的应用，这是有限差分方法的主要缺点。n n对于具有复杂形状的区域，可以采用曲线坐标系，并且使物面与坐标系内某一坐标线吻合。这种坐标系可以是正交的(如上述二例)，也可以是非正交的。网格的分布可以是均匀的，也可以是非均匀的，视问题的需要而定。此外网格还可以设计为自动调节的，即所谓自适应网格。第83页/共1
27、21页第八十四页，共121页。第84页/共121页第八十五页，共121页。第85页/共121页第八十六页，共121页。第86页/共121页第八十七页，共121页。有限有限(yuxin)差分方程的一般差分方程的一般形式形式n n以二维网格系统内的差分方程(fngchng)为例说明第87页/共121页第八十八页，共121页。第88页/共121页第八十九页，共121页。n n为了能由差分方程求出P点的函数值P，直观上看必须具备两个条件：一是求出所有各项的系数值和常数项，二是知道邻近各网格(wn)点上的值(包括时间坐标的网格(wn)。前一个问题的解决主要依靠把待解的微分方程离散化，后一个问题是靠人为地
28、给定在积分区域中的初始值。注意这里解差分方程需要的初始值与解微分方程需要的初始条件并非一回事。第89页/共121页第九十页，共121页。差分差分(ch fn)方程建立方程建立的方法的方法n n泰勒级数公式(gngsh)n n恰好在第三项之后截断级数，将二个方程相加及相减，我们得到n n ,第90页/共121页第九十一页，共121页。n n这个方法含有这样的假设(jish)：的变化多少有点象是x的一个多项式，从而高阶导数是不那么重要的。但是当存在指数形式的变化时，这种假设(jish)就可能导致人们不希望要的那种公式。泰勒级数公式的推导是比较直截了当的，但是缺乏弹性并且其中各项的物理意义难以理解*
29、。第91页/共121页第九十二页，共121页。第92页/共121页第九十三页，共121页。变分公式变分公式(gngsh)n n变分法证明：求解某些微分方程的问题等效于使一称之为泛函的相关量最小化。这种等效关系就是所谓的变分原理。如果相关因变量的网格点值使泛函最小，那么所得到的条件即给出所需要的离散化方程。变分公式非常普遍地用于应力(yngl)分析的有限元法，这时可以把变分公式与虚功原理联系起来。除了代数和概念上的复杂性之外，变分公式的主要缺点在于它的适用范围有限，因为我们所感兴趣的所有的微分方程都不存在变分原理。第93页/共121页第九十四页，共121页。控制控制(kngzh)容积公式容积公式
30、n n控制容积公式可以看成是加权余数(ysh)法的一种特殊形式。控制容积公式的基本思想是易于理解的，并且适于进行直接的物理解释。把计算域分成许多互不重叠的控制容积，并使每一个网格结点都由一个控制容积所包围。对每一个控制容积积分微分方程。应用表示网格结点之间变化的分段分布关系来计算所要求的积分。这样做的结果，我们就得到了一个包含有一组网格结点处的值的离散化方程。第94页/共121页第九十五页，共121页。n n控制容积公式的最吸引人的特征(tzhng)是：所得到的结果(解)将意味着任何一组的控制容积内，当然也就是在整个计算域内，诸如质量、动量以及能量这样一些物理量的积分守恒都可以精确地得到满足。
31、对于任意数目的网格结点，这一特征(tzhng)都存在(不只是限于网格结点数变得很大时的极限意义上)。因而，即便是粗网格的解也照样显示准确的积分平衡.第95页/共121页第九十六页，共121页。交错网格交错网格(wn)系统系统n n在这种交错网格系统中，除速度以外的各个变量，例如压力(yl)、密度、温度、组分浓度、比热、交换系数等均采用同一个网格系统；三个速度分量分别采用它独用的网格系统。如果压力(yl)和温度等变量的网格点为P，那么速度分量的网格点将位于网格P的与x方向垂直的界面的中心，v分量的网格点位于网格P的与y方向垂直的界面的中心，w分量的网格点位于网格P的与z方向垂直的界面的中心。因此
32、，在同一个三维积分区域内就有四组不同的网格交错排列。第96页/共121页第九十七页，共121页。第97页/共121页第九十八页，共121页。第98页/共121页第九十九页，共121页。第99页/共121页第一百页，共121页。第100页/共121页第一百零一页，共121页。差分方程差分方程(fngchng)的求解的求解n n三对角线矩阵法n n SIMPLE方法n n.估计速度分布；n n.用式(9-214)计算准速度分布；n n.计算压力方程(9-216)的系数，并解出压力分布；n n.把解出的压力分布作为p*，解动量方程得到u*和v*；n n.建立并求解(qi ji)压力校正方程，得到p；
33、n n.求出速度校正值u和v，校正速度分布，但不用p校正压力分布；n n.如果需要，求解(qi ji)其它变量的方程；n n.重复上述第2至第7步，直至收敛。第101页/共121页第一百零二页，共121页。煤粉火焰综合煤粉火焰综合(zngh)求解及示例求解及示例n n求解程序n n在给定一个合理的，收敛的流场，调用颗粒两相流的计算，首先只考虑动量藕合，而不考虑燃烧和反应，这个藕合通常是很微弱的，但是若在入口处煤与气之间存在着很大的速度差，则在流动与两相流动之间迭代几次，再加上颗粒反应，这时再调用气相湍流燃烧、修正气体温度、化学组成，混合物分子量及流体质量密度。在开始使用气相湍流燃烧程序时，可以
34、根据进气是预混的还是未混合的来选择(xunz)两种不同的方法。迭代直接开始去计算需要的平衡态解或化学动力学解。这样迭代计算，如果收敛且计算结果合理，则可将辐射换热加入到循环之中，最后，要利用各个模型的最精确的公式，得到满足全部方程及边界条件的最终稳定解。第102页/共121页第一百零三页，共121页。第103页/共121页第一百零四页，共121页。模型模型(mxng)计算结果示例计算结果示例n n杨伯翰大学(dxu)一维模型计算结果第104页/共121页第一百零五页，共121页。第105页/共121页第一百零六页，共121页。第106页/共121页第一百零七页，共121页。二维模型二维模型(m
35、xng)计算及结果计算及结果第107页/共121页第一百零八页，共121页。第108页/共121页第一百零九页，共121页。第109页/共121页第一百一十页，共121页。第110页/共121页第一百一十一页，共121页。第111页/共121页第一百一十二页，共121页。第112页/共121页第一百一十三页，共121页。三维计算结果实例三维计算结果实例(shl)n n颗粒运动的轨迹n n水分蒸发规律n n挥发份析出(xch)规律和着火时间、距离n n焦炭燃烬规律第113页/共121页第一百一十四页，共121页。第114页/共121页第一百一十五页，共121页。第115页/共121页第一百一十六页，共121页。第116页/共121页第一百一十七页，共121页。第117页/共121页第一百一十八页，共121页。第118页/共121页第一百一十九页，共121页。第119页/共121页第一百二十页，共121页。第120页/共121页第一百二十一页，共121页。