comsol旋转燃烧器的模拟操作实例要点

第14卷第6期2008年12月 燃 烧 科 学 与 技 术Journa l of Co m busti on Science and T echnologyV o.l 14N o .6D ec .2008微小S w iss -rol l 燃烧器的数值模拟李军伟,钟北京,王建华(清华大学航天航空学院,北京100084)摘 要:为了解微小Sw iss -ro ll 燃烧室的工作特点,对二维S w i ss -ro ll 燃烧器进行数值模拟,采用了CH 4/空气的多步反应机理,考虑了燃气对室壁的辐射,研究当量比和入口气流速度对燃烧特性和火焰稳定性的影响.研究结果表明,S w iss -ro ll 燃烧器能够在大的当量比范围内稳定工作.但上下极限并不对称,富燃时的极限比较小,比化学当量比略小,而富氧的极限比较大.对于相同的当量比,流速较小时,甲烷/空气火焰停留在燃烧器中心区的入口,甲烷全部参加反应;随着气体流速的增加,在燃烧器的中心形成回流区,扩大了燃烧器的富氧可燃极限,有助于火焰稳定,但是甲烷的转化率在减小.关键字:微小燃烧室;预混燃烧;数值模拟;甲烷;燃烧特性中图分类号:TK 16 文献标志码:A 文章编号:1006-8740(2008)06-0533-07Nu m erical Si m ulati on ofM icro Sw iss -Roll Co mbustorLI Jun-w e,i Z HONG Be-i ji n g ,WANG Jian-hua(Schoo l o f A ero space ,T si nghua U niversity ,Beijing 100084,Ch i na)Ab stract :T o understand the w ork i ng features of t he m icro Sw iss -roll combustor ,nu m erical si m ulations for a 2D S w i ss -ro ll co m bust o r w ere carried out ,w it h CH 4/a ir m ult-i step reacti on m echan i s m and rad i a ti on fro m gas to channe lw a ll taken i nto consideration .E ffects of equ i va l ence ra ti o and fl ow ve l oc ity on co m busti on characteristi cs and fl ame stability we re stud i ed .The resu lts show that the S w iss -ro ll co m bust o r can wo rk w it h i n a w ide range of equivalence rati o s .But the fue -l r i ch li m it is no t symm etr ica l to t he fue-l lean li m it .The for m er i s just l ess than 1,wh ile t he l a ter ism uch larger t han 1.Fo r the sa m e e -quiva l ence ra ti o ,if v eloc i ty is s ma l,l fla m e i s stab ilized at the en trance o f t he center reg ion w it h co m plete combusti on o fm e t hane .W ith the i ncrease o f flow ve l o city ,there appears a recircu l a ti ng zone i n t he center reg i on ,hence i ncreasi ng t he fue-l lean li m it and enhanc i ng fla m e stab ili zati on .But the rec irculati ng zone can decrease the conv ers i on ra ti o o fm ethane .K ey word s :m i cro co m bustor ;pre m i xed combusti on ;num erical si m ulation ;m ethane ;combusti on character istics收稿日期:2007-07-06.基金项目:国家自然科学基金资助项目(50376027).作者简介:李军伟(1978) ),男,博士后,davi d78lee @gma i .l co m.通讯作者:钟北京,zhongb @j ts i nghua .随着ME M S 的迅速发展,基于M E M S 技术的微型动力系统也受到了人们的关注.微型动力系统的一个重要特点是利用燃料或者推进剂在微燃烧器中燃烧,产生高温高压的燃气,驱动涡轮或者其他动力装置,从而把化学能直接或间接地转化为其他形式的能量,如动能、热能或电能.燃料或推进剂的能量密度大是其最大优点,因而微动力系统有着非常好的应用前景[1-3].基于燃烧的动力M E M S 中,为了减小微尺度燃烧器的热损失、扩展可燃极限,通过回收燃烧产物中的热能来预热未燃气体,产生过量焓火焰的技术受到了广泛关注[4].过量焓燃烧方式最开始由W einberg 和L l o yd 提出[5-7].在过量焓燃烧室中,燃烧产物和反应物在相邻的通道中逆向流动,燃烧产物通过通道壁面预热反应物,提高反应物的焓值,因此在贫燃、低热值以及大的热损失条件下,燃料仍然能够在过量焓燃烧器中维持燃烧,扩展了燃烧器的可燃极限.W ei n berg 也首次提出了二维Sw iss-ro ll(/瑞士面包圈0)燃烧器,与简单的直线型逆向换热器相比,Sw iss-ro ll结构具有更大的换热表面积,因为每两条已燃和未燃通道中间有一条未燃和已燃通道.因此,采用Sw iss-r o ll结构的微尺度燃烧器能够维持稳定燃烧.目前,对Sw iss-r o ll结构的换热燃烧器的研究中,只有Jones等人[8]对Sw iss-roll结构的换热燃烧器的定量熄火极限进行了模型研究,但是没有气体的耗散输运、未燃和已燃气体之间的换热速率以及有限速率的化学反应.为改进这种分析方法,Ronney[9]提出了U 型逆向流动热循环燃烧模型.这个模型极大地简化了复杂的Sw iss-ro ll过量焓火焰,并且考虑了壁面的换热以及有限化学反应,但这种模型不适合偏离熄火极限时的火焰.在微型燃烧器中,可燃极限内火焰的稳定及其位置的控制变得特别重要.因此,有必要了解在基本的可燃极限以下大的当量比范围内火焰的稳定及传播性质.Kuo等人[10]对丙烷/空气螺旋型逆流Sw iss-ro ll燃烧器进行了二维数值模拟研究,但仅采用了一步总包化学反应,过高地预测了燃烧的最高温度;Ahn等人[11]在螺旋型逆向流动Sw iss-ro ll燃烧器中进行了丙烷和空气混合物的燃烧实验,可是他们研究的Sw iss-ro ll燃烧器通道的尺寸比较大,大于标准状态下丙烷/空气的熄火距离;K i m等人[12]对三种不同几何形状的小型Sw iss-roll燃烧室的燃烧特性进行了实验,研究了不同的保温措施对火焰稳定,温度分布的影响,但他们并没有对燃烧器内部的燃烧特性进行更详细地分析.采用详细的化学反应机理,考虑燃烧器内部的流动和换热,通过数值模拟的方法研究甲烷/空气在微小型Sw i s s-ro ll燃烧器的燃烧特点,流量对燃烧位置的影响,当量比对甲烷转化率的影响.文中定义当量比本为实际与理论空气/甲烷质量的比值.1物理模型和数值计算模型Sw iss-ro ll燃烧器是一个直径为60mm的圆柱体,为了忽略上下端面边界散热的影响,假设燃烧器高度为无限高,着重研究燃烧器二维平面内的燃烧和流动情况.在圆柱体上按照螺旋线方程切割出了螺旋型通道,通道宽度0.65mm,两条通道之间的分割墙壁厚度1.0mm.为了防止火焰温度过高,烧毁燃烧器结构,另外两条通道之间的墙壁厚度较大,为6.1mm.未反应的甲烷/空气预混气体从外侧通道进入燃烧器,在中心区点火燃烧,燃烧后的尾气从内侧通道排出.燃烧器的二维截面和三维几何形状如图1所示.甲烷/空气的预混气体从燃烧器的入口进入,在中心大空间内燃烧.(a)二维截面(b)三维几何形状图1微小Sw iss-roll燃烧器计算模型(单位:mm)为了比较准确地模拟微小Sw iss-ro ll燃烧器的工作过程,对甲烷/空气预混气体采用二维稳态连续方程、动量方程、能量方程和组分方程,对燃烧器基体结构采用了二维稳态能量方程.由于通道很小,并且流体流速较低,雷诺数小于2000,属于层流流动范围,因此燃烧过程可以描述为预混气体的层流燃烧.为了考虑流动过程中燃气辐射的影响,辐射模型采用DO(dis-crete ord i n ates)模型[13].由于燃烧产物中CO2气体和水蒸气是主要的辐射组分,高温燃气的吸收系数通过W SGG M(w e i g h ted-sum-o-f g ray-gases m odel)模型计算[13].对这些守恒方程组采用有限体积法进行离散,采用欠松弛法进行隐式求解.求解过程中,采用S i m ple 方法离散压力-速度的耦合,动量方程、组分方程和能量方程采用一阶精度的迎风格式离散求解.在燃烧器内部的气流通道内采用较密的结构化网格,微细通道宽度方向上有15个网格点.固体壁面的结构比较复杂,采用三角形网格进行划分.计算网格见图2,共有2.1万个网格单元.#534#燃烧科学与技术第14卷第6期图2 计算网格本文的研究目的是了解微型Sw iss -ro ll 燃烧器的工作特点,因此化学反应模型采用预混燃烧的C H 4/O 2/N 2基元反应动力学机理,该机理包含35个可逆基元反应和17种组分[14],该反应模型包含的链式反应途径为CH 4y C H 3y C H 2O y HCO y CO +H 2y CO 2+H 2O式中:C H 3、C H 2O 、H CO 为中间产物,CO 和H 2为不完全氧化产物,C O 2和H 2O 为主产物.此外,自由基OH 、O 、H 和H 2O 2用来维持上述链式反应的传播过程.反应混合物的黏性系数、比热容和导热系数由各组分的相应数值通过质量分数加权平均计算.反应物的密度通过理想气体状态方程计算.各组分的比热容根据温度的分段拟合公式计算,各组分的黏性系数和导热系数根据气体分子运动学的方法计算[13].计算模型使用到以下边界条件:微细通道的入口采用质量入口边界条件,给定进气的质量流量以及各个组分的质量分数;流体和固体之间的壁面采用无滑移、没有法向组分扩散的边界条件,同时保证耦合壁面温度连续;通过壁面的热流密度连续.在微细通道的出口,给定出口的压力,其他参数根据远场条件外推;在整个燃烧器的最外层边界采用绝热边界条件,不考虑燃烧器向外界的散热.整个计算过程分两步进行.首先,不考虑化学反应,根据边界条件计算得到收敛的流场;然后,加入化学反应模型,并在中心区域设置高温点火源,待混合气体在中心区稳定燃烧后,关闭高温点火源,继续进行计算,直到最后收敛.为了确定可燃极限,首先计算得到可燃极限E R =1时稳定燃烧的流场,然后逐渐改变当量比,得到不同当量比时的流场和温度场,从而判断可燃极限及工作特点.2 计算结果与讨论2.1 可燃极限在甲烷流量为55.6m g /s 、163.1m g /s 和277.8m g /s 时,对燃烧器进行了数值模拟,得到了不同流量时的可燃极限,如图3所示,图中E R 为当量比.在3种流量下,微型Sw iss -ro ll 燃烧器的可燃极限都在0.8与2.4之间.富燃极限较小,而富氧极限较大.下面针对甲烷流量等于55.56m g /s 、不同当量比时微型燃烧器内的工作状况进行分析.图3 微型Sw iss -roll 燃烧器的可燃极限图4是不同当量比时燃烧器内的温度场云图.从图中可以看出,当量比为1时,甲烷和空气的燃烧温度最高,为1983K.随着当量比的升高,最高燃烧温度逐渐下降,在上极限(E R =2.4)时,最高燃烧温度只有1357K.如果继续增加空气流量,燃烧温度还要下降,最终散热量大于燃烧的放热量,火焰熄灭.从图中还可以看出火焰位置在燃烧器中的变化,火焰主要稳定在燃气的温度梯度和组分梯度最大的地方.当量比为1时,火焰面是一条直线,把气体通道完全封闭,甲烷和空气完全燃烧,燃烧温度最高,如图4(b)所示.当量比为0.8时,甲烷过量,空气量不足,火焰面向中心燃烧室内移动,而且火焰面也是一条直线,如图4(a)所示.由于甲烷流量不变,随着空气流量的增加,气流的速度增大.当量比为1.4和1.6时,火焰面仍然在中心燃烧室的入口下游,但是由于气流速度增加,火焰面不再是一条直线,而是一条中心凸起的曲线,如图4(c)和(d)所示.随着空气流量的进一步增加,气流速度变得更大,最终导致火焰被吹到中心燃烧室的下游.由于中心燃烧室存在一个回流区,所以火焰能够稳定在这个回流区内,通过吸附外层的甲烷和空气的混合物进行稳定地燃烧,如图4(e)和(f).在中心有回流区的#535#2008年12月 李军伟等:微小Sw iss -ro ll 燃烧器的数值模拟情况下,随着空气流量的再增加,空气带走的热量增大,中心燃烧区的面积减小,最终当损失掉的热量小于燃烧的放热量时,火焰不能维持,最终熄灭.通过上面的分析可以发现,Sw iss-r o ll燃烧器能够(a)E R=0.8(b)E R=1.(c)E R=1.4(d)E R=1.6(e)E R=1.8(f)E R=2.4图4不同当量比时燃烧器内的温度分布在大的当量比范围内稳定工作,富燃极限比较小,而富氧极限比较大.这是因为富燃时,空气量不足,甲烷燃烧不充分,燃烧放热量减小,所以容易熄火.而当富氧时,甲烷不变,空气量增加,火焰位置向下游移动,在燃烧器的中心形成一个回流区,火焰就发生在回流区的外边界,回流区能够稳定火焰,所以燃烧器能够在大的富氧极限下工作.图5是不同当量比时,C H4在燃烧器中心的质量分数云图.图中C H4质量分数梯度最大的地方就是火焰面,C H4含量从最大值突然减小到最小值.从图中可以看出,当量比为1和1.4时,预混气体在燃烧室入口完全燃烧,燃烧产物中没有剩余C H4.由于后者的气流速度比前者大,所以火焰面形成了曲面.当量比等于1.8时,气流速度增加,火焰面被吹向燃烧室下游.在燃烧室中心的底部凹陷处形成了高温回流区,回流区内是完全燃烧后的产物,回流区被甲烷和空气燃烧的火焰面包围.但是C H4气体没有完全燃烧,部分被卷吸到回流区参加燃烧,部分没有反应就直接流出燃烧#536#燃烧科学与技术第14卷第6期室.当量比增加到2.4时,回流区内预混气体完全燃烧的区域面积减小,火焰面的厚度增加,燃烧强度和最高燃烧温度减小.如果再增加空气流量,燃烧温度会进一步降低,散热量大于燃烧放热量,从而熄火.(a)ER=1(b)E R=114(c)ER=118(d)E R=214图5不同当量比时的C H4质量分数分布图6是计算得到的Sw iss-ro ll燃烧器出口尾气中各个组分的体积分数.从图中可以看出,在当量比为1时,C H4气体完全参加反应,燃烧产物主要是水蒸气和CO2气体.当量比小于1时,空气量不足,CH4没有被完全氧化,产生了部分氧化的产物,氢气和C O气体.随着当量比的增加,空气过量,产物中氧气分数逐渐增加.在当量比小于1.6时,气流速度较小,火焰没有被吹开,主要在燃烧室中心区入口,C H4被完全氧化,产物中没有C H4.但是当量比大于1.6时,气体速度增加,火焰被吹向燃烧室中心下游,部分甲烷参与反应生成水蒸气和CO2气体,其余的甲烷流出.随着当量比的增加,参与反应的甲烷含量减小,燃烧反应放热减小,最终火焰熄灭.根据甲烷在出口的物质的量的分数,计算不同甲烷的转化率,如图7所示.当量比等于0.8时,甲烷的转化率为98%.这是由于空气量不足,甲烷没有完全反应.当量比大于1.6时,由于空气流速增大,火焰被吹开,在燃烧室中心形成了回流区,甲烷部分参加反应,其余的没有反应就流出.随着当量比的增大,参加反应的甲烷量减小,转化率也减小.图6Sw iss-roll燃烧器出口尾气成分图7Sw iss-ro ll燃烧器C H4转化率#537# 2008年12月李军伟等:微小Sw iss-ro ll燃烧器的数值模拟3.2 流速的影响图8是当量比为1,不同入口流速时燃烧器内部的温度云图.总流量为0.255m /s 和0.51m /s 时,入口流速较小,甲烷和空气在入口燃烧,火焰停留在中心燃烧室的入口,甲烷反应比较完全.当总流量大于0151m /s 时,燃烧室内形成了回流区,回流区外侧被甲烷/空气的火焰面包围.随着流量的增大,燃烧高温区逐渐从中心区的下游向上游扩大,中心回流区的面积增大,燃烧放热量增大,最高燃烧温度升高.(a)0.255m /s (b )0.51m /s (c)1.01m/s(d)2.5m /s(e)4.85m /s图8 不同入口流速时燃烧器内部的温度云图根据出口气流的组分计算了甲烷气体在燃烧器中的转化率,如图9所示.从图中可以看出,随着流量的增加,入口气体流速增大,但是甲烷的转化率却逐渐下降.这是由于气流速度增加,中心燃烧室入口的火焰被吹开,部分甲烷气体没有反应就直接流出.图9 不同流量时甲烷的转化率和燃烧器入口流速图10是不同总流量时燃烧室中心的流线图.从图中可以看出,随着总流量的增加,燃烧室中心的回流区面积增大.入口流速为0.51m /s 时,预混气体在燃烧室的入口完全燃烧,燃烧产物直接从出口流出,只在燃烧室的底部形成了一个小的回流区.当入口流速增大后,火焰不能维持在燃烧室入口,向下游移动,于是就形成了一个回流区.回流区卷吸、加热未燃气体,使其在回流区的上边界燃烧,燃烧产物被吸入回流区内,最终从回流区的下游排出.另外一部分没有参加燃烧的气体沿燃烧室上部排出.另外,从图中还可以看出,随着入口流速的增加,回流区的面积增大.这是由于对于相同的当量比,流量增加,回流区的气体流速增大.#538#燃 烧 科 学 与 技 术 第14卷第6期(a)0.51m/s(b)2.50m/s(c)4.85m /s 图10 不同入口流速时燃烧室中心的流线图3 结 论(1)Sw iss -r o ll 燃烧器能够在大的当量比范围内稳定工作,而且富燃时的极限比较小,而富氧的极限比较大.甲烷流量是55.56m g /s 时,计算得到的当量比下极限为0.8,上极限为2.4.(2)当Sw iss -roll 燃烧器中气流速度较小时,甲烷/空气火焰停留在燃烧器中心区的入口,所有进入燃烧器的甲烷都参加反应;如果气流速度较大,火焰被吹开,在燃烧器的中心形成一个回流区,火焰就发生在回流区的外边界,只有部分甲烷气体流入回流区参加反应,其余的直接从回流区外侧流走,不参加反应.(3)对于相同的当量比,随着预混气体流速的增加,燃烧放热量增加,Sw iss -roll 燃烧器中心的回流区面积增大,更加有利于甲烷/空气火焰的稳定.(4)Sw iss -ro ll 燃烧器中心的回流区能够稳定火焰,并加热未反应的预混气体,使其在更大的可燃极限和流速下工作.但是回流区的存在也会导致甲烷的转化率降低,减少燃烧放热量,因此需要合理选择气体流速和当量比.参考文献:[1] Car l os F A.M i cro -pow er gene ration us i ng co m busti on :Issuesand approaches[C ]//Twen t y-N int h International Sympos i -u m on C o m bus tion .Sappo ro ,Japan ,2002.[2]Jacobson S A,Epste i n A H.An i n f o r m al survey o f pow er M E M S [C ]//T he International Sympos i u m on M icro-M e -chanical Eng ineering .T sukuba ,Japan ,2003.[3] 黄 俊,薛 宏,潘剑锋,等.微动力系统的若干研究动态和进展[J].世界科技研究与发展,2005,1(1):5-9.H uang Jun ,Xue H ong ,Pan Jianfeng ,et a.l R ecent trendsand dev elopment on m icro power syste m s[J].W orl d Science T echno logy Research and D evelopm ent ,2005,1(1):5-9(in Chi nese).[4] Ju Y,Cho i CW.A n analysis of sub -li m it fla m e dynam ics u -si ng opposite propagati ng fla m es i n m esosca le channe l[J].Com bustion and F la m e ,2003,133(4):483-493.[5] W e i nberg F J .Co m busti on te m pera t ures :The f uture [J].N at ure ,1971,233:239-241.[6] L l oyd S A,W ei nberg F J .L m i its to energy re l ease and utili za -ti on from che m i cal f uels[J].N at ure ,1975,257:367-370.[7] L loyd S A,W e i nberg F J .A recircu lati ng fl u i d i zed bedco m busto r for ex tended fl ow rang es [J ].Combustion and F la m e ,1976,27(3):391.[8] Jones A R,L loyd S A,W e i nbe rg F J .Combusti on i n heat ex changers[J].Proceedings of the Roy al Societ y of L ondon (Ser ies A ),1978,360:97-115.[9] R onney P D.Ana l ys i s of non -adiabatic heat -recircu l a ti ngco m busto rs[J].Co m bustion and F la m e ,2003,135(4):421-439.[10]K uo C H,R onney P D.N u m ericalm ode ling of non -ad iabaticheat -recircu l a ti ng combusto rs [J].P roceed i ng s of the Com-bustion Institute ,2007,31:3277-3284.[11]A hn J ,East wood C ,S itzk iL,et a.l G as -phase and cata lyticco m bustion i n heat -rec irculati ng burners[J].Procee d ings of t he Co m bustion Institute ,2005,30:2463-2472.[12]K i m N,K ato S ,T akuya K.F la m e stab ilizati on and em iss i onof s m a ll S w iss -ro ll combustors as heate rs [J].Combusti on and F la m e ,2005,141(3):229-240.[13]F l uent Inc .FLUE NT 6.1U ser c s Gui de [M ].L ebanon ,2003.[14]S m ooke M D.Reduce d K inetic M echanis m s and A symp t o ticApp rox i m ati ons for M ethane -A ir F la m es [M ].Berli n ,G e r -m any :Springer -V er l ag ,1991.#539#2008年12月 李军伟等:微小Sw iss -ro ll 燃烧器的数值模拟。

