Matrix acidization is an important technique used to enhance oil production at the tertiary recovery stage, but its numerical simulation has never been verified. From one of the earliest models, i.e., the two-scale model (Darcy framework), the Darcy–Brinkman–Forchheimer (DBF) framework is developed by adding the Brinkman term and Forchheimer term to the momentum conservation equation. However, in the momentum conservation equation of the DBF framework, porosity is placed outside of the time derivation term, which prevents a good description of the change in porosity. Thus, this work changes the expression so that the modified momentum conservation equation can satisfy Newton’s second law. This modified framework is called the improved DBF framework. Furthermore, based on the improved DBF framework, a thermal DBF framework is given by introducing an energy balance equation to the improved DBF framework. Both of these frameworks are verified by former works through numerical experiments and chemical experiments in labs. Parallelization to the complicated framework codes is also realized, and good scalability can be achieved.