A fully Coupled thermal-hydrodynamic–chemical numerical model for simulating gas hydrate-bearing sediments dissociation based on peridynamic differential operator method

During methane hydrate-bearing sediments (MHBS) extraction multiphase flow with heat transfer transitioning MHBS phase followed by pore dynamic evolution. Which involves a complex thermo-hydrodynamic (THC) coupled process. Effective theoretical and fully integrated numerical models are essential for evaluating the economic feasibility of gas hydrate production and comprehending the fundamental physical concepts. This study presents a novel approach to constructing a fully integrated THC numerical model by employing a peridynamic differential operator (PDDO). The dimensionless governing equations of Darcy’s law and the heat source approach are transformed into a non-local integral form. The Euler forward difference method is used for time integration. The accuracy and reliability of the newly developed PDDO THC model were tested by comparing it with experimental examples from other popular simulators, which provides an alternative to explicit modeling of the phase change process involving multiphase flow with heat transfer and may find comprehensive and useful applications to a variety of industrial and geophysical processes.