流体力学方程的间断有限元方法 - 计算物理.pdfVIP

计 算 物 理 第22卷第2期 V01．22．No．2 CHINESE OFCOMPUTA7nONALPHYSICS 2005年3月 JOURNAL MⅡ．．2005 [文章编号】 1001．246x(2005)02一O108．09 流体力学方程的间断有限元方法 蔚喜军1， 周 铁2 (1．北京应用物理与计算数学研究所计算物理实验室，北京100088；2．北京大学数学系，北京100871) [摘 要] 在二维区域三角形网格上应用一阶、二阶和三阶精度间断有限元方法，对流体力学方程和方程组进 行了数值模拟．计算结果与差分方法计算结果比较，认为间断有限元方法在求解复杂边界条件和区域问题上有 一定的优势． [关键词】 三角形网格；间断有限元方法；前台阶问题；双马赫反射问题 [中图分类号】 035；0241．82 [文献标识码】 A O 引言 间断有限元方法是1973年由Reed和硒ll[13首先提出，并应用于求解中子输运方程，但这种方法长期以 给出了部分收敛性理论证明后b。1【，这一方法才引起人们的注意，并逐渐开始应用于流体力学计算领 域㈨2’131． 间断有限元方法应用于流体力学计算领域，并受到人们的关注，主要因其保持了通常有限元方法的优点： 1)能够处理复杂的区域边界和具有复杂的边界条件问题，并获得与区域内部一致的计算精度；易于网格 5|． 加密和高精度处理边界条件，实现自适应计算u4’1 2)可以得到任意阶精度的格式，同时又具有很好的局部紧致性，构造的高阶格式不需要非常宽的模板 (stencil)(实际上这比较符合守恒律方程的解的性质，一般来说，模板太宽，会抹平激波)． 同时由于间断有限元方法吸收了差分方法的一些特点，能够显式求解．因此 3)容易实现并行算法¨5|． 4)具有很好的稳定性，满足z2稳定性和熵相容性． 间断有限元方法的缺点是程序设计比较复杂和计算量大．然而近十年来，随着计算机条件不断的改善， 各种大型向量计算机和并行计算机的相继问世，间断有限元方法能够较为容易地把求解二维问题的方法推 广到三维问题，实现自适应算法和并行算法．对于原来应用差分方法求解的流体力学问题，由于定解的区域 不规则或者边界条件难以处理，逐步开始应用间断有限元方法进行求解．目前间断有限元方法除应用于守恒 8’19]． 划中，把这一方法应用于求解辐射流体力学问题等n 本文在二维区域三角形网格上应用一阶、二阶和三阶间断有限元方法，对于Euler可压缩流体力学方程 和方程组进行了数值模拟，并与差分方法计算结果比较，认为间断有限元方法在求解复杂边界条件和复杂区 域问题上有一定的优势． 1间断有限元方法 考虑二维非定常流体力学方程组 以+(，(u))。+(g(U))，=0， [收稿日期】2003一ll一26 [基金项目]863高技术惯性约束聚变主题、国家自然科学基金(1047lOll)和中物院基金资助项目 [作者简介]蔚喜军(1959一)，男，内蒙古，研究员，博士，从事流体力学计算方法的研究． 第2期 蔚喜军等：流体力学方程的问断有限元方法 其中 u=(1D，lDn，lD￡)，E)’，，(U)=(ID“，IDu2+p，lD“秽，u(E+p))7， ID是密度，u，秽是戈，y方向的速度分量，p是压力；对于气体动力学方程 E=南+吉lD(“2+勘2) 是总能，y=1．4是绝热比．状态方程是p=(y一1)俨，这里e是内能． 间断有限元方法哺3求解方程组(1)，首先将求解区域力剖分成互不重叠的小区域0i：n=U0；，n；可以 是任意形状，如四边形、三角形网格等，这里我们采用三角形网格．假定ni和q(i≠．『)相交，则它