Help Needed: EOS3 Non-convergence in Compressed Air Injection into Aquifer
Help! TOUGH2 EOS3 simulation of compressed air injection into aquifer fails after few steps. Adjusting parameters/time steps didn't work. Seeking root causes & debugging tips!
1 reply
-
Looking at the output, the run stops because the P exceeds the allowable limit of 1000 bar in eleme 41717. This means that probably the assigned injection rate of 20 kg/s is too high for the rock petrophysical properties. Either reduce the injection rate or reconsider the petrophysical properties.
Element "41717" belongs to rock domain UNDER with very low permeability (1.E-20 m2). Injected air should not enter into the UNDER rock domain, as AQUIF(er) permeability is much higher (1.E-12 m2).
I guess you should reconsider the option about the interpolation of permeability at the interface. Using MOP(11)=2 should be more appropriate than using 0 as you did.
I suggest to activate the printout of phase changes ( MOP(5) ) for further diagnostic purposes.
Regards,
Alfredo
MOP(11) = 0 *** CHOOSES EVALUATION OF MOBILITY AND ABSOLUTE PERMEABILITY AT INTERFACES.
= 0: MOBILITIES ARE UPSTREAM WEIGHTED WITH WUP. (DEFAULT IS WUP = 1.0). PERMEABILITY IS UPSTREAM WEIGHTED.
= 1: MOBILITIES ARE AVERAGED BETWEEN ADJACENT ELEMENTS. PERMEABILITY IS UPSTREAM WEIGHTED.
= 2: MOBILITIES ARE UPSTREAM WEIGHTED WITH WUP. (DEFAULT IS WUP = 1.0). PERMEABILITY IS HARMONIC WEIGHTED.
= 3: MOBILITIES ARE AVERAGED BETWEEN ADJACENT ELEMENTS. PERMEABILITY IS HARMONIC WEIGHTED.
= 4: MOBILITY * PERMEABILITY PRODUCT IS HARMONIC WEIGHTED.