“…This method can handle arbitrarily complex geometries of porous media with great ease, but it was reported to produce incorrect outcomes, especially when the contact angle is less than 908, because of unphysical mass transfer along the boundary [Leclaire et al, 2016]. Several strategies [Ba et al, 2013[Ba et al, , 2015Liu et al, 2015a;Liu and Zhang, 2015;Leclaire et al, 2016] have been proposed to decrease spurious velocities around the contact line and enhance the accuracy of modeling contact angle. However, they either suffer from the difficulty of dealing with curved solid boundaries [Ba et al, 2013[Ba et al, , 2015Liu et al, 2015a;Liu and Zhang, 2015] or may lead to nonsymmetric results (which arise from the least squares interpolation scheme) for a physically symmetric problem [Leclaire et al, 2016].…”