comsol流体仿真(层流-传质耦合场)(对流扩散)

1、流体的对流扩散仿真报告日期2016-4-5 11:27:37目录1.全局定义1.1.参数 12.组件 12.1.定义2.2.几何 12.3.材料2.4.多孔介质反应流2.5.网格 13.研究 13.1.参数化扫描3.2.稳态3.3.求解器配置4.结果4.1.数据集4.2.绘图组1 全局定义日期Mar 18, 2016 4:27:52 PM全局设定名称H NewModel.mph路径E:COMSOLh_NewModel.mphCOMSOL 版本COMSOL 5.2 (建立: 220)单位系统SI使用的模块COMSOL MultiphysicsCAD Import ModuleChemical R

2、eaction Engineering Module1.1 参数 1参数名称表达式值描述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空气的平均摩尔质量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组件设定单位系统SI几何形参阶次automatic2.1 定义2.1.1 选择显式 1选择类型Explicit选择没有边显式 12.1.2 坐标系边界坐标系 1坐标系类型边界坐标系标记sys1坐标名称第一轴第二轴第三轴t1t2n2.2 几何 1几何 1单位长度单位µm角度单位deg几何统计描述值空间维度3域数5边界数50边数116端点数722.2.1 工作平面 1 (wp1)联合对象描述值联合对象开平面几何 (wp1)正方形 1 (sq1)位置描述值位置0, 0基准中心尺寸描述值边长22000.02.2.2 拉伸 1 (

6、ext1)设定描述值工作平面工作平面 1到面的距离距离 (µm)-h_in比例比例 xw比例 yw11位移位移 xw (µm)位移 yw (µm)00扭转角扭转角 (deg)02.2.3 工作平面 2 (wp2)联合对象描述值联合对象开平面几何 (wp2)多边形 1 (pol1)对象类型描述值类型实体坐标描述值数据源表格坐标xw (µm)yw (µm)45011000.045035001503500150-3500350.0-3500350.0-11000.0750-110000-3500950.0350065035

7、0065011000.045011000.0多边形 2 (pol2)对象类型描述值类型实体坐标描述值数据源表格坐标xw (µm)yw (µm)115011000.0175011000.01750-11000.01150-11000.0115011000.02.2.4 拉伸 2 (ext2)设定描述值工作平面工作平面 2到面的距离距离 (µm)h_channal到面的距离描述值反向开比例比例 xw比例 yw11位移位移 xw (µm)位移 yw (µm)00扭转角扭转角 (deg)02.2.5 镜像 1 (mir1)设定描述值保留输入对象开反射面

8、上的点描述值面内点0, 0, 0反射面法向矢量描述值平面法线1, 0, 02.2.6 工作平面 3 (wp3)联合对象描述值联合对象开平面几何 (wp3)贝氏多线段 1 (b1)多边形线段描述值控制点-1750, -1149.9999761580998, -1149.9999761580998, -1750, -1750, 11000.0, 11000.0, -11000.0, -11000.0, 11000.0度1, 1, 1, 1权重1, 1, 1, 1, 1, 1, 1, 1类型实体贝氏多线段 2 (b2)多边形线段描述值控制点-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权重1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1类型实体贝氏多线段 3 (b3)多边形线段描述值控制点1149.9999761580998, 1750, 1750, 1149.9999761580998, 1149.9999761580998, 11000.0, 11000.0, -11000.0, -11000.0, 11000.0度1, 1, 1, 1权重1, 1, 1, 1, 1, 1, 1, 1类型实体贝氏多线段 4 (b4)多边形线段描述值控制点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权重1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1类型实体2.2.7 channal (ext3)设定描述值工作平面工作平面 3到面的距离距离 (µm)-h_channal比例比例 xw比例 yw11位移位移 xw (µm)位移 yw (µm)00扭转角扭转角 (deg)02.3 材料2.3.1 pdmspdms选择几何实体层次域选择域 1材料参数名称值单位密度p_mediakg/m3动

13、力粘度u_pdmsPa*s基本 设定DescriptionValue密度p_media动力粘度u_pdms2.3.2 gasgas选择几何实体层次域选择域 2, 5材料参数名称值单位密度p_gaskg/m3动力粘度u_gasPa*s基本 设定DescriptionValue密度p_gas动力粘度u_gas2.3.3 mediamedia选择几何实体层次域选择域 34材料参数名称值单位密度p_mediakg/m3动力粘度u_mediaPa*s基本 设定DescriptionValue密度p_media动力粘度u_media2.4 多孔介质反应流使用的模块COMSOL MultiphysicsCh

