Jump to content
RockWare Support Forum

Pierre Bachaud

  • Content count

  • Joined

  • Last visited

Community Reputation

0 Neutral

About Pierre Bachaud

  • Rank
  1. 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) 10 1.653e-6 30 3.527e-6 100 4.382e-6 1000 4.644e-6 10000 4.672e-6 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 ! Pierre Bachaud
  2. Diffusion with ToughReact

    Dear colleagues, I am working on CO2 sequestration and am currently trying to simulate CO2/brine/caprocks interactions with Petrasim/ToughReact. I have two quick questions for wich I wasn't able to find an answer elsewhere : - I am using the ECO2N EOS and wonder if the parameter DIFUN (diffusion coefficient for aqueous species) can be enter negative in order to get rid of the tortuosity dependency. I tried to do so and it made the simulation crash. - I guess that it is not possible, but is there any way to modify the global geometry after having defined boundaries ? Thanks in advance for your help. Pierre