反应流模型的课件资料

反应流模型的课件资料第1页/共55页TemperatureinaGasFurnaceCO2MassFractionStreamFunction反应流系统应用FLUENT

化学反应模型可以应用于很大范围内的均匀和非均匀反应流。熔炉锅炉加热器燃气涡轮火箭发动机可以预估流场和混合特征温度场介质浓度微粒和污染物 （pollutants）第2页/共55页反应流模拟的难点湍流大多数工业产生的流体流动都是湍流。因为具有很大的时间和空间尺度，所以对非反应湍流和反应湍流，DNS方法并不适用。化学现实中的化学反应不能使用单一的反应方程来表示。十几种组分，可能有上百种反应。只有有限数量的燃料知道反应细节。反应时间尺度范围大湍流和化学间的交互作用化学反应速率的变化和湍流中物质混合程度有十分密切的关系。第3页/共55页燃烧中的化学动力学模拟简化化学模型(使用有限速度/涡耗散方法).考虑全局的化学反应原理。化学反应与湍流混合过程分离考虑（使用混合比方法）平衡化学PDF模型层流火焰模型部分预混模型模拟详细化学反应(占资源)需要特别小的时间步长来满足数值的稳定性和收敛性使用stiffsolver

将可以允许稍大的时间步长。第4页/共55页反应流模拟总览输运控制方程质量动量湍流能量化学组分离散相模型液滴/质点动力学异相化学反应液化作用蒸发作用反应模型快速化学反应有限速率化学反应预混部分预混非预混污染物模型(NOx,烟灰等.)辐射模型第5页/共55页反应系统定义反应系统可以分成两类同相反应反应物与产物同相。异相反应反应物与产物不同相。表面沉积的化学物质被认为与相同组分的气相物质不同。沉积率同时受化学动力学和组分扩散率影响。壁面反应是产生体相中化学物质的来源，而且决定了

表面的沉淀比例。表面反应：CVD,催化反应第6页/共55页同相反应非预混反应系统反应物从火焰的两边对流/扩散进入。可以简化成混合问题。湍流旋涡扰动层流火焰形状加强了混合。预混反应系统传播速率（火焰速度）取决于内部火焰结构。比非预混燃烧问题更难建模。湍流扰动层流火焰形状，加速火焰传播。部分预混反应系统反应系统包括非预混和预混流。FuelOxidizerReactorOutletFuel+OxidizerReactorOutletFuel+OxidizerReactorSecondaryFuelorOxidizerOutlet第7页/共55页反应流问题的特征反应流问题通过一系列无量纲变量来表示。两种最重要的无量纲变量是Reynolds数和Damköhler数。Reynolds数

ρ,U,L,μ

分别是密度，速度，特征长度和动力粘性系数

高雷诺数将表现为湍流Damköhler数高Damköhler数（Da>1）表示需要使用气相湍流燃烧模型。第8页/共55页反应流问题的特征马赫数Machnumber高马赫数(Ma>0.3)表示可压缩性表现明显。对于低马赫数流动，可以使用非预混模型。Boltzmann数高Boltzmann数表示对流是主要的传热形式。当Bo<10时，需要考虑辐射传热。Stefan-Boltzmannconstant5.672×10-8W/m2·K4Localspeedofsound(sonicvelocity)第9页/共55页*Rateclassificationnottrulyapplicablesincespeciesmassfractionisnotdetermined.CHEMISTRYFLOWCONFIGURATIONPremixedNon-PremixedPartiallyPremixedFastChemistryEddyDissipationModel(SpeciesTransport)PremixedCombustionModelNon-PremixedEquilibriumModelPartiallyPremixedModelReactionProgressVariable*MixtureFractionReactionProgressVariable+MixtureFractionFinite-RateChemistryLaminarFinite-RateModelEddy-DissipationConcept(EDC)ModelCompositionPDFTransportModelLaminarFlameletModelSteady/UnsteadyFLUENT中反应流模型概览第10页/共55页涡耗散使用1步或者2步总包（热扩散）反应原理.化学反应速率由湍流混合速率决定。有限速率项选择预混合火焰。非预混模型使用混合物比例和假定的PDF替代求解化学平衡中的物质输运和反应比例方程。层流小火焰模型用于中等的非平衡化学反应。预混模型使用反应程序.c文件建立化学反应。湍流火焰速度闭合--Zimont模型。部分预混合联合非预混和预混模型。两种模型都假设成立。涡耗散概念(EDC)为湍流火焰建立湍流/化学交换作用的模型。包括详细的化学反应。PDF传输模型比EDC模型更加严格的考虑湍流-化学之间的相互作用。考虑更充分，但也占用更多的计算机资源!同相（气相）燃烧模型第11页/共55页基本概念FLUENT中的反应流模型总览快速化学模型涡耗散模型(EDM)预混模型非预混模型部分预混模型有限速度化学模型层流有限速度模型涡耗散概念模型(EDC)组分PDF输运模型非预混火焰中的层流小火焰其它模型污染物行成模型(PollutantFormationModels)表面反应附录大纲第12页/共55页涡耗散模型(EDM)适用条件流域：

