摘要

鉴于求得解析解的困难,需要采用数值方法和算法得到大量科学和工程问题的计算模拟结果。这对求解饱水孔隙岩石中化学溶解面非稳定性问题而言也不例外。由于这类非稳定性问题可具有常规解和非常规解,很有必要探讨不同数值算法对模拟饱水孔隙岩石中化学溶解面非稳定性的影响。因此,本文考虑了与常用有限元相关的2种不同数值算法。在第1种数值算法中,孔隙率、孔隙流体压力和酸液浓度被选为基本变量。在第2种数值算法中,孔隙率、孔隙流体流速和酸液浓度被选为基本变量。相关的计算模拟结果表明:(1)第1种数值算法可真实地模拟非稳定化学溶解面在化学溶解系统中传播时的演化过程;(2)第2种数值算法不能模拟非稳定化学溶解面在化学溶解系统中传播时的演化过程;(3)过度的数学微分运算是导致第2种数值算法失效的主要原因。