多组份流动质量分数保极值原理算法.pdfVIP

第 31卷 第3期 计 算 物 理 Vo1．31．NO．3 2014年 5月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS May．2014 文章编号 ：1001—246X(2014)03—0292—15 多组份流动质量分数保极值原理算法 唐维军 ， 蒋 浪 ， 程军波 (1．北京应用物理与计算数学研究所 ，计算物理实验室 ，北京 100088；2．中国工程物理研究院研究生部 ，北京 100088) 摘 要 ：对基于质量分数的Mie—Gruneisen状态方程多流体组份模型提 出了新 的数值方法 ．该模型保持混合流体的 质量 、动量、和能量守恒 ，保持各组份分质量守恒 ，在多流体组份界面处保持压力和速度一致 ．该模型是拟守恒型方 程系统．对该模型系统的离散采用波传播算法．与直接对模型中所有守恒方程采用相 同算法不同的是，在处理分介 质质量守恒方程时，对波传播算法进行了修正，使之满足质量分数保极值原理 ．而不作修改的算法则不能保证质量 分数在 [0，1]范围．数值实验验证了该方法有效． 关键词 ：Mie—Gruneisen状态方程 ；质量分数保极值原理 ；波传播算法 中图分类号 ：0241．82 文献标志码 ：A O 引言 可压缩多介质复杂流动有很广泛的应用 ，包括爆轰 ，燃烧 ，惯性约束聚变 ，天体物理等诸多领域．研究这 些复杂流动的目的是要了解多介质流动的流体动力学机制．这涉及到多介质流动的数学建模及相关模型的 数值离散方法．一般而言，如不考虑粘性 、表面张力 、热力学效应等物理 因素，单纯使用欧拉方程描述单介质 流体无粘流动时，其数学模型由质量 、动量、能量的三个守恒律组成的完全守恒系统构成 ，且该完全守恒系统 可用近三十年来发展成熟 的各种守恒格式进行数值离散．但如流动涉及多介质时，仅用完全守恒模型来描述 却难 以实现．在一些非常特殊的情形 ，有些作者尝试用完全守恒系统来描述多介质流动 ，但涉及混合流 动时，相应的数值格式处理则变得相对复杂 ；并且在计算纯界面问题时 ，界面附近压力或速度的数值解可能 会出现非物理振荡．目前 比较公认的是采用拟守恒模型来描述多介质复杂流动．例如对理想气体状态方程一 维两流体组份流动 ，利用对纯接触问题界面压力和速度数值解要求保持一致这种思想 ，Abgrall 最早揭示模 型可由混合介质质量 、动量和能量守恒方程 ，以及与混合状态方程参数有关 的量的一个对流方程来描述．并 且 Abgrall 还确认该对流方程的离散格式不是独立 的，而与三个守恒方程的离散格式有关．这个事实显示 多介质流动的数学模型和数值格式与单介质流动情形相比要复杂．随后 Abgrall 思想被推广到刚性气体状 态方程 。 、vanderWaals气体状态方程 、密实介质状态方程 、以及 Mie—Gruneisen状态方程 ．此 外 ，还有众多的其它模型来描述多介质流动，如加入 levelset方程的扩展欧拉方程组模型 13]，Karni的原始变 量模型 和杂交模型 ¨ ，Allaire等的5方程模型 ，Saurel和 Abgrall的两相流7方程模型 ，使用体积分 数的4方程模型 ¨。 ，反扩散界面模型等 ，两相流 BGK模型 ．至于有关上述模型的数值格式 ，散见 上述各文及相关引文，此处不一一列举． 讨论 Mie—Gruneisen状态方程一维模型的数值离散格式．该模型基于 Abgrall 等效方程思想 ，且分介质 质量保持守恒的模型 的数值离散格式．与基于体积分数对流方程模型 相比，本文所用模型 不仅能保 持混合流体的质量、动量、和能量守恒 ，还保持各组份分质量守恒 ，且对纯接触 问题保持界面压力速度一致． 本文使用波传播算法 ’’离散系统 ，但不是象文[12]对所有守恒方程进行相同数值离散 ，而是在对分介质 质量守恒方程离散时，采用 Larrouturou…对完全守恒系统设计 的质量分数保极值原理算法 ，对波传播算法进 行修改．即对分介质质量守恒方程的数值通量进行修改，同时引进一个类似 CFL条件的时间步长约束．其 目 收稿 日期：2013—06—17；修 回 日期 ：2013—10