Leaderboard

Hello Zhidi, Thank you for attaching your script. To troubleshoot, I checked "Follow Output" in the Results pane and "explain on" in the Stepping dialog before running your model to get a better idea what is limiting the size of the time step in the simulation. Doing so, I noticed that the rate limiting constraint right before the model runs into the newton-raphson iteration error is due to the change in Q/K for reactant Illite. To test, I simplified your model by setting Illite to a reaction rate of 0 and was able to run your simulation to completion. I also noticed that the reaction rates for various minerals do vary quite a lot be Could you provide a bit of context regarding your fluids and simulation? I noticed that a lot of component concentrations are the same but with a few differences. Are these fluids you have created or measured? I noticed that you have added all minerals as a kinetic reaction except Halite as a simple mineral. In general, kinetic reactions are useful for reactions that are proceeding slowly but measurable over the time period of interest. I noticed that Quartz is dissolving so slowly that it's not really doing anything. If the mineral barely reacts then perhaps you do not need to consider it in your model. This can help to simplify things in terms of the kinetic reactions running through your system. We don't typically see simple minerals, which are incrementally added to the domain over the course of the simulation in a reactive transport simulation. If a mineral is in equilibrium with the initial fluid in the system, you should consider swapping it in for a component in your Initial pane. In general, it is good to check that the log Ks in the thermo database are suitable for the type of reactions occurring in your system. Hope this helps, Jia Wang Aqueous Solutions LLC