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

39. INITIALIZE PARAMETERS B.17

printf("\"Chloride failure (yr)\"\t\t %8.01f\n", chloride_failure_year);
perm_factor = MAX($k$_0/$3.0\cdot 10^{-18}$, 1.0) * MAX (head.value/10., 1.0);
if (crack_width.value == 0.0) {
   crack-factor = 1.0;
}
else {
  if (sample_length = crack-depth > 0.05 * sample_length)
     crack-factor = sample.length/(sample_length crack_depth);
  else crack_factor = 10.0;
}
$\Delta$T = 1.0/(64.0 * perm_factor * crack_factor * anion [Cl].n * do);
sulfate_rate_ms = Youngs.value SQR($\beta$.value) * anion [S04].c[FIRST_CELL] * CE.value *
     $D_{0}$/(roughness.value * $\gamma$ .value * (1 - ν.value)); /* m/sec */
sulfate_rate= sulfate_rate_ms (L/$D_{0}$); /* $\Delta X$/math>\Delta T</math> */
if (sulfate_rate > 0.0) sulfate failure.year = $3.171\cdot 10^{-08}$ * (MaxDepth L)/sulfate_rate_ms;
else sulfate failure_year= 1.0-10E10;
printf("\"Sulfate failure (yr)\"\t\t%8.0lf\n", sulfate_failure_year);
     /* calculate number of days until termination */
earliest failure_day = 365. MIN (sulfate failure_year, chloride_failure_year);
    /* determine intervals between print outs */
print_day_interval = MIN( earliest_failure_day, MaxDay.value)/PRINT_OUTS;
print.day
print_day_interval;
Pout ($k_{0}$/($\mu$* $D_{0}$)) * $\rho$ * g * head.value;
dP = ($k_{0}$/($\mu$* $D_{0}$)) * $\rho$ * g * (sample_length/NUM_CELLS); /* change per cell */
for (k = 0; k< NUM_SURFACES; k++) {
  $\xi '$[k] = 6;
  $\phi '^{n}<math>\phi$-1[k] = $\phi$;
  $\Psi$[k] = Pout + (real) (NUM_SURFACES — k − 1)/(NUM_SURFACES — 1);
}
(Initial estimate of pH 50)
}

42. The material parameters of the concrete must be determined based upon the information given by the user. Missing pieces of information must be calculated from established relationships. The . is_default portion of the material data determines whether the user specified the value.

(Establish material parameters of concrete 42) =

if ($D_{\infty }$.is_default != TRUE) { /* calculate porosity */
if (wc.is_default == TRUE) wc.value = (log10 (10000. $D_{\infty }$.value) +9.84)/6.0;
if ($k_{\infty }$.is_default == TRUE) 
      /* kinfty.value = 1.0E-18*(0.8904+.002525*exp(15.07* wc.value)); */
   $k_{\infty }$value pow(10.0, (5.0* wc.value - 21.0));
}
if ($ka_{\infty }$.is_default != TRUE) {
   if ($D_{\infty }$.is_default = TRUE) {
      $D_{\infty }$.value = 1.0-10E04 * pow(10.0, 6.0 * wc.value - 9.84);
      if (wc.is_default = TRUE) wc.value = (log10 (10000. $D_{\infty }$.value) +9.84)/6.0;
   }
}
if (wc.is_default != TRUE) {   /* calculate diffusivity and porosity */