#! File written Wed Sep 16 14:57:06 2026
#! MFiX version 24.3.1 by gunjan on Gunjan

### Run controls
  description      = 'fbsqp_ar4'
  run_name         = 'FB_SQP'
  units            = 'SI'
  run_type         = 'new'
  tstop            = 20.0
  dt               = 1.0000e-05
  dt_min           = 1.0000e-07
  dt_max           = 5.0000e-05
  dt_fac           = 0.9
  res_dt           = 0.005 #!MFIX-GUI eq{float(1.0/200)}
  batch_wallclock  = 1296000.0
  chk_batchq_end   = .False.
  drag_type        = 'SQP_DIFELICE_HOLZER_SOMMERFELD'
  enable_dmp_log   = .False.
  energy_eq        = .True.
  nodesi           = 2
  nodesj           = 2
  nodesk           = 1
  term_buffer      = 180.0
  write_dashboard  = .False.
  full_log         = .True.
  momentum_x_eq(0) = .True.
  momentum_y_eq(0) = .True.
  momentum_z_eq(0) = .True.
  project_version  = '247'
  species_eq(0)    = .False.
  species_eq(1)    = .False.

### Physical parameters
  gravity_x = 0.0
  gravity_y = -9.81
  gravity_z = 0.0

### Cartesian grid
  cartesian_grid = .False.
  use_stl        = .False.

### Numeric
  detect_stall      = .False.
  max_inlet_vel_fac = 3.0
  max_nit           = 100
  norm_g            = 0
  tol_resid         = 5.0000e-03
  ### Discretization
  discretize(1)     = 2  leq_it(1) = 50
  discretize(2)     = 2  leq_it(2) = 20
  discretize(3)     = 2
  discretize(4)     = 2
  discretize(5)     = 2
  discretize(6)     = 2
  discretize(7)     = 0
  fpe_trap_invalid  = .False.
  fpe_trap_overflow = .False.
  fpe_trap_zero     = .False.
  leq_tol(1)        = 0.0001
  tol_resid_scalar  = 0.001
  tol_resid_t       = 0.001
  tol_resid_th      = 0.001
  tol_resid_x       = 0.001

### Geometry
  coordinates = 'CARTESIAN'
  imax        = 7
  jmax        = 70
  kmax        = 3
  x_max       = 0.1
  x_min       = 0
  y_max       = 1.0
  y_min       = 0
  z_max       = 0.04
  z_min       = 0
  no_k        = .False.

#### Fluid
  mu_g0      = 1.8000e-05
  mw_avg     = 29.0
  c_pg0      = 1005.0
  mu_g_model = 'CONSTANT'

#### Solids
  mmax = 1
  # Solid 1: Particle
    solids_model(1)     = 'SQP'
    d_p0(1)             = 0.012510760454394617
    ro_s0(1)            = 1442.5
    nmax_s(1)           = 0
    c_ps0(1)            = 1800.0
    des_em(1)           = 0.0
    e_young(1)          = 1.0000e+06
    e_young_actual(1)   = 1.0000e+06
    k_s0(1)             = 0.3
    ks_model(1)         = 'MUSSER'
    sqp_a(1)            = 0.0015
    sqp_b(1)            = 0.0015
    sqp_c(1)            = 0.006
    sqp_m(1)            = 2.0
    sqp_n(1)            = 8.0
    sqp_q1(1)           = 1.0
    sqp_q2(1)           = 0.0
    sqp_q3(1)           = 0.0
    sqp_q4(1)           = 0.0
    v_poisson(1)        = 0.3
    v_poisson_actual(1) = 0.3


### Initial conditions
  # Initial condition 1: init1
    ic_x_e(1)            = 0.1 #!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)            = 1.0 #!MFIX-GUI eq{float(ymax)}
    ic_z_b(1)            = 0.0 #!MFIX-GUI eq{float(zmin)}
    ic_z_t(1)            = 0.04 #!MFIX-GUI eq{float(zmax)}
    ic_ep_g(1)           = 0.99055
    ic_t_g(1)            = 298.15
    ic_u_g(1)            = 0.0
    ic_v_g(1)            = 0.0
    ic_w_g(1)            = 0.0
    ic_gama_rg(1)        = 0.0
    ic_p_star(1)         = 0.0
    ic_t_rg(1)           = 298.15
    ic_ep_s(1,1)         = 0.00945
    ic_t_s(1,1)          = 1273.15
    ic_theta_m(1,1)      = 0.0
    ic_u_s(1,1)          = 0.0
    ic_v_s(1,1)          = 0.0
    ic_w_s(1,1)          = 0.0
    ic_des_lattice(1,1)  = 'HEXA'
    ic_gama_rs(1,1)      = 0.0
    ic_sqp_random_q(1,1) = .True.
    ic_t_rs(1,1)         = 298.15


