王国玉
,
张永学
,
曹树良
工程热物理学报
为了分析射流放水阀内部的非定常不可压粘性流场,在一般曲线坐标系下求解了以逆变速度为参数的不可压N-S方程和压力椭圆方程。计算采用了基于隐式SMAC格式的有限差分格式。在这种格式中,非定常流场的速度是通过使用Crank-Nicholson格式,并在每一时间步处以牛顿迭代的方法得到的;而压力椭圆方程是通过使用Tschebyscheff SLOR格式并交替改变计算方向而求解。通过使用交错网格和迎风差分,诸如改进的QUICK格式,计算过程中的误差和数值不稳定能够得到有效地控制。计算结果和实验观察得到的流动图谱对比表明,二者十分吻合。
关键词:
数值计算
,
不可压N-S方程
,
非稳定流动SMAC格式
,
射流放水阀
王辉
,
徐明海
工程热物理学报
采用非结构化网格有限容积法求解了不可压N-S方程组,对流项采用GAMMA格式,扩散项采用二阶中心差分格式建立离散方程,用SOAR算法处理压力与速度的耦合关系,得到了一种求解不可压N-S方程的非结构网格耦合求解器。通过方腔顶盖驱动流、后台阶绕流以及方腔自然对流等几个典型的算例,考察了求解器的计算精度及收敛特性,并与SIMPLE算法进行了比较,结果表明该求解器是有效可行的。
关键词:
非结构化网格
,
不可压N-S方程
,
耦合求解器
,
有限容积法
田瑞雪
,
葛永斌
,
吴文权
工程热物理学报
通过对二维不可压缩N-S方程的涡量-流函数方程组消去涡量而得到仅以流函数为求解变量的控制方程,从而使不可压N-S方程的求解个数减到最少.求解方法采用本文提出的二阶精度的九节点紧致差分格式,因此无须对靠近边界的网格点作特殊处理.为了加快迭代收敛速度,采用多重网格加速技术.数值实验结果验证了方法的精确性和可靠性.
关键词:
不可压N-S方程
,
流函数方法
,
多重网格