田瑞雪
,
葛永斌
,
吴文权
工程热物理学报
通过对二维不可压缩N-S方程的涡量-流函数方程组消去涡量而得到仅以流函数为求解变量的控制方程,从而使不可压N-S方程的求解个数减到最少.求解方法采用本文提出的二阶精度的九节点紧致差分格式,因此无须对靠近边界的网格点作特殊处理.为了加快迭代收敛速度,采用多重网格加速技术.数值实验结果验证了方法的精确性和可靠性.
关键词:
不可压N-S方程
,
流函数方法
,
多重网格
田瑞雪
,
葛永斌
,
吴文权
工程热物理学报
利用一阶偏导数项的四阶紧致差分算子,直接推导出了数值求解二维对流反应方程的一种新的高精度紧致差分格式.为了提高差分方程的求解效率,采用多重网格加速技术,建立了与之相适应的多重网格V循环算法.数值实验结果验证了本文方法的精确性和可靠性.
关键词:
对流反应方程
,
紧致差分格式
,
高精度
,
多重网格