Click here to donate to keep PhreeqcUsers open

KINETICS 2-9Fe(cr,alpha) -formula Fe 1 -m 0.7968 -m0 0.7968 -parms 2.45e-011 4.9e-011 6.5e-011 -tol 1e-008-steps 1 in 0 steps # seconds-step_divide 1-runge_kutta 3-bad_step_max 500-cvode true -cvode_steps 100-cvode_order 5INCREMENTAL_REACTIONS TrueRATES Fe(cr,alpha)-start 1 REM a0 is surface area (cm2), simple surface of cylindrical rebar 6 mm diameter 200 mm long 2 REM v is solution volume (cm3) (pore solution -> solution 1-10) 10 a0 = 37.68 20 v = 78.15 30 sr_iron = SR("Fe(cr,alpha)") 40 IF(M <= 0) THEN GOTO 90 50 area = (M/M0)^(2/3)*a0 60 IF (-LA("H+") >= 14) THEN rate_const = PARM(1) 70 IF (-LA("H+") < 14) AND (-LA("H+") > 12) THEN rate_const = PARM(2) 80 IF (-LA("H+") < 12) THEN rate_const = PARM(3) 90 rate = (area/v)*rate_const*(1-sr_iron)100 moles = rate*TIME110 SAVE moles-end