Fluent-表面化學(xué)反應(yīng)模擬

1、導(dǎo)入網(wǎng)格2定義求解器Soiuboi Diuon 'elhDctez dLoti CnnS-nbMtn-wi-寫 cL&enSalcdtabm fctw3fr5：jjn Cdj Abcmd8 w宦彳®包/做上m 口Pfister StJ京ivwjr-u-Ki ord 4ni_*l» !5也當(dāng) senate!,iwrsOrii ZflHt CCOdtWG = xndar ccr,drti-rf *31b JnterfacM 了的dmK Muir 二/¥ro 'Hsurfaced FLUENTFile Mesh Define Solve Adap

2、t Surface Display Report Parallel View Help3 X甌畿 中，a/，入喟 Problem SetupGeneraleneralModelsMaterialsPhasesCell Zone ConditionsBoundary ConditionsMesh InterfacesDynamic Mesh Reference ValuesSolutionSolution MethodsSolution Controls MonitorsSolution Initoalization Calculation ActivitiesRun CalculabonRes

3、ultsGraphics and AnimationsPlotsReportsMeshType® Pressure-Based C Density-BasedVelocity Formulation0 Absolute0 RelativeTime©SteadyO TransientScale| Check - Report Quality I Displ |SolverV GravityGravitational AccelerationUnits.x (m/s2) oY (m/s2) 0Z(m/s2)gflIMMeshModels3開啟能量方程ProbJem SetupG

4、eneral'lateriafcPhasesCeJI Zone ConditionsBoundary ConditionsMesh InterfacesDynamic Mesh Reference ValuesSolutiori5oiut)Dn MethodsSgJution ControlMonitorsSolution InitializationCalcJation ActivitiesRun Calculation收 suitsGraphics and ArtPretions Plots4 操作工況參數(shù) operating conditions1 操作壓力的介紹關(guān)于參考?jí)毫Φ脑O(shè)

5、定，首先需了解有關(guān)壓力的一些定義。ANSYSFLUENT中有以下幾個(gè)壓力，即Static Pressure （靜壓）、Dynamic Pressure （動(dòng)壓）與 Total Pressure （總壓）；Absolute Pressure （絕對(duì) 壓力）、Relative Pressure（參考?jí)毫Γ┡c Operating Pressure （操作壓力）。這些壓力間的關(guān)系為，Total Pressure 總壓）=Static Pressure 靜壓）+Dynamic Pressure 動(dòng)壓）；Absolute Pressure絕對(duì)壓力）=Operating Pressure操作壓力）+Gaug

