Jump to content
Geochemist's Workbench Support Forum

Fluid evaporation using "React"


Penguin623

Recommended Posts

Hi Brian,

 

I have modified "thermo_phrqpitz.dat" based on "thermo.dat", since I need Aluminum and Silicon to be added in the initial condition. In my calculation, I considered fluid evaporation and used "Flow-through configuration". However, after loading the modified database, the calculation gave me "Initial solution is too supersaturated". There is no further calculation. But, when I use "thermo.dat" to do the calculation, I could get a plot. Could you please let me know whether I can modify "thermo_phrqpitz.dat" and how? Or should I just use "thermo.dat" for fluid evaporation reaction path models if there is more element in my system?

 

I have attached the script and the modified database with this post. Please let me know your opinion.

 

Thank you very much!

thermo_phrqpitz_WXM.dat

Minerals-thermo.dat.rea

Link to comment
Share on other sites

Hi,

 

For this example, I would stick with the thermo.dat dataset. Your initial solution does not have a very high ionic strength and you are only evaporating about a third of the water. Running this example, the highest ionic strength reached is about 0.05 molal, so the B-dot equation (extended form of Debye-Huckel) works just fine. Thermo.dat also has the advantage of being parameterized from 0-300 C, unlike the phrqpitz dataset, which is defined at 25 C, though it has a limited (and largely untested) provision for calculations at other temperatures.

 

Unfortunately you cannot just add thermodynamic data for Al and Si to the phrqpitz dataset (which uses the virial activity model), without adding appropriate virial coefficients to the bottom section of the dataset. And if you are editing a database, be sure to follow the conventions outlined in the Appendix to the GWB Reference Manual. (Elements added in alphabetical order by name (Oxygen is always first, however), the list of Basis species must begin with H2O, etc.)

 

Hope this helps,

 

Brian Farrell

Aqueous Solutions

Link to comment
Share on other sites

Hi Brian,

 

Thank you very much for your explanation, What iss the highest ionic strength the B-dot equation (extended form of Debye-Huckel) can accomodate? If the total water loss for my system is 50%, can I still use "thermo.dat"?

 

Thank you very much!

Link to comment
Share on other sites

From Craig Bethke's Geochemical and Biogeochemical Reaction Modeling text, Section 8.1.2:

 

"The B-dot equation...is considered reasonably accurate in predicting the activities of Na+ and Cl- ions to concentrations as large as several molal, and of other species to ionic strenghts up to about 0.3 to 1 molal."

 

You should take a look at Chapter 8 of that text, or Section 7.4 of the GWB Essentials Guide for more on the implementation of activity coefficients and activity models in GWB. As for your specific problem, I would use thermo.dat.

 

Hope this helps,

Brian

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