Jump to content
Geochemist's Workbench Support Forum

Solution not converging, Basis fluid chemistry is unconstrained


 Share

Recommended Posts

Hello all,

I am attempting to model alteration mineral formation from an ultramafic composition, but cannot get a solution that converges. I have tried changing the parameters in the Config -> Iteration and Config -> Stepping settings, but it still won't converge. The error I get is:

N-R didn't converge after 999 its., maximum residual =    2.46e-15, Xi = 0.0000
Largest residual(s):
                       Resid     Resid/Totmol   Cbasis
---------------------------------------------------------
 Fe++               -4.405e-20     2.46e-15    1.074e-12
---------------------------------------------------------
 

My epsilon value is also 1e-11, so not sure why these lesser residuals count as not converging.

This is for a planetary application, so the actual fluid chemistry is not constrained. Can someone please give me some advice on what changes should be made to the Basis fluid chemistry to make the solution converge? Or are there any "standard" Basis fluid chemistry packages that I could use to attempt to achieve solution convergence?

The current React model is attached.

Thank you in advance for your help and advice, it is much appreciated.

All the best,

Jesse

Chassigny_Freshstart_USE_AllMins.rea

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

×
×
  • Create New...