Dear. all,

I performed a direct numerical simulation of multiphase flow in a porous medium (considering surface tension). And I want to create a similar demonstration of the results like in this video on youtube: https://www.youtube.com/watch?v=onmS612m9-E&ab_channel=FSEelearning

Here, in the video, the relative permeability depends on saturation. How to determine the relative permeability?

In Wikipedia(https://en.wikipedia.org/wiki/Relative_permeability), the relative permeability for phase 'i' is

k_{ri} = - mu_i q_i /(k*grad p_i)

where where q_i is the flux, p_i is the pressure drop, mu_i is the viscosity, k is the absolute permeability.

I calculate the viscosity, pressure gradient, as $f_i = \int_{\Omega_i}f_i \varphi_r dV/\int_{\Omega_i}\varphi_i dV$, f = (\mu, \grad p)

The flux $q_i$ is calculated on the outlet.

Also, it is not clear which component is assumed for the absolute permeability calculation?

BR,

Evgenii