6、e Pressure表壓）。其中，靜壓、動(dòng)壓和總壓是流體力學(xué)中關(guān)于壓力的概念。靜壓是測(cè)量到的壓力，動(dòng)壓是有關(guān)速度 動(dòng)能的壓力，是流動(dòng)速度能量的體現(xiàn)。而絕對(duì)壓力、操作壓力和表壓是 FLUENT入的壓力參考量，在 ANSYS FLUENT,所有設(shè)定的壓 力都默認(rèn)為表壓。這是考慮到計(jì)算精度的問題。2 操作壓力的設(shè)定設(shè)定操作壓力時(shí)需要注意的事項(xiàng)如下：對(duì)于不可壓縮理想氣體的流動(dòng)，操作壓力的設(shè)定直接影響流體密度的計(jì)算，因?yàn)閷?duì)于理想氣體而言，流動(dòng)的密度由理想氣體方程獲得，理想氣體方程中的壓力為操作壓力。對(duì)于低馬赫數(shù)的可壓縮流動(dòng)而言，相比絕對(duì)靜壓，總壓降是很小的，因此其計(jì)算精度很容易受到數(shù)值截?cái)嗾`差的影響。需

7、要采取措施來避免此誤差的形成，ANSYS FLUENT過采用表壓（由絕對(duì)壓力減去操作壓力）的形式來避免截?cái)嗾`差的形成，操作壓力一般等于流場(chǎng)中的平均總壓。對(duì)于高馬赫數(shù)可壓縮流動(dòng)的求解而言，因?yàn)榇藭r(shí)的壓力比低馬赫可壓縮流動(dòng)的大得多，所以求解過程中的截?cái)嗾`差的影響不大， 可以不設(shè)定表壓。由于ANSYS FLUENT所有需輸入的壓力都為表壓，因此此時(shí)可以將操作壓力設(shè)定為0（這樣可以最小化由于壓力脈動(dòng)而引起的誤差），使表壓與絕對(duì)壓力相等。如果密度設(shè)定為常數(shù)或者其值由通過溫度變化的函數(shù)獲得，操作壓力并沒有在計(jì)算密度的過 程中被使用。默認(rèn)的操作壓力為101325Pa操作壓力的設(shè)定主要基于兩點(diǎn)考慮，一是流動(dòng)馬

8、赫數(shù)的大小，二是密度計(jì)算方法。表格1操作壓力的推薦設(shè)置密度關(guān)系式馬赫數(shù)操作壓力理想氣體定律50.10或約等于流場(chǎng)的平均壓力理想氣體定律小于0.1約等于流場(chǎng)的平均壓力關(guān)于溫度的函數(shù)/、口壓縮不使用常數(shù)/、口壓縮不使用/、可壓縮的理想氣 體/、口壓縮約等于流場(chǎng)的平均壓力3 關(guān)于參考?jí)毫ξ恢玫脑O(shè)定對(duì)于不涉及任何壓力邊界條件的不可壓縮流動(dòng)，ANSYS FLUENT每次迭代后要調(diào)整表壓值。這個(gè) 過程通過使用參考?jí)毫ξ恢锰帲ɑ蛟撐恢酶浇┕?jié)點(diǎn)的壓力完成。因此，參考?jí)毫ξ恢锰幍谋韷簯?yīng)一 直為0。如果使用了壓力邊界條件，則不會(huì)使用到上述關(guān)系，因此參考?jí)毫ξ恢貌槐皇褂?。參考?jí)毫ξ恢媚J(rèn)為等于或接近（0, 0,

9、0）的節(jié)點(diǎn)中心位置。實(shí)際計(jì)算中可能需要設(shè)置參考?jí)?力位置到絕對(duì)靜壓已知的位置處。在Operating Conditions對(duì)話框中的Reference Pressure Locatio選項(xiàng) 組中設(shè)置新的參考?jí)毫ξ恢玫?x, y, z的坐標(biāo)即可。如果要考慮某一方向的加速度，如重力，可以勾選Gravity復(fù)選框。對(duì)于VOF計(jì)算，應(yīng)當(dāng)選擇 Specified Operating Density,并且在Operating Density下為最輕相設(shè)置 密度。這樣做排除了水力靜壓的積累，提高了 round-off精度為動(dòng)量平衡。同樣需要打開Implicit Body Force,部分平衡壓力梯度和動(dòng)量方

10、程中體積力，提高解的收斂性。Reference Pressure Location參考?jí)簭?qiáng)位置）應(yīng)是位于流體永遠(yuǎn)是100%的某一相（空氣）的區(qū)域， 光滑和快速收斂是其基本條件。單擊 Defines Operating Conditions 在 Operating Pressur仲輸入 10000 Pa,選中重力 Gravity, 在 Z 中輸入 9.81 m/s2, Operating Temperature入 303 K,點(diǎn)擊 OK確認(rèn)。5 定義多組分模型（1）在Model （模型）中選擇 Species Transport （組元輸運(yùn)）。（2）在 Reactions （反應(yīng)）中選擇 Vol

11、umetric Reactions （體積反應(yīng)）。（ 3）在Mixture Material （混合物材料）中選擇所計(jì)算問題中涉及到的反應(yīng)物，則Numberof Volumetric Species （體積組元數(shù)量）中 自動(dòng)顯示混合物中的組元數(shù)量。（4）在Turbulence-Chemistry Interaction （湍流與化學(xué)反應(yīng)相干模型）中根據(jù)需要選擇相應(yīng)的模型。如果選擇了 Eddy-Dissipation Concept （EDC）則可以進(jìn)一步修改 Volume Fraction Constant （體積濃度常 數(shù)）和Time Scale Constant（時(shí)間尺度常數(shù)）（6）如果想

12、完整計(jì)算多組分的擴(kuò)散或熱擴(kuò)散，就選中Full Multicomponent Diffusion （完整多組分?jǐn)U 散）和Thermal Diffusion （熱擴(kuò)散）選項(xiàng)。在上面的設(shè)置過程中,如果需要查看混合物中組分和化學(xué)反應(yīng)的相關(guān)設(shè)置，可以在 Species （組 元）面板中，點(diǎn)擊Mixture Material （混合物材料）右邊的 View （觀看）按鈕。如果計(jì)算中用到的混 合物模型是一種新的?M合物，則需要在Material （材料）面板中創(chuàng)建混合物，然后再將新定義的混合物選作計(jì)算用的混合物?；旌衔锏亩x過程包含組分選取、反應(yīng)模型設(shè)定、反應(yīng)機(jī)制設(shè)定等幾個(gè)步驟，下面逐一介紹。層流條件下，M

13、odel 模型中只能選擇Species Transport 組元運(yùn)輸一項(xiàng)。Turbulence-ChemistryInteraction （湍流與化學(xué)反應(yīng)相干模型）中也只有一項(xiàng)。ftcHEtn SedtpMu配低Mzd-Hs.勺S址PhasesCel Zan± 5F ±ie 必 C Jidlo yiWMi】世產(chǎn)£必 fbmc Mr+iRen-KT 3uesSouHonSdL/Xo MeItoce5lL 'Jun Cen 口出Mft-iitanfSdL-aan InrUimbor CakiJ口Ucn Actft -ti-s 即 CBfcUiaFC0回的sah

14、cs ao: Ama sort 詠即K0一1心二4中pjja 門Ff匚 Sipec m ModtlLfcpha&e - or Erww-C"而叫-L on in JrNadwCon Dff nmt EyEmce - EANSYS FLUENT UOptl j在 Models 中選中 Species Transport選按鈕，在 Reactions中選中 Volumetric 和 Wall suface復(fù)選 框，在 Wall Surface Reaction Options 中選中 Mass Deposition Source復(fù)選框，在 Options 中選擇 Inlet Di

15、ffusion, Full Multicomponent Diffusion , Thermal Diffusion 復(fù)選框，單擊 OK按鈕確認(rèn)。6設(shè)置材料1）添加神化氫arsine1,雙擊air2, 在 name 中輸入 arsine 和 chemical formula 處輸入 ash3; 比熱容 cp 選擇 kinetic-theory ； 導(dǎo) 熱系數(shù) thermal conductivity 選擇 kinetic-theory ；粘度 viscosity 選擇 kinetic-theory ；分子量 molecular weight 選擇 constant 為 77.95;標(biāo)準(zhǔn)狀態(tài)始

16、standard state enthalpy 為 0;標(biāo)準(zhǔn)狀態(tài)嫡 standard state entropy 為 130579.1;基準(zhǔn)溫度 reference temperature 為 298.15。no。3,點(diǎn)擊change/create,創(chuàng)建新物質(zhì)，在彈出的是否覆蓋選擇259.8。4,雙擊 arsine在 L-J characteristic Length#征長(zhǎng)度輸入 4.145 （埃米）；L-J Energy Paramete能量參數(shù)中輸入 單擊 Change/Create按鈕ar卜ar5Pv*aterp*r Fluid2）添加三乙基錢、甲基、氫氣、錢（固體）、神（固體）、錢、神

17、參數(shù)Ga（CH_3）_ 三乙式春3CH_3 甲優(yōu)H_2氫氣Ga_s 錢As_s江4Gat （周）As神（固）Nametmgch3ghydro-genga_sas_sgaasChemical Formulagach33ch3h2ga_sas_sgaasCP（specific heat）kinetic-theo rykinetic-theo rykinetic-theo ry520.64520.641006.431006.43Thermal conductivit ykinetic-theo rykinetic-theo rykinetic-theo ry0.01580.0158kinetic-t

18、heo rykinetic-theo ryViscositykinetic-theo rykinetic-theo rykinetic-theo ry2.125e-0 52.125e-0 5kinetic-theo rykinetic-theo ryMolecular weight114.83152.0269.7274.9269.7274.92Standard state enthalpy02.044e+0703117.713117.7100Standard state entropy130579.1257367.6130579.1154719.3154719.300reference tem

19、peratur e298.15298.15298.15298.15298.15298.15298.15L-J characterist ic Length5.683.7582.827-00L-J Energy Parameter398148.659.7-00Degree of freedom005-7編輯組信息1,修改組命名雙擊 mixture-temple , name 中輸入 gaas_deposition 單擊 change, 點(diǎn)擊 yes確認(rèn)。2,編輯組分信息在 mixture species 中組分選取首先進(jìn)入 Materials （材料）面板：Define->Material

20、s.在Materials （材料）面板上，先在Name （名稱）中為新的混合物確定一個(gè)名稱，然后在Material Type（材料類型）里選擇mixture （混合物）。如果有與目標(biāo)相近的混合物模型， 可以在下面的Mixture Material （混合物材料）中選擇一樣，比如 methane-air （甲烷一空氣）,然后在下面 Properties （性 能）中做詳細(xì)設(shè)置，即按順序設(shè)置組元、反應(yīng)類型、反應(yīng)機(jī)制等等：（1）點(diǎn)擊 Mixture Species （混合物組元）右邊的 Edit （編輯）按鈕進(jìn)入 Species （組元）面板， 如圖7-11所示。在 Mixture （混合物）下面有

21、4個(gè)框，即Available Materials （可用材料），Selected Species（已選組元），Selected Site SpecieS已選吸收組元）和Seleted Solid SpecieS已選固體組元）。Available Materials （可用材料）是指材料數(shù)據(jù)庫中可供選用的材料；Selected Species （已選組元）是指當(dāng)前混合物中已經(jīng)選中的組元；Selected Site Species已選吸收組元）是指在存在物面反應(yīng)的計(jì)算中，氣相混合物中即將通過反應(yīng)被物面吸收的組元；Selected Solid Species（已選固體組元）是指物面反應(yīng)計(jì)算中將從物面

22、進(jìn)入氣流的組元。顯然，如果不存在物面反應(yīng)，則不用考慮后面兩個(gè)方 框中的內(nèi)容，整個(gè)設(shè)置過程會(huì)大大簡(jiǎn)化。組元設(shè)置的中心任務(wù)是選擇混合物組元， 即設(shè)定Selected Species已選組元）的內(nèi)容。在Materials （材料）面板中，點(diǎn)擊 Database（數(shù)據(jù)庫）按鈕打開 Database Materials （數(shù)據(jù)庫中材料）面板， 拷貝所需的組元后，再回到Species（組元）面板。在拷貝之前，需要確認(rèn)Material Type （材料類型） 必須是fluid （流體）。這里不必?fù)?dān)心如何添加物面反應(yīng)中參與反應(yīng)的固體組元，因?yàn)閷?shí)際上這些組元 也會(huì)出現(xiàn)在fluid （流體）的列表中。在Selec

23、tedSpecies（已選組元）中，最后一項(xiàng)必須是質(zhì)量濃度最大的一個(gè)組元 。如果最后一項(xiàng) 不是質(zhì)量濃度最大的一項(xiàng)，可以先將質(zhì)量濃度最大的組元從方框中刪除， 然后再重新添加進(jìn)來，以保 證這個(gè)組元處于方框的底部。添加和刪除操作是通過選擇相應(yīng)的組元，然后點(diǎn)擊Add （添加）和Remove （刪除）按鈕實(shí)現(xiàn)的。本例中點(diǎn)擊mixture specie的單擊edit按鈕彈出species（組分）對(duì)話框，調(diào)整各類型組分Selected SpeciesSelected Site SpeciesSeleted Solid Speciesash3ga_sgaga(ch3)3as sasch3h2注：這里需要設(shè)置好

24、selected species的排歹1順序，在后面邊界條件設(shè)置中 velocity inlet中species 體現(xiàn)的三個(gè)組分，是按順序排列的前三個(gè)化學(xué)反應(yīng)設(shè)定組元設(shè)置完成后，就可以開始設(shè)置組元間的化學(xué)反應(yīng)。在 Materials （材料）面板中，Reactions （反應(yīng)）下拉列表中顯示的反應(yīng)類型取決于Species Modal （組元模型）面板中Turbulence-ChemistryInteraction （湍流一化學(xué)反應(yīng)相干）模型的設(shè)置如果設(shè)置的是 Laminar Finite-Rate（層流有限速率）模型或EDC模型，則反應(yīng)類型顯示為 finite-rate （有限速率）；如果設(shè)置

25、的是 Eddy-Dissipation （渦 耗散）模型，則反應(yīng)類型顯示為 eddy-dissipation（渦耗散）；如果設(shè)置的是Finite-Rate/Eddy-Dissipation （有限速率/渦耗散）模型，則反應(yīng)類型顯示為 finite-rate/eddy-dissipation （有限速率/渦應(yīng)散）。點(diǎn)擊 右端的Edit （編輯）按鈕進(jìn)入 Reactions,（反應(yīng)）面板（如圖 7-13所示），在Reactions （反應(yīng)）面 板中完成對(duì)化學(xué)反應(yīng)模型的設(shè)置。化學(xué)反應(yīng)的設(shè)置主要包括下列幾項(xiàng)內(nèi)容：1）在Total Number of Reactions （總的反應(yīng)數(shù)量）中設(shè)定總的反應(yīng)數(shù)

26、量；2）在Reaction Name （反應(yīng)名稱）中指定反應(yīng)名稱；3）在Reaction ID （反應(yīng)編號(hào)）中指定每個(gè)反應(yīng)的編號(hào)；4）在Reaction Type（反應(yīng)類型）中指定反應(yīng)類型，即指定反應(yīng)類型為 Volumetric （體積）、Wall Surface （物面）或 Particle Surface （顆粒表面）；5）在 Number of Reactants （反應(yīng)物數(shù)量）和 Number of Products （生成物數(shù)量）中指定反應(yīng)的 反應(yīng)物數(shù)量和生成物數(shù)量。然后在Species （組元）下拉菜單中選擇反應(yīng)物和生成物，并在 Stoich.Coefficient （反應(yīng)系數(shù)）中

27、設(shè)定組元在指定反應(yīng)方程中的系數(shù)。在 Rate Exponent （速率指數(shù)）中設(shè)定 組元的速率指數(shù)，即組元生成速率方程中指定反應(yīng)方程中摩爾濃度項(xiàng)的指數(shù)。6）如果在 Species Modal（組元模型）面板的 Turbulence-Chemistry Interaction （湍流反應(yīng)相干 模型）中選擇的是 laminar finite-rate （層流有限速率）、finite-rate/eddy-dissipation （有限速率/渦擴(kuò) 散）或EDC模型時(shí)，則需要根據(jù)反應(yīng)模型設(shè)置Arrhenius Rate （ Arrhenius速率）下面的選項(xiàng)。由于篇幅所限，這里不再詳述。7）根據(jù)實(shí)際反應(yīng)

28、過程，確定是否選擇Include Backward Reaction （包含逆向反應(yīng)）選項(xiàng)。8）如果使用 eddy-dissipation （渦擴(kuò)散）或 finite-rate/eddy-dissipation （有限速率/渦擴(kuò)散）模型， 則還需要設(shè)定Mixing Rate （混合速率）。9）重復(fù)2）8）步直到設(shè)置完所有反應(yīng)，然后點(diǎn)擊 OK按鈕完成全部設(shè)定過程。反應(yīng)模型設(shè)定在區(qū)域上定義反應(yīng)機(jī)制本例中，在Reaction旁單擊edit按鈕，在Reaction對(duì)話框，Number of Reactants中輸入2,化學(xué) 反應(yīng)輸入數(shù)據(jù)如下表。參數(shù)反應(yīng)一反應(yīng)二Reaction namegallium-

29、depArsenic-depReaction ID12Reaction typeWall surfaceWall surfaceNumber of reactants22Speciesash3, ga sgach33, as sStoich. Coefficientash3=1, ga s=1gach33=1, as s=1Rate exponentash3=1, ga s=1gach33=1, as s=1Arrhenius RatePEF=1e+06 AE=Q TE=0.5PEF=1e+12 AE=0, TE=0.5Number of products33Speciesga, as_s,

30、h2as, ga_s, ch3Stoich. Coefficientga=1, as_s=1, h2=1.5as=1, ga_s=1, ch3=3Rate exponentas_s=0, h2=0ga_s=0, ch3=0Create/Ed t M-ateielsOr dei MflLe,iafa byQ Namerhfflitral =ofi 隔FLUBMT Dataoase.User-Defir ed D*taM5«.noneMi:<tLre SpeqesRcozborfcdata rBa suf1 一 been added. inductivitynamesloaded.

31、uid.(qaas_depojs btcn Ndd StdLe EnLlidPEF=Pre-Exponential FactoAE=Activation Energ yTE=Temperature Exponent反應(yīng)機(jī)制設(shè)定FLUENT中提到的“反應(yīng)機(jī)制（reaction mechanisms） ”指的是局限在特定區(qū)域中的化學(xué)反應(yīng)?！胺?應(yīng)機(jī)制”中涉及的反應(yīng)是前面設(shè)定的化學(xué)反應(yīng)的子集。在Materials （材料）面板中，點(diǎn)擊Mechanisms（機(jī)制）旁邊的Edit （編輯）按鈕，可以打開 Reaction Mechanisms（反應(yīng)機(jī)制）面板，如圖 7-14所 示。具體步驟如下：1、在N

32、umber of Mechanisms （反應(yīng)機(jī)制數(shù)量）中設(shè)定反應(yīng)的數(shù)量。2、設(shè)定 Mechanisms ID （反應(yīng)機(jī)制編號(hào)）。3、設(shè)定 Name （名稱）。4、在Reaction Type （反應(yīng)類型）中設(shè)定反應(yīng)類型。在反應(yīng)類型確定后，屬于這個(gè)類型的反應(yīng)就 會(huì)出現(xiàn)在Reactions （反應(yīng)）列表中。5、選定反應(yīng)機(jī)制中包含的反應(yīng)。如果選定的反應(yīng)類型是Wall Surface（壁面反應(yīng)）并且其中包含吸收反應(yīng)，則還需要對(duì)吸收反應(yīng)進(jìn)行專門的設(shè)定， 即設(shè)定Number of Sites（吸收反應(yīng)數(shù)量）、Site Name （吸收反應(yīng)名稱）、Site Density （吸收密度），點(diǎn)擊Define

33、（定義）按鈕還可以進(jìn)一步選定吸收反應(yīng) 中的被吸收組元和被吸收組元的 Initial Site Coverage （初始吸收覆蓋率）等參數(shù)。本例中，在 mechanism 旁單擊 Edit按鈕彈出，reaction mechanisms對(duì)話框。在 Number of Mechanisms 中輸入 1, Name 中輸入 gaas-ald, Reaction Type 選擇 Wall Surface, Reactions 中選擇gallium-dep 和 arsenic-dep, Number of Sites 中輸入 1, Site Density 中輸入 1e-08。單擊 define 按鈕，

34、七 Site parameters對(duì)話框，Total Number of Site Species中輸入 2, Initial SiteCoverage中 ga_s輸入 0.7, as_s輸入 0.3。單擊 Apply 確認(rèn)。X H'EsJyb ID "口 gCflX即川F E ：j I nnpresstae-idcai-gsfi如11nh1%Ita ter l ai nd x ture - tenpldte iH?w tH vpertyi "'MesliMnd當(dāng)rf， H白， i nq 產(chǎn)中廣亡Id on iwch; wtHew praperty '

35、;"ih?rmal Dandiue Selecting cnnstnit iwtHod +n l+t*w |jrii|itrly MUlhcuhlLy1' hdh Stltcting coi>st*n>t rwthod foIla ter lai jir:ji -i<!*- i b -jibe Nm.UejaLr 什21Cm cd |Hdt> |iH*w prapertjji MStanOaird StittXlvcting Ggst-wt iwthnd for 4M Stvtir Fntmlpy韭 - drti mqui1PHiError i Ha

36、user-Oeff inedl I unctions Hwie been IsatedLError 仙ject; *Error i d e le tt- material i still in use tiy fl mid.Error IJbjBct： riijdtuire-tenpldite在 導(dǎo)熱性 thermal conductivity 選 擇 mass-weighted-mixing-law , 粘 度 viscosity 選 擇 mass-weighted-mixing-law。在計(jì)算多組元混合物時(shí)，可以在Thermal Conductivity （導(dǎo)熱系數(shù)）右邊的下拉列表中選擇m

37、ass-weighted-mixing-law （質(zhì)量加權(quán)混合律）定義導(dǎo)熱系數(shù)。如果混合物為理想氣體，則應(yīng)該選擇 ideal-gas-mixing-law （理想氣體混合律）。如果用理想氣體混合律計(jì)算混合物的導(dǎo)熱系數(shù)，則必須同時(shí)使用ideal-gas-mixing-law （理想氣體混合律）或者 mass-weighted-mixing-law （質(zhì)量加權(quán)混合律）計(jì)算粘性，因?yàn)橹挥羞@兩種方法可以用 于指定多組元混合物的粘性多組元混合物的導(dǎo)熱系數(shù)也可以用UDF （用戶自定義函數(shù)）來定義，其選項(xiàng)為user-defined （用戶定義）或user-defined-mixing-law （用戶定義混合

38、律），詳細(xì)設(shè)置方法請(qǐng)參見 UDF的相關(guān)章節(jié)在FLUENT中定義組元相關(guān)粘度的步驟如下：在Viscosity （粘度）右邊的下拉列表中選擇 mass-weighted-mixing-law （質(zhì)量加權(quán)混合律）。如果 在密度定義中使用了理想氣體定律，則選擇ideal-gas-mixing-law （理想氣體混合定律）。如果準(zhǔn)備用自定義函數(shù)定義粘度，則選擇 user-defined （用戶自定義），或選擇 user-defined-mixing-law （用戶自 定義混合律）。user-de巾ned-mixing-law （用戶自定義混合律）是采用 UDF來定義混合物的粘度，其基本定義方 法也是先分

39、別定義混合物中的各種組元，再用各組元的物性參數(shù)計(jì)算出混合物的粘度。user-defined（用戶自定義）則是以混合物粘度參數(shù)為對(duì)象進(jìn)行定義的。相關(guān)內(nèi)容請(qǐng)參見UDF相關(guān)章節(jié)。8邊界條件設(shè)定速度入口設(shè)定Velocity magnitude 為 0.02189; thermal 中 temperature 為 293; species中 Species Mass Fractions （組元質(zhì)量濃度）ash3為0.4, gach33為0.15, ch3為0。Boundary ConditionsZonedefadt-interior oudet1: Meshvelcotynnlet壓力出口設(shè)置 保持默

40、認(rèn)ounclry Cofidlliomi二已一三1.什一|-|二三口匚0.elosLvdntetMIT rjal-2MSWE * (TG&IMnEWn Therrm Raciaiiar 通君 | CFM Nulfipbl U弱 |Qu弊teawri出后scaR斤|皿七dcfl皿DkrtonSp"同MW25心即1田IEqjrfibnun Presiare 口51nbJition口婦.聞己 Pressure Soetilceta-TyaedmturFv 國中eqj¥fuMt .可Ta qet 附年 Fh' Rfite-6dt | Qyy profileaw F/

41、iinTeiis.”cjper4?ng candflcnt 3iEpa, Meah Pencdc Z:mi：n3nsB.LEnkCxd-Jx -7.壁面設(shè)置wall-1Thermal conditions 中選擇 temperature 單選按鈕，填入 473壁面設(shè)置wall-2Thermal cond幣ons中選擇temperature單選按鈕，填入 343Boundary ConditionsZoneAdjacent Cell ZonewbI-2defeult-intericr outet velodty-nlet wallZone PlameCarN壁面設(shè)置wall-4Wall moti

42、on 中選擇 moving wall 單選按鈕，在 motion 中選擇 absolute 和 rotational 按鈕，speed 輸入80。Thermal conditions中選擇temperature單選按鈕，填入1023。Species選項(xiàng)卡中保持默認(rèn)。Boundary Conditionswal 4v/al-5Zone loarreJone defeuh-iniefiof outtet vcbot/rict wal-1 wal-3 (Acjocent Cel Zone|Iidhfaren卻m lhermd| Ration | Speers CPM | MJfphase| UD5

43、V/dlHh|Eimr。-pararretOcpavM -lipWd MotonSstonar, Tai。gingwail( Tr arise tonal。Rotator alCccjmtsC Reteu * tc y acav U m AbxjluteShwt Corddon。No SipC Spccited Shear Specuiflf t/Coeftclent Mxa-wi Sbeaswal Rougm«sme2gH (enj p-|axi5taitRCi jJtk-ss rK>nt (ccnstantZone Hanev/al-4Adjacent CM Zon

44、1;OK | .三金口 help |壁面設(shè)置wall-5Wall motion 中選擇 moving wall 單選按鈕，在 motion 中選擇 absolute 和 rotational 按鈕，speed 輸入80。Thermal conditions中選擇temperature單選按鈕，填入 720。Species選項(xiàng)卡中保持默認(rèn)。defdulUnteftof oude: veloctyTe: ivdl-1 Zl2WTivdl S2 Wall“y Flore gl5Adjacent Cel Zone lludMuTcrtm | Therna Radianor eoes jDPM I Mu

45、lphase | UC6 v/al Hh- |Phasentre7al Motion。Statorar Wall9 Mo/rg,端MotionfcdiL.Rpat<» ro AdjT*ir cel Zsr 。A55*tcSoeec (rad/sm"2c膽XM5Trancloti&nai 。Rota5onaComponentsRotd-ton-AJCS&UnX (an)RotMbn-AxJs arectkrY(on)2(6)XpYP平5lPISheer CottonO No Sip5peafirdShroiSpecuhrity CoeffcentMa,a

46、ngon STessWal RoughncwconstantRojo：55 卜。巾(CTi)府pojghrwe conctrr econstant壁面設(shè)置wall-6Thermal conditions 中選擇 temperature 單選按鈕，填入 3039重新設(shè)置modelOptions中取消選擇inlet diffusion，點(diǎn)擊OK確認(rèn)。如果入口質(zhì)量很小的話就應(yīng)該關(guān)掉入口擴(kuò)散選項(xiàng)的，否則會(huì)導(dǎo)致質(zhì)量損失A oblkm 5ekipGeneralMa*北Cd Zone ConditMna Dcunoer Candildhns 114Kh nt=，世群 D/manic Fesn Referen

47、ce ValuerSdutonSolution M?thockSoljiiM Csh7dIb Monitors5&uUan Lf itidlizdljor CefciiLation Acnm. 艮 q CacjabonkcsjtsG-achcs art Arina lor* PbsReportsModelsModelshMtidioK - OffEbeqy - OnLannarRodetKn OffH-Cft.Tr:rjOx - off sox - off Soor-OfFOeroiried Detain:Rescmg diamel V 1*1® Phase - off SJ

48、irliiranni ft *E tousUcsOfr Bj-ETian Wal Fta -fl才p .AMedelOff Sceoes transport Man-FreriM=d CDrrbuEJfln Premixes 匚，m4£tcn iartaily remixed Cowustiar C3irmsiQn 1soic "BiarispgRcKtnns0 聞1/11±5E WallSLrtBre1 力 rtoa Suraca-10求解控制Solution methods 保持默認(rèn)ixtre Prcpcrie£Njmber of VLmebic G

49、peoes找M.K3。疵的Vj'all Surface Reaction Optans21 Heiat of Surface Rea曲m/I Mass Deresitign SwrzcAg9ries&venesB Factor 后OfticraI trilel： DiflLisonIZ fgon trergr source回 Fil Mufazonpi&nert LrftufiionMl Themall Di FusonIdPiivd 匚quiibr UT1 S-lfFCheiilstiv Sol/er| LHEMCIN'Cro iron pFJWTifln Iwignipf叫Kirber oF Sold Speoco三曲Numb£?r d Si x 部emo hTurbulence-C he ntslrf Intact on* Lenirwr Finite-R«teFinTp-Rrtp/Md/-risi