comsol流體仿真(層流-傳質(zhì)耦合場(chǎng))(對(duì)流擴(kuò)散)

1、流體的對(duì)流擴(kuò)散仿真報(bào)告日期2016-4-5 11:27:37目錄1.全局定義1.1.參數(shù) 12.組件 12.1.定義2.2.幾何 12.3.材料2.4.多孔介質(zhì)反應(yīng)流2.5.網(wǎng)格 13.研究 13.1.參數(shù)化掃描3.2.穩(wěn)態(tài)3.3.求解器配置4.結(jié)果4.1.數(shù)據(jù)集4.2.繪圖組1 全局定義日期Mar 18, 2016 4:27:52 PM全局設(shè)定名稱H NewModel.mph路徑E:COMSOLh_NewModel.mphCOMSOL 版本COMSOL 5.2 (建立: 220)單位系統(tǒng)SI使用的模塊COMSOL MultiphysicsCAD Import ModuleChemical R

2、eaction Engineering Module1.1 參數(shù) 1參數(shù)名稱表達(dá)式值描述D_media2e-9m2/s2E9 m²/sD_pdms4e-9m2/s4E9 m²/sD_gas2e-5m2/s2E5 m²/su_media1e-3Pa*s0.001 Pa·su_gas1e-51E5p_media1e3kg/m31000 kg/m³p_gas1kg/m31 kg/m³v_gas4m/s4 m/sh_in3mm0.003 mR_in8.31J/mo

3、l/K8.31 J/(mol·K)T_in310K310 KV_standerd1m31 m³P_standerd1e5Pa1E5 Pan0_airP_standerd*V_standerd/R_in/T_in38.818 molc0n0_air/V_standerd38.818 mol/m³h_channal0.16mm1.6E4 mMn_Air0.21*32 + 0.74*28 + 0.05*4429.64空氣的平均摩爾質(zhì)量VbO225.6cm3/mol2.56E5 m³

4、/molVbN231.2cm3/mol3.12E5 m³/molVbCO234cm3/mol3.4E5 m³/molalpha_N2(VbO2/VbN2)0.60.88808 rad威爾基-張公式alpha_CO2(VbO2/VbCO2)0.60.84344 radc10c0*0.218.1519 mol/m³c20c0*0.7428.726 mol/m³c30c0*0.051.9409 mol/m³u_pdmsu_media*1e5100 kg/(m·s

5、)2 組件 1日期Mar 18, 2016 3:13:08 PM組件設(shè)定單位系統(tǒng)SI幾何形參階次automatic2.1 定義2.1.1 選擇顯式 1選擇類型Explicit選擇沒(méi)有邊顯式 12.1.2 坐標(biāo)系邊界坐標(biāo)系 1坐標(biāo)系類型邊界坐標(biāo)系標(biāo)記sys1坐標(biāo)名稱第一軸第二軸第三軸t1t2n2.2 幾何 1幾何 1單位長(zhǎng)度單位µm角度單位deg幾何統(tǒng)計(jì)描述值空間維度3域數(shù)5邊界數(shù)50邊數(shù)116端點(diǎn)數(shù)722.2.1 工作平面 1 (wp1)聯(lián)合對(duì)象描述值聯(lián)合對(duì)象開(kāi)平面幾何 (wp1)正方形 1 (sq1)位置描述值位置0, 0基準(zhǔn)中心尺寸描述值邊長(zhǎng)22000.02.2.2 拉伸 1 (

6、ext1)設(shè)定描述值工作平面工作平面 1到面的距離距離 (µm)-h_in比例比例 xw比例 yw11位移位移 xw (µm)位移 yw (µm)00扭轉(zhuǎn)角扭轉(zhuǎn)角 (deg)02.2.3 工作平面 2 (wp2)聯(lián)合對(duì)象描述值聯(lián)合對(duì)象開(kāi)平面幾何 (wp2)多邊形 1 (pol1)對(duì)象類型描述值類型實(shí)體坐標(biāo)描述值數(shù)據(jù)源表格坐標(biāo)xw (µm)yw (µm)45011000.045035001503500150-3500350.0-3500350.0-11000.0750-110000-3500950.0350065035

