#! File written Wed Apr 6 15:07:21 2022 #! MFiX version 21.1.4 by a on 1234 ### Run Controls description = 'MFiX Simulation' run_name = 'yanzheng1' units = 'SI' run_type = 'new' tstop = 20.0 dt = 1.0000e-05 dt_min = 1.0000e-07 dt_max = 1.0000e-04 dt_fac = 0.9 res_dt = 0.1 chk_batchq_end = .False. drag_c1 = 0.8 drag_d1 = 2.65 drag_type = 'GIDASPOW' energy_eq = .False. nodesi = 1 nodesj = 4 nodesk = 1 full_log = .True. momentum_x_eq(0) = .True. momentum_y_eq(0) = .True. momentum_z_eq(0) = .True. project_version = '51' species_eq(0) = .False. species_eq(1) = .False. ### Physical parameters gravity_x = 0.0 gravity_y = -9.80665 gravity_z = 0.0 p_ref = 0.0 ### Cartesian grid cartesian_grid = .False. use_stl = .False. ### Numeric detect_stall = .True. max_inlet_vel_fac = 1.0 max_nit = 100 norm_g = 0.0 ### Discretization norm_s = 0.0 ur_fac(2) = 0.5 ### Geometry coordinates = 'CARTESIAN' imax = 18 jmax = 36 kmax = 4 x_max = 0.044 x_min = 0 y_max = 0.12 y_min = 0 z_max = 0.01 z_min = 0 #### Fluid mu_g0 = 1.8000e-05 mw_avg = 29.0 ro_g0 = 1.0 #### Solids mmax = 1 # Solid 1: particle solids_model(1) = 'TFM' d_p0(1) = 0.0012 ro_s0(1) = 1000.0 nmax_s(1) = 0 c_ps0(1) = 830.0 close_packed(1) = .True. k_s0(1) = 1.0 ks_model(1) = 'BAUER' ### Initial conditions # Initial condition 1: Background IC ic_x_e(1) = 0.044 #!MFIX-GUI eq{float(xmax)} ic_x_w(1) = 0.0 #!MFIX-GUI eq{float(xmin)} ic_y_s(1) = 0.0 #!MFIX-GUI eq{float(ymin)} ic_y_n(1) = 0.12 #!MFIX-GUI eq{float(ymax)} ic_z_b(1) = 0.0 #!MFIX-GUI eq{float(zmin)} ic_z_t(1) = 0.01 #!MFIX-GUI eq{float(zmax)} ic_des_fit_to_region(1) = .False. ic_ep_g(1) = 1.0 ic_t_g(1) = 293.15 ic_u_g(1) = 0.0 ic_v_g(1) = 0.0 ic_w_g(1) = 0.0 ic_p_star(1) = 0.0 ic_ep_s(1,1) = 0.0 ic_t_s(1,1) = 293.15 ic_theta_m(1,1) = 1.0000e-08 ic_u_s(1,1) = 0.0 ic_v_s(1,1) = 0.0 ic_w_s(1,1) = 0.0 # Initial condition 2: bed ic_x_e(2) = 0.044 #!MFIX-GUI eq{float(xmax)} ic_x_w(2) = 0.0 #!MFIX-GUI eq{float(xmin)} ic_y_s(2) = 0.0 #!MFIX-GUI eq{float(ymin)} ic_y_n(2) = 0.03 ic_z_b(2) = 0.0 #!MFIX-GUI eq{float(zmin)} ic_z_t(2) = 0.01 #!MFIX-GUI eq{float(zmax)} ic_des_fit_to_region(2) = .False. ic_ep_g(2) = 0.37 ic_t_g(2) = 293.15 ic_u_g(2) = 0.0 ic_v_g(2) = 0.0 ic_w_g(2) = 0.0 ic_p_star(2) = 0.0 ic_ep_s(2,1) = 0.63 ic_t_s(2,1) = 293.15 ic_theta_m(2,1) = 1.0000e-08 ic_u_s(2,1) = 0.0 ic_v_s(2,1) = 0.0 ic_w_s(2,1) = 0.0 # Boundary conditions # Boundary condition 1: inlet bc_type(1) = 'MI' bc_x_e(1) = 0.044 #!MFIX-GUI eq{float(xmax)} bc_x_w(1) = 0.0 #!MFIX-GUI eq{float(xmin)} bc_y_s(1) = 0.0 #!MFIX-GUI eq{float(ymin)} bc_y_n(1) = 0.0 #!MFIX-GUI eq{float(ymin)} bc_z_b(1) = 0.0 #!MFIX-GUI eq{float(zmin)} bc_z_t(1) = 0.01 #!MFIX-GUI eq{float(zmax)} bc_ep_g(1) = 1.0 bc_p_g(1) = 1.0132e+05 bc_t_g(1) = 293.15 bc_u_g(1) = 0.0 bc_v_g(1) = 0.9 bc_w_g(1) = 0.0 # Solid 1: particle bc_ep_s(1,1) = 0.0 bc_t_s(1,1) = 293.15 bc_u_s(1,1) = 0.0 bc_v_s(1,1) = 0.0 bc_w_s(1,1) = 0.0 bc_theta_m(1,1) = 1.0000e-08 # Boundary condition 2: outlet bc_type(2) = 'PO' bc_x_e(2) = 0.044 #!MFIX-GUI eq{float(xmax)} bc_x_w(2) = 0.0 #!MFIX-GUI eq{float(xmin)} bc_y_s(2) = 0.12 #!MFIX-GUI eq{float(ymax)} bc_y_n(2) = 0.12 #!MFIX-GUI eq{float(ymax)} bc_z_b(2) = 0.0 #!MFIX-GUI eq{float(zmin)} bc_z_t(2) = 0.01 #!MFIX-GUI eq{float(zmax)} bc_p_g(2) = 101325.0 # Boundary condition 3: wall1 bc_type(3) = 'PSW' bc_x_e(3) = 0.0 #!MFIX-GUI eq{float(xmin)} bc_x_w(3) = 0.0 #!MFIX-GUI eq{float(xmin)} bc_y_s(3) = 0.0 #!MFIX-GUI eq{float(ymin)} bc_y_n(3) = 0.12 #!MFIX-GUI eq{float(ymax)} bc_z_b(3) = 0.0 #!MFIX-GUI eq{float(zmin)} bc_z_t(3) = 0.01 #!MFIX-GUI eq{float(zmax)} bc_jj_ps(3) = 1 bc_uw_g(3) = 0.0 bc_vw_g(3) = 0.0 bc_ww_g(3) = 0.0 # Solid 1: particle bc_hw_s(3,1) = 0.0 bc_uw_s(3,1) = 0.0 bc_vw_s(3,1) = 0.0 bc_ww_s(3,1) = 0.0 # Boundary condition 4: wall2 bc_type(4) = 'PSW' bc_x_e(4) = 0.044 #!MFIX-GUI eq{float(xmax)} bc_x_w(4) = 0.044 #!MFIX-GUI eq{float(xmax)} bc_y_s(4) = 0.0 #!MFIX-GUI eq{float(ymin)} bc_y_n(4) = 0.12 #!MFIX-GUI eq{float(ymax)} bc_z_b(4) = 0.0 #!MFIX-GUI eq{float(zmin)} bc_z_t(4) = 0.01 #!MFIX-GUI eq{float(zmax)} bc_jj_ps(4) = 1 bc_uw_g(4) = 0.0 bc_vw_g(4) = 0.0 bc_ww_g(4) = 0.0 # Solid 1: particle bc_hw_s(4,1) = 0.0 bc_uw_s(4,1) = 0.0 bc_vw_s(4,1) = 0.0 bc_ww_s(4,1) = 0.0 # Boundary condition 5: wall3 bc_type(5) = 'PSW' bc_x_e(5) = 0.044 #!MFIX-GUI eq{float(xmax)} bc_x_w(5) = 0.0 #!MFIX-GUI eq{float(xmin)} bc_y_s(5) = 0.0 #!MFIX-GUI eq{float(ymin)} bc_y_n(5) = 0.12 #!MFIX-GUI eq{float(ymax)} bc_z_b(5) = 0.01 #!MFIX-GUI eq{float(zmax)} bc_z_t(5) = 0.01 #!MFIX-GUI eq{float(zmax)} bc_jj_ps(5) = 1 bc_uw_g(5) = 0.0 bc_vw_g(5) = 0.0 bc_ww_g(5) = 0.0 # Solid 1: particle bc_hw_s(5,1) = 0.0 bc_uw_s(5,1) = 0.0 bc_vw_s(5,1) = 0.0 bc_ww_s(5,1) = 0.0 # Boundary condition 6: wall4 bc_type(6) = 'PSW' bc_x_e(6) = 0.044 #!MFIX-GUI eq{float(xmax)} bc_x_w(6) = 0.0 #!MFIX-GUI eq{float(xmin)} bc_y_s(6) = 0.0 #!MFIX-GUI eq{float(ymin)} bc_y_n(6) = 0.12 #!MFIX-GUI eq{float(ymax)} bc_z_b(6) = 0.0 #!MFIX-GUI eq{float(zmin)} bc_z_t(6) = 0.0 #!MFIX-GUI eq{float(zmin)} bc_jj_ps(6) = 1 bc_uw_g(6) = 0.0 bc_vw_g(6) = 0.0 bc_ww_g(6) = 0.0 # Solid 1: particle bc_hw_s(6,1) = 0.0 bc_uw_s(6,1) = 0.0 bc_vw_s(6,1) = 0.0 bc_ww_s(6,1) = 0.0 # VTK outputs write_vtk_files = .True. time_dependent_filename = .True. vtu_dir = 'vtk' # VTK output 1: Background IC vtk_filebase(1) = 'Background_IC' vtk_x_e(1) = 0.044 #!MFIX-GUI eq{float(xmax)} vtk_x_w(1) = 0.0 #!MFIX-GUI eq{float(xmin)} vtk_y_s(1) = 0.0 #!MFIX-GUI eq{float(ymin)} vtk_y_n(1) = 0.12 #!MFIX-GUI eq{float(ymax)} vtk_z_b(1) = 0.0 #!MFIX-GUI eq{float(zmin)} vtk_z_t(1) = 0.01 #!MFIX-GUI eq{float(zmax)} vtk_data(1) = 'C' vtk_dt(1) = 0.01 vtk_nxs(1) = 0 vtk_nys(1) = 0 vtk_nzs(1) = 0 vtk_ep_g(1) = .True. vtk_p_g(1) = .True. vtk_p_star(1) = .True. vtk_u_g(1) = .True. vtk_v_g(1) = .True. vtk_vel_g(1) = .True. vtk_w_g(1) = .True. vtk_ep_s(1,1) = .True. vtk_theta_m(1,1) = .True. vtk_u_s(1,1) = .True. vtk_v_s(1,1) = .True. vtk_vel_s(1,1) = .True. vtk_w_s(1,1) = .True. # SPx outputs spx_dt(1) = 0.1 spx_dt(2) = 0.1 spx_dt(3) = 0.1 spx_dt(4) = 0.1 spx_dt(5) = 0.1 spx_dt(6) = 0.1 spx_dt(7) = 0.1 spx_dt(8) = 0.1 spx_dt(9) = 0.1 spx_dt(10) = 1.0 spx_dt(11) = 1.0 ### Residuals resid_string(1) = 'P0' resid_string(2) = 'U0' resid_string(3) = 'V0' resid_string(4) = 'P1' resid_string(5) = 'U1' resid_string(6) = 'V1' ### Discrete element model des_en_input(1) = 0.9 des_en_wall_input(1) = 0.9 kn = 1000 kn_w = 1000 mew = 0.1 mew_w = 0.1 ### Two-fluid model blending_function = 'NONE' c_e = 0.97 c_f = 0.1 e_w = 1.0 ep_star = 0.36 friction_model = 'SCHAEFFER' jenkins = .False. kt_type = 'GTSH' phi = 30.0 phi_w = 11.3 phip = 0.6 rdf_type = 'CARNAHAN_STARLING' #***************************************************** # Miscellaneous keys * #***************************************************** bc_jj_m = .False. # MFIX-GUI section #!MFIX-GUI regions_dict = {"order": ["Background IC", "bed", "inlet", "outlet", #!MFIX-GUI "wall1", "wall2", "wall3", "wall4"], "regions": {"Background IC": #!MFIX-GUI {"to": [{"__type__": "equation", "eq": "float(xmax)"}, {"__type__": #!MFIX-GUI "equation", "eq": "float(ymax)"}, {"__type__": "equation", "eq": #!MFIX-GUI "float(zmax)"}], "from": [{"__type__": "equation", "eq": #!MFIX-GUI "float(xmin)"}, {"__type__": "equation", "eq": "float(ymin)"}, #!MFIX-GUI {"__type__": "equation", "eq": "float(zmin)"}], "color": #!MFIX-GUI "#ffffff"}, "bed": {"to": [{"__type__": "equation", "eq": #!MFIX-GUI "float(xmax)"}, 0.03, {"__type__": "equation", "eq": #!MFIX-GUI "float(zmax)"}], "from": [{"__type__": "equation", "eq": #!MFIX-GUI "float(xmin)"}, {"__type__": "equation", "eq": "float(ymin)"}, #!MFIX-GUI {"__type__": "equation", "eq": "float(zmin)"}], "color": #!MFIX-GUI "#dbb370"}, "inlet": {"to": [{"__type__": "equation", "eq": #!MFIX-GUI "float(xmax)"}, {"__type__": "equation", "eq": "float(ymin)"}, #!MFIX-GUI {"__type__": "equation", "eq": "float(zmax)"}], "from": #!MFIX-GUI [{"__type__": "equation", "eq": "float(xmin)"}, {"__type__": #!MFIX-GUI "equation", "eq": "float(ymin)"}, {"__type__": "equation", "eq": #!MFIX-GUI "float(zmin)"}], "color": "#aa6ab6", "type": "XZ-plane"}, "outlet": #!MFIX-GUI {"to": [{"__type__": "equation", "eq": "float(xmax)"}, {"__type__": #!MFIX-GUI "equation", "eq": "float(ymax)"}, {"__type__": "equation", "eq": #!MFIX-GUI "float(zmax)"}], "from": [{"__type__": "equation", "eq": #!MFIX-GUI "float(xmin)"}, {"__type__": "equation", "eq": "float(ymax)"}, #!MFIX-GUI {"__type__": "equation", "eq": "float(zmin)"}], "color": "#c4b5d4", #!MFIX-GUI "type": "XZ-plane"}, "wall1": {"to": [{"__type__": "equation", #!MFIX-GUI "eq": "float(xmin)"}, {"__type__": "equation", "eq": #!MFIX-GUI "float(ymax)"}, {"__type__": "equation", "eq": "float(zmax)"}], #!MFIX-GUI "from": [{"__type__": "equation", "eq": "float(xmin)"}, #!MFIX-GUI {"__type__": "equation", "eq": "float(ymin)"}, {"__type__": #!MFIX-GUI "equation", "eq": "float(zmin)"}], "color": "#98cde1", "type": #!MFIX-GUI "YZ-plane"}, "wall2": {"to": [{"__type__": "equation", "eq": #!MFIX-GUI "float(xmax)"}, {"__type__": "equation", "eq": "float(ymax)"}, #!MFIX-GUI {"__type__": "equation", "eq": "float(zmax)"}], "from": #!MFIX-GUI [{"__type__": "equation", "eq": "float(xmax)"}, {"__type__": #!MFIX-GUI "equation", "eq": "float(ymin)"}, {"__type__": "equation", "eq": #!MFIX-GUI "float(zmin)"}], "color": "#c26ad6", "type": "YZ-plane"}, "wall3": #!MFIX-GUI {"to": [{"__type__": "equation", "eq": "float(xmax)"}, {"__type__": #!MFIX-GUI "equation", "eq": "float(ymax)"}, {"__type__": "equation", "eq": #!MFIX-GUI "float(zmax)"}], "from": [{"__type__": "equation", "eq": #!MFIX-GUI "float(xmin)"}, {"__type__": "equation", "eq": "float(ymin)"}, #!MFIX-GUI {"__type__": "equation", "eq": "float(zmax)"}], "color": "#d8c882", #!MFIX-GUI "type": "XY-plane"}, "wall4": {"to": [{"__type__": "equation", #!MFIX-GUI "eq": "float(xmax)"}, {"__type__": "equation", "eq": #!MFIX-GUI "float(ymax)"}, {"__type__": "equation", "eq": "float(zmin)"}], #!MFIX-GUI "from": [{"__type__": "equation", "eq": "float(xmin)"}, #!MFIX-GUI {"__type__": "equation", "eq": "float(ymin)"}, {"__type__": #!MFIX-GUI "equation", "eq": "float(zmin)"}], "color": "#e3ceb3", "type": #!MFIX-GUI "XY-plane"}}} #!MFIX-GUI author = a #!MFIX-GUI created_date = 2022-04-04 09:55 #!MFIX-GUI modified_time = 2022-04-05 09:47 #!MFIX-GUI project_version = 51 #!MFIX-GUI gui_version = 21.1.4 #!MFIX-GUI project_notes = "" #!MFIX-GUI modified_by = a #!MFIX-GUI geometry = {"geometry_dict": {}, "tree": {}} #!MFIX-GUI visual_props = {"mesh": {"color": "#f44328", "visible": true, #!MFIX-GUI "opacity": 1.0, "rep": "solid", "color_by": "Volume", "edge": #!MFIX-GUI "#7a2114"}, "background_mesh": {"color": "#64b6f7", "visible": #!MFIX-GUI true, "opacity": 1.0, "rep": "wire", "edge": "#325b7b"}, #!MFIX-GUI "geometry": {"color": "#e0e0e0", "visible": true, "opacity": 1.0, #!MFIX-GUI "rep": "wire", "edge": "#707070"}, "regions": {"color": "#e0e0e0", #!MFIX-GUI "visible": true, "opacity": 0.5, "rep": "solid", "edge": #!MFIX-GUI "#707070"}, "normals": {"color": "#0000e0", "visible": false, #!MFIX-GUI "scale": 0.1, "count": 1000, "edge": "#000070"}, "axes": {"color": #!MFIX-GUI "#0000e0", "visible": true, "edge": "#000070"}, "boundary": #!MFIX-GUI {"color": "#0000e0", "opacity": 1.0, "rep": "edges", "visible": #!MFIX-GUI false, "color_by": "bc_id", "edge": "#000070"}} #!MFIX-GUI ic_regions = [[[1], ["Background IC"]], [[2], ["bed"]]] #!MFIX-GUI graphics = {"order": ["Residuals", "DT"], "data": {"Residuals": #!MFIX-GUI {"plot": true, "vtk": false, "monitors": false, "histogram": #!MFIX-GUI false}, "DT": {"plot": true, "vtk": false, "monitors": false, #!MFIX-GUI "histogram": false}}} #!MFIX-GUI ui_state = {"mode": "modeler", "navigation": "Solids", "graphic_tab": 2} #!MFIX-GUI parameters = {"order": ["min", "max", "xmin", "xmax", "ymin", "ymax", #!MFIX-GUI "zmin", "zmax", "km", "m", "cm", "mm", "um", "mile", "yard", "ft", #!MFIX-GUI "ins"], "parameters": {"min": 0.0, "max": 0.0, "xmin": 0, "xmax": #!MFIX-GUI 0.044, "ymin": 0, "ymax": 0.12, "zmin": 0, "zmax": 0.01, "km": #!MFIX-GUI 1000.0, "m": 1.0, "cm": 0.01, "mm": 0.001, "um": {"__type__": #!MFIX-GUI "equation", "eq": "float(1e-06)"}, "mile": 1609.34, "yard": 0.9144, #!MFIX-GUI "ft": 0.3048, "ins": 0.0254}} #!MFIX-GUI solids_phase_name(1) = particle #!MFIX-GUI bc_regions = [[[1], ["inlet"]], [[2], ["outlet"]], [[3, 4, 5, 6], #!MFIX-GUI ["wall1", "wall2", "wall3", "wall4"]]] #!MFIX-GUI vtk_regions = [[[1], ["Background IC"]]] #!MFIX-GUI BUILD_DMP = 1 #!MFIX-GUI BUILD_INTERACTIVE = PYMFIX #!MFIX-GUI BUILD_PARALLEL = 1 #!MFIX-GUI BUILD_SMP = 0 #!MFIX-GUI BUILD_TYPE = RELWITHDEBINFO #!MFIX-GUI mfix_exe = /home/a/z/mfix/stress/GTSH/yanzheng1/mfixsolver #!MFIX-GUI queue_template = {"JOB_NAME": "${PROJECT_NAME}", "CORES": 40, #!MFIX-GUI "QUEUE": "general", "LONG": "", "MODULES": [], "template": "Joule"} #!MFIX-GUI OMP_NUM_THREADS = 1 #!MFIX-GUI submit_to_queue = 0 #!MFIX-GUI retained_keys = {"ro_g0": 1.0}