论文部分内容阅读
渗流现象是多孔介质中普遍客观存在的运动规律,其对水工建筑物稳定性、地下水资源利用、污染控制、石油天然气开采等工程问题均具有重要的影响,随着国内外在上述领域中工程需求的增加,渗流研究日益成为热点问题。本文在饱和、非饱和渗流理论基础上,结合非饱和渗透系数研究成果,对达西速度的表达式进行变形,建立了饱和一非饱和渗流数值模型。该模型中以测压管水头和饱和度为变量对渗流场进行描述。将渗流域分为饱和区与非饱和区,为避免研究复杂的多相流问题,可以认为非饱和区是与大气相通的,则测压管水头中的压力与大气压相等,使测压管水头等于位置水头,从而渗流控制方程在两个区域中分别退化为关于测压管水头的拉普拉斯方程和关于饱和度的对流方程。空间用伽辽金有限元方法,时间采用隐式差分方法对方程进行离散。分别对饱和区的拉普拉斯方程采用预处理共轭梯度法(PCG)、非饱和区的对流方程采用Gauss—Seidel方法进行耦合迭代求解,得到渗流场的测压管水头分布和饱和度分布,并引入VOF的概念,由饱和度的分布自动扑捉渗流自由面。为测定非恒定水位边界条件下的土坝渗流自由面和坝体内测压管水头的变化,在实验室内进行了模型土坝的非恒定渗流实验。通过对实验结果进行分析,初步认识了不透水地基上均质土坝在非恒定边界条件下的非稳定渗流的物理力学特性。推导了双介质土柱表层水入渗的下渗锋面解析表达式、稳定状态时测压管水头分布与饱和度分布的解析表达式,并给出了有周期性边界条件的有限潜水层内渗流自由面的解析表达式。用本文提出的饱和—非饱和渗流模型分别对几个解析解问题以及渗流实验进行了数值模拟,通过所得结果与解析解和实验的对比,证明了模型的合理性、可靠性和准确性。将饱和—非饱和渗流模型用于几个工程实例,其二维模拟结果同渗流电模拟实验结果进行了对比,证明本文的渗流模型完全能够准确计算坝内的渗流浸润线以及逸出面位置;三维数值模拟的结果在相近剖面处与二维结果作了比较,考虑三维网格剖分及剖面截取位置的差异,三维结果与二维结果吻合的较好,证明本文的饱和—非饱和渗流模型在三维计算中是可靠的且具有较高的实用价值。在串行程序的基础上,应用OpenMP编程接口进行了并行化改造,对渗流场实现了并行求解。通过相同位置剖面与二维结果的比较很好的证明了计算结果的可靠性和准确性;通过对不同规模问题的加速比的测定,证明了应用OpenMP改造的串行程序能够很好的实现渗流场并行求解,达到减小计算耗时、便于工程应用的目的。