7、0065011000.045011000.0多邊形 2 (pol2)對(duì)象類型描述值類型實(shí)體坐標(biāo)描述值數(shù)據(jù)源表格坐標(biāo)xw (µm)yw (µm)115011000.0175011000.01750-11000.01150-11000.0115011000.02.2.4 拉伸 2 (ext2)設(shè)定描述值工作平面工作平面 2到面的距離距離 (µm)h_channal到面的距離描述值反向開(kāi)比例比例 xw比例 yw11位移位移 xw (µm)位移 yw (µm)00扭轉(zhuǎn)角扭轉(zhuǎn)角 (deg)02.2.5 鏡像 1 (mir1)設(shè)定描述值保留輸入對(duì)象開(kāi)反射面

8、上的點(diǎn)描述值面內(nèi)點(diǎn)0, 0, 0反射面法向矢量描述值平面法線1, 0, 02.2.6 工作平面 3 (wp3)聯(lián)合對(duì)象描述值聯(lián)合對(duì)象開(kāi)平面幾何 (wp3)貝氏多線段 1 (b1)多邊形線段描述值控制點(diǎn)-1750, -1149.9999761580998, -1149.9999761580998, -1750, -1750, 11000.0, 11000.0, -11000.0, -11000.0, 11000.0度1, 1, 1, 1權(quán)重1, 1, 1, 1, 1, 1, 1, 1類型實(shí)體貝氏多線段 2 (b2)多邊形線段描述值控制點(diǎn)-750, -349.99999403954, -349.9

9、9999403954, -150.00000596046002, -150.00000596046002, -449.99998807906996, -449.99998807906996, -649.9999761581399, -649.9999761581399, -949.99998807907, -949.99998807907, -750, -750, -11000.0, -11000.0, -3500, -3500, 3500, 3500, 11000.0, 11000.0, 3500, 3500, -3500, -3500, -11000.0度1, 1, 1, 1, 1, 1,

10、 1, 1, 1, 1, 1, 1權(quán)重1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1類型實(shí)體貝氏多線段 3 (b3)多邊形線段描述值控制點(diǎn)1149.9999761580998, 1750, 1750, 1149.9999761580998, 1149.9999761580998, 11000.0, 11000.0, -11000.0, -11000.0, 11000.0度1, 1, 1, 1權(quán)重1, 1, 1, 1, 1, 1, 1, 1類型實(shí)體貝氏多線段 4 (b4)多邊形線段描述值控制點(diǎn)449.9

11、9998807906996, 649.9999761581399, 649.9999761581399, 949.99998807907, 949.99998807907, 750, 750, 349.99999403954, 349.99999403954, 150.00000596046002, 150.00000596046002, 449.99998807906996, 449.99998807906996, 11000.0, 11000.0, 3500, 3500, -3500, -3500, -11000.0, -11000.0, -3500, -3500, 3500, 3500,

12、 11000.0度1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1權(quán)重1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1類型實(shí)體2.2.7 channal (ext3)設(shè)定描述值工作平面工作平面 3到面的距離距離 (µm)-h_channal比例比例 xw比例 yw11位移位移 xw (µm)位移 yw (µm)00扭轉(zhuǎn)角扭轉(zhuǎn)角 (deg)02.3 材料2.3.1 pdmspdms選擇幾何實(shí)體層次域選擇域 1材料參數(shù)名稱值單位密度p_mediakg/m3動(dòng)

13、力粘度u_pdmsPa*s基本 設(shè)定DescriptionValue密度p_media動(dòng)力粘度u_pdms2.3.2 gasgas選擇幾何實(shí)體層次域選擇域 2, 5材料參數(shù)名稱值單位密度p_gaskg/m3動(dòng)力粘度u_gasPa*s基本 設(shè)定DescriptionValue密度p_gas動(dòng)力粘度u_gas2.3.3 mediamedia選擇幾何實(shí)體層次域選擇域 34材料參數(shù)名稱值單位密度p_mediakg/m3動(dòng)力粘度u_mediaPa*s基本 設(shè)定DescriptionValue密度p_media動(dòng)力粘度u_media2.4 多孔介質(zhì)反應(yīng)流使用的模塊COMSOL MultiphysicsCh