旋流燃烧器混合特性实验方案设计 如图所示的旋流燃烧器，由同轴的两根同心管道组成，中心管通燃料气，外层管道通助燃空气（带有旋流），当空气和燃料气喷入炉膛之后发生混合，并通过旋转射流的回流区卷吸炉膛内的高温烟气，因此射流中的气体由三种成分混合而成：燃料气、空气、炉膛内烟气。

为掌握燃烧器的燃烧特性，需要了解炉膛空间中各处的气体成分比例（假定暂不考虑化学反应引起的成分变化）。

一． 实验原理：由于不考虑化学反应，可考虑采用用温度场来模拟浓度场、（1） 与 均小于1，说明：动量交换过程不如热量和质量交换更强烈，∴温度和浓度混合边界层比速度边界层发展得快。

Pr 0.75,0.7~0.75t t Sc a D νν=≈=≈1~0.9t a Le D∴=≈Pr t a ν=t Sc D ν=（2）由于Le t =a /D ≈1，说明：温度和浓度边界层的发展十分相近，可以用传热过程的基本规律近似描写质量交换。

在上图中，用不同温度T 1=T 2实验，实测混合点xy 处的温度T xy （介于T 1和T 2之间，T 1>T xy > T 2）分布与浓度C xy 相似研究两股射流的混合实验，通过实验混合边界层中任一点浓度C 。

