mfix.dat is a generic input deck having parameters: theta1, theta2, theta3, theta4, theta5 and x1. Theta values are held constant (provided below), only x1 is varied.

theta1: Pressure linear scale factor (P0) - 1.6100105968e+01
theta2: Exponential factor (beta) - 2.0394023796e+00
theta3: Statistical weight - 1.0506421399e+01
theta4: Void fraction at maximum packing (ep_star) - 4.0121517743e-01
theta5: solids slip velocity factor (mppic_velfac_coeff) - 5.3145624742e-01
x1: Initial solids concentration - varied

1. Input_post_n119_Settling_11June2020.csv - Sampling locations from midpoints (of initial solids concentration) of OLH based design of experiments (n=120).

2. replace.py - To change ic_ep_g(2) based on ic_ep_s(2). This was not happening automatically through mfix+nodeworks. Input for ic_ep_s(2) is given through Input_Nodeworks_Settling.csv. Initially ic_ep_g(2) is set to 0.85 once the design of experiments is generated through mfix+nodeworks.

3. readAV_postmfix.in - Input deck for postmfix to extract volume fraction. The results are averaged along x and written at every y and z location. This is copied to all the sub-directories.

4. runjob - mfix slurm submission script. This is copied to all the sub-directories.

5. runjob.postmfix - postmfix slurm submission script. Input is an array denoting each sub-directory from which postmfix is called sequentially.
