! ! Bubbling Fluidized Bed Simulation ! ! Bubbling Fluidized Bed Simulation -- S. Benyahia SEP-19-2006 ! using a porous plate distributor with dP_plate about 30% dP_bed ! ! Run-control section RUN_NAME = 'R' DESCRIPTION = 'Fluid Bed with uniform gas distribution' RUN_TYPE = 'NEW' UNITS = 'cgs' TIME = 0. !start time TSTOP = 5.0 DT = 5.0E-4 !time step DT_MIN = 1.0E-8 !min time-step !! NORM_G = 0.0d0 NORM_S = 0.0d0 DISCRETIZE = 9*2 ! using superbee 2nd order scheme for all equations DEF_COR = .TRUE. ! with deffered correction ! MAX_NIT = 50 !! ENERGY_EQ = .FALSE. !do not solve energy eq SPECIES_EQ = .FALSE. .FALSE. !do not solve species eq ! GRANULAR_ENERGY = .TRUE. ! granular energy pde will be solved K_EPSILON = .F. FRICTION = .T. ! using frictional model of Srivastava-Sundaresan: SAVAGE = 2 ! Powder Tech (2003) vol. 129, pp. 72-85 ! ! Physical and Numerical Parameters ! ! ! Geometry Section # COORDINATES = 'Cartesian' XLENGTH = 40 IMAX = 124 YLENGTH = 78.44 JMAX = 124 NO_K = .TRUE. !2D, no k direction ZLENGTH = 1.0 ! Gas-phase Section MU_g0 = 1.8E-4 ! constant gas viscosity MW_avg = 29. MU_gmax = 1000. ! must define mu_gmax with k_epsilon P_REF = 1.01E6 ! using atmospheric reference pressure ! Solids-phase Section DRAG_TYPE = 'GIDASPOW' RO_s = 2.61 !solids density D_p0 = 0.05 !particle diameter e = 0.8 !restitution coefficient e_w = 0.7 !particle-wall restitution coefficient Phi = 30.0 !angle of internal friction EP_star = 0.4 !void fraction at packing PHIP = 0.001 !Specularite coef. for use with JJ BC phi_w = 11.31 ! angle of internal friction to use with Jenkins BC # # Initial Conditions Section # ! Bed Freeboard Plenum + IS IC_X_w = 0.0 0.0 0.0 IC_X_e = 40 40 40 IC_Y_s = 21.0 50.22 0.0 IC_Y_n = 50.22 78.44 21.0 IC_EP_g = 0.5 1.0 1.0 IC_U_g = 0.0 0.0 0.0 IC_V_g = @(60/0.5) 60 60.0 IC_U_s(1,1) = 0.0 0.0 0.0 IC_V_s(1,1) = 0.0 0.0 0.0 IC_P_star = 0.0 0.0 0.0 IC_T_g = 297. 297. 297. IC_Theta_m(1,1) = 10. IC_Theta_m(2,1) = 10. # # Boundary Conditions Section # ! Grid Exit BC_X_w = 16.0 0.0 BC_X_e = 24.0 40 BC_Y_s = 0.0 78.44 BC_Y_n = 0.0 78.44 BC_TYPE = 'MI' 'PO' BC_EP_g = 1.0 BC_U_g = 0.0 # BC_V_g = 300.0 !60cm/s superficial vel. BC_MASSFLOW_g = @(1.2D-03*300*8) BC_P_g = 2*1.01E6 BC_T_g = 2*297. # ############################################### # Wall BC's ############################################### BC_X_w(6) = 0. BC_X_e(6) = 0. BC_Y_s(6) = 0. BC_Y_n(6) = 78.44 BC_TYPE(6) = 'PSW' ! BC_Uw_g(6) = 0. BC_Vw_g(6) = 0. BC_Uw_s(6,1) = 0. BC_Vw_s(6,1) = 0. BC_Thetaw_m(6,1) = 0. BC_JJ_PS(6) = 1 ! use Johnson-Jackson BC for walls # ############################################### # BC_X_w(5) = 40.0 BC_X_e(5) = 40.0 BC_Y_s(5) = 0. BC_Y_n(5) = 78.44 BC_TYPE(5) = 'PSW' ! BC_Uw_g(5) = 0. BC_Vw_g(5) = 0. BC_Uw_s(5,1) = 0. BC_Vw_s(5,1) = 0. BC_Thetaw_m(5,1) = 0. BC_JJ_PS(5) = 1 ! ############################################### # semi-permeable internal surface of 1 cm thickness ############################################### IS_X_w(1) = 0. IS_X_e(1) = 40.0 IS_Y_s(1) = 20. IS_Y_n(1) = 21.0 IS_TYPE(1) = 'Y_SP' ! pressure drop along Y-flow-direction ! IS_PC(1,1) = 6.113151D-07 ! this is constant for Darcy's law (default 1D+32) IS_PC(1,2) = 0.0 ! inertia coefficient set to default zero value IS_VEL_S(1,1) = 0.0 # # ! Output Control ! ! TOL_DIVERGE = 1.E+08 OUT_DT = 10. !write text file R.OUT ! every 10 s RES_DT = 0.01 !write binary restart file ! R.RES every 0.01 s NLOG = 100 !write logfile R.LOG !every 100 time steps FULL_LOG = .TRUE. !display residuals on screen !SPX_DT values determine how often SPx files are written. Here BUB01.SP1, which !contains void fraction (EP_g), is written every 0.05 s, BUB01.SP2, which contains ! ! EP_g P_g U_g U_s ROP_s T_g X_g theta scalar rrate k_epsilon ! P_star V_g V_s T_s1 X_s ! W_g W_s T_s2 SPX_DT = 0.01 0.01 0.01 0.01 100. 100. 100. 0.01 100. 100. 0.01 RESID_STRING = 'P0' 'P1' 'U0' 'V0' 'U1' 'V1' 'G1' # end of mfix.dat file.