14、emical Reaction Engineering Module多孔介質(zhì)反應(yīng)流選擇幾何實(shí)體層次域選擇域 15EquationsSettingsDescriptionValue濃度線性計(jì)算邊界通量開(kāi)應(yīng)用平滑到邊界通量開(kāi)流體離散化P1 + P1使用分裂復(fù)數(shù)變量時(shí)的值類型實(shí)數(shù), 實(shí)數(shù), 實(shí)數(shù)變量名稱表達(dá)式單位描述選擇rfds.nxnx法矢，x 分量邊界 150rfds.nyny法矢，y 分量邊界 150rfds.nznz法矢，z 分量邊界 150rfds.nxmeshroot.nxmesh法矢（網(wǎng)格），x 分量邊界 150rfds.nymeshroot.nymesh法矢（網(wǎng)格），y 分量邊界 1

15、50rfds.nzmeshroot.nzmesh法矢（網(wǎng)格），z 分量邊界 150rfds.R_c10mol/(m3*s)總速率表達(dá)式域 15rfds.R_c20mol/(m3*s)總速率表達(dá)式域 15rfds.R_c30mol/(m3*s)總速率表達(dá)式域 15domflux.c1xrfds.dfluxx_c1mol/(m2*s)域通量 x 分量域 15domflux.c1yrfds.dfluxy_c1mol/(m2*s)域通量 y 分量域 15domflux.c1zrfds.dfluxz_c1mol/(m2*s)域通量 z 分量域 15domflux.c2xrfds.dfluxx_c2mol

16、/(m2*s)域通量 x 分量域 15domflux.c2yrfds.dfluxy_c2mol/(m2*s)域通量 y 分量域 15domflux.c2zrfds.dfluxz_c2mol/(m2*s)域通量 z 分量域 15domflux.c3xrfds.dfluxx_c3mol/(m2*s)域通量 x 分量域 15domflux.c3yrfds.dfluxy_c3mol/(m2*s)域通量 y 分量域 15domflux.c3zrfds.dfluxz_c3mol/(m2*s)域通量 z 分量域 15rfds.nxcroot.nxc/sqrt(root.nxc2+root.nyc2+root

17、.nzc2+eps)法矢，x 分量邊界 150rfds.nycroot.nyc/sqrt(root.nxc2+root.nyc2+root.nzc2+eps)法矢，y 分量邊界 150rfds.nzcroot.nzc/sqrt(root.nxc2+root.nyc2+root.nzc2+eps)法矢，z 分量邊界 150rfds.bndFlux_c1-dflux_spatial(c1)mol/(m2*s)邊界通量邊界 15, 7, 910, 12, 16, 19, 21, 24, 31, 34, 36, 39, 44, 4647, 4950rfds.bndFlux_c10.5*(uflux_s

18、patial(c1)-dflux_spatial(c1)mol/(m2*s)邊界通量邊界 6, 8, 11, 1315, 1718, 20, 2223, 2530, 3233, 35, 3738, 4043, 45, 48rfds.ntflux_c1rfds.bndFlux_c1+rfds.cfluxx_c1*rfds.nxc+rfds.cfluxy_c1*rfds.nyc+rfds.cfluxz_c1*rfds.nzcmol/(m2*s)法向總通量邊界 150rfds.ndflux_c1rfds.bndFlux_c1mol/(m2*s)法向擴(kuò)散通量邊界 150rfds.ncflux_c1rf

19、ds.cfluxx_c1*rfds.nxc+rfds.cfluxy_c1*rfds.nyc+rfds.cfluxz_c1*rfds.nzcmol/(m2*s)法向?qū)α魍窟吔?150rfds.bndFlux_c2-dflux_spatial(c2)mol/(m2*s)邊界通量邊界 15, 7, 910, 12, 16, 19, 21, 24, 31, 34, 36, 39, 44, 4647, 4950rfds.bndFlux_c20.5*(uflux_spatial(c2)-dflux_spatial(c2)mol/(m2*s)邊界通量邊界 6, 8, 11, 1315, 1718, 20,

20、 2223, 2530, 3233, 35, 3738, 4043, 45, 48rfds.ntflux_c2rfds.bndFlux_c2+rfds.cfluxx_c2*rfds.nxc+rfds.cfluxy_c2*rfds.nyc+rfds.cfluxz_c2*rfds.nzcmol/(m2*s)法向總通量邊界 150rfds.ndflux_c2rfds.bndFlux_c2mol/(m2*s)法向擴(kuò)散通量邊界 150rfds.ncflux_c2rfds.cfluxx_c2*rfds.nxc+rfds.cfluxy_c2*rfds.nyc+rfds.cfluxz_c2*rfds.nzcm

21、ol/(m2*s)法向?qū)α魍窟吔?150rfds.bndFlux_c3-dflux_spatial(c3)mol/(m2*s)邊界通量邊界 15, 7, 910, 12, 16, 19, 21, 24, 31, 34, 36, 39, 44, 4647, 4950rfds.bndFlux_c30.5*(uflux_spatial(c3)-dflux_spatial(c3)mol/(m2*s)邊界通量邊界 6, 8, 11, 1315, 1718, 20, 2223, 2530, 3233, 35, 3738, 4043, 45, 48rfds.ntflux_c3rfds.bndFlux_c3

22、+rfds.cfluxx_c3*rfds.nxc+rfds.cfluxy_c3*rfds.nyc+rfds.cfluxz_c3*rfds.nzcmol/(m2*s)法向總通量邊界 150rfds.ndflux_c3rfds.bndFlux_c3mol/(m2*s)法向擴(kuò)散通量邊界 150rfds.ncflux_c3rfds.cfluxx_c3*rfds.nxc+rfds.cfluxy_c3*rfds.nyc+rfds.cfluxz_c3*rfds.nzcmol/(m2*s)法向?qū)α魍窟吔?150rfds.f_eff11等效傳遞因子域 15rfds.prefP_standerdPa參考?jí)毫λ?/p>

