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

68. ADJUST PHYSICAL PARAMETERS 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.