Effect of contact angle and interfacial tension on relative permeability curves

[Matin-Bagheri]

Dear openpnm team
I am new to oopenpnm, and my question might be irrelevant, but since I didn't find any previous discussions on the matter, I'm asking here.

As I was working on a simple model to calculate relative permeability curves, it came to my notice, that using the InvasionPercolation algorithm with the example code provided on the website, the relative permeability results are independent of contact angle and IFT. Myu understanding is that the provided algorithm calculates the relative permeability based on the effective area of each phase at each saturation. Using the IP algorithm, the contact angle and IFT (if defined constant) have no influence on the pore invasion sequence and the phase occupancy. Thus, the relative permeability curves depend only on the pore structure, and not on the wettability. As long as theta > 90, the results are similar.

I tried to use Drainage algorithm for calculation of relative permeability, but so far I have only been able to shift the PC curve by changing theta and IFT. Even by applying the trapping, I didn't see any difference in the results. The only change was that half of the Kr curves disappeared. I would like to know if you have any suggestions for me to be able to capture the effect of contact angle and IFT on relative permeability curves.