地质流体的性质和动力学行为是当前地球科学研究的前沿领域.铜陵冬瓜山层控夕卡岩型铜矿床成矿作用中矿质输运-化学反应耦合过程的数学模型是一个非局部的耦合抛物方程组初边值问题.本文考虑这一数学模型的数值模拟,用降阶法对此耦合方程组建立了一个具有二阶精度的差分格式.用能量方法给出了差分方程问题解的先验估计式,并证明了差分格式的可解性、稳定性和收敛性,其收敛阶在L2范数下关于时间步长和空间步长均是二阶的.最后给出了数值例子,数值结果和理论分析结果是吻合的.
The mathematical model of the transport-chemical reaction coupling process of the Dongguashan skarn copper deposit in Tongling is a nonlocal and coupled system of parabolic equations. A difference scheme with the second order accuracy is derived for the medel by the method of the reduction of order. It is proved that the constructed scheme is uniquely solvable and convergent with the convergence of order two both in space and in time. The numerical examples show the validity of the theoretical results.