search

je-santos / MPLBM-UT

Library for performing multiphase simulations (based on the Shan-Chen model) in complicated geometries (i.e. porous media 3D images)
GNU General Public License v3.0
150 stars 57 forks source link

why first 2 phase and than 1 phase #36

Closed shan-zhong closed 3 years ago

shan-zhong commented 3 years ago

In the example of steady-state flow, why should we first simulate the 2-phase flow and then simulate the 1-phase flow? If I simulate a 2-phase flow with different saturation, do I only need to simulate the 1-phase flow ones? And for the result of the example, the krn curve starts from 1 at Sw=0 and end at 0 with Sw=1. At every point, the sum of krn and krw is equal to 1, is it true for reality?

Larry-Kong commented 3 years ago

I have the same question...

je-santos commented 3 years ago

Dear @shan-zhong and @Larry-Kong

Thank you for your interest in the simulator. I apologize for the late reply. I have updated the readme to make it a little more clear. Let me know if this is sufficient:

Run:

run2-phase.sh

to run the two-phase simulation. Once this simulation converges, it will output the position of both fluids (wetting and non-wetting) in the pore-space. We now would like to calculate the permeability of these fluid configurations to obtain one point of the relative permeability curve. This method assumes very low capillary numbers. To create these new domains we run:

create_geoms_4_kr.m

which builds the geometries to calculate the relative permeability points at the given saturation (in this case, 25%). These simulations are then carried out with: