Jump to content
Geochemist's Workbench Support Forum
Jesse Tarnas

Solution not converging, Basis fluid chemistry is unconstrained

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

Share this post


Link to post
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now

×