mfix.dat

########################################################################
#                                                                      #
# Author:  Sreekanth Pannala                            Date: 01-2010  #
# Modified: J. Musser                                   Date: 05-2014  #
#                                                                      #
#  Description:                                                        #
#  Freely falling particles under gravity that collide with a wall.    #
#  This is an extension of the regular freely-falling particle test    #
#  case to run parallel.                                               #
#                                                                      #
#  References:                                                         #
#  R. Garg, J. Galvin, T. Li, and S. Pannala, Documentation of open-   #
#     source MFIX-DEM software for gas-solids flows," from URL:        #
#     https://mfix.netl.doe.gov/documentation/dem_doc_2012-1.pdf,      #
#                                                                      #
########################################################################

  RUN_NAME = 'squeeze'
  DESCRIPTION = 'const vel'

#_______________________________________________________________________
# RUN CONTROL SECTION

  RUN_TYPE = 'NEW'
  UNITS = 'SI'

  TIME =  0.0
  TSTOP = 5 # 0.5s seems like a realistic dispense time, however squeeze time could be larger than this

  DT =     1.0E-4
  DT_MAX = 1.0E-3

  ENERGY_EQ = .F.

  SPECIES_EQ(0:1) =  .F.  .F.

  MOMENTUM_X_EQ(0:1) = .F.  .F.
  MOMENTUM_Y_EQ(0:1) = .F.  .F.
  MOMENTUM_Z_EQ(0:1) = .F.  .F.

  CALL_USR = .T.

#_______________________________________________________________________
# Physical Parameters
  GRAVITY_Y = -9.81

#_______________________________________________________________________
# NUMERICAL SECTION
 MAX_NIT = 50

#_______________________________________________________________________
# GEOMETRY SECTION

  COORDINATES = 'CARTESIAN'

   XLENGTH  = 0.02   IMAX = 30 # 21.05
   YLENGTH  = 0.001861   JMAX = 2 # 1.3 0.0012139
   ZLENGTH  = 0.0008035   KMAX = 1 # 2.1 0.0019967

# IMAX, JMAX and KMAX set such that NODESI, NODESJ and NODESK yield at least 3 cells per DMP division
  #XLENGTH  = 0.01   IMAX = 15
  #YLENGTH  = 0.0008   JMAX = 2
  #ZLENGTH  = 0.0028   KMAX = 12


#_______________________________________________________________________
# MATERIAL SECTION

! Gas Section
!---------------------------------------------------------------------//

  RO_g0 = 0.0
  MU_g0 = 0.0

! PARTICLE SECTION
!---------------------------------------------------------------------//
  MMAX = 1

  PARTICLES = 966                ! Number of particles
  gener_part_config = .False.
! Friction coefficients. (1)
  MEW =   0.5                  ! particle-particle
  MEW_W = 0.5                  ! particle-wall

  DES_COLL_MODEL = 'LSD' # is this required to be specified? (see dem doc, pg 11)

  DES_NEIGHBOR_SEARCH =  4     ! Grid based neighbor search
  NEIGHBOR_SEARCH_N   =  10     ! Steps between neighbor search

  DESGRIDSEARCH_IMAX = 50
  DESGRIDSEARCH_JMAX = 4
  DESGRIDSEARCH_KMAX = 2 # 3.5

  DES_EN_INPUT = 0.7
  DES_EN_WALL_INPUT = 0.7 # Particles must almost immediately stick the substrate
  DES_INTG_METHOD = 'ADAMS_BASHFORTH' # 'EULER'
  KN = 1.0E+5
  KT_FAC = @(2/7)
  KN_W = 1.0E+5
  KT_W_FAC = @(2/7)

!......................................................................!
! The following keywords are required inputs for the simulation. They  !
! are not specified here because they are passed as run time arguments.!
!                                                                      !
! Normal collision spring constant. (N/m)                              !
! KN =   xxxxx                 ! particle-particle                     !
! KN_W = xxxxx                 ! particle-wall                         !
!                                                                      !
! Restitution coefficient. (1)                                         !
! DES_EN_INPUT =      xxxxx    ! particle-particle                     !
! DES_EN_WALL_INPUT = xxxxx    ! particle-wall                         !
!......................................................................!


! Solids phase 1
!---------------------------------------------------------------------//
  SOLIDS_MODEL(1) = 'DEM'

  D_p0(1) =    100E-06   ! (m)
  RO_s0(1) =   9000.00   ! (kg/m)


# Drag force model
  #MODEL_B = .T.
  DRAG_TYPE = 'GIDASPOW' # GIDASPOW_PCF

#_______________________________________________________________________
# INITIAL CONDITIONS SECTION

  # Full domain	
  IC_X_w(1) = 0.0 # Initial solids volume fraction
  IC_X_e(1) = 0.02
  IC_Y_s(1) = 0.0
  IC_Y_n(1) = 1.861E-3 #1.2139E-3
  IC_Z_b(1) = 0.0
  IC_Z_t(1) = 0.8035E-3 #1.9967E-3
  IC_EP_g(1) = 1.0

 # IC_X_w(2) = 0.00 
 # IC_X_e(2) = 0.0016
 # IC_Y_s(2) = 0.0
 # IC_Y_n(2) = 0.0028
 # IC_Z_b(2) = 0.0
 # IC_Z_t(2) = 0.0008

  #IC_EP_S(2,1) = 0.05
  #IC_EP_g(2) = 0.95
  
 # IC_U_S(2,1) = 0.0
 # IC_V_S(2,1) = 0.0
 # IC_W_S(2,1) = 0.0
  #IC_EP_S(2,1) = 0.0	
 # IC_THETA_M(2,1) = 0.0
 # IC_DES_FIT_TO_REGION = .F.

	
  #_______________________________________________________________________