其中：C 1和C 2是被比拟的实际两股气流的浓度T 1和T 2是被比拟的实际两股气流的温度m 1和m 2是被比拟的实际两股气流的在空间中混合后的质量分数试验中采用热电偶，测量温度，取参考点为冰水混合物。

多只热电偶共用一台显示仪表的方式。

如图所示：221212xy xy T T C C T T C C --=--,1,2121,11121212122,2211122211p p p c c c p p p x x y xy y xy m m m c T m c T c T m C m C m m m m m m T C T C C C T ==⎧⎧+=+=⎪⎪+=−−−−−→+=⎨⎨⎪⎪+=+=⎩⎩二．实验系统的示意图：1➢x 2y三．实验步骤1．建立如图所示的系统；2．分别测定1.2处温度，并记录；3．沿x方向从出口处向外等距取8个点，沿y方向与同心管道中心线对称取5个点，记各点为x i y j分别测量温度，并记录4． 按 计算各点浓度； ,1,2121,11121212122,2211122211p p p c c c p p p x x y xy y xy m m m c T m c T c T m C m C m m m m m m T C T C C C T ==⎧⎧+=+=⎪⎪+=−−−−−→+=⎨⎨⎪⎪+=+=⎩⎩。

comsol旋转盘腔换热模拟
铜盘管在有冷空气流经的管道内输运热水。

