Is there a possibility to extent the TOUGH2 pressure range?
Hello everyone,
I am using TOUGH2 to simulate two phase flow in deep subsurface waste repositories. The pressures in the systems can reach quite high values due to gas production and high hydrostatic pressures.
Tough2 quits my simulations (convergence failure) with the messages that:
P= ... IS OUTSIDE THE RANGE (-0.10000E+06, 0.20000E+08)
I would however prefer to continue the simulation, since that pressue is not unrealistic. Is there a possibility to change those pressure boundaries?
If nexessary I could change the source code, but I havn't found where those boundaries are defined.
Thanks everyone in advance for their help.
Best regards, André.
3 replies
-
by the way...
I forgot to mention that I am using module EOS5
-
Hi André,
You would have to extend the pressure range in subroutine VISCO (file eos5.f) using a suitable function. I am not sure, however, whether this will allow you to continue the simulation, as you may run into another pressure limit somewhere in the code. I do not know the details of your application, but consider using a pore compressibility and potentially introduce permeability increases due to pathway dilation due to pressure buildup, both approximations of relevant processes that may reduce the maximum pressure in your system.
Stefan
-
Thank you for the information, Stefan. I will have a look into the VISCO subroutine.