23、域 15rfds.pArfds.pref+pPa絕對(duì)壓力域 152.4.1 pdms傳遞屬性pdms傳遞屬性選擇幾何實(shí)體層次域選擇域 1方程SettingsDescriptionValue電勢(shì)用戶定義電勢(shì)0擴(kuò)散系數(shù)用戶定義擴(kuò)散系數(shù)D_pdms, 0, 0, 0, D_pdms, 0, 0, 0, D_pdms材料pdms (mat1)密度來(lái)自材料動(dòng)力粘度來(lái)自材料參考長(zhǎng)度1參考長(zhǎng)度比例自動(dòng)混合長(zhǎng)度限制自動(dòng)流體材料域材料擴(kuò)散系數(shù)用戶定義擴(kuò)散系數(shù)D_pdms*alpha_N2, 0, 0, 0, D_pdms*alpha_N2, 0, 0, 0, D_pdms*alpha_N2擴(kuò)散系數(shù)用戶定義擴(kuò)散系

24、數(shù)D_pdms*alpha_CO2, 0, 0, 0, D_pdms*alpha_CO2, 0, 0, 0, D_pdms*alpha_CO2來(lái)自材料的屬性屬性材料屬性組密度pdms基本動(dòng)力粘度pdms基本變量名稱表達(dá)式單位描述選擇rfds.Dxx_c1D_pdms*rfds.f_effm2/s擴(kuò)散系數(shù)，xx 分量域 1rfds.Dyx_c10m2/s擴(kuò)散系數(shù)，yx 分量域 1rfds.Dzx_c10m2/s擴(kuò)散系數(shù)，zx 分量域 1rfds.Dxy_c10m2/s擴(kuò)散系數(shù)，xy 分量域 1rfds.Dyy_c1D_pdms*rfds.f_effm2/s擴(kuò)散系數(shù)，yy 分量域 1rfds.Dz

25、y_c10m2/s擴(kuò)散系數(shù)，zy 分量域 1rfds.Dxz_c10m2/s擴(kuò)散系數(shù)，xz 分量域 1rfds.Dyz_c10m2/s擴(kuò)散系數(shù)，yz 分量域 1rfds.Dzz_c1D_pdms*rfds.f_effm2/s擴(kuò)散系數(shù)，zz 分量域 1rfds.Dav_c1(rfds.Dxx_c1+rfds.Dyy_c1+rfds.Dzz_c1)/3m2/s平均擴(kuò)散系數(shù)域 1rfds.tfluxx_c1-rfds.Dxx_c1*c1x-rfds.Dxy_c1*c1y-rfds.Dxz_c1*c1z+rfds.cfluxx_c1mol/(m2*s)總通量，x 分量域 1rfds.tfluxy_c1

26、-rfds.Dyx_c1*c1x-rfds.Dyy_c1*c1y-rfds.Dyz_c1*c1z+rfds.cfluxy_c1mol/(m2*s)總通量，y 分量域 1rfds.tfluxz_c1-rfds.Dzx_c1*c1x-rfds.Dzy_c1*c1y-rfds.Dzz_c1*c1z+rfds.cfluxz_c1mol/(m2*s)總通量，z 分量域 1rfds.dfluxx_c1-rfds.Dxx_c1*c1x-rfds.Dxy_c1*c1y-rfds.Dxz_c1*c1zmol/(m2*s)擴(kuò)散通量，x 分量域 1rfds.dfluxy_c1-rfds.Dyx_c1*c1x-rfd