铜盘管呈螺旋缠绕，以便沿轴向插入圆形空气管道的内部。

冷空气在管道内流动，热水由盘管泵入。

使用共轭传热接口计算空气流场，以及空气与铜管的温度。

由于几何基本为轴对称，我们可以假定模型与气流完全轴对称，以简化模拟。

因此，我们可以使用二维轴对称共轭传热接口。

气流速度很高，所以使用了一个湍流模型；此例中，为k-epsilon模型。

我们可以假定管道内流动的水为经充分发展的流动。

同时我们也可以假定水温的变化很小，密度不变，所以平均速度保持恒定。

因此，我们不需要模拟水的流动，而是可以利用强制对流传热关系式模拟流体和管壁之间的传热。

对流热通量边界条件使用强制内部对流Nusselt数关系式来计
算水与铜盘管之间的传热，在铜盘管的所有内部边界上都应用了此边界条件。

输入项包括管道几何、流体类型、流体速度，以及流体温度。

除了流体温度之外，所有这些量在盘管每匝中保持不变。

基于COMSOL Multiphysics的二次燃烧炉设计“有限元的未来是多物理场耦合”，因为多物理场耦合能最大程度地体现分析对象的真实工况，在真实环境下，研究对象往往同时受多个场的作用，越是逼近真实的也就越是合理和可靠的。

