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

if ($D_{\infty }$.is_default = TRUE) $D_{\infty }$.value = $1.0\cdot 10^{-04}$ * pow(10.0, 6.0 * wc.value -9.84);
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));
}
$\vartheta _{0}$ = $D_{\infty }$.value/anion [Cl]. $D_{f}$;
if $\vartheta _{0}$ < DLIMIT )  $\vartheta _{0}$ = DLIMIT +0.0001;
if $\vartheta _{0}$ < DLIMIT + 0.07 * SQR(.18)) $\phi$ = sqrt(( $\vartheta _{0}$ - DLIMIT)/.07);
else $\phi$ = .17326 + sqrt(.03002 - (.05832 + DLIMIT - $\vartheta _{0}$)/1.87);
$\alpha$= (3.2 * wc.value)/1.36 - $\phi$ *(1 +3.2* wc.value)/1.36;
ƒ = 1.0/(1+ 3.2 * wc.value);
Vsample = 1000.0/$\phi$; 
molesCH = $\alpha$ + (0.61) * f * Vsample/$\rho _{CH}$;
F = 1.0/$\vartheta _{0}$;
$D_{0}$=$D_{\infty }$.value;
$k_{0}$=$k_{\infty }$.value;
litre [0] = 0.001* Vsample * $\phi$;
for (k = 1; k< NUM_SURFACES; k++) {
  sol_array[Ca] [OH].s[k] = molesCH;
  litre [k] = 0.001* Vsample * $\phi$;
}
for (k = 0; k< NUM_CELLS; k++) {
  $\xi$[k] = $\vartheta _{0}$
  $\phi ^{n}$[k] = 6;
}

This code is used in section 41. 43. The routine aci-211 estimates the strength, in psi, of the concrete in the individual com- putational elements. The estimate is based upon Table 5.3.4(a) of ACI 211.1-81 for air-entrained concrete. The strength can be adjusted to account for changes in porosity using: ${\textit {psi}}'\left[i\right]={\textit {psi}}_{0}{\frac {\ln {\left(\phi ^{n}\left[i\right]\right)}}{\ln {\left({\textit {phi}}_{0}\right)}}},$

(Function declarations 23) +=

void aci-211 (void);

44. void aci211()

{
  int i;
  real $\psi$;
  $\psi$ = pow(exp(wc.value - 3.575), -2.6817);
  for (i = 0; i < NUM_SURFACES; i++) strength [i] = $\psi$;
}

45. The function yyparse sets the concentration in the concrete by putting the user-specified value into the c[1] cell. The internal concentration is established by copying this value into the remaining cells. However, this procedure is not performed for OH and H.