Keishi

Hello Jia, Thank you very much for your suggestion. We did the lab-scale experiment, and the rate laws greatly represent the experimental data. So I believed that the rate laws used in the React script have no problems. (We also think that the dissolution of Periclase and Dolomite should be treated as an equilibrium reaction, but the kinetic calculation showed better fitting). On the other hand, I also confirmed that the dissolution rates of Periclase and Dolomite dropped quickly just after the calculation was initiated. Then I tried to replace the rate laws of Periclase and Dolomite with equilibrium reactions (Please see the attached React script "dolomite dissolution_2"). When I did this, the calculation did not converge again. I made some small modifications to fluid chemistry such as units, but this problem cannot be solved. I also tried another calculation where a different database was used (thermo_minteq.tdat) to see the effect of the database, but it still did not work well. I attach its script herein as well (Please see "dolomite dissolution_minteq). I am really sorry for your trouble but could you check if there is any problem with my scripts again? (actually, I start thinking that I need to give up this series of calculations if no radical error is found this time...) Best wishes, Keishi Dolomite dissolution_2.rea Dolomite dissolution_minteq.rea

Hello Jia, Thank you very much for your reply. I understand that you are still investigating this issue. I confirmed that the initial system converged when I input the pH information as an H+ concentration. Thank you again for your kind cooperation. Regarding the custom rate laws, I also suspect that it causes the calculation error. Therefore, this was checked by making the other script in React and it was confirmed that the calculation was done without any error. Please see the attached script, where dolomite dissolution and the formation of the secondary mineral were taken into account. Best regards, Keishi Dolomite dissolution.rea

Hello Jia, It has been a while since the last reply. How is going to solve this issue? Since this calculation is crucially important for my research, I do appreciate it if you let me know the current status of this issue. Is this the problem with my input script or with the algorism in GWB? I am really sorry for your trouble, but I am waiting for hearing from you. Sincerely, Keishi

Hello Jia, Thank you very much for your reply. Herein, I attached the custom thermo and surface dataset as per your comment. Actually, I have already tested and confirmed that these thermo and surface datasets can work properly in "React" calculation. I hope this helps you to find some mistakes in my script. Best wishes, Keishi Brucite_Dolomite_Pb.sdat thermo.com.V8.R6+.+Pb.tdat

Hello, I'm trying to model the dissolution of half-burned dolomite in column experiment with "X1t". However, I can't get the model to converge the solution composition of "Fluids" . I tried same composition in "React" calculation and I found that it worked. I attached a copy of my input script. Can I get some advice to solve this issue? Thank you very much in advance. Keishi Half-burned dolomite column.x1t