Jump to content
Geochemist's Workbench Support Forum

Recommended Posts

Hi all,

I am running a reactive transport model where I apply  a vertical and horizontal fluid to a reactive transport model that contains a variable abundance of hydroxyapatite. However, when the fluid it applied, the whole domain (approx. 800 x 20 m) becomes undersaturated w.r.t. hydroxyapatite almost immediately. When applying a fluid that is undersaturated I would expect this to happen, but over a much slower timescale.

Attached is the X2t file. Any advice would be appreciated.



Sam Bingham


Link to comment
Share on other sites

Hi Sam,

The inlet fluids don’t infiltrate very far into your domain over the timescale of your simulation, let alone by the first step in the plot file. You can plot “pore volumes displaced” or add a conservative tracer to verify this. And looking at the Q/K for Hydroxyapatite, I don’t think it ever gets undersaturated in most of the domain, it just approaches equilibrium – a log Q/K of ~0. The Hydroxyapatite is supersaturated in the initial porewater, and I think it’s just precipitating really quickly, and thus approaching equilibrium very quickly. The undersaturated inlet fluid doesn’t really come into play here.

I think it would be a lot easier to troubleshoot a scenario like this with a 1D model in X1t, or even a flush model in React.

Hope this helps,

Brian Farrell
Aqueous Solutions   

Link to comment
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.

Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

  • Create New...