12.1.8. Initial conditions (IC)¶

Each initial condition is specified over a rectangular region (or pie-shaped for cylindrical coordinates) that corresponds to the scalar numerical grid. These are 3D regions: X_w X_e, Y_s Y_n, and Z_t Z_b. The region is defined by the constant coordinates of each of the six faces, which may be specified as the physical coordinates or the cell indices. The physical coordinates are easier to specify than the cell indices. If cell sizes are not small enough to resolve a region specified using physical coordinates, MFIX will indicate this problem with an error message.

In cylindrical coordinates, when the theta direction crosses the 0 value, split that region into two regions: e.g., Split a region spanning 1.9 pi to 0.1 pi as 1.9 pi to 2 pi and 0 to 0.1 pi.

Initial condition regions may overlap. When an overlap occurs, MFIX uses the conditions specified for the higher IC number.

12.1.8.1. IC_DES_ABORT_ON_FAILED_SEEDING(IC, PHASE)¶

Data Type: LOGICAL

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Flag to abort if unsuccessful seeding occurs.

12.1.8.2. IC_DES_LATTICE(IC, PHASE)¶

Data Type: CHARACTER

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Type of initial lattice distribution (simple cubic or hexagonal) to the entire IC region.

Table 12.47 Valid values¶
Name	Default?	Description
`HEXA`	◉	Particles are distributed on a hexagonal close-packed lattice.
`CUBIC`		Particles are distributed on a cubical lattice.

12.1.8.3. IC_DES_NP(IC, PHASE)¶

Data Type: INTEGER

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial particle count within the entire IC region.

