/*---------------------------------------------------------------------------*\
Granular bed of Quartzite
\*---------------------------------------------------------------------------*/

FoamFile {
  version     2.0;
  format      ascii;
  class       dictionary;
  object      constantProperties;
}
// * * * * * *  Units * * * * * [kg m s K mol A cd] * * * * * * * * * * * * * //
// e.g. W: kg m^2 s^{-3}        [1 2 -3 0 0 0 0]

/****           Universal constants                                        ****/
R               R               [1 2 -2 -1 -1 0 0]      8.314471469;
sigmaPlanck     sigmaPlanck     [1 0 -3 -1 0 0 0]       5.6697e-8;

/***            MATERIAL PROPERTIES                                         ***/
/***            Initial densities and volume fractions                      ***/
// number of phases (index starts at 1)
nSolidPhases    1;
// intrinsic density of the quartzite
rhoI[1]         rhoI[1]         [1 -3 0 0 0 0 0]        2595;
// volume fraction of the quartzite grain = 1-0.4=0.6
epsI[1]         epsI[1]         [0 0 0 0 0 0 0]         0.6;

/***            Averaged momentum-conservation related data (Darcy's law)   ***/
// virgin material permeability (second order tensor)
K_v             K_v             [0 2 0 0 0 0 0]         (5.5e-7 0 0 0 5.5e-7 0 0 0 5.5e-7);
// char permeability (second order tensor)
K_c             K_c             [0 2 0 0 0 0 0]         (5.5e-7 0 0 0 5.5e-7 0 0 0 5.5e-7);
// virgin Forchheimer coefficient (second order tensor)
Beta_v          Beta_v          [0 -1 0 0 0 0 0]        (7.4e2 0 0 0 7.4e2 0 0 0 7.4e2);
// char Forchheimer coefficient (second order tensor)
Beta_c          Beta_c          [0 -1 0 0 0 0 0]        (7.4e2 0 0 0 7.4e2 0 0 0 7.4e2);
// virgin material porosity
eps_g_v         eps_g_v         [0 0 0 0 0 0 0]         0.4;
// char porosity
eps_g_c         eps_g_c         [0 0 0 0 0 0 0]         0.4;

/***  Convective heat exchange between solid and gas phases for 2 T model   ***/
// heat transfer coefficient
Hv0             Hv0             [1 -1 -3 -1 0 0 0]      3000;

// 1- Express the main directions (ijk) of the diagonal conductivity matrix in the basis of the mesh (xyz)
//              The transposed conductivity passage matrix (tP) expresses (i,j,k) coordinates in the basis (x,y,z)
//              where (i,j,k) is a direct orthonormal basis of the diagonal conductivity matrix
//              and (x,y,z) is the direct orthonormal basis of the mesh (in case of doubt, mesh directions can be seen in paraview)
//                                                      (i j k)                 ex. rotation a (in radians) around axis z
tP              tP              [0 0 0 0 0 0 0]         (1 0 0   // x           (cosa -sina 0
    0 1 0   // y            sina  cosa 0
    0 0 1); // z            0       0  1)
// 2 - Linear factors
kiCoef          kiCoef          [0 0 0 0 0 0 0]         1;       // to multiply column ki of the input files 'char' and 'virgin' by a linear factor: ki' = kiCoef*ki
kjCoef          kjCoef          [0 0 0 0 0 0 0]         1;       // idem for kj
kkCoef          kkCoef          [0 0 0 0 0 0 0]         1;       // idem for kk

/***               Mechanical Properties                                    ***/
nu_v    nu_v    [0 0 0 0 0 0 0]      0.17;    // virgin material Poisson's ratio
nu_c    nu_c    [0 0 0 0 0 0 0]      0.17;    // char Poisson's ration
E_v     E_v     [1 -1 -2 0 0 0 0]    42e9;    // virgin material Young modulus : 70e9 * 0.6
E_c     E_c     [1 -1 -2 0 0 0 0]    42e9;    // char Young modulus :            70e9 * 0.6
alpha_v alpha_v [0 0 0 -1 0 0 0]     6.6e-6;  // virgin thermal expansion coefficient
alpha_c alpha_c [0 0 0 -1 0 0 0]     6.6e-6;  // char thermal expansion coefficient
xi      xi      [0 0 0 0 0 0 0]      0.5;     // Jomaa Function (pyrolysis shrinkage)

// ************************************************************************* //
