论文部分内容阅读
土壤的结构是相当复杂的,包含于其中的固、液、气三者相互联系,相互制约,构成了一个矛盾的统一体,为生长于其中的植物提供必要的生活条件,是土壤肥力的物质基础.然而,近年来随着工农业生产的发展、工业“三废”大量排放、农药化肥大面积使用、核能利用带来的放射性废物日益增多、人口增长和城市化带来的垃圾不断增加,土壤的化学废物污染已经成为人们关注的又一重点,从而刺激了对于在土壤介质中的溶质运动情况的相关探索.对于溶质迁移过程的研究已经成为当代土壤物理学,特别是在预测杀虫剂、营养矿物质、重金属或其他溶质在土壤中的运动过程中充当着非常重要的工具.
在已有的土壤溶质运移模型中,对流-弥散模型是一类基本的运动方程,该模型具有坚实的物理基础,可以描述质量、热量的运输及反应扩散过程等众多物理现象,且形式简单.但是,由于土壤结构的复杂性,使得该模型不能用于野外大尺度范围的溶质运移问题.在聚合介质中,土壤是由一些列孔隙组成的,溶质通过这些孔隙的速度有快有慢,从而分为裂缝孔隙和岩石孔隙,岩石孔隙就伴随着不动水的存在.在[21]中研究了吸附性多孔介质中化学溶质的横向扩散的相关问题.在多孔介质中,流体分为不动区和动区,流动的水主要存在于大的孔隙中,而溶质在水中的流动也主要存在于这一区域.溶质的在水中迁移主要是基于对流以及纵向的扩散等原因.不动水主要存在于聚合体的内部以及不同聚合体的衔接处.溶质在两种流体区域之间的扩散相容是与它们的浓度差异成比例的,流动液体中的溶质与吸附在土壤颗粒表面上的溶质的动态平衡使得在动水区附近形成了动态的土壤区域,而静止土壤区域主要位于聚合体内的小孔隙附近或者在不动水区边界处.只有在化学溶质弥散渗透穿过不动水区域的流体外层后吸附才会发生.在对流-弥散模型基础上建立起来的非均衡溶质运移两区模型,由于考虑了土壤中不动水体的存在以及不动水体在溶质运移过程中所起的作用和溶质运移过程中的吸附、降解等因素的影响,因此能更准确的描述土壤中溶质的运移过程.
对于较复杂的现实问题,很难用解析的方法求解,而数值模拟是一种有效的求解方法.常用的数值计算方法有有限差分法和有限元法以及有限体积元法,但用标准的有限差分方法求解数学模型时常常失效,根本原因在于对流项的存在.由于对流为主的弥散方程具有双曲性质,中心差分格式虽然关于空间步长具有二阶精度,但会产生数值弥散和非物理力学特性的数值振荡,使数值模拟失真.特征线和经典的有限差分、有限元方法相结合所得到的特征差分、特征有限元方法,可以更好的反映出对流弥散方程的一阶双曲特性,这一类方法考虑沿着特征线方向的离散,利用了物理力学性质,可以有效的克服数值振荡,减少截断误差,大大提高计算精度.崔明老师和张德生老师曾对两区模型应用特征有限元方法进行数值分析.本文就是在以上思想的基础上,在第二章中将混合有限元、有限体积元相结合来求解,给出相应的误差估计并各处相应的一维实例.