两类流体问题的数学和数值分析

发布时间：2017-03-22 05:15

  本文关键词：两类流体问题的数学和数值分析，，由笔耕文化传播整理发布。

【摘要】：1973年,一篇题为解定常Stokes方程的协调和非协调元有限元方法的论文[23]中首次提到了Crouzeix-Raviart (C-R)有限元.这是一个定义在三角形/四面体上的P1有限元,自由度由形函数在剖分单元边或面上的中点处的取值构成.对计算的区域进行三角剖分,因此这种有限元的有限元空间中的函数在边的中点或面的中心是连续的,这也意味着此有限元空间并不是Sobolev空间H1(Ω)的子空间,而H1(Ω)空间是二阶椭圆边值问题的弱解空间.因此与自由度在剖分单元顶点上定义的协调P1有限元或Courant有限元[22]相比,C-R有限元经常被称作非协调有限元.在工程领域被称作Loof有限元[52].但是在数学领域人们是通过Crouzeix和Raviart (1973)写的那篇文章才开始深入了解C-R有限元.若想更加深入地了解C-R有限元,可以参考文献[16,17].本文中我们围绕着两种流体问题的数学和数值分析展开,流体速度都用的C-R有限元离散,压力用分片常数离散.第一类问题是不可压缩流体或微可压缩流体的Darcy-Forchheimer模型.另一种模型是自由流体、多孔介质流体和质量传递方程组成的耦合系统.多孔介质中流体运动的数学模型广泛应用于地下水、环境科学和油藏开发过程中的油水的运移等领域[5,7,63].该模型主要基于流体的质量、动量和能量守恒律.工程上一般比较关心压力、速度、温度和浓度等物理量,而关于地质、流体描述的许多参数例如重力、黏度、压缩系数、密度、孔隙度、渗透率、相对渗透率等对物理过程都起着重要作用,因此模型中围绕着以上物理量并且遵循以上各种守恒律引入了这些参数.由于介质的非均质性、各向异性导致引入的参数会剧烈地变化,因此许多物理量都是在平均意义下的,引入的模型及其相应的近似形式都有一定的适用范围.基于一些合理的假设,流体流动的数学模型可以简化,但仍表现为依赖于时间的强耦合非线性偏微分方程组.由于该微分方程组结构复杂,只在特定情况下才有解析解存在,因此提出能够保持系统物理性质的高精度的、有效的数值格式成为科学与工程中的迫切需求.Darcy定律描述了Newton流体压力梯度和速度呈现线性关系,已由多年的实验数据证实.这个线性关系只在合理的物理假设条件下才成立,例如流体流动地很缓慢,所有的惯性项可以忽略时[5]. Darcy定律的理论推导可以参考[58,78], Darcy流体对于油气采收和预防地下水污染有着重要意义.Forchheimer(1901)[37]在沙包中做流体实验时观察到当Reynolds数比较大(大约Re1)时,Darcy定律不再充分成立且会出现压力梯度和速度间的非线性关系.实际上,他通过对一大组实验数据进行观察发现,这个非线性关系是二次的[5].虽然目前还存在着对Darcy-Forchheimer方程泛函性的争议[6],但是非线性性已经从实验上、数值上[46]和理论上[55,69]被证实.Darcy-Forchheimer方程仍被用为模拟多孔介质中的高速流体的流动,尤其是在天然气井附近[3]. Forchheimer模型的推导可以参考[2,18,39,48,69,78],Forchheimer模型数学理论研究可以参考[4,47,33,74], Forchheimer方程是单调非线性非退化的,类似的问题有p-Laplacian问题、拟Newtonian流和Ladyzhenskaya流体问题,处理这一类单调非线性算子的技巧和方法可以参考[30,31,34,35,44,72].近年来,Darcy-Forchheimer模型的数值离散方法已有一些,[43,53]中引入了primal非协调混合有限元方法,[71]中给出了协调混合有限元方法,这种方法使得我们构造的迭代格式中有限元系数矩阵更加稀疏,从而节省了内存空间和CPU运行时间Park [62]提出了微可压缩流体Darcy-Forchheimer模型的对偶混合元方法Girault和Wheeler [43]验证了Darcy-Forchheimer问题弱解的存在唯一性.他们提出一种primal混合元离散格式：速度用分片常数元离散、压力使用不连续的C-R线性元离散(参考[23]).他们提出了交替方向迭代算法来解非线性方程,并对迭代算法和有限元格式的收敛性都给予证明Lopez等[53]使用[43]中给出的格式进行数值实验.对离散后的非线性方程组提出了Newton迭代方法并与[43]中的交替方向算法对比.进一步,在[53]中也提出了另外一种有限元,这种有限元使得对压力的近似更加光滑,压力使用连续的P1 Lagrange有限元逼近,速度还是用不连续的分片常数元近似Salas等[71]对这种数值方法给出了详细的理论研究并提出了另外两种数值方法：两种方法中压力还是用连续的P1 La-grange有限元逼近,速度的近似空间变成一次多项式,一种是速度用连续的P1 Lagrange有限元离散,另一种方法速度用不连续的C-R线性元离散.但是提出的这两种新方法到目前为止还没给出数学和数值研究.本文目标之一在于提出不可压缩质量守恒方程与Darcy-Forchheimer方程耦合问题的一种新的对偶离散格式.速度和压力分别由不连续的C-R元和分片常数元离散,这种混合元经常用来处理Darcy-Stokes模型.在这种情形下,消去速度,该模型表现为压力的二阶椭圆方程.进一步地,假设流体微可压缩,我们对微可压缩质量守恒方程与Darcy-Forchheimer方程耦合问题使用相同的离散格式,此时消去速度,模型表现为压力的抛物方程.到目前为止,已经有很多关于耦合的自由流和多孔介质流问题的文献.这类耦合问题的数学模型是由自由流体区域的Stokes方程,多孔介质区域的Darcy方程再加上一些合适的界面条件组成.这些界面条件包括Beavers-Joseph-Saffman条件[8,70],通量连续条件,力的平衡条件.这个问题在数学和数值分析方面也很具有挑战性：两区域内方程的解具有不同的正则性,交界面上流体的切向速度不连续,变分形式中积分项在交界面上要比在区域内部少一维,要保证不降低解的正则性和收敛阶难度不小.对于这个模型,已有不少文献[24,38,51,54,66]提出了稳定的和收敛的数值格式.当今社会面临的一个很严重的问题就是由地下储藏设备泄漏、化学药品泄漏还有各种人类活动导致的地表水和地下水污染.由Stokes方程、Darcy方程和质量传递方程组成的耦合系统可以用来描述水中泄漏的污染物的传播并且评估污染风险.这个模型仅在[76]中给出系统的研究,其中流体的粘性系数假设与溶质的浓度无关,这样的假设实际已经解耦了流动方程和浓度方程.我们接下来研究的就是二维区域上由耦合的Darcy-Stokes方程和对流扩散方程组成的系统,其中对流扩散方程可以用来模拟流体中溶质的运动,并且此时流体的粘性系数是依赖于溶质的浓度.这个耦合问题实际含有两层耦合含义：一是两区域间的耦合,在不同的区域上有不同的流动方程、不同的扩散系数和不同的源汇项,只在交界面上进行物理量的传递；二是流动方程和质量传递方程之间的耦合,通过流速和浓度彼此相互影响.因此这双重耦合会导致整个系统异常复杂.对两个区域上的方程都用混合形式离散的数值算法又可分为两类：一类是在不同区域用不同的有限元离散；另一种是在两区域使用相同的有限元离散.使用同一有限元的优势在于不论是在理论分析还是在程序实现中处理交界面条件更加方便,同时也使得编写程序代码时可以较少考虑单元所在区域,从而编写效率更高.我们也用同一元的思想对这个系统提出了一种稳定化混合元方法.在整个Stokes和Darcy区域对流场压力和速度分别采用分片常数和C-R有限元空间来逼近.这里使用C-R元是因为它具有与分片常数压力组合易于满足inf-sup条件、能保持分片单元质量守恒、二维和三维情形都容易实现等好处.然而,Mardal等人[54]指出对Darcy方程使用C-R元时离散格式并不收敛.众所周知,C-R元不满足离散的Korn's不等式,因此Hansbo和Larson [45]通过使用一个罚项[16]加罚速度在单元边界上的跳跃来满足离散的Korn's不等式.对于质量传递方程中浓度的离散,我们使用经典的Lagrange有限元.单独区域内(特别是渗流区)的耦合流动和传质问题的工作,也就是多孔介质区域内混溶驱动问题已有很多文献[9,73,49,36,19,27,28,68,20,47,32,29].本文的组织结构如下：在第一章中,我们简要地介绍了多孔介质中Darcy-Forchheimer律的数学模型,基于流体的物理特性给出质量守恒方程,并根据组分质量守恒推导出多孔介质流体中流体浓度的对流扩散方程.接着介绍了本文中常用的一些符号,包括函数空间及其范数定义,并给出后面章节理论推导过程中常用的几个引理.在第二章中,我们提出了稳定的对偶混合有限元格式来离散不可压缩Darcy-Forchheimer流体方程.速度和压力分别由非协调的Crouzeix-Raviart有限元和分片常数元近似.我们验证了inf-sup条件并由这个条件成立和非线性算子的单调、强制、半连续性证明了离散解的存在唯一性.给出了速度L2和L3范数,压力L2范数的先验误差估计.最后我们用数值实验验证了理论分析的正确性,并对[59]中提出的数值格式和我们提出的格式做了数值比较,从某些方面说明我们的方法具有优越性.在第三章中,我们使用了跟第二章相同的稳定的对偶混合有限元格式来离散微可压缩Darcy-Forchheimer流体方程.通过引入Darcy-Forchheimer速度和压力的投影来推导出半离散数值格式和全离散数值格式的先验误差估计.在第四章中,我们考虑了一个由Stokes方程、Darcy方程以及质量传递方程组成的耦合系统,这个系统用来描述带传质过程的自由流和多孔介质流耦合问题模型,系统中流体的粘性系数是与浓度相关的.我们给出一种全离散有限元格式,从而把流体流动方程和浓度方程分离开来.在整个区域对速度和压力分别使用Crouzeix-Raviart有限元和分片常数元近似并给出一个类似于[16]中的稳定项加罚速度在单元边界上的跳跃,对于浓度方程,我们用协调的Lagrange有限元来离散.通过对方程弱解做一些合理的正则性假设,我们给出了离散格式速度、压力和浓度的先验误差估计.
【关键词】：Darcy-Forchheimer方程 Stokes和Darcy方程 Crouzeix-Raviart元 流动和迁移 存在性 误差估计
【学位授予单位】：山东大学
【学位级别】：博士
【学位授予年份】：2015
【分类号】：O241.82
【目录】：

