两相可压缩溶混流驱动问题的特征有限元方式.pdf

二维可压溶混流驱动问题的特征有限元方法 妙356 795 致谢 从我一九九三年进入山东涝：学学习到【今j高酗鞠贡±．毕业邂汶已经1i七 个年头了．七年中我不仅在山东大学浓厚的学习氛圈中学习到文化知识， 学习妻瞻津澈学问，而且学习劐良P阿做^，如伺．处世．在即柏锆毋罔l芟， 告射辛i盼嘲朐瑚舛蛹翮亲如兄期篚l妹的同学之即，请允{焖璐踣一患谢 谢l 首先，捌酗嘲我的导J幌暖搁喇夷授、粱睛掀．程既斓n彰峻韵 划髓嘲l和严i村瞩崛本文得以蒯a帙撼燃捌财喔蒯勖g启迪使 我获益良多． 同时，我翘廖琳粕蔷【让先生．我正唇盎抛的}寸论班上开烟较深入的了 解专碰呐内涵，学习重蝴入豫树嘲啪硼嗍喇湘治学彷法· 其次，剩幽嘲制乎平苦师。王申株封民置啪期吼猢敏薨狮， 鲁统超老师，孙文涛翱雨辱各位老师和同学，是他们在生活和学习上给我 提供了诸多的关明讲明日． 籀l页 二维可压溶混流驱动问题的特征有限元方法 中文摘要 本文共分两苹。 第’用特，征有限元方法处理两相完全可酗溶混流马嘞问题，j跫j：治理选择 误差发展不等式。得到了解的最优P0；Ⅳ1(Q)№甑是差估计。共分三节。 第一节是引言。介绍了问题及-前人-_的研究成／果，简述本文与前人工作的不同之 处及所得到的结论。 ／‘每虑平面E有界域内两相可压缩可涸溶驱动问题： ∥(c)罂+V．“+d(c-．即一gDVc．即=g “．而=D．Vc．元=0G，f)∈(蓖】×J) pb，o)=风G)cG，o)=C00) 其中Q为平面馏2上单位厚度有界域，边界为鼢，亓为其边界外法向．时间域 J=lo，丁j·／ 第二节给出了方程的特征有限元全离散隘武及椭圆投影。 ／／首先引入一些记号，给出对系数的假定和正则性条件。构造时空有限元空间 埘。，％，引入特征方向委=古(妒昙+e“·V)，然后给出了方程的变分形式： Vv，z∈H1Q)，求0，c)满足 引入全离散特征有限元格式，． 第3页 二维可压溶混流驱动问题的特征有限元方法 Plh采用格式： 己知∽～，硝一’)E％×tlh，求，=广时，∽，硝)专％X ∞(c：。_，肼，v)+(口仅。如：，Vv)一pG≯b(c：“)蹦～．耳≯：～，v) Vv∈甩。 一(g协vcf‘．1蹦～，v)=0”‘，v) “：=一口(c≯扔： (≯!争n--n-I，z]+慨乳：，耽)+∞G∥p∥一+G妊％‘即“力 VzE％ =b”’p∥一％n-’)z) 最J舌引入椭圆船：亭：J一／lvl^；歹：J斗％满足： ㈨妒0一声l≮～)+(p(c_一gore)·V0一歹lV)+五0一芦，V)=0 (DV(c一巩耽)+∥(c一手，z)=0忱∈％，fly∈nh，j∈J 式中Ⅳ，^为足够大常数，以保证左端双线性形式的强制椭圆性／ ／ 第昙节进衔受差f占计a ／首先对初值进行估计，在此基础t进行误差估计t棚＆lh肌不等式得到 r∽H1(Q))唰自膨砉果： 定理：削分参数满足如下条件： j≥2，七≥2，△，=口G，l△f=o(趣X^；=D魄X艇=D协，) 全离散格式(2．2)●}茸(1l埔啊擞间的误差满足： 苫f帜。一c。r《2+l|％幻一所圳：)△，≤M江2+^≯+酵) M是不依赖于址，^的常数，估计对rU；H’Q)扣颧绸毗的，故结果是完