摘要

地质流体的性质和动力学行为是当前地球科学研究的前沿领域.铜陵冬瓜山层控夕卡岩型铜矿床成矿作用中矿质输运-化学反应耦合过程的数学模型是一个非局部的耦合抛物方程组初边值问题.本文考虑这一数学模型的数值模拟,用降阶法对此耦合方程组建立了一个具有二阶精度的差分格式.用能量方法给出了差分方程问题解的先验估计式,并证明了差分格式的可解性、稳定性和收敛性,其收敛阶在L2范数下关于时间步长和空间步长均是二阶的.最后给出了数值例子,数值结果和理论分析结果是吻合的.