27、s.Dyy_c1*c1y-rfds.Dyz_c1*c1zmol/(m2*s)擴(kuò)散通量，y 分量域 1rfds.dfluxz_c1-rfds.Dzx_c1*c1x-rfds.Dzy_c1*c1y-rfds.Dzz_c1*c1zmol/(m2*s)擴(kuò)散通量，z 分量域 1rfds.gradx_c1c1xmol/m4濃度梯度，x 分量域 1rfds.grady_c1c1ymol/m4濃度梯度，y 分量域 1rfds.gradz_c1c1zmol/m4濃度梯度，z 分量域 1rfds.dfluxMag_c1sqrt(rfds.dfluxx_c12+rfds.dfluxy_c12+rfds.dfluxz

28、_c12)mol/(m2*s)擴(kuò)散通量大小域 1rfds.tfluxMag_c1sqrt(rfds.tfluxx_c12+rfds.tfluxy_c12+rfds.tfluxz_c12)mol/(m2*s)總通量大小域 1rfds.Dxx_c2D_pdms*alpha_N2*rfds.f_effm2/s擴(kuò)散系數(shù)，xx 分量域 1rfds.Dyx_c20m2/s擴(kuò)散系數(shù)，yx 分量域 1rfds.Dzx_c20m2/s擴(kuò)散系數(shù)，zx 分量域 1rfds.Dxy_c20m2/s擴(kuò)散系數(shù)，xy 分量域 1rfds.Dyy_c2D_pdms*alpha_N2*rfds.f_effm2/s擴(kuò)散系數(shù)，yy

29、 分量域 1rfds.Dzy_c20m2/s擴(kuò)散系數(shù)，zy 分量域 1rfds.Dxz_c20m2/s擴(kuò)散系數(shù)，xz 分量域 1rfds.Dyz_c20m2/s擴(kuò)散系數(shù)，yz 分量域 1rfds.Dzz_c2D_pdms*alpha_N2*rfds.f_effm2/s擴(kuò)散系數(shù)，zz 分量域 1rfds.Dav_c2(rfds.Dxx_c2+rfds.Dyy_c2+rfds.Dzz_c2)/3m2/s平均擴(kuò)散系數(shù)域 1rfds.tfluxx_c2-rfds.Dxx_c2*c2x-rfds.Dxy_c2*c2y-rfds.Dxz_c2*c2z+rfds.cfluxx_c2mol/(m2*s)總通量

30、，x 分量域 1rfds.tfluxy_c2-rfds.Dyx_c2*c2x-rfds.Dyy_c2*c2y-rfds.Dyz_c2*c2z+rfds.cfluxy_c2mol/(m2*s)總通量，y 分量域 1rfds.tfluxz_c2-rfds.Dzx_c2*c2x-rfds.Dzy_c2*c2y-rfds.Dzz_c2*c2z+rfds.cfluxz_c2mol/(m2*s)總通量，z 分量域 1rfds.dfluxx_c2-rfds.Dxx_c2*c2x-rfds.Dxy_c2*c2y-rfds.Dxz_c2*c2zmol/(m2*s)擴(kuò)散通量，x 分量域 1rfds.dfluxy_

31、c2-rfds.Dyx_c2*c2x-rfds.Dyy_c2*c2y-rfds.Dyz_c2*c2zmol/(m2*s)擴(kuò)散通量，y 分量域 1rfds.dfluxz_c2-rfds.Dzx_c2*c2x-rfds.Dzy_c2*c2y-rfds.Dzz_c2*c2zmol/(m2*s)擴(kuò)散通量，z 分量域 1rfds.gradx_c2c2xmol/m4濃度梯度，x 分量域 1rfds.grady_c2c2ymol/m4濃度梯度，y 分量域 1rfds.gradz_c2c2zmol/m4濃度梯度，z 分量域 1rfds.dfluxMag_c2sqrt(rfds.dfluxx_c22+rfds.

32、dfluxy_c22+rfds.dfluxz_c22)mol/(m2*s)擴(kuò)散通量大小域 1rfds.tfluxMag_c2sqrt(rfds.tfluxx_c22+rfds.tfluxy_c22+rfds.tfluxz_c22)mol/(m2*s)總通量大小域 1rfds.Dxx_c3D_pdms*alpha_CO2*rfds.f_effm2/s擴(kuò)散系數(shù)，xx 分量域 1rfds.Dyx_c30m2/s擴(kuò)散系數(shù)，yx 分量域 1rfds.Dzx_c30m2/s擴(kuò)散系數(shù)，zx 分量域 1rfds.Dxy_c30m2/s擴(kuò)散系數(shù)，xy 分量域 1rfds.Dyy_c3D_pdms*alpha_C

