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

68. Adjust physical parameters.

Due to dissolution/precipitation.

⟨Adjust physical parameters 68⟩ =
for ( k = FIRST.CELL; k < NUM_CELLS; k++)
$\phi ^{n}$[k]=0.5*($\phi '^{n}$[k] +$\phi ^{n}$[k+1]);    /*interpolate phi' values */
$\vartheta$=DLIMIT+0.007*SQR($\phi '^{n}$[k])
if ($\phi ^{n}$[k]>0.180) $\vartheta$ += 1.8 * SQR ($\phi ^{n}$[k]-.180);
if ($\phi ^{n}$[k]>$\phi$)    /* leaching */
  $\xi$=$\vartheta _{0}$+5.0*($\vartheta$-$\vartheta _{0}$);
else    /* precipitation */
  $\xi$[k]=$\vartheta$

}
$\xi '$[FIRST_CELL]-($\Delta X$[FIRST_CELL]/(1.0 + $\Delta X$[FIRST_CELL])) * ($\xi$[FIRST_CELL + 1]-$\xi$[FIRST_CELL]
for ( k = FIRST.CELL +1; k < NUM_SURFACES - 1; k++) $\xi$[k] = 0.5 * ($\xi$[k]+$\xi$[k-1]);
$\xi '$[NUM_SURFACES - 1] = 1.5 * $\xi$[NUM_CELLS - 1] - 0.5 * $\xi$[NUM_CELLS - 2];

This code is used in section 9.