COMSOL Multiphysics通过耦合变量来求解各个物理场的反演方程和积分方程，耦合变量作为各个独立场的公用变量，进行同时求解，形成一个完整的总刚矩阵，从而实现与这类直接推导耦合方程类似的求解过程。

COMSOL Multiphysics在当今科学及工程领域已得到了大量成功的应用，赢得了广泛的好评。

荷兰Hi End线材名厂Siltech在25周年庆典上推出39对限量Pantheon喇叭庆祝，为了得到最优化的音响效果，特意采用COMSOL Multiphysics模拟设计出独有的内部构造。

瑞典Outokumpu Copper R&D (Västerås, Sweden)的工程师Jonas Fjellstedt通过COMSOL Multiphysics来获得提高生产金属棒效率的方法。

借助COMSOL Multiphysics的帮助，Jonas Fjellstedt成功的将原来的加工速度提高40％，这意味着公司无需新建生产线就可以增加产量，从而节约大量成本。

欧洲空中客车公司采购了100多套COMSOL Multiphysics用于先进飞机的设计，确保新研制的飞机在技术上的竞争力。

而在高校里，COMSOL Multiphysics则作为优秀的教学辅助工具，以其强大的开放性，得到了日益广泛的应用和好评。

Rowan大学K.K.Bhatia教授正在指导学生使用COMSOL Multiphysics对汽车引擎进行模拟 意大利Catania大学工业与机械工程系Giuliano C.博士等人的一个研究成果向工程界展示了COMSOL Multiphysics在计算化学燃烧，设计燃烧炉等方面的优异性能。

