Some questions about simulation of dynamic displacement process

[dreamofzhouxin]

Dear Dr JE-santos:
Thanks for your great work helping me do my research a lot! Also I have met some problems that need your instruction:
(1) In file ShanChen.cpp, the value of con_theta is always positive due to abs, how can I simulate a process of wetting phase invading non-wetting phase? When I set Gads_f1_s1 > 0 such as 0.2, the results of rho seem to be divergent, however that for Gads_f1_s1= -0.2 is ok, why?


(2) From previous discussion, I know the values of rho1 and rho2 should be unchangeable for computational convergence, so how can I simulate the displacement process of two fluids with different density like oil and water? In other words, how can I set the parameters in file input.yml to compute the real flow process?
(3) In file ShanChen.cpp, the value of delta_rho that related to deltaP is dependent on Gads_f1_s1 and min_radius that is suggested to be 3, thus the deltaP is constant under the same wetting condition even for different fracture geometry. This may cause no occurrence of residual capture of wetting phase under drainage condition, which is inconsistent with facts. How can I simulate the process of residual capture?
(4) Under drainage condition, the wetting phase saturation decreases dramatically when Pc reaches a critical value, which leads to there is a straight line segment in kr curves. Can I save some rho.dat files at unconverted state for the critical Pc step to compute more values of Sw and kr? or can you give me other advises?

(5) In file create_geom_for_rel_perm.py in folder mplbm_utils, whether the computation of pore_vol has problem or not, the result of pore_vol seems to be the volume of solid portion? About tol and rho1_mask, why you choose these values, and can you give me some indications that how to determine the interface between the two phase in rho.dat file and calculate its areas?

Thanks for your help!

With kind regards,

zhouxin