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
153 stars 59 forks source link

Some questions about simulation of dynamic displacement process #87

Closed dreamofzhouxin closed 1 year ago

dreamofzhouxin commented 1 year ago

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

xiangWu-WW commented 4 months ago

Hi，Zhouxin，I have also been doing some wetting work to displace non-wetting phase recently, which requires setting different density viscosity ratios. I wonder if your problem has been solved yet. Looking forward to your reply and suggestions.