component1 湍流的参数设置？反应方程的各种参数设置？:流体性质设置入口边界条件设置出口边界条件设置网格的设置component2:建模部分矩形1矩形2chamfer1mesh control edge 1Form Composite Edges 1(cme 1)That concludes the geometry for the reacting jet. Now define a coupling variable that can be used to apply the outlet conditions from the previous model to the inlet of the current。

Linear Extrusion 1 (linext1)REACTING FLOW(RSPF)入口边界条件1设置物质浓度入口1边界条件物质浓度入口2边界条件出口边界条件反应项的定义Use the tabulated heat capacities to create interpolation functions, one for each species.the heat transfer interface.初始值的设置热源项设置网格设置：Edit Physics-Induced Sequence.求解设置：求解1直接求解求解2设置Now solve the reacting isothermal jet. The complicated reactions reduire adjustment of the CFL-number controller parameters.。

旋流燃烧器三维等温湍流流动的数值模拟Artit Ridluan a,Smith Eiamsa-ard a,Pongjet Promvonge ba Department of Mechanical Engineering, Faculty of Engineering, Mahanakorn Universityof Technology, Bangkok 10530, Thailandb Department of Mechanical Engineering, Faculty of Engineering, King Mongkut's Instituteof Technology Ladkrabang, Bangkok 10520, Thailand摘要本文对旋流燃烧器中的强烈旋转湍流流动进行了数值模拟。