# Boundary conditions
  # Boundary condition 1: Bottom_inlet
    bc_type(1) = 'MI'
    bc_x_e(1)  = 0.1 #!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.04 #!MFIX-GUI eq{float(zmax)}
    bc_ep_g(1) = 1.0
    bc_p_g(1)  = 1.0132e+05
    bc_t_g(1)  = 298.15
    bc_u_g(1)  = 0.0
    bc_v_g(1)  = 2.2
    bc_w_g(1)  = 0.0
    # Solid 1: Particle
      bc_ep_s(1,1) = 0.0
      bc_t_s(1,1)  = 1273.15
      bc_u_s(1,1)  = 0.0
      bc_v_s(1,1)  = 0.0
      bc_w_s(1,1)  = 0.0

  # Boundary condition 2: Top_outlet
    bc_type(2)            = 'PO'
    bc_x_e(2)             = 0.1 #!MFIX-GUI eq{float(xmax)}
    bc_x_w(2)             = 0.0 #!MFIX-GUI eq{float(xmin)}
    bc_y_s(2)             = 1.0 #!MFIX-GUI eq{float(ymax)}
    bc_y_n(2)             = 1.0 #!MFIX-GUI eq{float(ymax)}
    bc_z_b(2)             = 0.0 #!MFIX-GUI eq{float(zmin)}
    bc_z_t(2)             = 0.04 #!MFIX-GUI eq{float(zmax)}
    bc_p_g(2)             = 101325.0
    bc_po_apply_to_des(2) = .True.

  # Boundary condition 3: Leftwall
    bc_type(3)   = 'NSW'
    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)    = 1.0 #!MFIX-GUI eq{float(ymax)}
    bc_z_b(3)    = 0.0 #!MFIX-GUI eq{float(zmin)}
    bc_z_t(3)    = 0.04 #!MFIX-GUI eq{float(zmax)}
    bc_c_t_g(3)  = 0.0
    bc_hw_t_g(3) = 14000.0
    bc_jj_ps(3)  = 0
    bc_tw_g(3)   = 298.15
    bc_uw_g(3)   = 0.0
    bc_vw_g(3)   = 0.0
    bc_ww_g(3)   = 0.0
    # Solid 1: Particle
      bc_c_t_s(3,1) = 0.0
      bc_tw_s(3,1)  = 298.15

  # Boundary condition 4: Rightwall
    bc_type(4)   = 'NSW'
    bc_x_e(4)    = 0.1 #!MFIX-GUI eq{float(xmax)}
    bc_x_w(4)    = 0.1 #!MFIX-GUI eq{float(xmax)}
    bc_y_s(4)    = 0.0 #!MFIX-GUI eq{float(ymin)}
    bc_y_n(4)    = 1.0 #!MFIX-GUI eq{float(ymax)}
    bc_z_b(4)    = 0.0 #!MFIX-GUI eq{float(zmin)}
    bc_z_t(4)    = 0.04 #!MFIX-GUI eq{float(zmax)}
    bc_c_t_g(4)  = 0.0
    bc_hw_t_g(4) = 14000.0
    bc_jj_ps(4)  = 0
    bc_tw_g(4)   = 298.15
    bc_uw_g(4)   = 0.0
    bc_vw_g(4)   = 0.0
    bc_ww_g(4)   = 0.0
    # Solid 1: Particle
      bc_c_t_s(4,1) = 0.0
      bc_tw_s(4,1)  = 298.15

  # Boundary condition 5: Frontwall
    bc_type(5)   = 'NSW'
    bc_x_e(5)    = 0.1 #!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)    = 1.0 #!MFIX-GUI eq{float(ymax)}
    bc_z_b(5)    = 0.0 #!MFIX-GUI eq{float(zmin)}
    bc_z_t(5)    = 0.0 #!MFIX-GUI eq{float(zmin)}
    bc_c_t_g(5)  = 0.0
    bc_hw_t_g(5) = 14000.0
    bc_jj_ps(5)  = 0
    bc_tw_g(5)   = 298.15
    bc_uw_g(5)   = 0.0
    bc_vw_g(5)   = 0.0
    bc_ww_g(5)   = 0.0
    # Solid 1: Particle
      bc_c_t_s(5,1) = 0.0
      bc_tw_s(5,1)  = 298.15

  # Boundary condition 6: Backwall
    bc_type(6)   = 'NSW'
    bc_x_e(6)    = 0.1 #!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)    = 1.0 #!MFIX-GUI eq{float(ymax)}
    bc_z_b(6)    = 0.04 #!MFIX-GUI eq{float(zmax)}
    bc_z_t(6)    = 0.04 #!MFIX-GUI eq{float(zmax)}
    bc_c_t_g(6)  = 0.0
    bc_hw_t_g(6) = 14000.0
    bc_jj_ps(6)  = 0
    bc_tw_g(6)   = 298.15
    bc_uw_g(6)   = 0.0
    bc_vw_g(6)   = 0.0
    bc_ww_g(6)   = 0.0
    # Solid 1: Particle
      bc_c_t_s(6,1) = 0.0
      bc_tw_s(6,1)  = 298.15


