多孔介质渗流问题守恒特征线数值方法及理论

多孔介质渗流问题守恒特征线数值方法及理论

作者:师大云端图书馆 时间:2015-07-24 分类:期刊论文 喜欢:4290
师大云端图书馆

【摘要】多孔介质中的流体运动方程广泛应用于地下水、环境科学和油藏模拟等领域[11,62].模型主要包括了流体的流动和质量的转移,体现着流体本身的质量守恒、能量守恒等物理性质.其中重力、粘度、毛细管力、密度等对该物理过程都起着重要作用,而工程上考虑的重要参数一般包括流体速度、压力、温度以及浓度等.数学上该模型由描述多孔介质中的速度方程即Darcy定律和描述组分混溶传质过程的浓度方程构成.基于一些合理的假设条件可以将方程简化,但仍表现为依赖于时间的非线性耦合问题.对此模型研究保持物理性质的高效数值模拟格式有着重要的实际意义和价值.渗流的速度方程通常由Darcy定律描述[26,45],它体现了流体的平均速度与压力梯度的线性关系.由于压力和速度的重要性,对不可压缩条件下的质量守恒方程和Darcy方程构成的模型一般采用的是混合元离散格式,这样即保证速度和压力的精度,又能保证局部的质量守恒.经典的混合元空间Raviart-Thomas、Brezzi-Douglas-Marini等保证了速度法向分量的连续性,并且都给出了解的存在唯一性及最优的误差估计[17,18,19,65].大量的稳定化技术也可以用来求解Darcy问题.例如Masud和Hughes[60]加入的稳定项使得连续的速度压力都是有效的,当然还有加一个平方项或者引入Bubble函数等稳定化技术.有时需要速度是连续的,也就是切向方向速度也是连续的,Arbogast和Wheeler[6]给出了一类连续速度通量的逼近格式,虽然损失了散度范数的最优收敛率,但还是得到了最优的L2模误差估计.Darcy流与Stokes流耦合问题在实际中也有着广泛应用.即在多孔介质的基础上耦合一个自由流区域的Stokes问题,两个子问题内边界由Beavers-Joseph-Saffman条件确定[72],这两个子问题的正则性要求不一样,并且内边界上切线速度是不连续的Layton等人[56]详细介绍了该模型,并引入Lagrange乘子证明了弱解的存在唯一性Yotov[66]给出了一类求解技术,在两个子区域采用不同的离散空间,即在Stokes区域上用DG求解,而Darcy区域用混合元技术求解.这样不利于处理内边界条件和编程.而相同网格剖分下的同一种混合元离散无疑是高效的,Burman[21]给出非协调的Crouzeix-Raviart有限元空间逼近格式.芮和张[69]对其改进,引入了更简单的罚项,保证稳定的基础上还有局部质量守恒性Arbogast和Brunson[3]还将[6]提出的连续速度空间用到该耦合问题中,并给出了最优的L2误差估计.对于奇异扰动问题即Brinkman方程,Mardal[59]给出了该问题弱形式解的存在唯一性.提出了一种绝对稳定的有限元格式,并对Brinkman问题给出了系统的理论分析和最优的误差估计.在此基础上后面陆续有很多的离散格式[61,83].而组分的混溶传质过程最终可以由一个关于浓度的对流扩散方程所表示.当对流占优时,由于数值振荡和物理弥散的效应,传统的有限元和有限差分在进行模拟时效果并不理想,此时特征线技术很好的解决了这一难题.该技术从数学上将对流扩散问题转化一个等价的易处理的扩散问题,体现出很好的优越性,并且对于时间步长不再有很强的限制,用来模拟大时间步长的实际问题时效果明显.该思想早期是由[43]提出的向前追踪特征线方法,但是这样的技术破坏了原始空间剖分,存在很大的局限性Douglas和Russell[54]1982年提出了向后追踪的修正特征线技术(MMOC),给出基于有限元和有限差分方法下的离散格式,得到有限元离散下最优的H1与L2模估计.此后该技术迅速得到广泛应用,Russell[71]将其运用到了不可压缩的混溶驱动模型中去,压力方程采取有限元离散,并给出了最优阶的误差估计Ewing等[35,36]也采用MMOC技术求解浓度方程,但对压力方程采用更加有效的混合元技术,使得速度、压力和浓度均能达到最优估计.此外在模拟对流扩散问题时,一个重要的性质就是要保证质量守恒,也就是组分的总量在不考虑源汇项的情况下随时间是平稳不变的,上述技术并不能满足这一点.芮[68]提出了对流扩散问题的质量守恒格式,此时的速度有连续性的要求Celia[22]给出了ELLAM技术保证了质量守恒,但是计算有些复杂,Wang[77,78]给出了ELLAM格式在对流扩散问题中的收敛性及最优估计,并给出了一些实际模拟.对于更为复杂的可压缩混溶驱动的模型Douglas和Roberts[53]给出可压缩混溶驱动的数学模型,并给出了基于有限元和混合元方法的半离散格式的误差估计,此后陈[23]给出了混溶驱动的全离散数值格式.袁和程也分别给出了基于MMOC的有限元误差估计[24,85].韩和吴[47]对于Stokes方程提出了基于交错网格的混合元离散技术,即速度分量和压力使用三套不同的网格剖分,数值模拟更加方便.我们在此基础上加以修改得到对Darcy方程适用的连续速度逼近.参照对流扩散问题的守恒格式[68],得到渗流中的不可压缩混溶驱动模型的质量守恒的特征线数值格式及误差分析,随后给出了可压缩混溶驱动模型的全离散数值格式及误差分析,并分别给出数值实验进行验证.本文的组织结构如下:在第一章中,介绍了多孔介质渗流问题的数学模型.基于流体的物理性质给出Darcy定律及质量守恒方程.结合状态方程给出了混溶驱动中组分传质的质量守恒方程,得到所谓的浓度方程.同样也给出了可压缩流体的浓度方程.给出了常用的函数空间的记号和范数定义,最后给出了本文理论推导需要的一些不等式.在第二章中,对于Darcy司题,给出了基于交错网格剖分下的连续速度逼近.借助RT空间,得到很好的插值性质.给出了数值格式的误差分析.虽然得不到散度范数的误差,但还是给出了L2误差估计.数值算例中与Arbogast和Wheeler[6]中提到的混合元空间进行了比较,数值结果接近,但是用的自由度相对较少,这样可以减少时间复杂度.在第三章中,首先用连续元求解稳态的Dracy-Stokes耦合问题,给出了全离散格式及误差分析,数值实验验证了收敛性.第二部分考虑了扰动的Darcy-Stokes问题,即Brinkman模型,误差分析得到L2误差估计,数值实验验证了未知量受扰动因子影响下数值逼近的收敛率.在第四章中,给出了多孔介质中不可压缩混溶驱动的守恒特征线数值格式.运用[68]提出的守恒特征线格式,得到混溶驱动的守恒特征线技术(MCC),理论证明了数值格式的质量守恒性.利用外插技术,将速度方程和浓度方程解耦,在一些归纳假设的条件下给出离散格式的误差分析和最优的L2模误差估计.数值算例验证了该逼近格式的质量守恒性,未知量的收敛阶也是与理论一致的.最后给出了一个实际问题的模拟,验证了数值格式的有效性.在第五章中,讨论了多孔介质可压缩混溶驱动模型的特征线数值格式.该模型为较强的非线性耦合问题,利用投影算子首先给出未知量的初始值,然后基于一些归纳假设进行误差分析,对速度方程和浓度方程对应的误差方程分别进行估计,再整合到一起最终得到最优的误差估计.最后用数值实验验证了理论分析.
【作者】李新栋;
【导师】芮洪兴;
【作者基本信息】山东大学,计算数学,2014,博士
【关键词】混合元;特征线;质量守恒;交错网格;混溶驱动;误差估计;