在当前的工作中，对于三维等温旋流燃烧流动的综合科学研究使用了三种一阶湍流模型：标准k-ε湍流模型，RNG湍流模型，SST k-ω模型；二阶湍流模型，雷诺应力模型（RSM）与二阶数值差分格式。

计算结果表明RSM在测量旋流流动影响上要优于其他湍流模型。

旋流燃烧器流动的数值模拟结果为旋流燃烧器的相关设计与运行参数，包括轴向与径向速度，压力场，湍流动能的流动特点进行了描述。

1. 引言在过去的数十年中，由于旋流流动在工业上的应用，比如熔炉，燃气涡轮燃烧器，旋风燃烧炉，旋风燃烧器，灰尘收集等，它已受到相当大的关注。

一种研发用于加热小中型锅炉的新型不结渣旋流燃烧器可以在干燥的粉状及泥浆的情况下点燃含煤燃料。

在对涡流室及低温燃烧中气固两相的强烈旋转流动与燃烧的基本理解基础上，人们提出了旋流燃烧器的概念。

旋流燃烧器不仅融合了旋风燃烧器，多级燃烧器，旋风燃烧炉，粉碎性燃煤燃烧器，流化床燃烧器的优点，而且消除了他们本身的大部分缺点。

旋流燃烧器是由Nieh和Fu开发研究以求在高效，清洁方面增加煤的利用。

有一种燃烧器与Nieh和Fu的相似，但以稻壳作为替代燃料，Promvonge and Silapabanleng 对这种燃烧器进行了实验研究。

comsol操纵技能comsol操作技巧下面是本人在利用comsol解决实际问题时碰到的一些问题，通过上网查询，以及自己想出的方法进行解决，很多是网络是无法直接查询到，希望和大家一起分享，也许其中的某条正是你下载冥思苦想要解决的问题，希望能够帮到你。

上网查找的部分如有侵权，请告之删除，谢谢！（一）利用comsol的计算源程序，来建立新的循环计算如果需要利用其它计算后的数值代入到comsol进行后计算，这就不可避免的要利用comsol的源程序进行后学的连续计算，这里主要需解决在次计算利用上次计算初始值的问题，下面两段就是有无利用上次计算结果作为初始值的程序：%正常的求解% Solve problemfem.sol=femtime(fem, ...'solcomp',{'T'}, ...'outcomp',{'T'}, ...'blocksize','auto', ...'tlist',[colon(0,0.1,1)], ...'estrat',1, ...'tout','tlist', ...'linsolver','gmres', ...'itrestart',100, ...'prefuntype','right', ...'prefun','gmg', ...'prepar',{'presmooth','ssor','presmoothpar',{'iter',3,'relax',0.8} ,' postsmooth','ssor','postsmoothpar',{'iter',3,'relax',0.8},'csolver', 'pardiso'}, ...'mcase',[0 1]);%将上次求解结果作为下次计算的初始值的求解% Mapping current solution to extended mesh init = asseminit(fem,'init',fem0.sol,'xmesh',fem0.xmesh,'blocksize','a uto');%（此处是关键）% Solve problemfem.sol=femtime(fem, ...'init',init, ... %（此处是关键）'solcomp',{'T'}, ...'outcomp',{'T'}, ...'blocksize','auto', ...'tlist',[colon(0,0.1,1)], ...'estrat',1, ...'tout','tlist', ...'linsolver','gmres', ...'itrestart',100, ...'prefuntype','right', ...'prefun','gmg', ...'prepar',{'presmooth','ssor','presmoothpar',{'iter',3,'relax',0.8} ,' postsmooth','ssor','postsmoothpar',{'iter',3,'relax',0.8},'csolver', 'pardiso'}, ...'mcase',[0 1]);（二）对常数进行赋值下面一段程序是comsol中常数的赋值语句：% Constantsfem.const = {'U0','3.9', ...'I0','35', ...'It1','35', ...'It2','35', ...'It3','35', ...'It4','35', ...'Vt1','4', ...'Vt2','4', ...'Vt3','4', ...'Vt4','4', ...'OpenVt1','4.1', ...'OpenVt2','4.1', ...'OpenVt3','4.1', ...'OpenVt4','4.1'};我们可以通过以下语句对其中的35,4,4.1进行修改It1=700;It2=700;It3=700;It4=700;fem.const(1,6)={It1};fem.const(1,8)={It2};fem.const(1,10)={It3};fem.const(1,12)={It4};这样就可以实现常量变成变量，进行计算了（三）启动COMSOL with MATLAB 2007b时，MATLAB 和COMSOL Multiphysics用户界面都没有出现该问题多出现在Windows上使用COMSOL with MATLAB 2007b，MATLAB 2007a及更早的版本没有这种问题。

