多孔介质中化学—热—水力—力学耦合分析与混合元方法

近年来,非饱和多孔介质的热—水力—力学耦合行为的研究引起了工程界的广泛的注意。人们已经投入了很大的精力来进行多孔介质瞬态THM耦合行为的定量研究。 本论文致力于这个课题的三方面的工作:1)建立了求解多孔介质中化学—热—水力—力学耦合行为的数学模型;2)发展了这个数学模型初值、边值问题的混合元方法的数值解法;3)建立多孔介质中的化学—热—水力—力学（CTHM）耦合本构模型。 这个数学模型由一组偏微分方程组成:包括孔隙水、孔隙气的质量守恒方程;混溶于孔隙流体中的化学污染物的质量守恒方程;热量守恒方程和多孔介质混合物的总体动量守恒方程。为了模拟混溶污染物在非饱和多孔介质中的传输,模型中包含了污染物传输的六种控制机制:对流、分子弥散、机械逸散、吸附、降解、不动水效应。 采用有限元法求解数学模型的初值、边值问题。考虑到化学污染物质量守恒方程在物理上的对流—扩散的特性和数学上的非自伴随和双曲线的特性,采用交错算法进行求解,即:对化学污染物的质量守恒方程和其余的控制方程分别进行离散和求解。 化学污染物的质量守恒方程采用依赖于时间的对流—扩散方程的隐式特征线Galerkin方法进行离散和求解,其中污染物的浓度作为基本未知量。该方法的基本思想是对浓度在时域中进行针对物质粒子（Lagrangain）而不是空间点（Eulerain）的离散。 非饱和多孔介质中热—水力—力学（THM）耦合行为的有限元模拟通常要求解u-pw-pa-T混合形式的方程。混合方程的基本未知数u,pw,pa,T分别为固相位移、孔隙水压力、孔隙气压力和温度。位移场、水压力和气压力场及温度场分属于不同的函数空间,因此应该对他们采用不同的有限元插值近似,以进行满足自然边界条件的耦合控制方程弱形式的空间离散。 流体和固体力学u-p混合形式的有限元公式研究表明,u-p有限元插值的函数空间必须满足Babuska-Brezzi条件或更简明的Zienkiewicz-Taylor分片试验。这些限制条件判定了应用方便的位移（或速度）和压力采用等阶插值的低阶形函数的单元在不可压缩条件下压力场会产生虚假的数值振荡。 与u-p形式的混合公式相比,非饱和多孔介质的u-pw-pa-T公式更加复杂,求解以u,pw,pa,T为基本未知数的整体半离散方程、高斯积分点上耦合的非线性本构方程迭代及计算一致性切线模量矩阵所需要的时间要长很多。因此就更加迫切需要发展低阶高精度粗网格的有限元方法避免虚假的数值振荡,同时在满足数值计算精度基础上又能降低计算费用。 本论文基于一点积分混合应变元在固体中的成功应用,提出了饱和、非饱和多孔介质中THM耦合问题的稳定的一点积分混合应变元。本文提出的偏微分方程组的Galerkin弱形式基于固体力学中的胡海昌-Washizu三变量广义变分原理,将之推广到由固体骨架和不混溶的两种孔隙流体组成的三相非饱和多孔介质与温度场的相互作用的耦合场分析中。