动网格在非饱和饱和界面数值模拟中的应用研究进展

为了探讨结构和非结构动网格技术在地下水非饱和-饱和数值模拟领域未来的发展趋势, 总结了非饱和-饱和耦合数值模拟研究现状, 介绍了动网格技术原理及运动边界结构和非结构网格的变形方法, 综述了动网格技术在非饱和-饱和分界面的应用现状及存在的不足, 探讨了相关研究的未来发展趋势。综合分析表明: 结构动网格和非结构动网格均存在其固有优缺点, 结构、非结构混合网格以及多种动边界处理方法的结合使用在非饱和-饱和耦合数值模拟研究中具有重要的应用价值。在模拟潜水面的变动时, 可将多种网格变形方法结合使用, 当潜水面位置和形状变动较小时, 采用弹簧法更新网格; 当潜水面位置变化较大但形状变化较小时, 采用重叠结构动网格技术或铺层法更新网格; 当潜水面形状变动较大时, 则采用网格重构法更新网格, 从而更精确地模拟非饱和-饱和分界面的变化和移动。相关研究为场地非饱和带土壤与饱和带地下水协同防治工作提供了科学指导。      

基于FEFLOW的三维土壤-地下水耦合铬污染数值模拟研究

`土壤-地下水耦合数值模拟是定量刻画水流和溶质运移的主要手段。现有大范围场地尺度的研究受到数据采集难度及模拟计算量的限制,多是将土壤和地下水分成两个系统,这种方式不利于模型之间的计算反馈,易出现计算误差,因此将土壤和地下水作为整体系统研究具有重要意义。为精确刻画实际场地土壤-地下水系统中污染物迁移规律,揭示变饱和反应溶质迁移模型的参数敏感性,以某铬污染场地为研究对象,基于现场试验及前人研究所获数据,采用Galerkin有限元法建立三维土壤-地下水模型,定量描述六价铬在土壤-地下水中的迁移规律。在此基础上,通过改变补给条件,研究潜水面在土壤-地下水系统中的波动。并讨论阻滞系数和反应常数对溶质运移的影响。结果表明:在土壤中,污染物最大水平迁移距离为场地东南侧300 m;地下水中污染晕最大分布面积约为1.632 km2;垂向上土壤中的六价铬仅需15.6 h即可下渗至潜水面,第6天贯穿含水层。当潜水面随着补给量变化而波动时,地下水中六价铬会随水流进入土壤,影响土壤中污染分布。对溶质运移参数的讨论显示,当反应常数由0增大至10?6 s-1时,迁移出场区边界时地下水中污染物浓度约减少2000 mg/L,较难迁移至涟水河。基于FEFLOW的数值模型,能够解决各系统之间交互性差的问题,提供较为精确的模拟结果。