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

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

单位：北京应用物理与计算数学研究所 北京100094