湍流(高雷诺数)化学： 快速化学反应(高Da)构造：

预混/非预混/部分预混案例：气体反应煤燃烧限制条件当混合和运动时间尺度大小相似时(Da~1)，结果不可靠。不能预测中间组分和离解效果。不能模拟动力学细节现象，比如点火，熄火和低Da数流。求解物质输运方程，湍流混合控制反应速率。第13页/共55页预混燃烧模型适用条件：流域： 湍流(高雷诺数)化学： 快速化学反应(高Da数)构造： 只有预混气流案例：预混合反应流系统贫氧预混燃气涡轮燃烧室限制条件不能模拟化学动力学中的细节（比如点火，熄火和低Da数流动）。使用反应进展变量追踪火焰前沿轨迹(Zimont模型)。Fuel+OxidizerReactorOutlet第14页/共55页非预混平衡模型适用条件流域： 湍流(高雷诺数)化学： 快速化学反应(高Da),假定平衡环境构造： 只有非预混合案例气体反应（熔炉，火炉）。对于需要模拟中间产物浓度和离解效应的问题，那么采用这种模型是个较好的选择，不需要知道详细的化学反应速率。限制条件当混合和运动时间尺度大小相似时(Da~1)，结果不可靠。不能模拟化学动力学细节现象（比如点火，熄火，低Da数流动）。求解输运方程中的混合分数和混合分数变化。(不是单独的物资方程)。FuelOxidizerReactorOutlet第15页/共55页部分预混燃烧模型适用条件流域： 湍流(高雷诺数)化学： 平衡或者适度的非平衡（火焰）构造： 只有部分预混合案例有冷却孔的燃气涡轮燃烧室。同时有预混和非预混流的系统。限制条件当混合和运动时间尺度大小相似时(Da~1)，结果不可靠。不能模拟化学动力学细节现象（比如点火，熄火，低Da数流动）。部分预混模型中，结合了反应进度变量和混合分数方法。输运方程需要求解反应进度变量，混合分数和混合分数变化。Fuel+OxidizerReactorSecondaryFuelorOxidizerOutlet第16页/共55页大纲基本概念FLUENT中的反应流模型总览快速化学模型涡耗散模型(EDM)预混模型非预混模型部分预混模型有限速度化学模型层流有限速度模型涡耗散概念模型(EDC)组分PDF输运模型非预混火焰中的层流小火焰其它模型污染物行成模型(PollutantFormationModels)表面反应附录第17页/共55页详细化学反应模型详细化学反应的控制方程一般难以求解：上十种物质上百种反应反应时间尺度间的差别较大详细化学动力学建模：火焰点燃和熄灭污染物(NOx,CO,UHCs)慢化学反应(非平衡)液体间反应可用的模型：层流有限速率涡耗散概念模型(EDC)PDF输运动力学模型(需要附加的license)CHEMKIN-格式的反应机理和热物理性质，可以直接输入。FLUENT

使用In-SituAdaptiveTabulation(ISAT)算法加速计算（可以用于层流，EDC,PDF输运模型）

。第18页/共55页层流有限速率模型适用条件流域： 层流(低雷诺数)化学： 有限速率化学反应构造： 预混/非预混/部分预混案例层流扩散火焰。煤气炉中的层流预混合火焰。限制条件需要占用较多的CPU资源。默认使用ISAT算法进行计算加速。一般不适合湍流燃烧模型。层流有限速率模型需要使用刚性化学求解器（stiffchemistrysolver）(pressure和density求解器都有效)。第19页/共55页涡耗散概念模型(EDC)适用条件：流域：

