Mie-Griineisen状态方程下多介质守恒型欧拉方程组的数值计算.pdf

计 算 物 理 第28卷第6期 V01．28。No．6 2011年11月 CHINESEJOURNALOFCOMPUTATl0NALPHYSICS Nov．，2011 文章编号：1001-246X(2011)06-0803-07 方程组的数值计算 吴宗铎， 宗 智 (大连理工大学船舶3-程学院，工业装备结构分析国家重点实验室．辽宁大连116024) 程中，使原方程结构得到简化．同时引入质量分数以表示流场中不同介质的变化，并对通量作近似处理，再利 结果表明，通过本方法既能使方程组满足守恒，又能较好地处理多种介质形成的间断面． 关键词：多介质；质量分数；Mie-Grflneisen；守恒 中图分类号：035 文献标识码：A O 引言 近年来，多介质的可压缩流动中的强间断面已成为大家关注的一个焦点，很多工程应用如爆炸、侵彻、冲 击等都涉及到强间断面的处理．强间断面给数值计算带来了极大的麻烦，如果介质的状态方程的结构形式比 Shyue…通过在Euler方程组中引入一定数量以体积分数作变量的非守恒偏微分方程，将方程组转化为准守 程的基础上，再对变量加入限制器，完成了一些复杂界面的模拟．Lagouti爸re¨1在方程组里同时加入变量为质 量分数和体积分数的偏微分方程，使之能处理一般的状态方程．但体积分数的引入虽然能准确捕捉界面位 置，但在界面方程为非守恒形式．如果要使方程组满足全守恒，质量分数的效果要强于体积分数． 使得它能比较方便地处理材料的初始参数带来影响，同时在处理高强度压缩、碰撞、爆轰等问题上也能取得 不错的效果．但其结构复杂，因此使计算也变得复杂． 分数，对流场中的不同介质构建质量守恒方程，使整个方程组满足守恒．再通过对变量做简化处理，将Mie． Grttneisen方程中的复杂参数引入到变量中，使方程的形式变得简单易于计算．这样，不同介质的界面处，流 场的状态转化为多种介质以不同质量分数共存的状态，消除了强间断面带来的数值振荡，较好地处理了 Mie．Grflneisen状态方程下多介质的强间断面． 1控制方程及计算方法 1．1控制方程 以二维直角坐标系为例，守恒型欧拉方程可以写为 收稿日期：2010—12—29；修回日期：2011一∞一02 基金项目：创新研究群体科学基金和973项目(2010CB832700)资助项目 作者简介：吴宗铎(1984一)，男，江西，博士生，从事计算流体力学激波研究，E-mail：wuzengduoO@yahoo．corn．cn 万方数据 804 计 算 物 理 第28卷 型at+鲁+鲁．o’a髫 a_y (1) 其中U为守恒变量，F，，F，分别为x,y方向的通量．表达形式为 pⅡ 1 移 2 1 m U= F，= 掣+p}， F，= + puv I p ∥眇“¨+ 、●●●●●●●Lr●●●●J (pE+p)nJ D p” 和P，E建立关系