TY - JOUR
T1 - Mesh-free simulation of two-phase fluid flow in porous media based on the shock-fitting method
AU - Lee, En Jui
AU - Wang, Wei
AU - Chen, Po
AU - Jiao, Zunsheng
AU - Gong, Yanbin
AU - Mu, Dawei
AU - Liao, Wu Yu
N1 - Publisher Copyright:
© 2022 Elsevier B.V.
PY - 2022/8
Y1 - 2022/8
N2 - The physical interactions between two immiscible fluids in porous media may lead to the formation of a shock front at the interface between the two fluid phases. Numerical methods for modeling fluid flow with shocks can be classified into two categories: the shock-capturing and the shock-fitting methods. The shock-fitting method can in general provide more accurate solution of the shock front than the shock-capturing method, but suffers from the extra complexity in accommodating the moving shock front inside the background mesh. In this study, we explore the possibility of integrating the mesh-free methods for solving partial-differential equations with the shock-fitting method. In this integrated method, which we call “mesh-free shock-fitting (MFSF)”, the nodes needed for mesh-free calculations of spatial derivatives are generated on the fly to adapt to the moving shock front. We demonstrate the implementation of MFSF using numerical experiments in one and two spatial dimensions. Preliminary results show that MFSF has the potential to provide more accurate solutions with lower computational cost than conventional shock-capturing methods and can simplify many of the operations in shock-fitting methods with unstructured meshes.
AB - The physical interactions between two immiscible fluids in porous media may lead to the formation of a shock front at the interface between the two fluid phases. Numerical methods for modeling fluid flow with shocks can be classified into two categories: the shock-capturing and the shock-fitting methods. The shock-fitting method can in general provide more accurate solution of the shock front than the shock-capturing method, but suffers from the extra complexity in accommodating the moving shock front inside the background mesh. In this study, we explore the possibility of integrating the mesh-free methods for solving partial-differential equations with the shock-fitting method. In this integrated method, which we call “mesh-free shock-fitting (MFSF)”, the nodes needed for mesh-free calculations of spatial derivatives are generated on the fly to adapt to the moving shock front. We demonstrate the implementation of MFSF using numerical experiments in one and two spatial dimensions. Preliminary results show that MFSF has the potential to provide more accurate solutions with lower computational cost than conventional shock-capturing methods and can simplify many of the operations in shock-fitting methods with unstructured meshes.
UR - http://www.scopus.com/inward/record.url?scp=85130887820&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85130887820&partnerID=8YFLogxK
U2 - 10.1016/j.petrol.2022.110637
DO - 10.1016/j.petrol.2022.110637
M3 - Article
AN - SCOPUS:85130887820
SN - 0920-4105
VL - 215
JO - Journal of Petroleum Science and Engineering
JF - Journal of Petroleum Science and Engineering
M1 - 110637
ER -