混合物状态方程的计算?

摘要：

多介质流体动力学过程的数值模拟往往涉及混合物状态方程的计算。做图法和Newton法是混合物状态方程计算常采用的方法，前者虽直观精度却差，后者计算效率高却只具有局部收敛性，当解与其初始猜测值相差较远时Newton法不一定能够获得收敛解。为此，本文给出一种具有大范围收敛性的嵌入算法(imbedding method)求解混合物状态方程，其基本思想是通过引入嵌入参数，将待解的混合物状态方程和易解的混合物状态方程线性组合，构成嵌入方程组，当嵌入参数从0连续地变化到1时，嵌入方程组的解由易解的混合物状态方程的解连续地变化为待解的混合物状态方程的解。嵌入方程组可由Newton法迭代求解，也可转化为以嵌入参数为自变量的常微分方程组，从而易于由成熟的计算方法如梯形法等进行求解。进一步利用热力学基本关系， Maxwell形式的微分方程描述了压力和温度随嵌入参数的演化速率与应变速率和组分质量分数演化速率的关系。对铅锡混合物热力学量的计算表明了本文算法的有效性。

关键词：混合物、状态方程、迭代法、Maxwell方程

分类号：O1 ;TP3

资助基金：中国工程物理研究院科学技术发展基金2013A0201010;国家自然科学基金批准号:11272064资助的课题.* Project supported by the Science and Technology Development Foundation of China Academy of Engineering PhysicsGrant 2013A0201010;the National Natural Science Foundation of ChinaGrant 11272064

在线出版日期：2015-04-09（万方平台首次上网日期，不代表论文的发表时间）

页数：共1页

页码：064702-1-064702-6