In this work, a simplified wetting boundary scheme in the phase-field lattice Boltzmann model is developed for wetting phenomena on curved boundaries. The proposed scheme combines the advantages of the fluid-solid interaction scheme and geometric scheme-easy to implement (no need to interpolate the values of parameters exactly on solid boundaries and find proper characteristic vectors), the value of contact angle can be directly prescribed, and no unphysical spurious mass layer-and avoids mass leakage. Different from previous works, the values of the order parameter gradient on fluid boundary nodes are directly determined according to the geometric formulation rather than indirectly regulated through the order parameters on ghost solid nodes (i.e., ghost contact-line region). For this purpose, two numerical approaches to evaluate the order parameter gradient on fluid boundary nodes are utilized, one with the prevalent isotropic central scheme and the other with a local gradient scheme that utilizes the distribution functions. The simplified wetting boundary schemes with both numerical approaches are validated and compared through several numerical simulations. The results demonstrate that the proposed model has good ability and satisfactory accuracy to simulate wetting phenomena on curved boundaries under large density ratios.
Read full abstract