• 中文摘要10-15
• 英文摘要15-22
• 第一章 预备知识22-30
• 1.1 Darcy-Forchheimer模型22-23
• 1.2 质量守恒定律23-26
• 1.3 基本符号和引理26-30
• 第二章 不可压缩Darcy-Forchheimer流模型的稳定Crouzeix-Raviart元方法30-54
• 2.1 引言30-32
• 2.2 弱形式32-34
• 2.3 有限元离散34-38
• 2.4 误差估计38-41
• 2.5 数值算例41-54
• 第三章 微可压缩Darcy-Forchheimer流模型的稳定Crouzeix-Raviart元方法54-64
• 3.1 引言54-55
• 3.2 半离散有限元格式及误差分析55-57
• 3.3 全离散有限元格式及误差分析57-59
• 3.4 数值算例59-64
• 第四章 耦合自由流和多孔介质流的传质过程的稳定非协调有限元方法64-84
• 4.1 引言64-67
• 4.2 弱形式67-69
• 4.3 有限元离散69-74
• 4.4 误差估计74-84
• 参考文献84-92
• 致谢92-94
• 攻读博士学位期间完成论文情况94-96
• 作者简介96-97
• 附件97

【相似文献】

中国期刊全文数据库 前10条

1 ;Darcy's Law Expressed by Chemical Index[J];Journal of China University of Geosciences;2000年04期

