有限体积元法是求解偏微分方程的主要方法之一,其优点是不仅能够更灵活地处理复杂的几何区域和边界条件,而且能够保持某些物理量的守恒性,因此近些年得到了较大的发展,但是现有的研究结果几乎都集中在常规的单介质问题,也就是非界面问题.然而,在许多实际问题中,所遇到的模型是具有多介质、多物理性质的界面问题,这类问题的数值求解具有更大的挑战,同时也更具有应用价值.在本文中,我们针对三类界面问题,研究了其相关的有限体积元法.首先考虑的界面问题是以压力函数和速度函数为未知函数的Darcy流模型,该模型的特点是计算区域被分解为有限个非重叠的子区域,在任意两个相邻子区域间的交界面上扩散系数矩阵可能是间断的.这种区域结构常用于带有断层的多孔介质流的模拟,另外为了在局部上获得高精度的梯度近似,也会用到这种区域结构.对于这类问题的数值求解,在每个子区域上需要独立地定义网格剖分,那么在相邻子区域间的交界面的两侧,网格节点是非匹配的.在这种非匹配的网格上,构造数值格式并进行相应的收敛性分析是有很大难度的.本文在混合有限体积元法的框架下,研究了该问题的数值求解.接下来考虑的界面问题是辐射三温能量方程组.辐射三温能量方程组属于非平衡辐射扩散方程组,在惯性约束聚变、磁约束聚变、天体物理、高超声速等领域中有广泛的应用.三温指的是电子、离子、光子的温度函数,该方程组用来描述辐射能在多介质物理区域内的传输过程以及电子、离子、光子这三者之间的能量交换过程.由于这类方程组的系数是强间断、强非线性及高度耦合的,因此相关的数值求解非常困难.首先要求得到的数值解具有良好的单调性(保正性)和守恒性.此外由于计算量较大,因此如何提高计算效率也是大家需要考虑的问题之一.本文将用于求解二阶椭圆方程的有限体积元法,推广到辐射三温能量方程组,研究了两种有限体积元格式并给出了一种网格自适应算法.最后考虑的是浸入界面问题.这类问题的求解区域被一条充分光滑的曲线划分为两个非重叠的子区域,这条曲线被称为界面,方程的扩散系数在每个子区域上是连续的,但在界面上是间断的.在做网格剖分时,不考虑界面的存在,在求解区域内独立地进行网格剖分,这时会出现界面穿过剖分单元的情形,此时如果采用传统的有限元法来进行数值求解,所得的计算结果无法达到最优阶的数值精度.为了改善这一缺陷,一些学者构造了一种新的有限元空间,即浸入界面有限元空间.若采用浸入界面有限元空间作为试探函数空间,所得的浸入界面有限元法的收敛阶可以接近最优.虽然这类方法在实际应用中发展迅速,但是相应的收敛性分析非常困难,至今仍有许多问题亟待解决.本文将浸入界面有限元空间应用到有限体积元法框架下,研究了一种带有惩罚项的浸入界面有限体积元法.本文前一部分共分五章,最后一部分是结论,前五章的内容包括：第一章是绪论,先简单介绍有限体积元法,然后介绍本文研究内容的背景和发展状况,即非匹配网格上Darcy流的混合法,非平衡辐射扩散方程以及浸入界面问题的数值解法.第二章考虑的是非匹配网格上的Darcy流模型,即前面提到的第一类界面问题.在每个子区域上,采用的是三角形网格剖分,其相应的对偶剖分为重心对偶剖分.在每个子区域网格上,选取最低阶的Raviart-Thomas空间来近似速度函数和压力函数,并按照标准的混合有限体积元法来离散原方程.由于网格的不匹配,近似速度空间在界面上不再满足法向流连续,并且变分方程中会涉及压力函数在界面上的迹.在本章中,我们在界面上引入线性Mortar元空间来近似压力函数的迹,并添加了一个用于提高近似速度流连续性的界面条件.这样所得的数值格式称为Mortar元混合有限体积元法,我们从理论和数值实验两个方面证明了该格式按L2范数具有最优的收敛阶.第三章考虑的仍然是非匹配网格上的Darcy流模型.在每个子区域上,采用的网格剖分、对偶剖分以及近似速度函数空间和近似压力函数空间与前一章中的相同,并仍然按照标准的混合有限体积元法来离散原方程.但是在界面上,本章采用双重拉格朗乘子空间来近似压力函数的迹.相对于每个子区域,其上的网格剖分在界面上会诱导出一个一维的网格,在这个网格上引入一个分片常数函数空间来近似压力函数的迹,这样便得到一个双重取值的拉格朗日乘子空间.此外,我们添加一个Robin型界面条件来增强近似函数在界面处的连续性.这样所得的数值格式称为非Mortar元混合有限体积元法.同样,我们从理论和数值实验两个方面证明了该格式按L2范数具有最优的收敛阶.第四章考虑的是辐射三温能量方程组,即前面提到的第二类界面问题.对于这类问题的数值求解,需要克服的主要问题有两个,分别是单调性和守恒性.本章从辐射三温能量方程组的守恒形式出发,采用合理的数值积分公式和近似方法来处理非线性项和间断系数.借助不同的积分公式,构造了两种守恒的有限体积元格式.由单调性分析和数值实验来看,第一种格式在许多网格上是单调的,我们推导出了相应的网格限制条件.然而,第二种格式不可能保持单调性并且在数值模拟的一开始便会产生大量的负数温度,这与实际问题不相符.因此,第二种格式通常被认为是不可用的.但是我们设计了两种后处理技术来克服这一问题,包括全局修补技术和截断法.数值结果表明,这两种后处理技术都具有较好的计算效果.最后,我们设计了一个基于残量型后验误差估计的自适应算法,使网格能够灵活地局部加细和粗化,这在很大程度上提高了计算效率.第五章针对浸入界面问题,即前面提到的第三类界面问题,构造了一种新的浸入界面有限体积元法.称被界面穿过的网格单元为界面单元,称不被界面穿过的网格单元为非界面单元.浸入界面有限元空间的构造方式为：在非界面单元上选取以节点为自由度的多项式,在界面单元上利用界面条件重新构造一个分片多项式.当选取浸入界面有限元空间作为试探函数空间时,相对应的有限体积元法被称为浸入界面有限体积元法.本章在已有的浸入界面有限体积元法的基础上,将原有的格式作了修正,在界面分划和与界面相交的边上,添加了两个惩罚项,分别用来限制函数值和法向流在界面上的跳跃.数值结果表明,这样所得的数值格式具有较好的稳定性,即使当扩散系数在界面上的跳跃较大时,数值解仍然能保持最佳的收敛性.通过严格的理论分析,我们论证了修正后的浸入界面有限体积元法的稳定性,进而得到该方法解的存在性与唯一性.