12.1.8.4. IC_DES_RAND(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Random factor applied to particle positions within the lattice distribution.

12.1.8.5. IC_DES_RAND_FACTOR_X(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Random factor in x-direction (multiplier of IC_DES_RAND)

12.1.8.6. IC_DES_RAND_FACTOR_Y(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Random factor in y-direction (multiplier of IC_DES_RAND)

12.1.8.7. IC_DES_RAND_FACTOR_Z(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Random factor in z-direction (multiplier of IC_DES_RAND)

12.1.8.8. IC_DES_SM(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial solids mass within the entire IC region.

12.1.8.9. IC_DES_SPACE_FACTOR_X(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Spacing factor in x-direction (multiplier of IC_DES_SPACING).

12.1.8.10. IC_DES_SPACE_FACTOR_Y(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Spacing factor in y-direction (multiplier of IC_DES_SPACING).

12.1.8.11. IC_DES_SPACE_FACTOR_Z(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Spacing factor in z-direction (multiplier of IC_DES_SPACING).

12.1.8.12. IC_DES_SPACING(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Spacing between particle within initial lattice distribution (expressed as a fraction of particle diameter) to the entire IC region.

12.1.8.13. IC_E_TURB_G(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Initial value of Epsilon in K-Epsilon equation.

12.1.8.14. IC_EP_G(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Initial void fraction in the IC region.

12.1.8.15. IC_EP_S(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial solids volume fraction of solids phase in the IC region. This may be specified in place of IC_ROP_s.

12.1.8.16. IC_GAMA_RG(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Gas phase radiation coefficient in the IC region. Modify file rdtn2.inc to change the source term.

12.1.8.17. IC_GAMA_RS(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Solids phase radiation coefficient in the IC region. Modify file energy_mod.f to change the source term.

12.1.8.18. IC_GSP_RANDOM_Q(IC, PHASE)¶

Data Type: LOGICAL

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initialize GSP with random orientation (quaternion) in the IC region.

12.1.8.19. IC_I_E(IC)¶

Data Type: INTEGER

$1 \leq I C \leq 500$

I index of the east-most wall.

12.1.8.20. IC_I_W(IC)¶

Data Type: INTEGER

$1 \leq I C \leq 500$

I index of the west-most wall.

12.1.8.21. IC_J_N(IC)¶

Data Type: INTEGER

$1 \leq I C \leq 500$

J index of the north-most wall.

12.1.8.22. IC_J_S(IC)¶

Data Type: INTEGER

$1 \leq I C \leq 500$

J index of the south-most wall.

12.1.8.23. IC_K_B(IC)¶

Data Type: INTEGER

$1 \leq I C \leq 500$

K index of the bottom-most wall.

12.1.8.24. IC_K_T(IC)¶

Data Type: INTEGER

$1 \leq I C \leq 500$

K index of the top-most wall.

12.1.8.25. IC_K_TURB_G(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Initial value of K in K-Epsilon equation.

12.1.8.26. IC_L_SCALE(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Turbulence length scale in the IC region.

12.1.8.27. IC_OVERRIDE_PART_INPUT(IC, PHASE)¶

Data Type: LOGICAL

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial conditions take precedence over data in particle_inpupt.dat file for this region and solids phase.

12.1.8.28. IC_P_G(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Initial gas pressure in the IC region. If this quantity is not specified, MFiX will set up a hydrostatic pressure profile, which varies only in the y-direction.

12.1.8.29. IC_P_STAR(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Initial solids pressure in the IC region. Usually, this value is specified as zero.

12.1.8.30. IC_PIC_CONST_STATWT(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Flag to specify the initial constant statistical weight for computational particles/parcels. Actual number of parcels will be automatically computed.

12.1.8.31. IC_PSD_MAX_DP(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial maximum particle size of particle size distribution of solids phase in the IC region.

12.1.8.32. IC_PSD_MEAN_DP(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial mean value of particle size distribution of solids phase in the IC region. For a log-normal psd, The mean diameter (µ) and standard deviation (𝜎) will be used to compute the log-normal parameters µ_ln and 𝜎_ln : µ_ln = ln(µ²/√(µ² + 𝜎²)) and 𝜎_ln = √(ln(1 + (𝜎/µ)²)))

12.1.8.33. IC_PSD_MIN_DP(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial minimum particle size of particle size distribution of solids phase in the IC region.

12.1.8.34. IC_PSD_STDEV(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial standard deviation of particle size distribution of solids phase in the IC region. For a log-normal psd, The mean diameter (µ) and standard deviation (𝜎) will be used to compute the log-normal parameters µ_ln and 𝜎_ln : µ_ln = ln(µ²/√(µ² + 𝜎²)) and 𝜎_ln = √(ln(1 + (𝜎/µ)²)))

12.1.8.35. IC_PSD_TYPE(IC, PHASE)¶

Data Type: CHARACTER

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial particle size distribution type of solids phase in the IC region.

Table 12.48 Valid values¶
Name	Default?	Description
`MONO`	◉	Uniform particle size.
`NORMAL`		Normal (Gaussian) particle size distribution.
`LOG_NORMAL`		Log-normal particle size distribution.
`CUSTOM`		Particle size distribution specified in external file.

12.1.8.36. IC_ROP_S(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial bulk density (rop_s = ro_s x ep_s) of solids phase in the IC region. Users need to specify this IC only for polydisperse flow (MMAX > 1). Users must make sure that summation of ( IC_ROP_s(ic,m) / RO_s(m) ) over all solids phases is equal to ( 1.0 - IC_EP_g(ic)).

12.1.8.37. IC_SCALAR(IC, SCALAR EQ.)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq S c a l a r E q . \leq 100$

Initial value of Scalar n, assigned to scalar equation.

12.1.8.38. IC_SQP_RANDOM_Q(IC, PHASE)¶

Data Type: LOGICAL

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initialize SQP with random orientation (quaternion) in the IC region.

12.1.8.39. IC_T_G(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Initial gas phase temperature in the IC region.

12.1.8.40. IC_T_RG(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Gas phase radiation temperature in the IC region.

12.1.8.41. IC_T_RS(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Solids phase radiation temperature in the IC region.

12.1.8.42. IC_T_S(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial solids phase temperature in the IC region.

12.1.8.43. IC_THETA_M(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial solids phase granular temperature in the IC region.

12.1.8.44. IC_TYPE(IC)¶

Data Type: CHARACTER

$1 \leq I C \leq 500$

Type of initial condition. Mainly used in restart runs to overwrite values read from the .RES file by specifying it as _PATCH_. The user needs to be careful when using the _PATCH_ option, since the values from the .RES file are overwritten and no error checking is done for the patched values.

12.1.8.45. IC_U_G(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Initial x-axial component of gas velocity in the IC region.

12.1.8.46. IC_U_S(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial x-axial component of solids-phase velocity in the IC region.

12.1.8.47. IC_V_G(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Initial y-axial component of gas velocity in the IC region.

12.1.8.48. IC_V_S(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial y-axial component of solids-phase velocity in the IC region.

12.1.8.49. IC_W_G(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Initial z-axial component of gas velocity in the IC region.

12.1.8.50. IC_W_S(IC, PHASE)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$

Initial z-axial component of solids-phase velocity in the IC region.

12.1.8.51. IC_X_E(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

X coordinate of the east face.

12.1.8.52. IC_X_G(IC, SPECIES)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq S p e c i e s \leq 100$

Initial mass fraction of gas species.

12.1.8.53. IC_X_S(IC, PHASE, SPECIES)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$
$1 \leq P h a s e \leq 10$
$1 \leq S p e c i e s \leq 100$

Initial mass fraction of solids species.

12.1.8.54. IC_X_W(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

X coordinate of the west face.

12.1.8.55. IC_Y_N(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Y coordinate of the north face.

12.1.8.56. IC_Y_S(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Y coordinate of the south face.

12.1.8.57. IC_Z_B(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Z coordinate of the bottom face.

12.1.8.58. IC_Z_T(IC)¶

Data Type: DOUBLE PRECISION

$1 \leq I C \leq 500$

Z coordinate of the top face.