MFI 教学分析和总结.docxVIP

MFIX 教学 1、首先来看看 mfix.dat 文件，它是 mfix 的工况设置文件，里面包含了物理问题的描述，几何描述，网格设置，物性参数，数学模型以及计算过程中的数据输出。它相当于fluent 中的cas 文件。 下面是文件的内容： # #和!表示注释行 # CFB model test simulation # M.Syamlal Date: 04/24/95 # Run time for F90 allocatable arrays version on an Octane is 1.8 days # and the disk space required is 85 Mb. Decrease TSTOP to shorten the run. # Increase SPX_DT values to reduce storage. # RUN CONTROL SECTION ! 下面是 case 执行控制段 ! 等号左边为关键字，右边是赋的值 RUN_NAME = CFB3D ! 计算结果输出文件的名字DESCRIPTION = CFB modeling test - 301 -- 3D UNITS = CGS ! 采用 CGS 单位制质量、长度、时间、能量分别用g, cm, s, cal RUN_TYPE = new ! 新运行一个工况，非重启一个工况 TIME = 0.0 ! 计算初始时间 TSTOP = 5.0 ! 计算终止时间 DT = 20.0e-4 ! 第一步时间步长；不写则为稳态计算；DT_FAC = 1 设定定时间步长计算 ENERGY_EQ = .F. ! 不求解能量方程 SPECIES_EQ = 2*.F. ! 不求解气相和固相的组分方程CALL_USR = .F. ! 不调用用户自定义程序 NORM_s = 0.0d0 ! 归一化连续性方程残差的参数；默认时用第一步计算的残差NORM_g = 0.0d0 # Physical parameters ! 物性段 e = 0.8 ! coefficient of restitution Phi = 0.0 ! angle of internal friction # GEOMETRY SECTION ! 网格 COORDINATES = CYLINDRICAL ! 选坐标系 XLENGTH = 20. ! x 方向的长度IMAX = 16 ! x 方向的网格数 YLENGTH = 955. JMAX = 191 ZLENGTH = @( 2*pi) ! @表示计算后面的式子KMAX = 12 # GAS SECTION ! 气相物性 MW_avg = 29.0 ! 分子量MU_g0 = 1.8e-4 ! 气相运动粘度 # PARTICLE SECTION ! 固相物性 MMAX = 1 ! 只有一个固相D_p0 = 0.012 ! 粒 径 RO_s = 2.6 ! 密 度EP_star = 0.4 ! 堆积空隙率 ur_fac(2) = 0.8 !solids pr. ! 欠松弛因子，2 为颗粒分数 # INITIAL CONDITIONS SECTION ! 初始条件 ! Bed IC_x_w = 0.0 ! x,y,z 的起止点 IC_x_e = 20.0 IC_y_s = 0.0 IC_y_n =955.0 IC_z_b = 0.0 IC_z_t =@( 2*pi) IC_ep_g = 0.9 ! 初始堆积空隙率IC_P_star = 0.0 ! 颗粒压力 IC_U_g = 0.0 ! 气相 x,y,z 速度 IC_V_g IC_W_g = 420.0 = 0.0 IC_U_s(1,1) = 0.0 ! 固相 x,y,z 速度 IC_V_s(1,1) = 0.0 ! 第一个下标为初始条件下标，第二个为固相下标IC_W_s(1,1) = 0.0 IC_T_g = 294.0 ! 气温 # BOUNDARY CONDITIONS SECTION ! 边界条件 ! ! air inlet BC_x_w(1) = 0.0 BC_x_e(1) = 20.0 BC_y_s(1) = 0.0 BC_y_n(1) = 0.0 BC_z_b(1) = 0.0 BC_z_t(1) = @( 2*pi) BC_TYPE(1) = MI BC_ep_g(1) = 0.9 BC_P_g(1) = 1151424.5 BC_T_g(1) = 294.0 BC_U_g(1) = 0.0 BC_V_g(1) = 420.0 BC_W_g(1) = 0.0 BC_U_s(1,1) = 0.0 BC_MASSFLOW_s(1,1)=