# VTK outputs
  write_vtk_files         = .True.
  time_dependent_filename = .True.
  vtu_dir                 = 'VTK'
  # VTK output 1: Entire domain
    vtk_filebase(1) = 'CELL'
    vtk_x_e(1)      = 0.1 #!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)      = 1.0 #!MFIX-GUI eq{float(ymax)}
    vtk_z_b(1)      = 0.0 #!MFIX-GUI eq{float(zmin)}
    vtk_z_t(1)      = 0.04 #!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_c_pg(1)     = .True.
    vtk_ep_g(1)     = .True.
    vtk_mu_g(1)     = .True.
    vtk_p_g(1)      = .True.
    vtk_ro_g(1)     = .True.
    vtk_t_g(1)      = .True.
    vtk_u_g(1)      = .True.
    vtk_v_g(1)      = .True.
    vtk_vel_g(1)    = .True.
    vtk_w_g(1)      = .True.

  # VTK output 2: Entire domain
    vtk_filebase(2)            = 'SQP'
    vtk_x_e(2)                 = 0.1 #!MFIX-GUI eq{float(xmax)}
    vtk_x_w(2)                 = 0.0 #!MFIX-GUI eq{float(xmin)}
    vtk_y_s(2)                 = 0.0 #!MFIX-GUI eq{float(ymin)}
    vtk_y_n(2)                 = 1.0 #!MFIX-GUI eq{float(ymax)}
    vtk_z_b(2)                 = 0.0 #!MFIX-GUI eq{float(zmin)}
    vtk_z_t(2)                 = 0.04 #!MFIX-GUI eq{float(zmax)}
    vtk_data(2)                = 'P'
    vtk_dt(2)                  = 0.01
    vtk_nxs(2)                 = 0
    vtk_nys(2)                 = 0
    vtk_nzs(2)                 = 0
    vtk_select_mode(2)         = 'C'
    vtk_part_angular_vel(2)    = .True.
    vtk_part_c_ps(2)           = .True.
    vtk_part_col_force(2)      = .True.
    vtk_part_coordination(2)   = .True.
    vtk_part_density(2)        = .True.
    vtk_part_diameter(2)       = .True.
    vtk_part_drag_coeff(2)     = .True.
    vtk_part_orientation(2)    = .True.
    vtk_part_overlap_stats(2)  = .True.
    vtk_part_residence_time(2) = .True.
    vtk_part_temp(2)           = .True.
    vtk_part_vel(2)            = .True.
    vtk_part_volume(2)         = .True.


# SPx outputs
  spx_dt(1)  = 100.0
  spx_dt(2)  = 100.0
  spx_dt(3)  = 100.0
  spx_dt(4)  = 100.0
  spx_dt(5)  = 100.0
  spx_dt(6)  = 100.0
  spx_dt(7)  = 100.0
  spx_dt(8)  = 100.0
  spx_dt(9)  = 100.0
  spx_dt(10) = 100.0
  spx_dt(11) = 100.0

