LAUSR

Matrix acidization simulation is a challenging task in the study of flows in porous media, due to the changing porosity in the procedure. The improved DBF framework is one model to do this simulation, and its numerical scheme discretises the mass and momentum conservation equations together to form a pressure-velocity linear system. However, this linear system can only be solved by direct solvers to solve for pressure and velocity simultaneously, since zeros appear in the diagonal of the coefficient matrix. Considering the large-scale attribute of matrix acidization simulation, the solving time of direct solvers is not intolerant. Thus, a decoupled scheme is proposed in this work to decouple the coupled pressure-velocity linear system into two independent linear systems: one is to solve for pressure, and the other one is to solve for velocity. Both of the new linear systems can be solved by parallel and iterative solvers, which guarantees the large-scale simulation can be finished in a reasonable time period. A numerical experiment is carried out to demonstrate the correctness of the decoupled scheme and its higher computing efficiency.