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\cdot10^{-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 */ if (sulfate_rate > 0.0) sulfate failure.year = $$3.171\cdot10^{-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$$\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 */