nonequilibrium transport model in rate module

The HP1 code incorporates modules simulating
(1) transient water flow in variably-saturated media,
(2) transport of multiple components, and
(3) mixed equilibrium/kinetic geochemical reactions.
Post Reply
Bowen Luo
Posts: 2
Joined: Fri Apr 28, 2023 6:02 am
Location: China

nonequilibrium transport model in rate module

Post by Bowen Luo » Thu Jul 06, 2023 2:06 pm

hi,
I am using the phreeqc module of HYDRUS-1D. I want to write two formulas of nonequilibrium transport model (two-site model) in rate module. I've written a piece of code that works, but I'm not sure it does, exactly as the formulas say. So I am hoping you could help me correct it! Here are two formula which are shown as jpg and code.

Code: Select all

Rates
Sea
-start
20 beta = parm(1)
30 rate  = TOT("Cd")^ beta
40 moles = rate * time
50 PUT(rate, 10)
200 SAVE moles
-end

Seb
############################
-start
20 rho = parm(1)
25 f = parm(2)
26 kd = parm(3)
30 rate1 = GET(10)
40 rate2  = f * rho * kd * rate1
50 moles = - rate2 * time
200 SAVE moles
-end

Sk
############################
-start
20 rho = parm(1)
25 f = parm(2)
26 kd = parm(3)
27 beta = parm(4)
28 w = parm(5)
30 rate = rho * w * (1 - f) * kd * TOT("Cd") ^ beta - rho * w * kin("Sk")
50 moles = - rate * time
200 SAVE moles
-end
Attachments
2.jpg
2.jpg (8.74 KiB) Viewed 69812 times
1.jpg
1.jpg (5.62 KiB) Viewed 69812 times

Post Reply