基于传递矩阵的两相介质波动问题时域显式有限元方法稳定性研究

饱和两相介质波动问题的时域显式有限元方法具有条件稳定性。这类方法的稳定性可通过方法中动力反应时域递推计算格式传递矩阵传递因子的数值进行评估，传递因子的数值越小，方法的稳定性越好。针对饱和两相介质波动问题的两种时域显式有限元方法，基于其动力反应时域递推计算格式的传递矩阵的性质，在两相介质体系物理参数取值不同的条件下，对两种方法稳定性性质的差异进行了对比研究。结果表明：渗透系数的取值对两种时域显式有限元方法稳定性性质的差异具有较为显著的影响。当渗透系数取值较大时，两种方法的稳定性性质没有明显的差异；当渗透系数的取值较小时，两种方法的稳定性性质将表现出明显的差异，方法1的稳定性要优于方法2。