湍流(高雷诺数)化学： 有限速率化学反应构造： 预混/非预混/部分预混案例：湍流反应流中的有限速率化学反应CO的缓慢燃烧NOx形成限制条件：需要占用很多的CPU资源。默认使用ISAT算法进行计算加速。 （附录中有详细说明）。EDC模型是涡耗散模型在湍流中模拟有限速率化学反应的扩展。第20页/共55页PDF输运模型适用条件：流域： 湍流(高雷诺数)化学： 有限速率化学反应构造： 预混/非预混/部分预混案例：湍流反应/燃烧流，预混有限速率现象，例如CO和NOx的形成液体间的反应限制条件PDF输运模型使用了蒙特卡洛方法，十分占用计算机资源。默认使用ISAT算法进行计算加速。湍流火焰加入了有限速率化学反应，使用了复合概率密度函数 （PDF）输运模型。比EDC模型更严格使用蒙特卡洛方法第21页/共55页定常和非定常小火焰模型适用条件：流域： 湍流(高雷诺数)化学： 由于空气动力学约束形成的适度 非平衡构造： 只有非预混合案例：预测喷流火焰内燃机（只有柴油机非定常火焰子模型）液体间的反应(非定常火焰子模型)限制条件定常方法不能模拟详细化学动力学中的实际现象(比如点火，熄火和低Da数流)。层流火焰模型中，解算出输运方程中的混合分数变化。为考虑非平衡效果，计算了梯度扩散（火焰应变比例）。Flameletstructure第22页/共55页大纲基本概念FLUENT中的反应流模型总览快速化学模型涡耗散模型(EDM)预混模型非预混模型部分预混模型有限速度化学模型层流有限速度模型涡耗散概念模型(EDC)组分PDF输运模型非预混火焰中的层流小火焰其它模型污染物行成模型(PollutantFormationModels)表面反应附录第23页/共55页污染物形成模型NOx

形成模型(定性预测NOx

形成趋势).FLUENT

包括三种计算NOx产物的反应模型ThermalNOx

PromptNOxFuelNOxNOx

再燃烧模型选择Non-CatalyticReduction(SNCR)模型烟灰形成模型一步和两步模型烟灰影响辐射吸收(在Soot

面板中选择Soot-Radiation选项)SOx

形成模型求解SO2,H2S的附加方程,以及SO3

附加方程可选。通常，SOx预测是个后处理过程。第24页/共55页表面反应化学物质沉淀在表面，和气体中相同的化学物质有个明显的分界。每个表面吸附（或者点）物质都要求解点平衡方程。可以考虑详细的表面反应机理(任意数量的反应步和任意数量的气体相或/和点物质)。CHEMKIN格式的表面化学反应机理可以导入。表面反应可以发生在壁面或者多孔介质上。不同的表面反应机理可以运用在不同的表面。案例催化反应CVD第25页/共55页Outline基本概念FLUENT中的反应流模型总览快速化学模型涡耗散模型(EDM)预混模型非预混模型部分预混模型有限速度化学模型层流有限速度模型涡耗散概念模型(EDC)组分PDF输运模型非预混火焰中的层流小火焰其它模型离散相模型(DPM)污染物行成模型(PollutantFormationModels)表面反应附录第26页/共55页附录第27页/共55页涡耗散模型(EDM)估计平均反应率（物质k），Rk模型。使用eddy-breakupmodel模拟反应率假设化学反应比湍流混合化学物质进入火焰的速度快得多。(Da>>1)燃烧完全通过湍流混合比例控制，湍流混合比例是指大涡存在时间标量比例,k/ε。总包反应（一步或两步）每种物质都需要解算雷诺（时间）平均物质质量分数方程。有限比例/涡耗散选项反应率是Arrhenius比例和涡耗散比例的最小值。第28页/共55页非预混模型湍流中无限快速反应，当油和氧化剂混合（有限的混合）的时候，瞬间达到化学反应平衡。物质和焓输运方程合并成为单一的混合分数守恒方程，f，它满足下列假设条件：物质扩散系数相等Speciesdiffusioncoefficientsareequal每种物质的Lewis数必须一致相等中等马赫数混合分数可以直接和物质质量分数，混合密度，混合温度相关。这些关系称作瞬间（不是时间平均的）数据。使用概率密度函数把瞬间数据转变成为时间平均数据（湍流反应流模拟预估时间平均属性）。第29页/共55页复合概率密度函数（PDF）输运模型用来模拟湍流焰中的有限速率化学反应。比EDC模型更严格可以从PDF中通过下式计算出平均反应率：PDF输运方程有十分高的维度，而且不能用有限体积方法求解，是通过蒙特卡洛方法解算的。复合PDF输运模型第30页/共55页In-SituAdaptiveTabulation(ISAT)算法初衷当使用详细化学反应模型（EDC，综合PDF，层流有限比例）时，每个单元（或者每个粒子）在每步迭代都需要求解刚性反应系统（ODE的集合）需要一种高效的方法获得只关心部分的数据。In-situadaptivetabulation(ISAT)方法可以加快反应表格的生成和数据的查看。优点避免重复的数值积分，节省了计算机CPU的计算时间。初始迭代比较慢（表格正在建立），但是计算在加速运转中。适用条件ISAT可以使用于PDF输运,EDC,和层流有限比例模型第31页/共55页NOx