# BOUNDARY CONDITIONS SECTION: assuming (0,0,0) starts at a corner
# North wall
  BC_X_w(1) = 0.0
  BC_X_e(1) = 0.02
  BC_Y_s(1) = 1.861E-3 #0.0012139
  BC_Y_n(1) = 1.861E-3 #0.0012139
  BC_Z_b(1) = 0.0
  BC_Z_t(1) = 0.8035E-3 #0.0019967
 
  BC_TYPE(1) = 'NSW'

  #BC_EP_S(1,1) = 0.2 #; # Solids vol fraction at inlet
  #BC_U_S(1,1) = 0.0
  #BC_V_S(1,1) = -0.01 #; # m^3/s, equal to 60 cc/min
  #BC_W_S(1,1) = 0.0

# West face
  BC_X_w(2) = 0.0
  BC_X_e(2) = 0.0
  BC_Y_s(2) = 0.0
  BC_Y_n(2) = 1.861E-3 #0.0012139 # y_max
  BC_Z_b(2) = 0.000
  BC_Z_t(2) = 0.8035E-3 #0.0019967 # z_max
 
  BC_TYPE(2) = 'NSW'

# East face
  BC_X_w(3) = 0.02
  BC_X_e(3) = 0.02
  BC_Y_s(3) = 0.0
  BC_Y_n(3) = 1.861E-3 #0.0012139
  BC_Z_b(3) = 0.0
  BC_Z_t(3) = 0.8035E-3 # 0.0019967
 
  BC_TYPE(3) = 'NSW'

# South face
  BC_X_w(4) = 0.0
  BC_X_e(4) = 0.02
  BC_Y_s(4) = 0.0
  BC_Y_n(4) = 0.0
  BC_Z_b(4) = 0.0
  BC_Z_t(4) = 0.8035E-3 #0.0019967
 
  BC_TYPE(4) = 'NSW'

# Back face
  BC_X_w(5) = 0.0
  BC_X_e(5) = 0.02
  BC_Y_s(5) = 0.0
  BC_Y_n(5) = 1.861E-3 #0.0012139
  BC_Z_b(5) = 0.0
  BC_Z_t(5) = 0.0
 
  BC_TYPE(5) = 'NSW'

# Front face
  BC_X_w(6) = 0.0
  BC_X_e(6) = 0.02
  BC_Y_s(6) = 0.0
  BC_Y_n(6) = 1.861E-3 #0.0012139
  BC_Z_b(6) = 0.8035E-3 #0.0019967
  BC_Z_t(6) = 0.8035E-3 #0.0019967
 
  BC_TYPE(6) = 'NSW'

# Periodic BCs
  CYCLIC_X = .F.
  CYCLIC_Y = .F.
  CYCLIC_Z = .T.
  
  CYCLIC_X_PD = .F.
  CYCLIC_Y_PD = .F.
  CYCLIC_Z_PD = .F.		
# ---------------------------------------------------- #

# Side walls (possible with a different cft of restitution?)
 # IS_X_w(1) = 0.0025
 # IS_X_e(1) = 0.0025
 # IS_Y_s(1) = 0.000
 # IS_Y_n(1) = 0.002
 # IS_Z_b(1) = 0.00
 # IS_Z_t(1) = 0.001
 # IS_TYPE(1) = 'IMPERMEABLE' 

#_______________________________________________________________________
# USER OUTPUT CONTROL

#_______________________________________________________________________
# OUTPUT CONTROL SECTION

  RES_DT = 1.0d-2  ! interval to update restart (.RES) file

  FULL_LOG = .T.  ! display residuals on screen
  NLOG = 250      ! time steps between updates to (.LOG) file

! Interval at which .SPX files are written
  SPX_DT(1:9) = 9*1.0d3

  GROUP_RESID = .T.

  PRINT_DES_DATA = .F.


#_______________________________________________________________________
# DMP SETUP

  NODESI = 1
  NODESJ = 1
  NODESK = 1

  ENABLE_DMP_LOG = .T.
  CHK_BATCHQ_END = .T.
  
#------------Wall Time-----------

  BATCH_WALLCLOCK = 432000.0 # 5 days
  TERM_BUFFER = 300.0  
  
  
#------------VTK FILES-----------
write_vtk_files = .True.
time_dependent_filename = .True.
vtk_data(1) = 'P'
vtk_x_w(1) = 0.00
vtk_x_e(1) = 0.02
vtk_y_s(1) = 0.00	
vtk_y_n(1) = 1.861E-3 #0.0012139
vtk_z_b(1) = 0.0
vtk_z_t(1) = 0.8035E-3 #0.0019967
vtk_filebase(1) = 'Particle'
vtk_dt(1) = 1E-4
vtk_nxs(1) = 0
vtk_nys(1) = 0
vtk_nzs(1) = 0
vtk_part_diameter(1) = .True.
vtk_part_vel(1) = .True.
vtk_part_id(1) = .True.