### Residuals
  group_resid     = .True.
  resid_string(1) = 'P1'
  resid_string(2) = 'U1'
  resid_string(3) = 'V1'

### Discrete element model
  des_coll_model            = 'HERTZIAN'
  des_conv_corr             = 'GUNN'
  des_en_input(1)           = 0.9
  des_en_wall_input(1)      = 0.9
  des_en_wall_input(2)      = 0.5
  des_epg_clip              = 0.1
  des_et_input(1)           = 0.5
  des_et_wall_input(1)      = 0.5
  des_et_wall_input(2)      = 0.5
  des_etat_fac              = 0.5
  des_etat_w_fac            = 0.5
  des_explicitly_coupled    = .False.
  des_interp_mean_fields    = .True.
  des_interp_on             = .True.
  des_interp_scheme         = 'DPVM_SATELLITE'
  des_neighbor_search       = 4
  desgridsearch_imax        = 5
  desgridsearch_jmax        = 30
  desgridsearch_kmax        = 2
  ew_young                  = 1.0000e+06
  ew_young_actual           = 1.0000e+06
  kt_fac                    = 0.28571428571429 #!MFIX-GUI eq{float(2/7)}
  kt_w_fac                  = 0.28571428571429 #!MFIX-GUI eq{float(2/7)}
  mew                       = 0.4
  mew_w                     = 0.3
  neighbor_search_n         = 20
  neighbor_search_rad_ratio = 1.0
  nfactor                   = 0
  particles                 = 2782
  vw_poisson                = 0.3
  vw_poisson_actual         = 0.3

### Two-fluid model
  friction_model = 'SCHAEFFER'
  phi_w          = 11.31

### Output control
nlog           = 25
write_part_out = .True.


