Help Needed: EOS3 Non-convergence in Compressed Air Injection into Aquifer

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.