矩形区域两相渗流驱动问题有限体积元数值计算和分析-计算数学专业毕业论文.docxVIP

山东大学硬士学位论文审文祷要 山东大学硬士学位论文 审文祷要 本文采用有限体积元方法她麓多孔介质中两褶不可隘缩流俸的运翻+黟孔介质 中的两栩对汛瓣动问题是由一缄非线性偏微分方程来描述的，分为可混涔与不漫溶 两释情形。冀凝煎方法煮狠多，漱一般翁Gaierkin肖黻元法羚，Do锚秽托sf2】对其 运霪簿瑟畜袋麓分方洼，Russeli‘矧运露蒋荭煮陵元方法，麓零嚣最凌翁投藏静。 本文在上这王捧蕊基礁一曼，运麓有藩露强元方浚，裘短影篷壤主寒楚璎这个麓题。 对于可漏溶问题半离散的非弥激情形得到了三2的最优估计，全离散的弥散情形得 到了嚣1翡皴傀信诗，对予不溅溶瓣邀金盔鼗褥澎褥到了嚣1懿最霞傣计。舂蔽体 积元方法具有襁强的物理特性，计算简单，有利于进行数值运算【lo]。 文零懿绻濑效下： 鼙一章一胃浞溶不可基嚣静两糟渗滚驱凌阏瑟。 篥一节一孕；富。萼；入哥瀛漆不可羼臻鼹挺渗淡聪淤阀题静褪理鹜祭帮茭数学 模型的微分方程系统(，)。 ， ￡，，邮、 {f8》可，毪=-V r(誊笋翼V≯)=g 嚣∈轮、0t《T 溺{ 彝《 段划 {(喾番转)≥+嚣t翠s～V，《DVs》一肇～sj譬茹∈ll，0￡茎T 、 V§ 第=节～肖限体积元格式的戆立。本文对区域进行殴边形剐分和对偶剖分， 教斌掭黼羧变瓣鸯分片双一次黼数室游，蕺捡骏瑟数滋藤为努片常数函数枣瘸，缮 剜了压力方橼与饱和度方程的半离散及金离散有限体积元格式： (2．3) 嚣(；臻，镪》=溆‰j 敦‰芒该 (2．8)(≯警，坝)+(rAVXh，％)十(廖(誊，DV，V‰)+(《，弼)=《鸭‰) (2．9)3h(x，0)=岛 (2．Io}(毋威碟，％)+e嚣酽V霹，冁)+fD(EU”)V霹，V蝴)+(譬嚣，珊j一向哥，‰j (2．11)鲤*岛A 第三节～灏备知识和引理。本文在矩形嗣格剖分上证明了有限元掇式中阿积 毡逮0，双袋浚澎式嚣知；；，嚣转稷8A汐^写专戆黯黎淫茨强鬟慈。莠葶}入了美 山东大学硕士学位论文力函数P和饱和度函数s在试探函数空间上的椭圆投影、插值多项式及其相关的性 山东大学硕士学位论文 力函数P和饱和度函数s在试探函数空间上的椭圆投影、插值多项式及其相关的性 质。 第四节 误差估计。就矩形网格情形建立收敛阶的估计，分析了不含弥散项 的半离散格式和含有弥散项的全离散格式的收敛性，分别得到了L2和H1的最优 估计。 利用 (4．3) 11『7《-S K(h2+【1s～s^1{o)≤Ⅳ(^2+11tlu) 分别可以得到 半离散格式的误差估计： 定理』』设D=D(。．，)．s和＆分别为问题(，)和半离散有限体积元格式 (2．8)(2．9)的解，则存在常数K，当h充分小的时候，有： (4 17) 忪一S^f|o茎Kh2 全离散格式的误差估计： 定理J 2设D=D(x，“)，(p：s)和(R，鼠)分别为问题(，)和全离散有限体积 元格式(2 3)．(2．10)．(2．11)的解，(P，S)∈L。(日3)nHl(日3)n日2(日1)，空间和时 间步长满足： At。=o(h)，At：=o((At。)；)，Atp=o((At。)÷)那么存在常数K， 当h充分小的时候，有： (4 48) sup iIp～BPlll+sup IIs一．霸||1≤K(h+At。) ”2 n 第五节 数值算例分析。我们对区域进行正方形的网格剖分，利用有限体积 元方法求解一个初边值问题。根据这个算例，我们发现对区域进行矩形剖分可以简 化计算，并且所得出的精确的数值结果进一步验证了有限体积元的优点，既保持有 限差分方法的计算简单性，又兼有有限元方法的精确性。 第二章一不混溶不可压缩的两相渗流驱动问题。 第一节一引言和有限体积元格式的建立。引入不混溶不可压缩两相渗流驱动 问题的物理背景和其数学模型的微分方程系统(』J)。 l(n)V·札=一V·(≈(z．s)Vp)=q z∈Q，0t≤T 、 ’l(b) 击(z)篆+6(s)u Vs—V-(D(s)Vs)=，(z、s)z∈Q，0t≤T 山东大学硕士学位论文 山东大学硕士学位论文 对应于第一章的试探函数空间和检验函数空间，得到了压力方程和饱和度方程 的全离散有限体积元格式： (1．4) B(Sh；R，Yh)=(g，Vh) Vv／。∈K (1．5)(妒dc簖，Vh)+(6(器。)Eu“V岛，Vh)+(D(器～)vsZ，Vv^)=(f(z，霹一1)，Vh) (1．6)s2=Sob 第二节一误差估计。分析了全离散格式，得到了H1模的最优估计． 定理2』设D=D(x，s)，(P，5)和(R，Sh)分别为问题(II)和(1．4)，(1．5)，(1．6) 的解，(p，s)∈LOO(日3)r]H1(日3)nH2(日1)，空间和时间步长满足：At。=o(^)，△￡：=D((△￡。)；