【参考文献】
[1]刘金涛.薄膜太阳能电池板自动装箱设备的设计与分析[D].兰州理工大学,机械制造及其自动化,2013,硕士.
[2]贾德东.α-二亚胺溴化镍配合物合成及乙烯催化性能的研究[D].齐鲁工业大学,高分子化学与物理,2014,硕士.
[3]郑舟舟.论猥亵犯罪的立法完善[D].吉林大学,法律,2014,硕士.
[4]丁立,方攸同,蒋毅.高速列车牵引电机的同步虚拟主轴控制[J].机电工程,2014,10:1311-1315.
[5]王兴抗.“宣恩民俗文化”在初中泥塑教学中的应用与实践[D].华中师范大学,学科教学,2013,硕士.
[6]王维琴,盖玲,王剑平.高压脉冲电场预处理对甘薯干燥的影响[J].农业机械学报.2005(08)
[7]吴俊.基于32位浮点DSP的变压器保护的研制[D].浙江大学,电力系统及其自动化,2004,硕士.
[8]程民杰.棕榈油替代鱼油对半滑舌鳎生长、生理生化和肌肉营养品质影响的研究[D].天津农学院,水产养殖,2014,硕士.
[9]张倩倩.认知视角下的个体量词教学研究[D].吉林大学,语言学及应用语言学,2014,硕士.
[10]汤淼.合肥供水集团有限公司发展战略研究[D].安徽大学,工商管理,2014,硕士.
[11]晏梅珍.江西渝州科技职业学院学生作业管理系统的分析与设计[D].云南大学,软件工程,2012,硕士.
[12]张永龙,兰惠清,符泽第.埋地成品油管道泄漏强度模型与试验验证[J].安全与环境学报,2015,02:134-137.
[13]张春华,李迪,陈生隘.焊接车间管理中融合知识发现的探讨和实践[J].现代制造工程,2001,12:7-8.
[14]吕鹏辉,陆科,叶茂,成英,钟永恒.我国新能源汽车及其轻量化材料科技成果计量分析[J].现代情报,2011,07:93-97.
[15]程百花.基于不同材料摩擦副的旋转式能量回收装置性能研究[D].天津大学,化学工程,2013,硕士.
[16]郭素琴.河北省G县农村空巢家庭养老现状及社会工作介入[D].内蒙古师范大学,社会工作(专业学位),2013,硕士.
[17]闭俏丽.广西瑶族地区中小学教育特色发展的问题及对策研究[D].西南大学,教育学原理,2013,硕士.
[18]商金花.无线MIMO系统球检测算法的研究[D].南京邮电大学,电子与通信工程(专业学位),2013,硕士.
[19]许燕,邓自立.广义系统Wiener状态滤波新算法[J].控制与决策,2003,03:328-331.
[20]贺继胜.可再生能源微电网的建模与控制[D].广东工业大学,电气工程,2014,硕士.
[21]鲍艳辉.d波超导体隧道结中准粒子输运特性的理论研究[D].东北师范大学,2014.
[22]孙瑞玲.基于石墨烯复合材料及石墨烯量子点的制备并用于一些酚类环境污染物的分析[D].南昌大学,2014.
[23]于军香.桑螟寄生性天敌资源及其优势种群的生物学特性研究[D].苏州大学,特种经济动物饲养,2003,硕士.
[24]杨明,孙志挥,吉根林.快速挖掘全局频繁项目集[J].计算机研究与发展,2003,04:620-626.
[25]孙自金.对我国世居高原竞走运动员选材的研究[D].北京体育大学,体育教育训练学,2013,硕士.
[26]卢小佳.多普勒天气雷达信号处理的研究[D].安徽大学,通信与信息系统,2014,硕士.
[27]肖皖林.别墅群冷热电三联供系统配置与性能评价的研究[D].华中科技大学,供热、供燃气、通风及空调工程,2013,硕士.
[28]徐芃.原产国与制造国相分离对消费者购买意愿影响的实证研究[D].南京师范大学,企业管理学,2012,硕士.
[29]杨清鹏.基于中学生心理特点的作文教学策略研究[D].西北师范大学,教育,2013,硕士.
[30]王悦.2012总统大选研究及中法关系发展前景展望[D].西安外国语大学,法语语言文学,2013,硕士.
[31]张纬.芯片堆叠中散热分析方法研究[D].西安电子科技大学,计算机应用技术,2012,硕士.
[32]谢婷.中国中等职业教育经费保障机制研究[D].云南财经大学,财政学,2014,硕士.
[33]徐立敏.女性意识的嬗变[D].郑州大学,英语语言文学,2013,硕士.
[34]张萍.概率意义下量子系统可观性的研究[D].郑州大学,应用数学,2013,硕士.
[35]宋何煜.Timeless对肝癌细胞HepG_2增殖与凋亡影响的研究[D].北京交通大学,2014.
[36]马月萍.基于FPGA与DDS的磁共振成像射频脉冲发生器的研制[D].北京化工大学,控制工程(专业学位),2013,硕士.
[37]秦岚.轮状病毒NSP1抑制天然免疫机制及其对宿主蛋白调控的研究[D].北京协和医学院,病原生物学,2010,博士.
[38]李华.Smad2与Smad3参与TGF-β诱导的后发性白内障发生的研究[D].天津医科大学,眼科学,2014,博士.
[39]张俊俊,张振华,郭超,李杰,邓小清.正三角锯齿型石墨烯的电子输运特性[J].物理化学学报,2012,07:1701-1706.
[40]汪萌萌.农村高中生偶像崇拜现象研究[D].南京师范大学,教育学原理,2012,硕士.
[41]谭佳音.企业并购中支付方式及交易价格问题研究[D].沈阳工业大学,会计学,2003,硕士.
[42]刘森林.社会工作伦理视角下救助站对“职业跑站者”的救助管理问题研究[D].中国社会科学院研究生院,社会政策(专业学位),2014,硕士.
[43]边丽匣.含氮杂环化合物的电化学氧化与降解过程中废水的毒性变化[D].浙江工商大学,环境工程(专业学位),2014,硕士.
[44]王建.从《南方人物周刊》封面人物看公共人物形象的建构[D].兰州大学,新闻学,2013,硕士.
[45]袁泉,杨振坤,柴洁,王冬雁.早龄期H型钢混凝土梁受力性能试验研究及有限元分析[J].建筑结构学报,2014,03:193-200.
[46]程冉冉.WDM光网络中基于节能的分段保护算法研究[D].大连海事大学,信息与通信工程,2013,硕士.
[47]王华.Dobby CAD系统的研究与开发[D].浙江大学,控制理论与控制工程,2004,硕士.
[48]李铖铖.共轭亚油酸对肝细胞胰岛素抵抗的影响[D].大连工业大学,食品科学,2010,硕士.
[49]赵哲.基于关键链技术的EPC进度管理研究[D].河北工程大学,企业管理,2012,硕士.
[50]刘开志.纳米ZnO/微晶白云母复合材料的制备与表征[D].成都理工大学,矿物材料学,2013,硕士.

相关推荐
更多