# MFIX-GUI section
#!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.1, "ymin": 0, "ymax": 1.0, "zmin": 0, "zmax": 0.04, "km": 1000.0,
#!MFIX-GUI   "m": 1.0, "cm": 0.01, "mm": 0.001, "um": {"__type__": "equation",
#!MFIX-GUI   "eq": "float(1e-06)"}, "mile": 1609.34, "yard": 0.9144, "ft":
#!MFIX-GUI   0.3048, "ins": 0.0254}}
#!MFIX-GUI gui_version = 24.3.1
#!MFIX-GUI project_notes = ""
#!MFIX-GUI author = tutorial
#!MFIX-GUI geometry = {"geometry_dict": {}, "tree": {}}
#!MFIX-GUI regions_dict = {"order": ["Bottom_inlet", "Top_outlet", "Entire
#!MFIX-GUI   domain", "init1", "Leftwall", "Rightwall", "Frontwall",
#!MFIX-GUI   "Backwall"], "regions": {"Bottom_inlet": {"to": [{"__type__":
#!MFIX-GUI   "equation", "eq": "float(xmax)"}, {"__type__": "equation", "eq":
#!MFIX-GUI   "float(ymin)"}, {"__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": "#6b70c5", "type":
#!MFIX-GUI   "XZ-plane"}, "Top_outlet": {"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(xmin)"}, {"__type__":
#!MFIX-GUI   "equation", "eq": "float(ymax)"}, {"__type__": "equation", "eq":
#!MFIX-GUI   "float(zmin)"}], "color": "#98c2da", "type": "XZ-plane"}, "Entire
#!MFIX-GUI   domain": {"to": [{"__type__": "equation", "eq": "float(xmax)"},
#!MFIX-GUI   {"__type__": "equation", "eq": "float(ymax)"}, {"__type__":
#!MFIX-GUI   "equation", "eq": "float(zmax)"}], "from": [{"__type__":
#!MFIX-GUI   "equation", "eq": "float(xmin)"}, {"__type__": "equation", "eq":
#!MFIX-GUI   "float(ymin)"}, {"__type__": "equation", "eq": "float(zmin)"}],
#!MFIX-GUI   "color": "#ac999d"}, "init1": {"to": [{"__type__": "equation",
#!MFIX-GUI   "eq": "float(xmax)"}, {"__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": "#996edd"}, "Leftwall":
#!MFIX-GUI   {"to": [{"__type__": "equation", "eq": "float(xmin)"}, {"__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": "#ae9287",
#!MFIX-GUI   "type": "YZ-plane"}, "Rightwall": {"to": [{"__type__": "equation",
#!MFIX-GUI   "eq": "float(xmax)"}, {"__type__": "equation", "eq":
#!MFIX-GUI   "float(ymax)"}, {"__type__": "equation", "eq": "float(zmax)"}],
#!MFIX-GUI   "from": [{"__type__": "equation", "eq": "float(xmax)"},
#!MFIX-GUI   {"__type__": "equation", "eq": "float(ymin)"}, {"__type__":
#!MFIX-GUI   "equation", "eq": "float(zmin)"}], "color": "#c67088", "type":
#!MFIX-GUI   "YZ-plane"}, "Frontwall": {"to": [{"__type__": "equation", "eq":
#!MFIX-GUI   "float(xmax)"}, {"__type__": "equation", "eq": "float(ymax)"},
#!MFIX-GUI   {"__type__": "equation", "eq": "float(zmin)"}], "from":
#!MFIX-GUI   [{"__type__": "equation", "eq": "float(xmin)"}, {"__type__":
#!MFIX-GUI   "equation", "eq": "float(ymin)"}, {"__type__": "equation", "eq":
#!MFIX-GUI   "float(zmin)"}], "color": "#ddb39e", "type": "XY-plane"},
#!MFIX-GUI   "Backwall": {"to": [{"__type__": "equation", "eq": "float(xmax)"},
#!MFIX-GUI   {"__type__": "equation", "eq": "float(ymax)"}, {"__type__":
#!MFIX-GUI   "equation", "eq": "float(zmax)"}], "from": [{"__type__":
#!MFIX-GUI   "equation", "eq": "float(xmin)"}, {"__type__": "equation", "eq":
#!MFIX-GUI   "float(ymin)"}, {"__type__": "equation", "eq": "float(zmax)"}],
#!MFIX-GUI   "color": "#ab65df", "type": "XY-plane"}}}
#!MFIX-GUI bc_regions = [[[1], ["Bottom_inlet"]], [[2], ["Top_outlet"]], [[3],
#!MFIX-GUI   ["Leftwall"]], [[4], ["Rightwall"]], [[5], ["Frontwall"]], [[6],
#!MFIX-GUI   ["Backwall"]]]
#!MFIX-GUI ic_regions = [[[1], ["init1"]]]
#!MFIX-GUI vtk_regions = [[[1], ["Entire domain"]], [[2], ["Entire domain"]]]
#!MFIX-GUI OMP_NUM_THREADS = 1
#!MFIX-GUI submit_to_queue = 0
#!MFIX-GUI created_date = 2019-09-05 14:16
#!MFIX-GUI modified_time = 2026-09-16 14:51
#!MFIX-GUI modified_by = gaox|highx|cgw|dietikej|vikrant|gunjan|user
#!MFIX-GUI retained_keys = {"k_s0(1)": 1.0, "des_usr_var_size": 9, "ro_g0":
#!MFIX-GUI   1.2005, "ic_ep_s(1,1)": 0.02, "ic_des_np(1,1)": 2800, "dt_fac":
#!MFIX-GUI   0.9, "dt_min": 1e-07, "dt_max": 5e-05, "part_out_zero_vel": true,
#!MFIX-GUI   "c_pg0": 1005.0, "k_g0": 0.0257, "particles": 2782, "mu_g0":
#!MFIX-GUI   1.8e-05, "mw_avg": 29.0}
#!MFIX-GUI BUILD_DMP = 0
#!MFIX-GUI BUILD_INTERACTIVE = PYMFIX
#!MFIX-GUI BUILD_PARALLEL = 1
#!MFIX-GUI BUILD_SMP = 0
#!MFIX-GUI BUILD_TYPE = RELWITHDEBINFO
#!MFIX-GUI queue_template = {"JOB_NAME": "FB_SQP", "CORES": 4, "QUEUE":
#!MFIX-GUI   "general", "LONG": "", "template": "Joule 3"}
#!MFIX-GUI mfix_exe = /home/gunjan/miniforge3/envs/mfix-24.3.1/bin/mfixsolver_dmp
#!MFIX-GUI mpirun_flags = --use-hwthread-cpus
#!MFIX-GUI project_version = 247
#!MFIX-GUI solids_phase_name(1) = Particle