模型NOx

由氮的氧化物组成，对环境有害FLUENT

有三种计算NOx产生的机理。ThermalNOx

高温时十分有效PromptNOx在富燃区有效FuelNOx燃料中氮的含量高而且温度也相对较低时使用User-definedfunctionNOx

再燃NO可以通过和富含燃料区域中的碳氢化合物反应而减少。Non-CatalyticReduction(SNCR)模型。NOx

形成的定性预测第32页/共55页示例富燃，预混火焰EDCPDF传输层流有限速率IFRF旋涡煤粉火焰详细化学反应模型DPM辐射GELM-1600燃气涡轮燃烧室非预混，

层流小火焰温度和燃料NOx

预测D形西瓜火焰（SandiaFlameD）非预混湍流扩散火焰EDC/PDF传输比较催化燃烧表面反应多孔介质第33页/共55页案例学习:富燃预混合火焰富燃预混合预蒸发（LPP）丙烷/空气火焰减少燃烧污染物(CO和NOx)的扩散轴对称火焰管几何量 d=10.25cm,L=42.92cm丙烷/空气混合物 V=25m/s,T=800K火焰稳定段 51%openareaPropane/air7,056quadelements*D.N.Anderson,NASALewisResearchCenter,NASA-TM-X-71592,March1975.第34页/共55页富燃预混合火焰湍流化学相互作用弱完全预混合火焰区薄化学反应正好在火焰阵面之后发生丙烷详细的化学机理NO和N2O形成的Kundu机理17种物质和23步反应有限速率化学燃烧模型层流有限速率模型完全忽略湍流-化学间的相互影响。EDC模型可能导致NOx的低预测—反应发生非常快，而且NOx只有很少的时间累积复合PDF输运正确模拟湍流-化学相互作用*K.P.Kunduetal.,NASALewisResearchCenter,AIAAPaperNo.98-3986,July1998第35页/共55页PDFmodelLaminarFinite-RateModelEDCmodel富燃预混合火焰–温度云图第36页/共55页

IFRF旋涡煤粉火焰IFRF工业生产规模的电炉1997年，Peters和Weber进行的模拟仿真

，“建立2.4兆瓦的旋涡煤粉火焰数学模型,”CombustionScienceandTechnology,122,131[Ref.1]Exhaust7CoolingloopsMeasurementLocationsz=0.25m,z=0.85mInletSwirlNumber=0.923FlowRate=2684kg/hrCoalflowrate=263kg/hrCoal+transportair=421kg/hr1/4Geometry第37页/共55页网格三维周期的四分之一几何模型非结构六面体网格自适应前7万单元对靠近入口区域自适应后有26万单元最大的网格扭曲度是0.53第38页/共55页气相燃烧模拟ContoursofMeanTemperature(K)第39页/共55页离散相GottelbornhvBb煤起始(重量%,干燥)挥发物55.0,固定碳36.7,灰8.3最终(重量%,daf)80.36C,5.08H,1.45N,0.94S,12.17O.较低的生热值(LCV,MJ/kgdaf)挥发物32.3,烧焦32.9Rosin-Rammler大小分布最小1μm,最大300μm,平均45μm,带宽1.36单比例devolatization模型A=2×105s-1E=7.4×107J/kmolKinetics/diffusion-limited表面燃烧随机轨道模型每步DPM迭代21600轨迹,10个粒子大小每步DPM迭代中进行25种气相迭代第40页/共55页温度/组分结果平均温度(K)平均NO(ppm,dry)第41页/共55页GELM-1600燃气涡轮燃烧室CourtesyNovaResearchandTechnologyCorp.,Calgary,Canada非预混合,天然气12.8MW,19:1压力比(全负荷)环形燃烧室,18个喷嘴SwirlvanesFuelinletnozzlesDilutionairinlets第42页/共55页气相燃烧模拟ContoursofMeanTemperature(K)第43页/共55页NOx

模拟ThermalandfuelNOxContoursofMeanNO(ppm,wet)第44页/共55页NOx

模拟NO出口流量vs.燃烧室负荷第45页/共55页D形西瓜火焰-PDF输运试验非预混合甲烷，

湍流扩散火焰D形西瓜火焰几何2352个矩形网格单元的轴对称网格湍流标准k–ε，