Page:4SIGHT manual- a computer program for modelling degradation of underground low level waste concrete vaults (IA 4sightmanualcomp5612snyd).pdf/113

39. INITIALIZE PARAMETERS B.16

real = 0.16; /* porosity */
real F; /* formation factor */
real $\alpha$; /* degree of hydration */
real f; /* initial solids fraction */
real Vsample; /* volume of element */
real $\rho _{CH}$ = 33.1;     /* density of CH */
real molesCH = 13.0;    /* moles CH/ltr soln. */
real $\vartheta _{0}$= 0.00281;    /* initial 1/F */
real $\vartheta _{0}$    /* current 1/F */
real print day; /* next day to print interm. results */
real print day_interval; /* days between interm. results */
real earliest failure_day; /* earliest structural failure */

41. initialize.parameters(). This routine initializes the scaling parameters and physical parameters.

void initialize parameters()
{
 int i, j, k;
 int n, m;
 int iter;
 real moles_cation, moles_anion, charge, z, dz, z, critical_Cl_concentration, T_crit,
     concen_ratio, darcy_velocity, perm_time, chloride_time, perm_factor,
     /* modify DT to account for high permeability */
  perm_depth, A, crack_factor;    /* modify DT to account for cracks */
  for (k = FIRST_CELL; & < NUM_SURFACES; k++)
    $\Delta {X}$[k] = 1.0;
  ⟨Establish material parameters of concrete 42⟩
  aci_211();
  ⟨Propagate ion concentration in concrete 45⟩
  ⟨Print parameters to stdout 47⟩
  ⟨Initialize ion 'eta' parameter 48⟩
  L=sample_length/NUM_CELLS; /* the universal length scale */
  MaxDepth MarDepth/L; /* dimensionless */
  ⟨Calculate crack and joint adjustment to permeability 49⟩
  critical_Cl_concentration 0.000400* Vsample (1.0) * 2.5/35.4;
  if(anion [Cl].c[FIRST_CELL] > critical_Cl_concentration) {
       /* calculate T for [Cl_crit]/[Cl(x=0)] = erfc(X/2 sqrt(DT)) */
    concen_ratio= inv_erfc (critical_Cl_concentration/anion [Cl].c[0]);
    T.crit SQR((rebar_depth/L)/(2.0* concen_ratio));
    chloride_failure_year = $3.171\cdot 10^{-08}$ * T_crit * SQR(L)/</math>D_0</math>;
    if(head.value> 0.0) {    /* correct diffusion estimate to include darcy penetration */
      darcy velocity = ($k_{0}$/$\mu$) + ($\rho$*g*head_value)/sample_length;
      A = 1.0/(4.0 <math>D_0<?math> SQR (concen_ratio));
      perm_depth = rebar_depth +1.0/(2.0* A* darcy-velocity) - sqrt(rebar_depth/(A *
          darcy-velocity) + 1.0/(4.0 SQR(A darcy_velocity)));
      perm_time = perm_depth/darcy_velocity;
      chloride.failure-year = 3.17·10E-08 * perm_time;
    }
  }
  else chloride_failure-year=<math>1.0\cdot10^{+10}<>math>;