计算物理 ›› 2015, Vol. 32 ›› Issue (5): 505-513.

• 论文 •    下一篇

分片光滑物态方程的能量方程非线性迭代解法

杭旭登1, 李双贵1, 杨容1, 袁光伟1,2   

  1. 1. 北京应用物理与计算数学研究所, 北京市海淀区丰豪东路2号 100094;
    2. 计算物理实验室, 北京市海淀区花园路6号 100088
  • 收稿日期:2014-08-04 修回日期:2014-11-29 出版日期:2015-09-25 发布日期:2015-09-25
  • 通讯作者: 袁光伟,E-mail:yuan_guangwei@iapcm.ac.cn
  • 作者简介:杭旭登(1976-),男,江苏无锡,博士,研究员,从事辐射输运数值模拟研究,E-mail:hang_xudeng@iapcm.ac.cn
  • 基金资助:
    国家自然科学基金(11171036,11001023);中物院科学技术发展基金(2014A0202009)资助项目

A Nonlinear Iterative Method for Energy Equations with Piecewise Smooth EOS

HANG Xudeng1, LI Shuanggui1, YANG Rong1, YUAN Guangwei1,2   

  1. 1. Institute of Applied Physics and Computational Physics, No.2 Fenghao Donglu, Haidian, Beijing 100094, China;
    2. Laboratory of Computational Physics, No.6 Huayuan Road, Haidian, Beijing 100088, China
  • Received:2014-08-04 Revised:2014-11-29 Online:2015-09-25 Published:2015-09-25

摘要: 实际应用中的物态方程由分片光滑曲面拼接而成,拼接处存在间断.隐式求解相应的能量方程时,经常出现迭代收敛慢的情况和非物理解.本文通过构造对应的新的非线性问题,提出一种非线性迭代算法.该算法适用于求解有间断的分片光滑物态方程的非线性能量方程,其中引入一个度量能量变化的参数用于自动判断跳段是否发生,在求解时无需事先知道物态方程间断的位置,且能精确计算物态方程间断带来的能量盈亏,用于评估物态方程间断对能量的影响.典型算例验证了新算法具有稳定的收敛性,并给出符合物理规律的解.

关键词: 间断状态方程, 能量方程, 非线性迭代算法

Abstract: In practical applications, equation of states (EOS) consists of several piecewise smooth surfaces, which leads to discontinuity at interface. As a traditional nonlinear iterative algorithm is applied to an energy equation with discontinuous EOS, it may lead to slow convergence and unphysical solutions. To overcome the difficulties, a nonlinear problem is designed, and a nonlinear iterative algorithm is proposed to solve the problem. The algorithm is fit for energy equations with discontinuous EOS of piecewise smooth functions. A parameter of energy change is defined in the algorithm so that it is unnecessary to know discontinuity position in advance. The algorithm calculates precisely net gain or leakage of energy, which can be used to assess influence of discontinuity in EOS. Typical numerical experiments verify that the algorithm converges stably, and gives physical solutions.

Key words: discontinuous EOS, energy equation, nonlinear iterative algorithm

中图分类号: