New reaction sandbox example for mineral precipitation-dissolution

117 views
Skip to first unread message

Hammond, Glenn E

unread,
Jul 29, 2022, 6:03:03 PM7/29/22
to pflotra...@googlegroups.com, pflotran-dev (pflotran-dev@googlegroups.com)
PFLOTRAN users,

Due to the large number of PFLOTRAN users/developers who are developing custom mineral precipitation-dissolution rate expressions in a reaction sandbox, I have created an example sandbox for Calcite precipitation-dissolution. The example contains two rate expressions. The code for the first rate expression is very similar to what is implemented in reaction_mineral.F90:RKineticMineral() (without all the prefactors, Temkin constant, etc.). And the first leverages the stoichiometries and logKs from the database. The second has all the stoichiometries and logK hardwired. Both produce identical results, just slightly different formulation.

The following pull request tracks the change: https://bitbucket.org/pflotran/pflotran/pull-requests/635

Anyone is welcome to provide feedback.

Glenn

Glenn

unread,
Jul 29, 2022, 7:54:03 PM7/29/22
to pflotran-users
In response to a direct question from Peter:

> Glenn: when I try it the rate comes up zero in the .tec output file.

> …Peter

We have to go about printing the rates in a different way through the auxiliary_data array.  I pushed a commit that adds the output of the calcite rates through the Sandbox.  See the following

https://bitbucket.org/pflotran/pflotran/commits/01841b13380c4dca1f6326bad81d0215c0aa202d

Glenn

Peter Lichtner

unread,
Jul 29, 2022, 7:59:45 PM7/29/22
to pflotra...@googlegroups.com
Did you check—I still am getting zero for the calcite rate. 

...Peter <iPhone>

On Jul 29, 2022, at 17:54, Glenn <ghamm...@gmail.com> wrote:

In response to a direct question from Peter:
--
You received this message because you are subscribed to the Google Groups "pflotran-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pflotran-user...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/pflotran-users/677f3ad6-6096-45db-8b68-413d98c34a1cn%40googlegroups.com.

Hammond, Glenn E

unread,
Jul 29, 2022, 8:01:42 PM7/29/22
to pflotra...@googlegroups.com

Did you add AUXILIARY to the CHEMISTRY OUTPUT block?

 

https://bitbucket.org/pflotran/pflotran/commits/01841b13380c4dca1f6326bad81d0215c0aa202d#Lregression_tests/default/reaction_sandbox/reaction_sandbox_calcite.inT85

 

Glenn

 

From: pflotra...@googlegroups.com <pflotra...@googlegroups.com> On Behalf Of Peter Lichtner
Sent: Friday, July 29, 2022 5:00 PM
To: pflotra...@googlegroups.com
Subject: Re: [pflotran-users: 6776] Re: New reaction sandbox example for mineral precipitation-dissolution

 

Check twice before you click! This email originated from outside PNNL.

 

Hammond, Glenn E

unread,
Jul 29, 2022, 8:03:17 PM7/29/22
to pflotra...@googlegroups.com

The “Calcite Rate” will be zero, but the “Calcite Sandbox Rate” should be non-zero for most of the cells.

Reply all
Reply to author
Forward
0 new messages