33、O2*rfds.f_effm2/s擴(kuò)散系數(shù)，yy 分量域 1rfds.Dzy_c30m2/s擴(kuò)散系數(shù)，zy 分量域 1rfds.Dxz_c30m2/s擴(kuò)散系數(shù)，xz 分量域 1rfds.Dyz_c30m2/s擴(kuò)散系數(shù)，yz 分量域 1rfds.Dzz_c3D_pdms*alpha_CO2*rfds.f_effm2/s擴(kuò)散系數(shù)，zz 分量域 1rfds.Dav_c3(rfds.Dxx_c3+rfds.Dyy_c3+rfds.Dzz_c3)/3m2/s平均擴(kuò)散系數(shù)域 1rfds.tfluxx_c3-rfds.Dxx_c3*c3x-rfds.Dxy_c3*c3y-rfds.Dxz_c3*c3z+rf

34、ds.cfluxx_c3mol/(m2*s)總通量，x 分量域 1rfds.tfluxy_c3-rfds.Dyx_c3*c3x-rfds.Dyy_c3*c3y-rfds.Dyz_c3*c3z+rfds.cfluxy_c3mol/(m2*s)總通量，y 分量域 1rfds.tfluxz_c3-rfds.Dzx_c3*c3x-rfds.Dzy_c3*c3y-rfds.Dzz_c3*c3z+rfds.cfluxz_c3mol/(m2*s)總通量，z 分量域 1rfds.dfluxx_c3-rfds.Dxx_c3*c3x-rfds.Dxy_c3*c3y-rfds.Dxz_c3*c3zmol/(m2*s

35、)擴(kuò)散通量，x 分量域 1rfds.dfluxy_c3-rfds.Dyx_c3*c3x-rfds.Dyy_c3*c3y-rfds.Dyz_c3*c3zmol/(m2*s)擴(kuò)散通量，y 分量域 1rfds.dfluxz_c3-rfds.Dzx_c3*c3x-rfds.Dzy_c3*c3y-rfds.Dzz_c3*c3zmol/(m2*s)擴(kuò)散通量，z 分量域 1rfds.gradx_c3c3xmol/m4濃度梯度，x 分量域 1rfds.grady_c3c3ymol/m4濃度梯度，y 分量域 1rfds.gradz_c3c3zmol/m4濃度梯度，z 分量域 1rfds.dfluxMag_c3s

36、qrt(rfds.dfluxx_c32+rfds.dfluxy_c32+rfds.dfluxz_c32)mol/(m2*s)擴(kuò)散通量大小域 1rfds.tfluxMag_c3sqrt(rfds.tfluxx_c32+rfds.tfluxy_c32+rfds.tfluxz_c32)mol/(m2*s)總通量大小域 1rfds.umodel.input.u1m/s速度場(chǎng)，x 分量域 1rfds.vmodel.input.u2m/s速度場(chǎng)，y 分量域 1rfds.wmodel.input.u3m/s速度場(chǎng)，z 分量域 1rfds.cfluxx_c1c1*model.input.u1mol/(m2*s

37、)對(duì)流通量，x 分量域 1rfds.cfluxy_c1c1*model.input.u2mol/(m2*s)對(duì)流通量，y 分量域 1rfds.cfluxz_c1c1*model.input.u3mol/(m2*s)對(duì)流通量，z 分量域 1rfds.cfluxMag_c1sqrt(rfds.cfluxx_c12+rfds.cfluxy_c12+rfds.cfluxz_c12)mol/(m2*s)對(duì)流通量大小域 1rfds.cbf_c10mol/(m2*s)對(duì)流邊界通量邊界 16, 8, 1115, 1718, 20, 2230, 3233, 35, 3743, 45, 4850rfds.cflu

