Numerical modelling and analysis of reactive flow and wormhole formation in fractured carbonate rocks

Piyang Liu, Jun Yao, Gary Douglas Couples, Zhaoqin Huang, Hai Sun, Jingsheng Ma

Research output: Contribution to journalArticlepeer-review

52 Citations (Scopus)
165 Downloads (Pure)


Combining the two-scale continuum model and the discrete fracture network model, a continuum-based model is developed that calculates the reactive flow of acid in carbonate rock with a complex fracture network. The locations of fractures in the model are explicitly defined and the method can capture complex geometric relationships. The governing equations are discretized by the finite-volume method, where the diffusion term and convection term are discretized using the two-point flux approximation (TPFA) scheme and the upwind scheme, respectively. The physical domain is discretized by Delaunay triangulation. To keep the robustness and efficiency of the TPFA scheme, the optimization algorithm is used to move the centroid node of the control volume to improve the orthogonality of the grids. Numerical simulations of reactive flow in 2D fractured porous media, in cases with simple and complex fracture arrays, under linear and radial flow conditions, are presented. In particular, a sensitivity analysis of the dissolution process with respect to the presence of fractures, fracture aperture, fracture distribution, and acid injection rate, is conducted.
Original languageEnglish
Pages (from-to)143-157
Number of pages15
JournalChemical Engineering Science
Early online date17 Jun 2017
Publication statusPublished - Nov 2017


Dive into the research topics of 'Numerical modelling and analysis of reactive flow and wormhole formation in fractured carbonate rocks'. Together they form a unique fingerprint.

Cite this