2 四旭飞;陈占清;缪协兴;岳建华;;利用瞬态法提取岩样非Darcy流渗透特性(英文)[J];湖南科技大学学报(自然科学版);2006年03期

3 唐平,黄先伍,李天珍;岩样非Darcy流的渗流特性试验研究[J];湘潭师范学院学报(自然科学版);2002年02期

4 许友生,刘慈群,林机;Darcy渗流定律的微观界定及其应用[J];应用数学和力学;2004年03期

5 张茂省;孙萍萍;朱立峰;裴赢;毕俊擘;;原位“Darcy实验”方法研究与实践[J];地质通报;2013年06期

6 许友生,刘慈群,林机;MICROCOSMIC BOUND THEOREM OF DAYCY'S LAW AND ITS APPLICATION[J];Applied Mathematics and Mechanics(English Edition);2004年03期

7 郑健;南北;;承压水完整井Darcy流和非Darcy流水头分布计算方法的分析比较[J];地质学刊;2010年04期

8 孙明贵;黄先伍;李天珍;雷光宇;茅献彪;;石灰岩应力-应变全过程的非Darcy流渗透特性[J];岩石力学与工程学报;2006年03期

9 涂洪亮;;多孔介质中一类Darcy-Forchheimer流的结构稳定性[J];黄冈师范学院学报;2007年03期

10 刘忠玉;纠永志;乐金朝;孙丽云;;基于非Darcy渗流的饱和黏土一维非线性固结分析[J];岩石力学与工程学报;2010年11期

中国博士学位论文全文数据库 前5条

1 王颜;两类流体问题的数学和数值分析[D];山东大学;2015年

2 潘浩;多孔介质中Darcy-Forchheimer渗流数值计算[D];山东大学;2012年

3 刘玉;水沙混合物非Darcy裂隙渗流的试验研究[D];中国矿业大学;2014年

4 沈宇晶;Navier-Stokes/Darcy多区域耦合问题的多重网格方法[D];浙江大学;2014年

5 董李娜;Darcy-Stokes问题的一致收敛有限单元及离散的de Rham复形[D];郑州大学;2014年

中国硕士学位论文全文数据库 前5条

1 赵丹汇;一般Darcy-Forchheimer问题的块中心有限差分法[D];山东大学;2014年

2 张杰;基于非Darcy渗流的变荷载作用下饱和黏土一维固结分析[D];郑州大学;2011年

3 韩小妹;低渗透岩石非Darcy渗流实验研究[D];清华大学;2004年

4 孟俊玲;奇异摄动Darcy-Stokes问题的有限元分析[D];郑州大学;2012年

5 沈漪;多孔介质中二相流的基于动态网络的Darcy尺度建模[D];苏州大学;2007年

  本文关键词：两类流体问题的数学和数值分析，由笔耕文化传播整理发布。

本文编号：260943