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

39. INITIALIZE PARAMETERS B.16 real = 0.16; /* porosity */ real F; /* formation factor */ real $$\alpha$$; /* degree of hydration */ real f; /* initial solids fraction */ real Vsample; /* volume of element */ real $$\rho_{CH}$$ = 33.1;    /* density of CH */ real molesCH = 13.0;   /* moles CH/ltr soln. */ real $$\vartheta_0$$= 0.00281;   /* initial 1/F */ real $$\vartheta_0$$   /* current 1/F */ real print day; /* next day to print interm. results */ real print day_interval; /* days between interm. results */ real earliest failure_day; /* earliest structural failure */ 41. initialize.parameters. This routine initializes the scaling parameters and physical parameters.

void initialize parameters { int i, j, k;  int n, m;  int iter; real moles_cation, moles_anion, charge, z, dz, z, critical_Cl_concentration, T_crit, concen_ratio, darcy_velocity, perm_time, chloride_time, perm_factor, /* modify DT to account for high permeability */ perm_depth, A, crack_factor;   /* modify DT to account for cracks */ for (k = FIRST_CELL; & < NUM_SURFACES; k++) $$\Delta{X}$$[k] = 1.0; ⟨Establish material parameters of concrete 42⟩ aci_211; ⟨Propagate ion concentration in concrete 45⟩ ⟨Print parameters to stdout 47⟩ ⟨Initialize ion 'eta' parameter 48⟩ L=sample_length/NUM_CELLS; /* the universal length scale */ MaxDepth MarDepth/L; /* dimensionless */ ⟨Calculate crack and joint adjustment to permeability 49⟩ critical_Cl_concentration 0.000400* Vsample (1.0) * 2.5/35.4; if(anion [Cl].c[FIRST_CELL] > critical_Cl_concentration) { /* calculate T for [Cl_crit]/[Cl(x=0)] = erfc(X/2 sqrt(DT)) */ concen_ratio= inv_erfc (critical_Cl_concentration/anion [Cl].c[0]); T.crit SQR((rebar_depth/L)/(2.0* concen_ratio)); chloride_failure_year = $$3.171 \cdot 10^{-08}$$ * T_crit * SQR(L)/ D_0 ; if(head.value> 0.0) {   /* correct diffusion estimate to include darcy penetration */ darcy velocity = ($$k_0$$/$$\mu$$) + ($$\rho$$*g*head_value)/sample_length; A = 1.0/(4.0 D_0<?math> SQR (concen_ratio)); perm_depth = rebar_depth +1.0/(2.0* A* darcy-velocity) - sqrt(rebar_depth/(A * darcy-velocity) + 1.0/(4.0 SQR(A darcy_velocity))); perm_time = perm_depth/darcy_velocity; chloride.failure-year = 3.17·10E-08 * perm_time; }  }   else chloride_failure-year=1.0\cdot10^{+10}<>math>;