Hi Renchao and everyone,
I was working on the same project I have posted in this topic and faced another problem during solving the model. I don’t want to create another topic in the discourse so I wanted to ask the question in here.
I am trying to model the bacteria transport within in a sand column (1-D). There are two mechanisms during the transport, bacteria decay in time and attachment of the bacteria on sand grains.
I’ve got two components of the bacteria, suspended (unattached) and attached.
change in attached bacteria wrt to time is => kattSuspended - kdSuspended in which katt is attachment rate (per second) and kd is the decay rate (per second).
change in suspended bacteria wrt to time is => -kattSuspended - kdAttached
Initially, both suspended and attached concentrations are zero.
When I try to run the simulation (attached below), simulation stuck at “Phreeqc: Executing chemical calculation.”. I waited for a long time but couldn’t pass this part.
However, if I remove the kd*Attached component of the rate, simulation can be completed, but it is wrong of course.
I am working on it but couldn’t find a solution. Any kind of help would be greatly appreciated.
Only_Bacteria_v2.zip (55.2 KB)