Jump to content
Geochemist's Workbench Support Forum

BCF Crystal Growth Theory


dirk
 Share

Recommended Posts

I am interested in writing a simple script that incorporates the BCF crystal growth theory equation but when I run React I get fluctuating precipitation/dissolution for the precipitating mineral. The result is fluctuating positive then negative mineral mass for siderite (with initial mass = 0). This occurs when the mineral mass is very low 1e-186 but this low number is the nature of the nonlinear BCF. When left to run, the mass keeps going from positive to negative.

Is there some way to write in the script something to prevent the negative numbers? It appears to be a complete transition from positive to negative to positive, ie. the next logical number in the series during the reaction can end up being the opposite sign.

 

Any help is much appreciated.

 

Attached is the script for siderite.

 

Dirk

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