14、emical Reaction Engineering Module多孔介质反应流选择几何实体层次域选择域 15EquationsSettingsDescriptionValue浓度线性计算边界通量开应用平滑到边界通量开流体离散化P1 + P1使用分裂复数变量时的值类型实数, 实数, 实数变量名称表达式单位描述选择rfds.nxnx法矢，x 分量边界 150rfds.nyny法矢，y 分量边界 150rfds.nznz法矢，z 分量边界 150rfds.nxmeshroot.nxmesh法矢（网格），x 分量边界 150rfds.nymeshroot.nymesh法矢（网格），y 分量边界 1

15、50rfds.nzmeshroot.nzmesh法矢（网格），z 分量边界 150rfds.R_c10mol/(m3*s)总速率表达式域 15rfds.R_c20mol/(m3*s)总速率表达式域 15rfds.R_c30mol/(m3*s)总速率表达式域 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)法向扩散通量边界 150rfds.ncflux_c1rf

19、ds.cfluxx_c1*rfds.nxc+rfds.cfluxy_c1*rfds.nyc+rfds.cfluxz_c1*rfds.nzcmol/(m2*s)法向对流通量边界 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)法向扩散通量边界 150rfds.ncflux_c2rfds.cfluxx_c2*rfds.nxc+rfds.cfluxy_c2*rfds.nyc+rfds.cfluxz_c2*rfds.nzcm

21、ol/(m2*s)法向对流通量边界 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)法向扩散通量边界 150rfds.ncflux_c3rfds.cfluxx_c3*rfds.nxc+rfds.cfluxy_c3*rfds.nyc+rfds.cfluxz_c3*rfds.nzcmol/(m2*s)法向对流通量边界 150rfds.f_eff11等效传递因子域 15rfds.prefP_standerdPa参考压力水平

23、域 15rfds.pArfds.pref+pPa绝对压力域 152.4.1 pdms传递属性pdms传递属性选择几何实体层次域选择域 1方程SettingsDescriptionValue电势用户定义电势0扩散系数用户定义扩散系数D_pdms, 0, 0, 0, D_pdms, 0, 0, 0, D_pdms材料pdms (mat1)密度来自材料动力粘度来自材料参考长度1参考长度比例自动混合长度限制自动流体材料域材料扩散系数用户定义扩散系数D_pdms*alpha_N2, 0, 0, 0, D_pdms*alpha_N2, 0, 0, 0, D_pdms*alpha_N2扩散系数用户定义扩散系

24、数D_pdms*alpha_CO2, 0, 0, 0, D_pdms*alpha_CO2, 0, 0, 0, D_pdms*alpha_CO2来自材料的属性属性材料属性组密度pdms基本动力粘度pdms基本变量名称表达式单位描述选择rfds.Dxx_c1D_pdms*rfds.f_effm2/s扩散系数，xx 分量域 1rfds.Dyx_c10m2/s扩散系数，yx 分量域 1rfds.Dzx_c10m2/s扩散系数，zx 分量域 1rfds.Dxy_c10m2/s扩散系数，xy 分量域 1rfds.Dyy_c1D_pdms*rfds.f_effm2/s扩散系数，yy 分量域 1rfds.Dz

25、y_c10m2/s扩散系数，zy 分量域 1rfds.Dxz_c10m2/s扩散系数，xz 分量域 1rfds.Dyz_c10m2/s扩散系数，yz 分量域 1rfds.Dzz_c1D_pdms*rfds.f_effm2/s扩散系数，zz 分量域 1rfds.Dav_c1(rfds.Dxx_c1+rfds.Dyy_c1+rfds.Dzz_c1)/3m2/s平均扩散系数域 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)扩散通量，x 分量域 1rfds.dfluxy_c1-rfds.Dyx_c1*c1x-rfd

27、s.Dyy_c1*c1y-rfds.Dyz_c1*c1zmol/(m2*s)扩散通量，y 分量域 1rfds.dfluxz_c1-rfds.Dzx_c1*c1x-rfds.Dzy_c1*c1y-rfds.Dzz_c1*c1zmol/(m2*s)扩散通量，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)扩散通量大小域 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扩散系数，xx 分量域 1rfds.Dyx_c20m2/s扩散系数，yx 分量域 1rfds.Dzx_c20m2/s扩散系数，zx 分量域 1rfds.Dxy_c20m2/s扩散系数，xy 分量域 1rfds.Dyy_c2D_pdms*alpha_N2*rfds.f_effm2/s扩散系数，yy

29、 分量域 1rfds.Dzy_c20m2/s扩散系数，zy 分量域 1rfds.Dxz_c20m2/s扩散系数，xz 分量域 1rfds.Dyz_c20m2/s扩散系数，yz 分量域 1rfds.Dzz_c2D_pdms*alpha_N2*rfds.f_effm2/s扩散系数，zz 分量域 1rfds.Dav_c2(rfds.Dxx_c2+rfds.Dyy_c2+rfds.Dzz_c2)/3m2/s平均扩散系数域 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)扩散通量，x 分量域 1rfds.dfluxy_

