/*---------------------------------------------------------------------------*\
Constant property directory. Update as needed.
\*---------------------------------------------------------------------------*/

FoamFile {
  version     2.0;
  format      ascii;
  class       dictionary;
  location    "constant/porousMat/materialProperties";
  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                          ***/
nSolidPhases    2;                                              // number of phases (index starts at 1)
rhoI[1]         rhoI[1]         [1 -3 0 0 0 0 0]        1600;   // intrinsic density of the fibers
epsI[1]         epsI[1]         [0 0 0 0 0 0 0]         0.1;    // volume fraction of the fibrous preform
rhoI[2]         rhoI[2]         [1 -3 0 0 0 0 0]        1200;   // initial intrinsic density of the virgin matrix
epsI[2]         epsI[2]         [0 0 0 0 0 0 0]         0.1;    // volume fraction of the virgin matrix


/***            Averaged momentum-conservation related data (Darcy's law)       ***/
K_v             K_v             [0 2 0 0 0 0 0]         (1.6e-11 0 0 0 1.6e-11 0 0 0 1.6e-11); // virgin material permeability (second order tensor)
K_c             K_c             [0 2 0 0 0 0 0]         (2e-11 0 0 0 2e-11 0 0 0 2e-11);       // char permeability (second order tensor)
eps_g_v         eps_g_v         [0 0 0 0 0 0 0]         0.8;    // virgin material porosity
eps_g_c         eps_g_c         [0 0 0 0 0 0 0]         0.85;   // char porosity

/***            Anisotropic conductivity parameters: main directions and linear factors         ***/
//              kxyz = tP & kijk' & P
//              where, kijk'            =       (ki' 0  0
//                                               0  kj' 0
//                                               0  0  kk') (in direct orthonormal basis, i, j, k)
//
// 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


