mfix.dat is a generic input deck having parameters: theta1, theta2, theta3, theta4, theta5 and x1. These are used to generate the design of experiments using nodeworks.

theta1: Pressure linear scale factor (P0)
theta2: Exponential factor (beta)
theta3: Statistical weight
theta4: Void fraction at maximum packing (ep_star)
theta5: solids slip velocity factor (mppic_velfac_coeff)

1. Input_OLH_n120_i6_11June2020.csv - Design of experiments from Aytekin based on OLH sampling.

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.

The following command could be used for copying to all subdirectories:
echo 1000* | xargs -n 1 cp readAV_postmfix.in .

4. runjob - mfix slurm submission script.

5. runjob.postmfix - postmfix slurm submission script. Input is an array denoting each sub-directory from which postmfix is called sequentially. Once the script is executed, the results are written to VF.dat file under subfolders 1000000-1000119.

6. Height_stat_postmfix_final.py - Postprocessing script to extract y1: Settling shock location, y2: Filling shock location and y3: Void fraction at the first cell from the bottom.

------------------------------------------------------------------------------

-- Avinash Vaidheeswaran
-- 11/03/22