38、xx_c2c2*model.input.u1mol/(m2*s)對(duì)流通量，x 分量域 1rfds.cfluxy_c2c2*model.input.u2mol/(m2*s)對(duì)流通量，y 分量域 1rfds.cfluxz_c2c2*model.input.u3mol/(m2*s)對(duì)流通量，z 分量域 1rfds.cfluxMag_c2sqrt(rfds.cfluxx_c22+rfds.cfluxy_c22+rfds.cfluxz_c22)mol/(m2*s)對(duì)流通量大小域 1rfds.cbf_c20mol/(m2*s)對(duì)流邊界通量邊界 16, 8, 1115, 1718, 20, 2230, 32

39、33, 35, 3743, 45, 4850rfds.cfluxx_c3c3*model.input.u1mol/(m2*s)對(duì)流通量，x 分量域 1rfds.cfluxy_c3c3*model.input.u2mol/(m2*s)對(duì)流通量，y 分量域 1rfds.cfluxz_c3c3*model.input.u3mol/(m2*s)對(duì)流通量，z 分量域 1rfds.cfluxMag_c3sqrt(rfds.cfluxx_c32+rfds.cfluxy_c32+rfds.cfluxz_c32)mol/(m2*s)對(duì)流通量大小域 1rfds.cbf_c30mol/(m2*s)對(duì)流邊界通量邊界

40、16, 8, 1115, 1718, 20, 2230, 3233, 35, 3743, 45, 4850rfds.helemhm單元尺寸域 1rfds.Res_c1rfds.u*c1x+rfds.v*c1y+rfds.w*c1z-rfds.R_c1mol/(m3*s)方程殘差域 1rfds.Res_c2rfds.u*c2x+rfds.v*c2y+rfds.w*c2z-rfds.R_c2mol/(m3*s)方程殘差域 1rfds.Res_c3rfds.u*c3x+rfds.v*c3y+rfds.w*c3z-rfds.R_c3mol/(m3*s)方程殘差域 1rfds.rhomaterial.r

41、hokg/m3密度域 1rfds.mumaterial.muPa*s動(dòng)力粘度域 1rfds.divuux+vy+wz1/s速度場(chǎng)散度域 1rfds.srsqrt(0.5*(4*ux2+2*(uy+vx)2+2*(uz+wx)2+4*vy2+2*(vz+wy)2+4*wz2)+eps)1/s剪切速率域 1rfds.Fx0N/m3體積力，x 分量域 1rfds.Fy0N/m3體積力，y 分量域 1rfds.Fz0N/m3體積力，z 分量域 1rfds.Usqrt(u2+v2+w2)m/s速度大小域 1rfds.vorticityxwy-vz1/s渦流強(qiáng)度場(chǎng)，x 分量域 1rfds.vorticit

42、yy-wx+uz1/s渦流強(qiáng)度場(chǎng)，y 分量域 1rfds.vorticityzvx-uy1/s渦流強(qiáng)度場(chǎng)，z 分量域 1rfds.vort_magnsqrt(rfds.vorticityx2+rfds.vorticityy2+rfds.vorticityz2)1/s渦流大小域 1rfds.cellRe0.25*rfds.rho*sqrt(emetric(u,v,w)/emetric2)/rfds.mu1單元雷諾數(shù)域 1rfds.nurfds.mu/rfds.rhom2/s動(dòng)力學(xué)粘度域 1rfds.betaTd(rfds.rho,p)/rfds.rho1/Pa等溫壓縮系數(shù)域 1rfds.T_st

43、ressx2*rfds.mu*ux*rfds.nxmesh+rfds.mu*(uy+vx)*rfds.nymesh+rfds.mu*(uz+wx)*rfds.nzmesh-2*rfds.divu*rfds.mu*rfds.nxmesh/3-p*rfds.nxmeshN/m2總應(yīng)力，x 分量邊界 16, 8, 1115, 1718, 20, 2230, 3233, 35, 3743, 45, 4850rfds.T_stressyrfds.mu*(vx+uy)*rfds.nxmesh+2*rfds.mu*vy*rfds.nymesh+rfds.mu*(vz+wy)*rfds.nzmesh-2*rf