31、c2-rfds.Dyx_c2*c2x-rfds.Dyy_c2*c2y-rfds.Dyz_c2*c2zmol/(m2*s)扩散通量，y 分量域 1rfds.dfluxz_c2-rfds.Dzx_c2*c2x-rfds.Dzy_c2*c2y-rfds.Dzz_c2*c2zmol/(m2*s)扩散通量，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)扩散通量大小域 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扩散系数，xx 分量域 1rfds.Dyx_c30m2/s扩散系数，yx 分量域 1rfds.Dzx_c30m2/s扩散系数，zx 分量域 1rfds.Dxy_c30m2/s扩散系数，xy 分量域 1rfds.Dyy_c3D_pdms*alpha_C

33、O2*rfds.f_effm2/s扩散系数，yy 分量域 1rfds.Dzy_c30m2/s扩散系数，zy 分量域 1rfds.Dxz_c30m2/s扩散系数，xz 分量域 1rfds.Dyz_c30m2/s扩散系数，yz 分量域 1rfds.Dzz_c3D_pdms*alpha_CO2*rfds.f_effm2/s扩散系数，zz 分量域 1rfds.Dav_c3(rfds.Dxx_c3+rfds.Dyy_c3+rfds.Dzz_c3)/3m2/s平均扩散系数域 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、)扩散通量，x 分量域 1rfds.dfluxy_c3-rfds.Dyx_c3*c3x-rfds.Dyy_c3*c3y-rfds.Dyz_c3*c3zmol/(m2*s)扩散通量，y 分量域 1rfds.dfluxz_c3-rfds.Dzx_c3*c3x-rfds.Dzy_c3*c3y-rfds.Dzz_c3*c3zmol/(m2*s)扩散通量，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)扩散通量大小域 1rfds.tfluxMag_c3sqrt(rfds.tfluxx_c32+rfds.tfluxy_c32+rfds.tfluxz_c32)mol/(m2*s)总通量大小域 1rfds.umodel.input.u1m/s速度场，x 分量域 1rfds.vmodel.input.u2m/s速度场，y 分量域 1rfds.wmodel.input.u3m/s速度场，z 分量域 1rfds.cfluxx_c1c1*model.input.u1mol/(m2*s

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

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

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

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动力粘度域 1rfds.divuux+vy+wz1/s速度场散度域 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涡流强度场，x 分量域 1rfds.vorticit

42、yy-wx+uz1/s涡流强度场，y 分量域 1rfds.vorticityzvx-uy1/s涡流强度场，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单元雷诺数域 1rfds.nurfds.mu/rfds.rhom2/s动力学粘度域 1rfds.betaTd(rfds.rho,p)/rfds.rho1/Pa等温压缩系数域 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总应力，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 分量边界 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总应力，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粘性应力，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 分量边界 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粘性应力，z 分量边界 16, 8, 1115, 1718, 20, 2230, 3233, 35, 3743, 45, 4850rfds.upwind_helpxum/s迎风项，x 分量域 1rfds.upwind_helpyvm/s迎风项

47、，y 分量域 1rfds.upwind_helpzwm/s迎风项，z 分量域 1rfds.K_stress_tensorxxrfds.mu*(2*ux-2*rfds.divu/3)N/m2粘性应力张量，xx 分量域 1rfds.K_stress_tensoryxrfds.mu*(vx+uy)N/m2粘性应力张量，yx 分量域 1rfds.K_stress_tensorzxrfds.mu*(wx+uz)N/m2粘性应力张量，zx 分量域 1rfds.K_stress_tensorxyrfds.mu*(uy+vx)N/m2粘性应力张量，xy 分量域 1rfds.K_stress_tensoryyr

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

49、sor_testxxrfds.mu*(2*test(ux)-2*(test(ux)+test(vy)+test(wz)/3)N/m2粘性应力张量试函数，xx 分量域 1rfds.K_stress_tensor_testyxrfds.mu*(test(vx)+test(uy)N/m2粘性应力张量试函数，yx 分量域 1rfds.K_stress_tensor_testzxrfds.mu*(test(wx)+test(uz)N/m2粘性应力张量试函数，zx 分量域 1rfds.K_stress_tensor_testxyrfds.mu*(test(uy)+test(vx)N/m2粘性应力张量试函数，xy 分量域 1rfds.K_stress_tensor_testyyrfds.mu*(2*test(vy)-2*(test(ux)+test(vy)+test(wz)/3)N/m2粘性应力张量试函数，yy 分量域 1rfds.K_stress_tensor_testzyrfds.mu*(test(wy)+test(vz)N/m2粘性应力张量试函数，zy 分量域 1rfds.K_stress_tensor_testxzrf