Dear Petrasim/ToughReact users,
I'm currently having some difficulties having consistent results during batch reactions calculations. It seems that the final concentrations are depending on the defined max time step (DELTMX), which makes totally no sense for me.
For example, I did some tests with a 1 cubic meter-cell, with a porosity of 20% and a mineralogic composition of 40% of calcite, 20% of ankerite, 10% of illite, 10% of kaolinite, 10% of pyrite, and 10% of quartz. Calcite is set to equilibrium, while kinetic constraints are defined for the other minerals.
Here are the SiO2(aq) final concentrations obtained after a 10 days simulation for different max time steps :
Max time step(s) SiO2(aq)(mol/kgw)
Of course, except for DELTMX, no parameter is changed during these tests.
Are these discrepancies expected? Is this a solver problem? I'd really appreciate some help here, because I really haven't any idea where to look at...
Thanks in advance to you all !