论文部分内容阅读
多孔介质中复杂渗流系统刻画了诸如具复杂地质结构的油气资源运移、地下水污染传输等流体的溶混驱动过程,在实际工程实践中具有重要的应用。因此,对多孔介质中复杂渗流系统建立准确高效的数值模拟方法与完善的数值分析理论体系,对深刻揭示实际问题的物理机理、指导科学工程实践具有重要的理论价值与应用前景。但由于复杂渗流系统解的整体光滑性低及渗流区域形状的不规则,难以建立恰当的数值模拟方法与系统严谨的数值分析理论。本文基于Lagrange型有限元空间和Crouzeix-Raviart型有限元空间对流动界面的良好辨识能力与逼近性质,并结合间断Galergin有限元法、有限体积元方法保持物理量守恒、网格剖分灵活等显著优势,分别对不同地质结构下具有齐次界面和非齐次界面的一类复杂渗流系统(二阶椭圆界面渗流系统)建立相应的数值离散方法,并建立严谨的数值分析理论体系。据此,以实现对复杂渗流系统的高效数值模拟,为工程实践提供科学坚实的决策依据,并弥补当前理论研究成果的缺失,完善该类问题的理论分析体系。论文主要内容如下:1、在界面单元上构造依赖于界面的Lagrange型分片线性多项式空间,而在非界面单元上采用协调或非协调的线性有限元空间,结合间断体积元方法,对具齐次界面条件的各向同性二阶椭圆界面渗流系统建立间断浸入界面体积元离散格式,证明了离散格式解的存在唯一性,以及离散解的最优能量模与次最优L2模误差估计。2、采用水平集(level set)技术,化非齐次界面条件为齐次界面条件,将具非齐次界面条件的各向同性二阶椭圆界面渗流系统转化为具齐次界面条件的渗流系统。考虑到依赖于界面的Lagrange型分片线性有限元空间强烈的非协调性,我们在单元变分形式的基础上,在界面单元边界引入惩罚项,构造部分惩罚的浸入界面有限元离散格式,并利用Lax-Milgram引理证明格式的可解性。进一步,利用单元边界惩罚项良好的控制功能,建立部分惩罚浸入界面有限元格式的最优阶能量模和L2模误差分析理论。3、基于Crouzeix-Raviart型浸入界面有限元空间,对具非齐次界面条件的各向同性二阶椭圆界面渗流系统,构造部分惩罚的浸入界面Crouzeix-Raviart有限元离散格式,证明格式解的存在唯一性,并给出了离散格式的最优阶能量模估计。4、构造依赖于界面的Lagrange型浸入界面有限元空间,结合间断体积元方法,对具齐次界面条件的各向异性二阶椭圆界面渗流系统建立间断浸入界面体积元离散格式,证明了离散格式解的存在唯一性,以及离散解的最优能量模与次最优L2模误差估计。5、设计数值算例来验证理论分析的正确性.