不同旋流强度下燃烧器燃烧特性的数值模拟$芮文明1刘宇2王爱华1张炎1田赛1(1.东北大学国家环境保护生态工业重点实验室，2.上海鼎研智能科技有限公司)摘要为研究不同旋流强度下天然气旋流式纯氧燃烧器的燃烧特性，采用计算流体力学 (C F D)商业软件Fluent对燃烧器在三种叶片角度（15°、3〇°、45°)下的天然气纯氧燃烧进行数值模拟，获得了在不同旋流强度下，燃烧室内的温度场和生成物浓度分布。

结果表明：当轴向叶片旋流器叶片角度为30°时，燃烧室内温度分布比较均匀，没有局部高温区，温度梯 度小，燃烧效果较好；在燃烧器出口 1.3m处，燃烧产物《)2剧烈产生，质量分数达到了 85%,在燃烧室后部，C02浓度逐渐降低趋于平衡；在燃烧器出口 1.5m处，燃烧产物H20 的质量分数达到42%,在燃烧室后部，H20的含量下降到40%并趋于稳定。

关键词旋流强度纯氧燃烧数值模拟温度场浓度分布Numerical simulation of the combustion characteristicsof different swirl intensity burnerRui Wenming1Liu Yu2Wang Aihua1Zhang Yan1Tian Sai1 (1. Northeastern University, 2. ShangHai Dingyan Intelligent Technology Co. , Ltd.)A b s t r a c t In order to study the combustion characteristics of the oxygen gas swirling burners at differ­ent swirl intensities, the burner with three blade angle (15°,30° , 45°)und e r different swirl intensity of oxygen gas combustion b y utilizing computational fluid d y n a m i c s (C F D)commercial software Fluentw a s simulated. T h e temperature field a n d the resultant concentration distribution in the combustionc h a m b e r at different swirl intensity wer e obtained. T h e results s h o w w h e n the blade angle of 30 ° ,the combustion c h a m b e r of axial vane swirl burner temperature is m o r e evenly distributed, no local high temperature zone, temperature gradient is small, combustion effect is better ；In the combustor exit of 1. 3m,combustion products C02 produce violently, m ole m a s s fraction is 85% ,then, in the b ack of combustor, C02 concentration decreased gradually be c o m e s m o r e b a l a n c e d；In the combustor exit of 1. 5m，H20 concentration distribution of the combustion products, the quality score of 42% ,in the bac k of combustor ,the concentration of H20 decreased to 40%a n d stabilized.K e y w o r d s swirl intensity oxygen combustion numerical simulation temperature field c o n c e n­tration distribution我国是能源消耗大国，但能源利用率比发达 国家低很多，c o2的排放量仅次于美国，占全球 总排放量的11%[1]。

使用COMSOL分析各类旋转机械在模拟旋转机械时，可以通过研究振动对机器性能的影响来有效避免机器故障。

为了实现这一目标，一种方法是使用新的“转子动力学模块”，它是 COMSOL Multiphysics® 软件“结构力学模块”的扩展模块。

在本文中，我们将介绍“转子动力学模块”，带领你了解它的实用特征和功能，助你改进旋转机械设计流程。

转子动力学建模有哪些用途？首先，我们简要介绍一下转子动力学建模。

如之前发布的一篇文章所述，旋转机械具有广泛的应用，从航空航天技术到发电，遍及众多行业，转子动力学分析十分有助于加强旋转机械的功能，提高安全性。

举例来说，假设你要确保一台发电机（一种旋转机械）避免由设计欠佳导致的不稳定、破坏性共振和故障问题。

这时你可以执行转子动力学分析，研究影响发电机物理特性的振动现象，以及由发电机的旋转和结构引起的振动加剧。

发电机（左）及其三维模型（右）。

借助仿真软件，你能够提高转子动力学研究的准确性和简易性。

现在，加上“转子动力学模块”，这一过程会变得更加方便灵活。

“转子动力学模块”会帮助你设置正确的设计参数，分析共振、应力、应变以及横向和扭转振动效应对旋转机械的的影响，由此使响应保持在可接受的运行限制范围内。

此外，你能进一步了解固定和移动的转子组件如何影响产品设计，并计算临界速度、固有频率和振型。

在下一节中，我们将深入探讨一些具体的优点和功能。

为什么要使用“转子动力学模块”？“转子动力学模块”的突出优势之一是其出色的灵活性。

你可以轻松地自定义仿真分析，方便地研究旋转装配或整个结构的特定组成部分。

上述操作的第二项可通过“转子动力学模块”的实心转子接口来实现，在有限元建模中，此接口使用三维 CAD 几何表示转子和实体单元。

你可以通过研究旋转装配中的所有组件来生成最精确的结果。

在分析转子的应力和变形时，你不需要去模拟整个系统，这样会提高仿真的精度。

为了计算整个域中的应力分布和变形场分布，你必须将转子模拟为实心单元。