《FLUENT中组分输运及化学反应燃烧模拟.ppt》由会员分享，可在线阅读，更多相关《FLUENT中组分输运及化学反应燃烧模拟.ppt（32页珍藏版）》请在三一文库上搜索。
1、Temperature in a gas furnace,FLUENT中组分输运及化学反应（燃烧）模拟,热科学与能源工程系 2003年10月,哲嗓耳镣沿灸农椿孪颓凡蛀吭周陵毡壕剩厢屉锥谓还砸陀熏拳谚枢授阻饱FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,概要,应用 燃烧模拟简介 化学动力学 气相燃烧模型 稀疏相燃烧模型 污染物排放模拟 燃烧数值模拟步骤介绍,沙吹扩藐哮砾吝爽墒遇疹骗闽粳鲍胁糕狮滇崔宝洁腆健仔衷垣错牲托兔芜FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,燃烧模拟,广泛应用与均相和非均相燃烧过程模拟 燃烧炉 锅炉
2、加热器 燃气轮机 火箭发动机 求解内容 流场流动特性及其混合特性 温度场 组分浓度场 颗粒和污染物排放,Temperature in a gas furnace,CO2 mass fraction,Stream function,炊湾百助蕉殷皿蕾锹迪岛哗秋帘召锗坚葡蒂捌脾羔狗借骗漳构酿州蛋鹅板FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,燃烧模型概要,稀疏相模型,液滴/颗粒动力学 非均相反应 液化 蒸发,输运控制方程,质量 动量 (湍流) 能量 化学组分,燃烧模型,预混 局部预混 非预混燃烧,污染物模型,辐射换热模型,徊旗槛曼惧宅藏蚜尼褥试隶朵畔嘛硅拔喀视凛
3、塌把狡厦想侍反茬我少葵难FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,气相燃烧 有限速率模型 (Magnussen model) 守恒标量的 PDF模型 (一个或两个混合分数) 层流火焰面（小火焰）模型 (V5) Zimont model (V5) 稀疏相模型 湍流颗粒弥散 随机轨道模型（Stochastic tracking） 颗粒云团模型（Particle cloud model） (V5) 粉煤与喷油燃烧子模型 辐射模型: DTRM, P-1, Rosseland 和Discrete Ordinates (V5) 湍流模型: k-, RNG k-, R
4、SM, Realizable k- (V5) and LES (V5) 污染物排放模型: NOx with reburn chemistry (V5) and soot,FLUENT提供的燃烧模型,冒粗逢他裹典砒涕员锑蛾旗凛籍燥轰列氨欢盯淬傅忱割序舜猴俱蔷翟迷睬FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,模拟燃烧过程的化学反应动力学,难点与挑战 多数实际的燃烧过程是湍流 化学反应速率高度非线性; 湍流化学反应高度耦合，相互作用很重要。 真实化学反应机理包含数十个组分, 数百个基元反应，并且方程组极具刚性 (基元化学反应时间尺度相差大) 实际处理方法 简化化
5、学反应机理 有限速率燃烧模型 考虑湍流及其混合、弱化反应化学 混合分数模型 平衡化学的 PDF模型 层流火焰面模型 进展变量模型 Zimont 模型,唐唬砸免搬铺坡桐镣激钉沦犊躯驰耪伸泥渤占君钢呼刀钎础梗灌蔗易找扰FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,有限速率模型,用总包机理反应描述化学反应过程. 求解化学组分输运方程. 求解当地时间平均的各个组分的质量分数, mj. 组分 j的源项 (产生或消耗）是机理中所有k个反应的净反应速率 : Rjk (第k 个化学反应生成或消耗的j 组分)是根据 Arrhenius速率公式、混合或涡旋破碎（EBU）速率的小
6、值。. 混合速率与涡旋寿命相关, k /. 物理意义是湍流涡旋是决定化学反应的首要因素。对于非预混燃烧，湍流涡旋决定了组分混合；对于预混燃烧湍流决定了热输运（高温加热低温）。即：化学反应决定于湍流混合组分（非预混燃烧）和热量（预混燃烧）的速率。,丁虚姚倍扭时氧焦仿郡靡桔祁兔剂金款拦裁恳棒鲁鲸绚沮孔吨榆食福监源FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,Fluent燃烧模型,有限速率模型 求解组分的质量分数输运方程，化学反应机理由用户自己定义。 非预混燃烧模型 该模型中并不求解单个组分的输运方程，而是求解一个或者两个守恒标量（混合分数）的输运方程 预混燃烧模
7、型 模拟完全混合的燃烧问题。充分混合的燃烧物和产物被火焰前锋分隔，求解出的化学反应进展变量来描述该火焰前峰的位置 部分预混燃烧模型 该模型用来处理系统中同时具有非预混和充分预混的情况。该方法同时求解了混合分数和反应进展变量,牺亡夺曲弊悔拼峭哟票孙仓歉役纤县走谤郸戈霹头颧嫉歪氯矮带毡锦宣顿FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,有限化学反应速率模型设置,要求: 给出组分及其物性 给出化学反应及其反应速率在内的化学反应动力学数据 FLUENT V5 在mixture material database里面提供了数据 对于常用的燃料，数据库都会给定机理，组分
8、物性等信息. 如果用户需要给定个性化机理，则: 生成新的混合物. 改变已有混合物的物性/化学反应.,祟息莽讽竹辟伯何走颁骤祟臻球堵缚追囊伙萍雷冉偶朴斟壁食瑚迭炽每售FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,有限速率模型小节,优点: 可以应用于nonpremixed, partially premixed和premixed combustion 简单、直观 应用广泛 缺点: 不适合混合速率与化学反应动力学时间尺度相当时候的化学反应 (要求 Da 1). 没有严格考虑湍流化学反应之间的相互作用问题 不能考虑中间产物或组分、不能考虑分裂影响. 模型常数不确定,
9、 特别是用于计算多个化学反应的时候尤为如此，模型常数通用性较差。.,摄软慑樊点寄杜卉甫杖假孰白烤物些轧暑绥映碰蔓芋午怂娜糠讥鳖啸紧烤FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,守恒标量 (混合物分数) 模型: PDF 模型,只适应用于非预混 (扩散) 火焰燃烧 假定化学反应过程受混合速率控制 满足局部化学平衡. 控制体（计算单元）组分、物性决定于燃料和氧化剂在该处的混合程度. 化学反应机理不明确. 用化学平衡计算来处理化学反应 (prePDF). 只求解混合物分数及其方差的输运方程, 无需求解组分的输运方程. 可以严格考虑湍流与化学反应的相互作用,哟掳吞炙
10、舵悠碘袁修杉狡锑粟像震肄寞册蕾暗磨狠决月葬贯憨窒秒毛拴慕FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,混合分数定义,混合分数, f, 写成元素的质量分数形式: 其中， Zk 是元素k的质量分数 ；下标 F 和O 表示燃料和氧化剂进口流处的值。 对于简单的 fuel/oxidizer系统, 混合物分数代表计算控制体里的燃料质量分数. 混合物分数是守恒标量（conserved scalar）: 组分输运方程中没有化学反应源项.,曳狱鹊痰汕霞掸奸倘崭缺廊纷粉歹徘披潍谤豌症诧会旧裸拦豆蒜淹悉癌虫FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应
11、燃烧模拟,可以用单个混合物分数模拟的燃烧系统,Fuel/air 扩散火焰: 多氧化剂入口的扩散火焰: 多燃料进口的扩散火焰:,60% CH4 40% CO,21% O2 79% N2,f = 1,f = 0,35% O2 65% N2,60% CH4 40% CO,35% O2 65% N2,f = 1,f = 0,f = 0,60% CH4 20% CO 10% C3H8 10% CO2,21% O2 79% N2,f = 1,f = 0,f = 1,60% CH4 20% CO 10% C3H8 10% CO2,优伟淮舞硕毁凑葱歌肚霉购户辉夺唆跃嘱猿倒洲蛆专呛交泄蔗痕滴拥粘渠FLUENT中
12、组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,系统化学平衡假设,化学反应很快到达平衡. 可以考虑中间组分.,排遂却牵弱令司徊谩阁疑址部问槛借沧控篱虱遂嗓股诊封枉奏但潦位廉蹄FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,PDF 模拟Turbulence-Chemistry相互作用,Fluctuating mixture fraction is completely defined by its probability density function (PDF). p(V), the PDF, represents fraction of
13、 sampling time when variable, V, takes a value between V and V + V. p(f) can be used to compute time-averaged values of variables that depend on the mixture fraction, f: Species mole fractions Temperature, density,漆潍糟箱涧阀削棕岿官非瘟良迭芭嗽厨抡暑韵寓区脓佳镰有荧陌郧售祈他FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,PDF Model Fle
14、xibility,非绝热系统: In real problems, with heat loss or gain, local thermo-chemical state must be related to mixture fraction, f, and enthalpy, h. Average quantities now evaluated as a function of mixture fraction, enthalpy (normalized heat loss/gain), and the PDF, p(f). 第二守恒标量: FLUENT用第二守恒标量可以模拟: Two f
15、uel streams with different compositions and single oxidizer stream (visa versa) Nonreacting stream in addition to a fuel and an oxidizer Co-firing a gaseous fuel with another gaseous, liquid, or coal fuel Firing single coal with two off-gases (volatiles and char burnout products) tracked separately,
16、拒醉率罩扮瞬扬心周鹏贸庶澜晒督斋价启乘猖房卯晾僵腔渗豁卤住怜嗅冷FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,混合分数/PDF模型小节,优点: 可以计算中间组分. 考虑分裂影响. 考虑湍流化学反应之间作用. 无需求解组分输运方程（特别是多组分），简化计算量 性能好，经济 缺点: 系统必须满足（靠近）局部平衡. 不能用于可压速或非湍流流动. 不能用于预混燃烧.,兄掖圾谎埔瘫绦农衷鼠盛研赌渊镑废蚤今卢羞赚剧椿峦铝饭把察兔馋按幂FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,层流火焰面模型,用混合分数和标量耗散率来求解（绝热）温度
17、、密度和组分等量。 对于混合分数 PDF 模型 (绝热), 热化学状态只是 f 的函数 c 与当地应变率有关,把混合分数 PDF扩展到模拟中度化学非平衡燃烧模拟中 用层流拉伸火焰系综来模拟湍流火焰, 对撞扩散火焰,灯江囊淘熏米碰荤巢眩婴众几勋傍盟令意庞夜白局旺醛热巳脓订铃泄素捧FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,层流火焰面模型(2),用指定概率密度函数（PDF） P(f,c) 的方法来决定层流火焰面系综。假定的概率包括计算为：Pf (f) Pc (c), 其中， Pf (f) 用Beta 函数， Pc (c) 用delta函数 只适合绝热系统（FL
18、UENT V5） Import strained flame calculations prePDF or Sandias OPPDIF code Single or multiple flamelets Single:user specified strain, a Multiple: strained flamelet library, 0 a aextinction a=0 equilibrium a= aextinction is the maximum strain rate before flame extinguishes Possible to model local exti
19、nction pockets (e.g. lifted flames),桑典辟刚贾敖伺用社淖功悦赎喉腾迫碧浙絮猖腑贷柴捆姨腮阐守搓备狡雁FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,预混燃烧的Zimont模型,用单个过程变量来模拟热化学过程, 平均反应速率, 湍流火焰传播速度, Ut, 根据贫燃预混燃烧推导得到，并考虑 预混燃料当量比 湍流引起的火焰前锋皱折和增厚 湍流拉伸引起的火焰前锋淬熄 分子扩散 适合绝热燃烧, 非绝热燃烧必须求解焓方程,柏旱掳歧搜惊怔检境指棕涧侠娶央薯黄炙浅轻芋拐钦莎路盈霖棠勾垃圃吁FLUENT中组分输运及化学反应燃烧模拟FLUENT
20、中组分输运及化学反应燃烧模拟,稀疏相模型,在拉格朗日坐标系下，求解颗粒、液滴、气泡的轨道. 并与连续相（气相）进行热、质量和动量的耦合求解。 稀疏相体积分数必须 10% 质量分数可以比较高 不考虑颗粒之间作用、颗粒破碎 模拟湍流弥撒 随机轨道（Stochastic tracking） 颗粒云团模型（Particle cloud (V5)） Rosin-Rammler or linear size distribution 非定常流动中的颗粒轨道 (V5) 模拟颗粒分离, 喷雾烘干, 液体燃料或煤粉的燃烧过程。.,降葫夺绅汀友扫框陨冗枉缕桃穿蝉申射栗庇缄炭蘑稍哆植柏穷匀婉馆朴橡FLUENT中组分
21、输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,Monte-Carlo方法模拟湍流颗粒弥散 (discrete random walks) 颗粒运动计算中考虑气体的平均速度及随机湍流脉动速度的影响。 每个轨道包含了一群具有相同特性的颗粒，如相同的初始直径，密度等. 考虑湍流弥散尤为重要 更复合真实物理过程，但计算量更大。 可以通过光滑源项、消除与气相的耦合来强化计算稳定性；,颗粒弥散: 随机轨道模型,Coal particle tracks in an industrial boiler,惊峪市衙碉涨克委渴呆应缝旋当贿覆裴梢实资杰妈酉奔畅含栗志渊碘具蔡FLUENT中组分输运及化
22、学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,颗粒弥散: 颗粒云团模型,用平均速度决定颗粒的平均轨道 假定该平均轨道为 3D多变量的 Gaussian 分布，计算颗粒偏离的范围（3倍标准偏差） Rigorously accounts for inertial and drift velocities A particle cloud is required for each particle type (e.g. initial d,r etc.) Particles can escape, reflect or trap (release volatiles) at walls E
23、liminates (single cloud) or reduces (few clouds) stochastic tracking Decreased computational expense Increased stability since distributed source terms in gas phase BUT decreased accuracy since Gas phase properties (e.g. temperature) are averaged within cloud Poor prediction of large recirculation z
24、ones,抑性醉赴呛精侩梅耸刨韵锯俗改瓮绰卓耍为钝颁瓢官疤滩连凹茁藩胰绅楷FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,非定常流动的颗粒轨道模拟,Each particle advanced in time along with the flow For coupled flows using implicit time stepping, sub-iterations for the particle tracking are performed within each time step For non-coupled flows or coupled
25、flows with explicit time stepping, particles are advanced at the end of each time step,丁担程觅圾遣沾篓歧守薄踞究艳瑰直缴胺曼竟填穴添笼申阴缝朝仕书深菊FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,煤/油燃烧模型,变化模型可以模拟煤或油燃烧l Droplet 油燃烧模拟 Combusting particle 煤燃烧模拟 Several devolatilization and char burnout models provided. Note: These models
26、 control the rate of evolution of the fuel off-gas from coal/oil particles. Reactions in the gas (continuous) phase are modeled with the PDF or finite rate combustion model.,辱忆蟹茂罗鸥巧战竹炒卡陷屠蓉辰脆掐惮删民谦疏锯掏殿拖球灿状资弱帛FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,NOx Models,NOx 中含量最多的是一氧化氮 (NO). 引起烟雾 引起酸雨 导致臭氧损耗 NOx
27、 排放有三个机理: 热机理（Thermal Nox） - Zeldovich mechanism (氧化空气中的 N) 高温燃烧最明显 快速机理（Prompt NOx ）- De Soete, Williams提出的经验机理. 通常该机理对Nox排放量影响较小 富燃燃烧区域该机理有明显作用 燃料 NOx 机理- De Soete, Williams等人提出的经验机理. 在煤燃烧中， 燃料中含N高，并且燃烧温度不是很高时起作用。 NOx 再燃 (V5) 在燃料富足的燃烧区域，NO会与碳氢化合物反应，使得NO 减少,篷速豁瘦径燃愧护犬梁欺炽召等舰光毖剐艾防氮绰洪蝶怔姨码泅锡殊烁悸FLUENT中组分
28、输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,烟（Soot） 模型,两种烟模型: One-step 模型 (Khan and Greeves) 1个输运方程，烟质量分数 Two-Step 模型(Tesner) Transport equations for radical nuclei and soot mass fraction concentrations Soot 生成采用经验公式 其中, C, pf, 和 F 是模型常数, 燃料分压力和当量比。 采用Magnussen模型计算 Soot 燃烧 (消耗) Soot 对辐射的影响 激活Soot面板上的Soot-Radiat
29、ion 选项,智倚返蓉科哟泛籽语疤莹紫今牵钾审鹿醇著先豆碾阵帜颅与絮梭览读勺混FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,燃烧模拟步骤及求解策略,打开求解器（2D或3D） 检查物理模型适用性 网格分辨率需求 (resolve shear layers) 求解参数、收敛条件设置 边界条件 燃烧经常对进口边界条件十分敏感 正确的进口标量和速度的分布 很难估计壁面热条件；尽可能优先给定边界温度，很难给定对流换热或者辐射换热量 初始条件 稳态解跟初始值无关，但给的初始值不好，由于各个方程之间不和谐和输运方程的非线性，导致解分歧或发散不收敛。先求解冷态流动、接着气体
30、燃烧、颗粒燃烧、再考虑辐射； 对于强旋流动，渐渐增加旋流速度；,昨辊详毙知邓逆炉卤媚堕构钥廷镣属祥坷哦懂粤寐十浚叛剐硒反数滞薄蚕FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,燃烧模拟步骤及求解策略（2）,低松弛因子 松弛因子改变对结果影响是高度非线性的 采用混合分数PDF模型时，密度采用低松弛 (0.5) 考虑浮力流动中，速度采用低松弛因子 高速流动中压力采用低松弛因子 一旦解较为稳定, 尽可能让松弛因子提高到系统默认值 (温度， 0.9 ； P是1；旋流速度，组分（混合分数）也是1或接近1。 离散化 开始用一阶精度, 收敛后用二阶精度提高求解精度 对于三角
31、形或四面体网格，二级精度离散尤为关键 离散相模拟 增加稳定性, 增加随机轨道数目 (或者采用颗粒云团模型) 低松弛因子，增加每次耦合计算中颗粒相计算后气相叠代次数,练舒精确傲骂霍顷凿螟汞凋前试衡傣箕匀眠嗜誉娜旅躁砍瓶默夯挫党哩准FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,燃烧模拟步骤及求解策略(3),Magnussen 模型 默认设置 finite rate/eddy-dissipation (Arrhenius/Magnussen) 对于非预混 (扩散) 火焰，关掉有限速率模型 对于预混火焰，应该采用Arrhenius项，使得反应物不至于过早燃烧。 需要
32、高温点火，保证反应进行（ initialization/patch） 考虑比热 Cp随温度变化，以便计算的温度不至于过高而不切合实际。 混合分数 PDF 模型 如果符合该模型假设，可以选用该模型 建表中要保证有足够多的点以便保证插值精度，同时要不影响计算时间。 采用beta函数 PDF,乖幕装脊拨备抒慷刽锭甸堪哇监耐钻俺恩衍械蛹禾场善厉衷憎锻堵帕役妥FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,燃烧模拟步骤及求解策略(4),湍流 开始用 标准k-e 模型 收敛后转变为RNG k-e , Realizable k-e 或 RSM ，可以与实验结果比较、也可以比
33、较结果受湍流模型的影响情况 收敛判断 一般变量残差小于 10-3 ，温度、 P-1 和组分的残差要小于10-6 The mass and energy flux reports must balance 监视和检查所关心的变量 (如：出口平均温度) 检查变量的等值线图是否光滑、合理及稳定不变。,泊固磐抠躇娥明叉昨回纳元扁挝乘竟柒眶伟银帕慢途椎混均警曹籍琉条能FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,小节,FLUENT V5 用于模拟燃烧的优点在于： 一组较好的燃烧模型 方便、易用 数据库包含了机理与物性参数 网格与求解方案可以自由调节 FLUENT可以求解多数的工程实际燃烧问题. 留意选择的物理模型是否合适于实际需要模拟的燃烧过程.,泅混挪钻擦北主句桃讣慈云社芹廷篙扬里氨泰幸印呆逝箱挡说暖尊蛾让偿FLUENT中组分输运及化学反应燃烧模拟FLUENT中组分输运及化学反应燃烧模拟,