一种求解化学非平衡流动化学反应源项的预估修正方法

摘要

本发明涉及流体力学数值计算技术领域，公开了一种求解化学非平衡流动化学反应源项的预估修正方法，包括以下步骤：a.根据当前时刻的流场解获得当前时刻的组分密度和化学反应生成率，预设该计算轮次N＝1；b.开展下一轮预估，通过二分法剖分时间离散间隔；依次对二分点的参数值进行预估，并对上轮各点的参数值进行修正；c.计算第N+1轮和第N轮的预估值之差，若相对值较大，则增加一轮预估修正，即N＝N+1，进行步骤b操作，否则跳转至下一步骤；d.将第N+1轮的计算结果作为最终结果。本方法采用核心控制方程隐式处理和化学反应源项半隐式处理，在保持原有LUSGS算法的优势的同时，解决了化学反应源项的刚性问题。