gravity-capillary equilibrium in EOS3
I am using EOS3 (TOUGH3-eos3 windows executable) to simulate a gravity capillary equilibrium along a 1-D column but I keep running to convergence failure issue.
I was able to successfully simulate this problem in EOS7 but for some reasons that I cannot figure out, it’s not converging in eos3.
I was hoping if you could please take a look at my input and let me know what is wrong.
Attached are my input and output.
Thanks and best regards,
I took your input file without any modification and run it on my windows 10 pc (with the pc executable in the package we released), I got an error message:
$$$$$$$$$$ LIQUID PHASE DISAPPEARS AT ELEMENT *A11 1* $$$$$ SG = 0.105000E+02
!!!!!!!!!!! ERRONEOUS DATA INITIALIZATION !!!!!!!!!! STOP EXECUTION----------
+++++++++ CANNOT FIND PARAMETERS AT ELEMENT *A11 1 * XX(M) = 0.000000E+00
Then I notice you set MOP(19)=1, which used a different set of initialization for primary variables. I re-set MOP(19)=0. It runs all the way to the end without any problem.
My question to you is: when you run it on pc, did you use the executable in the package you purchased or did you make one yourself? if you made one, what compiler did you use?
I attach the input and output file from my run.
Thank you, Yingqi