SO2 Kinetics problem

[justinninja1]

Hello,

I'm designing a model based on the kinetics of SO2 into water with various metal ions acting as a catalyst.

Since I am new to PhreeqC, I am starting off very basically and am just focusing on SO2 gas in pure water solution.

My aim is to plot the formation of Sulphate ions, SO42-, as a function of time. To do this, I understand I need to use the rates and kinetics keywords, and I have designed a statement which allows my SO2 to convert to SO3, but the problem is my kinetics stop the reaction there:

 RATES
Suox
 -start
 
   10 k = 1.4e5
   20 Sx = TOT("SuO3--")
   30 IF (Sx <= 0)then goto 200
   #30 k = 1.4e5
   40 rate = k*Sx
   100 moles = rate*TIME
   200 SAVE moles
 -end

 KINETICS
 Suox
        -formula Su      -1.0   O   -1.0    Sf +1.0
        -steps 100 400 #3100 10800 21600 5.04e4 8.64e4 1.728e5 1.728e5 1.728e5 1.728e5
        -step_divide 1e-4





Su is tetra-valent Sulphur (S+4) and Sf is hexa-valent Sulphur (S+6). I've removed the reactions which link S+4 with S+6 and am attempting to use the above statement to form my sulphate but the problem is that no SO4 will form:

Kinetics 1.   

   Time step: 100 seconds  (Incremented time: 100 seconds)

   Rate name       Delta Moles Total Moles   Reactant        Coefficient

   Suox             -0.000e+00   1.000e+00   O                        -1
                                             Sf                        1

This continues no matter what the time is, the SO3 ions are simply not being used up and I have no idea why this is no the case. I can see from the speciation at equilibrium that that have been formed with the reaction with water.

Any help would be appreciated, thank you.