44、ds.divu*rfds.mu*rfds.nymesh/3-p*rfds.nymeshN/m2總應(yīng)力，y 分量邊界 16, 8, 1115, 1718, 20, 2230, 3233, 35, 3743, 45, 4850rfds.T_stresszrfds.mu*(wx+uz)*rfds.nxmesh+rfds.mu*(wy+vz)*rfds.nymesh+2*rfds.mu*wz*rfds.nzmesh-2*rfds.divu*rfds.mu*rfds.nzmesh/3-p*rfds.nzmeshN/m2總應(yīng)力，z 分量邊界 16, 8, 1115, 1718, 20, 2230, 323

45、3, 35, 3743, 45, 4850rfds.K_stressxrfds.mu*(2*ux*rfds.nxmesh+(uy+vx)*rfds.nymesh+(uz+wx)*rfds.nzmesh-2*rfds.divu*rfds.nxmesh/3)N/m2粘性應(yīng)力，x 分量邊界 16, 8, 1115, 1718, 20, 2230, 3233, 35, 3743, 45, 4850rfds.K_stressyrfds.mu*(vx+uy)*rfds.nxmesh+2*vy*rfds.nymesh+(vz+wy)*rfds.nzmesh-2*rfds.divu*rfds.nymesh/3

46、)N/m2粘性應(yīng)力，y 分量邊界 16, 8, 1115, 1718, 20, 2230, 3233, 35, 3743, 45, 4850rfds.K_stresszrfds.mu*(wx+uz)*rfds.nxmesh+(wy+vz)*rfds.nymesh+2*wz*rfds.nzmesh-2*rfds.divu*rfds.nzmesh/3)N/m2粘性應(yīng)力，z 分量邊界 16, 8, 1115, 1718, 20, 2230, 3233, 35, 3743, 45, 4850rfds.upwind_helpxum/s迎風(fēng)項(xiàng)，x 分量域 1rfds.upwind_helpyvm/s迎風(fēng)項(xiàng)

47、，y 分量域 1rfds.upwind_helpzwm/s迎風(fēng)項(xiàng)，z 分量域 1rfds.K_stress_tensorxxrfds.mu*(2*ux-2*rfds.divu/3)N/m2粘性應(yīng)力張量，xx 分量域 1rfds.K_stress_tensoryxrfds.mu*(vx+uy)N/m2粘性應(yīng)力張量，yx 分量域 1rfds.K_stress_tensorzxrfds.mu*(wx+uz)N/m2粘性應(yīng)力張量，zx 分量域 1rfds.K_stress_tensorxyrfds.mu*(uy+vx)N/m2粘性應(yīng)力張量，xy 分量域 1rfds.K_stress_tensoryyr

48、fds.mu*(2*vy-2*rfds.divu/3)N/m2粘性應(yīng)力張量，yy 分量域 1rfds.K_stress_tensorzyrfds.mu*(wy+vz)N/m2粘性應(yīng)力張量，zy 分量域 1rfds.K_stress_tensorxzrfds.mu*(uz+wx)N/m2粘性應(yīng)力張量，xz 分量域 1rfds.K_stress_tensoryzrfds.mu*(vz+wy)N/m2粘性應(yīng)力張量，yz 分量域 1rfds.K_stress_tensorzzrfds.mu*(2*wz-2*rfds.divu/3)N/m2粘性應(yīng)力張量，zz 分量域 1rfds.K_stress_ten

49、sor_testxxrfds.mu*(2*test(ux)-2*(test(ux)+test(vy)+test(wz)/3)N/m2粘性應(yīng)力張量試函數(shù)，xx 分量域 1rfds.K_stress_tensor_testyxrfds.mu*(test(vx)+test(uy)N/m2粘性應(yīng)力張量試函數(shù)，yx 分量域 1rfds.K_stress_tensor_testzxrfds.mu*(test(wx)+test(uz)N/m2粘性應(yīng)力張量試函數(shù)，zx 分量域 1rfds.K_stress_tensor_testxyrfds.mu*(test(uy)+test(vx)N/m2粘性應(yīng)力張量試函數(shù)，xy 分量域 1rfds.K_stress_tensor_testyyrfds.mu*(2*test(vy)-2*(test(ux)+test(vy)+test(wz)/3)N/m2粘性應(yīng)力張量試函數(shù)，yy 分量域 1rfds.K_stress_tensor_testzyrfds.mu*(test(wy)+test(vz)N/m2粘性應(yīng)力張量試函數(shù)，zy 分量域 1rfds.K_stress_tensor_testxzrf