摘要

潜蚀的模拟需要计算土中的渗流过程以及土中细粒随孔隙水流的运移过程。常规有限元法在求解纯运移方程时存在稳定性问题,为此提出了一个采用有限元法求解渗流方程、采用有限体积法求解细粒运移方程的交替方法。鉴于常规有限元法不能给出满足有限体积法输入要求的流速场,给出了一个基于单元不平衡流量总体再平衡修正单元边界流速的算法。利用该算法对有限元计算的单元边界流速进行局部守恒处理,即可实现在同一有限元网格上利用有限体积法求解细粒运移方程,因此可以方便地和现有的有限元计算程序结合。算例验证表明,所提出的局部守恒处理算法及交替解法计算效率高且具有可接受的精度,是模拟潜蚀问题的一条简单实用途径。

全文