Assistance with Capillary Pressure in Black Oil Simulation

[IGWEONU CHIBUZOR]

Good day,

I am currently working on a blackoil simulation which is 100% oil saturated, and I am looking to inject water into the reservoir. I would like to add capillary pressure to my simulation, and here’s the snippet of the code I’m working with:

"

pRef = 5000*psia;
s0 = [0, 1, 0];

fluid = initSimpleADIFluid('phases' , 'WOG', ...
    'mu' , [   1,  5, 0.2] * centi*poise     , ...
    'rho', [1000, 700, 250] * kilogram/meter^3, ...
    'c',   [1e-8, 1e-5, 1e-3] / barsa, ...
    'n'  , [2, 2, 2], ...
    'pRef' , pRef);

srw = 0.1; %residual water saturation          

 

pe = 5 * kilo * Pascal;
pcOW = @(sw) pe * sw.^(-1/2);
fluid.pcOW = @(sw) pcOW(max((sw - srw) / (1 - srw), 1e-5));"

The issue I am encountering is that my simulation is not converging. I’ve searched through the MRST folders for any implementation of capillary pressure using this approach, but I couldn’t find any relevant examples.

Could you please help me with this issue or point me to any examples that could guide me?

Thank you very much.

Best regards,

Chibuzor