Invasion percolation relative permeability code (2 phase): Different inlet boundaries for the non-wetting phase (oxygen).

[Alrnas]

Hi There,

I have been testing the relative permeability example for my problem, which is finding the relative permeabilities for water/oxygen in an electrolyzer electrode. I saw in that example a small hypothetical pressure difference was assumed (inlet 1 and outlet 0) for both air and water, and for each phase the entire side of the pore network is assigned to calculate the inlet_flow for that phase:
Snwp_num=10
flow_in = pn.pores('left')
flow_out = pn.pores('right')
max_seq = np.max([np.max(ip['pore.invasion_sequence']),
np.max(ip['throat.invasion_sequence'])])
start = max_seq//Snwp_num
stop = max_seq
step = max_seq//Snwp_num
Snwparr = []
relperm_nwp = []
relperm_wp = []
So in the code, how do I assign the inlet boundary condition for different oxygen content generated at the inlet at different cell operating condition to calculate the relative permeability for water/oxygen? Which sections of the code should be adjusted? Thank you