Click here to donate to keep PhreeqcUsers open
Welcome,
Guest
. Please
login
or
register
.
Did you miss your
activation email
?
1 Hour
1 Day
1 Week
1 Month
Forever
Login with username, password and session length
Forum Home
Login
Register
PhreeqcUsers Discussion Forum
»
USER_PUNCH/USER_GRAPH
»
BASIC functions
»
How to extract Debye-Hückel parameters for all species from the database
« previous
next »
Print
Pages: [
1
]
Go Down
Author
Topic: How to extract Debye-Hückel parameters for all species from the database (Read 1963 times)
md.muniruzzaman
Frequent Contributor
Posts: 16
How to extract Debye-Hückel parameters for all species from the database
«
on:
March 10, 2021, 01:15:54 PM »
Hello everybody,
I am wondering if there is a way to extract Debye-Hückel (A, B) parameters for all the species in the solution. I know that I can always use the basic functions DH_A and DH_B to extarct these parameters for selected species, which is somewhat a manual approach and I need to do species by speices basis. So, is there any smarter way to extract these values automatically for all the species present in a SOLUTION?
Any suggestions or tips would be highly appreciated.
Cheers.
Logged
dlparkhurst
Top Contributor
Posts: 3667
Re: How to extract Debye-Hückel parameters for all species from the database
«
Reply #1 on:
March 10, 2021, 04:22:50 PM »
Unfortunately, no, there is no function to extract a0 (the ion size parameter) or bdot (the ionic strength factor) for individual species. The DH_A and DH_B are temperature dependent constants from the Debye Huckel equation that are the same for all aqueous species; a0 and bdot are ion specific parameters.
-DH_A*z^2*mu^.5/(1 + DH_B*a0*mu^0.5) + bdot*mu
Logged
md.muniruzzaman
Frequent Contributor
Posts: 16
Re: How to extract Debye-Hückel parameters for all species from the database
«
Reply #2 on:
March 10, 2021, 07:11:28 PM »
Hi David,
First of all, sorry for the mix-up between the parameters (smaller and capilal A and B )! But it seems you understood my point. Thanks for the clarification! So, it means if I want to extract these parameters (a0 or bdot), the only possible way would be to manually grab them from the database, right?
Thanks again!
Logged
dlparkhurst
Top Contributor
Posts: 3667
Re: How to extract Debye-Hückel parameters for all species from the database
«
Reply #3 on:
March 10, 2021, 07:20:44 PM »
Yes, I think you would have to do it manually.
Note that llnl.dat uses the same a0 and bdot parameters for all charged species, and one bdot for selected uncharged species. Also, Pitzer databases do not use a0 or bdot, nor does sit.dat.
Logged
dlparkhurst
Top Contributor
Posts: 3667
Re: How to extract Debye-Hückel parameters for all species from the database
«
Reply #4 on:
March 11, 2021, 03:08:36 PM »
I was writing other Basic functions, so I added DH_A0("H+") and DH_BDOT("H+"). No guarantees on when the next version is released, but we are working on it.
Logged
md.muniruzzaman
Frequent Contributor
Posts: 16
Re: How to extract Debye-Hückel parameters for all species from the database
«
Reply #5 on:
March 15, 2021, 02:03:35 PM »
Hi David,
Thanks a lot again for the clarification, and for considering these basic functions in the new release. It would surely be helpful for my purpose of implementating the activity gradient flux term in Nernst-Planck based diffusive flux calculations. In this direction, I wonder if it would also be possible to have "Get" methods to extarct these parameters in PhreeqcRM.
Regards,
Munir
Logged
dlparkhurst
Top Contributor
Posts: 3667
Re: How to extract Debye-Hückel parameters for all species from the database
«
Reply #6 on:
March 15, 2021, 02:48:26 PM »
There is a GetLog10Gammas method, and the next version will have a GetLog10Molalities. I'm kind of sick of writing PhreeqcRM methods.
If that is not enough, you could always use the new functions and SELECTED_OUTPUT/USER_PUNCH to extract the values in a one-time run to get ready for the main calculation. The list of species is available from GetSpeciesNames.
Logged
md.muniruzzaman
Frequent Contributor
Posts: 16
Re: How to extract Debye-Hückel parameters for all species from the database
«
Reply #7 on:
March 16, 2021, 10:47:16 AM »
Hi David,
No worries! It was only a suggestion for some future steps..... The selected_output method in connection with the new function would be perfectly sufficient :)
Logged
Print
Pages: [
1
]
Go Up
« previous
next »
PhreeqcUsers Discussion Forum
»
USER_PUNCH/USER_GRAPH
»
BASIC functions
»
How to extract Debye-Hückel parameters for all species from the database