fluent物质输送与有限速率化学反应.doc
《fluent物质输送与有限速率化学反应.doc》由会员分享,可在线阅读,更多相关《fluent物质输送与有限速率化学反应.doc(30页珍藏版)》请在咨信网上搜索。
1、第十三章 物质输送和有限速率化学反应FLUENT可以通过求解描述每种组成物质的对流、扩散和反应源的守恒方程来模拟混合和输运,可以模拟多种同时发生的化学反应,反应可以是发生在大量相(容积反应)中,和/或是壁面、微粒的表面。包括反应或不包括反应的物质输运模拟能力,以及当使用这一模型时的输入将在本章中叙述。注意你可能还希望使用混合物成分的方法(对非预混系统,在14章介绍)、反应进程变量的方法(对预混系统,在15章介绍),或部分预混方法(在16章介绍)来模拟你的反应系统。见12章FLUENT中反应模拟方法的概述。本章中的分为以下章节:l 13.1 容积反应l 13.2 壁面表面反应和化学蒸汽沉积l 1
2、3.3 微粒表面反应l 13.4 无反应物质输运131 容积反应与容积反应有关的物质输运和有限速率化学反应方面的信息在以下小节中给出:l 13.1.1 理论l 13.1.2 模拟物质输运和反应的用户输入概述l 13.1.3 使能物质输运和反应,并选择混合物材料l 13.1.4 混合物和构成物质的属性定义l 13.1.5 定义物质的边界条件l 13.1.6 定义化学物质的其他源项l 13.1.7 化学混合和有限速率化学反应的求解过程l 13.1.8 物质计算的后处理l 13.1.9 从CHEMKIN导入一个化学反应机理13.1.1 理论物质输运方程当你选择解化学物质的守恒方程时,FLUENT通过
3、第种物质的对流扩散方程预估每种物质的质量分数,Yi。守恒方程采用以下的通用形式: (13.1-1)其中是化学反应的净产生速率(在本节稍后解释),为离散相及用户定义的源项导致的额外产生速率。在系统中出现N种物质时,需要解N-1个这种形式的方程。由于质量分数的和必须为1,第N种物质的分数通过1减去N-1个已解得的质量分数得到。为了使数值误差最小,第N种物质必须选择质量分数最大的物质,比如氧化物是空气时的N2。层流中的质量扩散在方程13.1-1中,是物质的扩散通量,由浓度梯度产生。缺省时,FLUENT使用稀释近似,这样扩散通量可记为: (13.1.2)这里是混合物中第种物质的扩散系数。对于确定的层流
4、流动,稀释近似可能是不能接受的,需要完整的多组分扩散。在这些例子中,可以解Maxwell-Stefan方程,详细情况见7.7.2节。湍流中的质量扩散在湍流中,FLUENT以如下形式计算质量扩散: (13.1.3)其中是湍流施密特数,(缺省设置值为0.7)。注意,湍流扩散一般淹没层流扩散,在湍流中指定详细的层流性质是不允许的。能量方程中的物质输送处理在许多多组分混合流动中,物质扩散导致了焓的传递。这种扩散对于焓场有重要影响,不能被忽略。特别是,当所有物质的Lewis数 (13.1-4)远离1时,忽略这一项会导致严重的误差。FLUENT缺省地包含这一项。在方程13.1-4中,k为热导率。进口处的扩
5、散在FLUENT的非耦合求解器中,入口的物质净输送量由对流量和扩散量组成,对耦合解算器,只包括对流部分。对流部分由你指定的物质浓度确定。扩散部分依赖于计算得到的物质浓度场。因此,扩散部分(从而使净输送量)不预先指定。见13.1.5节有关指定入口净输送量的信息。反应建模的一般有限速率形式反应速率作为源项在方程13.1-1中出现,在FLUENT中根据以下三种模型中的一个计算:层流有限速率模型:忽略湍流脉动的影响,反应速率根据Arrhenius公式确定。涡耗散模型:认为反应速率由湍流控制,因此避开了代价高昂的Arrhenius化学动力学计算。涡耗散概念(EDC)模型:细致的Arrhenius化学动力
6、学在湍流火焰中合并。注意详尽的化学动力学计算代价高昂。通用有限速率对于范围很广的应用,包括层流或湍流反应系统,预混、非预混、部分预混燃烧系统都适用。层流有限速率模型层流有限速率模型使用Arrhenius公式计算化学源项,忽略湍流脉动的影响。这一模型对于层流火焰是准确的,但在湍流火焰中Arrhenius化学动力学的高度非线性,这一模型一般不精确。对于化学反应相对缓慢、湍流脉动较小的燃烧,如超音速火焰可能是可以接受的。化学物质的化学反应净源项通过有其参加的NR个化学反应的Arrhenius反应源的和计算得到。其中是第种物质的分子量,为第种物质在第r个反应中的产生/分解速率。反应可能发生在连续相反应
7、的连续相之间,或是在表面沉积的壁面处,或是发生在一种连续相物质的演化中。考虑以如下形式写出的第r个反应: (13.1-6)其中系统中化学物质数目;反应r中反应物i的化学计量系数;反应r中生成物i的化学计量系数;第i种物质的符号;反应r的正向速率常数;反应r的逆向速率常数;方程13.1-6对于可逆和不可逆反应(FLUENT中缺省为不可逆)都适用。对于不可逆反应,逆向速率常数简单地被忽略。方程13.1-6中的和是针对系统中的所有物质,但只有作为反应物或生成物出现的物质才有非零的化学计量系数。因此,不涉及到的物质将从方程中清除。反应r中物质i的产生/分解摩尔速度以如下公式给出: (13.1-7)其中
8、:反应r 的化学物质数目;反应r中每种反应物或生成物j的摩尔浓度;反应r中每种反应物或生成物j的正向反应速度指数;反应r中每种反应物或生成物j的逆向反应速度指数;见13.1.4节有关输入整体正向反应(不可逆)和单元反应(可逆)的化学计量系数和速率指数方面的内容。表示第三体对反应速率的净影响。这一项由下式给出:其中为第r个反应中第j种物质的第三体影响。在缺省状态,FLUENT在反应速率计算中不包括第三体影响。但是当你有它们的数据时,你可以选择包括第三体影响。反应r的前向速率常数通过Arrhenius公式计算:其中,指数前因子(恒定单位); 温度指数(无量纲);反应活化能();气体常数()你(或者
9、数据库)可以在FLUENT的问题定义中提供,并可选择提供。如果反应是可逆的,逆向反应常数可以根据以下关系从正向反应常数计算:其中为平衡常数,从下式计算:其中表示大气压力(101325Pa)。指数函数中的项表示Gibbs自由能的变化,其各部分按下式计算:其中和是标准状态的熵和标准状态的焓(生成热)。这些值在FLUENT中作为混合物材料的属性指定。压力独立反应FLUENT可以用以下三种方法之一来表示压力独立反应(或压力下降)反应的速率表达式。“压力下降”反应是发生在Arrhenius 高压和低压限制之间的反应,因而不仅仅依赖于温度。有三种方法表示在“fall-off”区域的速率表达式,最简单的是L
10、indemann140形式。还有其它良种相关的方法,Troe方法77和SRI方法230,它们提供了更精确的”fall-off”区域表达形式。Arrhenius速率参数对于高压和低压限制都是需要的。两个限制的速率系数融合以产生光滑的压力独立表达式。在Arrhenius形式中,高压限制和低压限制的参数如下:在任意压力下,净反应速率常数为:其中定义为:为溶液气体的浓度,可以包括第三体效率。如果方程13.1-16函数F为1,则是Lindemann形式。FLUENT提供了两种其他形式来表述F,称为Troe方法和SRI方法。在Troe方法中,F按下式给出:其中,参数做为输入确定。在SRI方法中,缝合函数F
11、近似为除了低压限制表达式中的三个Arrhenius参数以外,你还需要提供F表达式中的a, b, c, d, e。!化学动力学机理中有很高的非线性并且形成了一组强烈耦合的方程。求解过程指导见13.1.7节。如果你有一个CHEMKIN形式的化学反应机理112,你可以将这一机理导入FLUENT,如13.1.9节。涡耗散模型大部分燃料快速燃烧。整体反应速率由湍流混合控制。在非预混火焰中,湍流缓慢地通过对流/混合燃料和氧化剂进入反应区,在反应区它们快速地燃烧。在预混火焰中,湍流对流/混合冷的反应物和热的生成物进入反应区,在反应区迅速地发生反应。在这些情况下,燃烧称为混合限制的,复杂,常常是未知的化学反应
12、动力学速率可以安全地忽略掉。FLUENT提供了湍流-化学反应相互作用模型,基于Magnussen 和Hjertager149的工作,称为涡耗散模型。反应r中物质的产生速率由下面两个表达式中较小的一个给出:在方程13.1-25和13.1-26中,化学反应速率由大涡混合时间尺度控制,如同Splading227的涡破碎模型一样。只要湍流出现(),燃烧即可进行,不需要点火源来启动燃烧。这通常对于非预混火焰是可接受的,但在预混火焰中,反应物一进入计算区域(火焰稳定器上游)就开始燃烧。为了修正这一点,FLUENT提供了有限速率/涡耗散模型,其中Arrhenius(方程13.1-7)和涡耗散(方程13.1-
13、25和13.1-26)反应速率都进行计算。净反应速率取两个速率中较小的。实际上,Arrhenius反应速率作为一种动力学开关,阻止反应在火焰稳定器之前发生。一旦火焰被点燃,涡耗散速率通常会小于Arrhenius反应速率,并且反应是混合限制的。!尽管FLUENT允许采用涡耗散模型和有限速率/涡耗散模型的多步反应机理(反应数2),但可能会产生不正确的结果。原因是多步反应机理基于Arrhenius速率,每个反应的都不一样。在涡耗散模型中,每个反应都有同样的湍流速率,因而模型只能用于单步(反应物产物)或是双步(反应物中间产物,中间产物产物)整体反应。模型不能预测化学动力学控制的物质,如活性物质。为合并
14、湍流流动中的多步化学动力学机理,使用EDC模型(下面介绍)。!涡耗散模型需要产物来启动反应(见方程13.1-26)。当你初始化求解的时候,FLUENT设置产物的质量比率为0.01,通常足够启动反应。但是,如果你首先聚合一个混合解,其中所有的产物质量比率都为0,你可能必须在反应区域中补入产物以启动反应。详细内容见13.1.7节。LES的涡耗散模型当使用LES湍流模型时,湍流混合速率(方程13.1-25和13.1-26中的)被亚网格尺度混合速率替代。计算为:涡-耗散-概念(EDC)模型涡-耗散-概念(EDC)模型是涡耗散模型的扩展,以在湍流流动中包括详细的化学反应机理148。它假定反应发生在小的湍
15、流结构中,称为良好尺度。良好尺度的容积比率按下式模拟80:其中*表示良好尺度数量,容积比率常数=2.1377;运动粘度认为物质在好的结构中,经过一个时间尺度后开始反应。其中为时间尺度常数,等于0.4082在FLUENT中,良好尺度中的燃烧视为发生在定压反应器中,初始条件取为单元中当前的物质和温度。反应经过时间尺度后开始进行,由方程13.1-7的Arrhenius速率控制,并且用普通微分方程求解器CVODE进行数值积分45。经过一个时间的反应后物质状态记为物质的守恒方程13.1-1中的源项计算公式为:EDC模型能在湍流反应流动中合并详细的化学反应机理。但是,典型的机理具有不同的刚性,它们的数值积
16、分计算开销很大。因而,只有在快速化学反应假定无效的情况下才能使用这一模型,例如在快速熄灭火焰中缓慢的CO烧尽、在选择性非催化还原中的NO转化。推荐使用双精度求解器以避免刚性机理中固有的大指数前因子和活化能产生的舍入误差。见13.1.7节获得使用EDC模型求解的指导。13.1.2 物质输送和反应模拟输入概览设定涉及物质输送和反应问题的基本步骤如下,每一步的详细执行过程见13.1.3-13.1.5节。有关设定和求解的附加信息在13.1.6-13.1.8节中提供。1 选定物质输送和容积反应,指定混合物材料。见13.1.3节(混合物材料概念在下面解释)。2 如果你还要模拟壁面或微粒表面反应,则要打开壁
17、面和/或微粒表面反应。细节见13.2和13.3节。3 检查和/或定义混合物的属性。(见13.1.4节)。混合物属性包括:l 混合物中的物质l 反应l 其他物理属性(如粘度、比热)4 检查和/或设置混合物中单个物质的属性(见13.1.4节)5 设置物质边界条件(见13.1.5节)在很多情况下,当你选择混合物材料是,求解器从材料数据库中得到物质性质、反应等,因而你将不需要修改任何物理属性。但有一些性质可能在数据库中没有定义。如果有任何性质需要设置时,你将被警告,这样你可以指定这些性质的适当值。你还可能希望检查数据库中这些性质的值,以确定它们对你的应用是否正确。修改已存在的混合物材料或从最开始创建一
18、个新的材料的详细内容见13.1.4节。混合物材料的修改包括以下方面:l 物质的添加和删除l 改变化学反应l 修改混合物的其他材料属性l 修改混合物本构物质的材料属性如果你在求解一个反应流问题,你常常希望将混合物的比热定义为组成的函数,将每种物质的比热定义为温度的函数。你还可能对其它一些性质希望也做这样的定义。缺省状态下,将使用恒定的属性,但对一些物质的性质,在数据库中存在一个温度的分段多项式函数可供你使用。如果你知道更多适合于你的问题的函数的话,你还可以选择指定一个不同的温度依赖函数。混合物材料在FLUENT中提出混合物材料的概念以方便物质输送和反应流动的设置。混合物材料可以认为是一组物质和一
19、列控制它们相互作用的规律。混合物材料带有以下性质:l 一列本构物质,相对于“流体”物质l 一列混合定律,指示如果希望得到组分依赖的属性,混合属性(密度,粘度,比热等)如何从单个物质得到l 如果希望属性不依赖组分,直接指定混合物属性l 其它与耽搁物质无关的材料属性(如吸收和辐射系数)l 一组反应,包括反应类型(有限速率,涡耗散等)和化学计量和速率常数混合物材料和流体材料都储存在FLUENT的材料数据库中。包括许多常见的混合物材料(如甲烷-空气,丙烷-空气)。通常,在数据库中定义了一步/两步反应机理和大量混合物及其构成物质的属性。当你指定了你希望使用哪种混合物材料后,适当的混合物材料,流体材料和属
20、性将被装载到求解器中。如果缺少任何所选材料(或构成流体材料)必须的属性,求解器将通知你需要指定它。另外,你可以选择修改任何预定义的属性。见7.1.2节了解有关FLUENT数据库属性数据源的信息。例如,如果你计划模拟一种甲烷-空气的燃烧,你不需要明确指定反应中涉及的物质和反应本身。只需要简单地选择甲烷-空气作为使用的混合物材料,相关的物质(CH4,O2,CO2,H2O和N2)和反应数据将从数据库装入求解器。然后你可以检查物质、反应和其它属性并定义其它任何缺少的属性,和/或修改任何你希望使用不同值或函数的属性。通常你希望定义一个与组分、温度相关的比热,还可能希望将其它属性定义为温度和/或组分的函数
21、。混合物材料的使用给你提供了一种灵活性,可以使用大量预定义混合物中的一种,修改这些混合物,或是创建你自己的混合物材料。自定义混合物材料在Materials面板中进行,在13.1.4节中讲述。13.1.3 选定物质输送和反应,并选择混合物材料物质输送和容积反应的问题设置总物质模型(Species Model)面板开始(图13.1.1)。1 在Model下,选择Species Transport.2 在Reaction下,选择Volumetric reactions3 在Mixture Properties下的Mixture Material下拉列表中选择在你的问题中希望使用的混合物材料下拉列表中
22、将包括所有在当前数据库中定义的混合物。为检查一种混合物材料的属性,选择它,并点击View按纽。如果你所希望使用的混合物不在列表中,选择混合物模板(mixture-template)材料,并参阅13.1.4节了解设置你自己的混合物属性的详细内容。如果有一种混合物材料和你你所希望使用的混合物相似,可以选择这一材料并参阅13.1.4节修改已存在材料性质的详细信息。当选择Mixture Material时,混合物中的Number of Volumetric Species将在面板中显示,表达你的信息。!注意如果你在已经选定物质输送后,重新打开Species Model面板时,只有你的案例中可得到的混合
23、物材料才会显示在列表中。你可以通过从数据库中拷贝在你的案例中增加更多混合物材料,如7.1.2节所述,或是创建一个新的混合物,如7.1.2节和13.1.4节所述。正如在13.1.2节中提到的,物质输送的模拟参数和反应(如果有关)将自动从数据库中装入。如果缺少任何信息,当你点击Species Model面板中的OK按纽后将被告知缺少什么。如果你希望检查或修改混合物材料的任何属性,你将使用Materials面板,如13.1.4节所述。4 选择湍流-化学反应相互作用模型,可以使用四种模型:层流有限速率:只计算Arrhenius速率(见方程13.1-7),并忽略湍流-化学反应相互作用。涡耗散模型(针对湍
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- fluent 物质 输送 有限 速率 化学反应
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【w****g】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【w****g】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。