Jump to content
Geochemist's Workbench Support Forum

X1t convergence problem


dirk

Recommended Posts

Sorry posted this on the archive thread...

 

Hi, I am running a simulation of a flow through column experiment and slowly working up to the full rtm. I run the initial solution through for 5 days then the contaminated solution. calcite reacts based on kinetic constraints, precipitation of gypsum, gibbsite, Fe(OH)3 at equilibrium all else suppressed, adsorption using the D&M model and the FeOH database (decreased site densities of w and s sites by an order of magnitude). after adding adsorption I cannot get it to finish. current error is "Step is constrained by rate of increase in length of reaction step at node 0." after 8 days of 10 day run. Any ideas?

 

Dirk

Link to comment
Share on other sites

Ok, sorted some out. Turns out the contaminated fluid was undersaturated with respect to Fe(OH)3 so by day 8 it has exhausted the Fe(OH)3 and that made convergence not possible for some reason that seems to be related to increase in the length of the reaction step in the first cell. So using goethite instead, seems to be working.... So I am wondering about the why of it. Fe(OH)3 is the adsorbant and is the breakdown because there is no adsorbant left , other minerals (calcite for instance) do go to 0 in the first cell, so does that mean that there has to be some adsorbant left in a cell in order for simulation to go forward? 

Link to comment
Share on other sites

Hello Dirk,

I am glad you got your model up and running. It sounds like your convergence issues might be linked to your model running out of minerals. When you run out of mineral(s) or some other reactant in your system, the model will take smaller time steps to try to calculate the numerical solution, which might explain why you received the error “Step is constrained by rate of increase in length of reaction step at node 0”. In this case, you were trying to sorb solutes onto Fe(OH)3 when it is running out, which is probably why your model couldn’t converge. However, I can’t be sure unless I take a look at your input file but that would be my guess. If you would like someone to help take a closer look, please attach your input and thermo file.

Hope this helps,

Jia Wang

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.

Guest
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.

Loading...
×
×
  • Create New...