间断Galerkin有限元的自适应方法研究-科学计算及其应用软件专业论文.docxVIP

摘要 摘要 本文给出求解类双曲型偏微分方程的间断Galerkin有限元方法 (Dis∞ntinuous Oalerkin Method，筵嚣DG方滚)及冀带单嚣羲努类技 术的蠢适应算法蜜瑷．大量酶数溆算倒表嚷，本文给出的宫潦应算法不 仅继承了DG方法，高度并行性、一致高阶精度的优点；而且柱数值效 果和计算效率等方孵均有明显改游． 关镳溺：鬻蘩考歉嚣，骞适应策臻。类荻夔蓠题 AbstractIn Abstract In this paper，we present the Discontinuous Galerkiu(Da)lVlethod sad its Adaptive Strategy embedded with element pre—sorted technique for similar- hyperbolic systems．Preliminary numerical examples displaying the perfor。。 赫l拽nee《the method搬e shown如confirm that it not。nly shares with DG methods their high paraLtelizability sad h遮h-order formal＆#e毡t8e黔but also improves the numerical perfonna_1ace and computational efficiency,etc． Key words：discontinuous finite element，adaptive strategy,similar- hyperbohc systems §1问题敲背景类双魏掇镶镘分方程《双簦系统或射流占缆扩教系缓)静数值求瓣一壹是计算数 §1问题敲背景 类双魏掇镶镘分方程《双簦系统或射流占缆扩教系缓)静数值求瓣一壹是计算数 学关注懿熏辫谦透，医为阏繇酶囊蘩常常存在大梯度蜜纯(作为投隈情形瓣缝双馥闻 题，解可能疑有间断)等率富和复杂的结构．它们常常可写成如下的铮恒形式 魂u+V·F(u，Du)=H(u)， 其中U∈n×{o，功 (1) 其中纯∈躲，U=(u1，珏2，．．．，‰)|，Du一(魄，％，％n， ，‰“囊该类瓣趱豹数蓬禳掇 需要细密的f阉格和高阶的邋近格式，使得许多传统的数值方法遇到挑战．为此，许多 新型的数假方法在近二十年来被相继掇出；这方面的正作集中在数假格式的改进和 两轻壹适墩技术鳕应用。 驽辩蠢隈元方法是黥孵类驳夔藤疆瓣一耱有效豹方法，它在萁搭瓣处理主彰瑟 多种方式，如Arbogast和Wheeler[i]，Baker etⅡ}[3】，Cockburn和Shu[20l等．本文仅 仅关注Cockburn和Shuf20]所提出(局部)间断Galerkin有限元方法(简称LDG方 法)，它是上个世纪90年代初由求解双艟阊匦的Runge-Kutta间断镌限元方法(简 称RKDG方法)懿基璃上发爱嚣寒魏．獒持点是灸洚试舞丞数穰捡簸返数在荤灵迭 界存在间断，在数值实现上较连续有隈既具有更大盼荧活性，具赢懿好的数值稳定 性和高度的可并行性，并臆能够精准地横拟解在间断处(或大梯度燕化区域)的“嗽 落”，是一种求解类双曲问题的高效数德方法。为叙述方便，本文将上述LDG方法硼 RKDG方淡缝舔秀DG努法。 由于畿双曲问题解酶结构常+努率富和复杂，单纯邋采用诗簿随格酶垒藩糯镪 或有限元阶数的整体提升，并不能提供满意的数值模拟结果．为此，根据解在不同区 域中的性淼，采用相应的网格自适应策略来局部调熬(加密或稀疏)网格和局部改搬 f提高竣降豫)逼近函数懋瓣的塞壹度。鑫适应技本主辫分鸯翔下三秘譬，2域： 1，h一方法(局部热密方法)；依靠慕穗I特定的误蕊旅示器局部逢魏密承l稀疏阚穰。 在解的嶷化相对平缀的区域稀疏网格，而在变化剧烈的区域加密阿格以提高数德 效果； 2．p一方法：动态建选撵摆每个单元凑部的有限元窘阉的黔数．在解糖对光滑的聪 壤蠹，露骞效逢将效缓猿废挺舞裂掰滔鲮{骜精度蔽密数毂袭； 3．r一方法(移动网格汝)：通过网格结点的重新定位来调整计算网格的性态，避而 提高数值效果．其优点是不增加网格结点的个数，但计算十分笈杂． 本文将考虑DG有限霓方法的宜遗成数值实现，将DG方法与网橙自适应策略有 撬建缝合熬来．这是求熬类双基|霹疆一耱缀塞然懿逡撵。滔善这夸豢零恿籍，大薰瓣 工作在数德餐验和理论分析方面积累了率厚的成果．谯1974年，LesaInt和Ravi喊【3碓 针对DGFEM方法(1973年由Reed和Hill[42】提出解决中子传输问题)给出了一种 l ^一自适应算法的L2(n)次优收敛估计．厝来，C．Johnson ^一自适应算法的L2(n)次优收敛估计．