二维非定常不可压涡量速度NS方程组的高精度紧致差分方法.pdfVIP

中国工程热物理学会 传热传质学 学术论文会议 编号：083290 二维非定常不可压涡量-速度 N-S 方程组 的高精度紧致差分方法 1 1 2 杨国锋 葛永斌 田振夫 (1.宁夏大学应用数学与力学研究所, 银川,750021；2.复旦大学力学与工程科学系,上海,200433 ) (Email:ygf0411@163.com Tel 摘要：提出了数值求解二维非定常不可压涡量-速度变量 N-S 方程组的高精度隐式紧致差分格式. 其空间为四阶精度，时间为二阶精度，并且是无条件稳定的.为了验证本文方法的精确性和可靠性进行 了数值实验，数值实验结果与精确解或文献中的结果吻合的非常好. 关键词：非定常不可压 N-S 方程组；涡量-速度方法；高精度紧致差分格式； 有限差分方法 0.前言 对不可压 Navier-Stokes(N-S)方程组数值求解方法的研究一直是计算流体力学领域 一个重要的研究方向. 由于采用原始变量的控制方程存在强的压力和速度的耦合问题,尽 管也发展起来了较为成熟和完善的 SIMPLE 类等算法解决这一问题,但这些方法往往很 难被推广到高精度.为此一些学者提出了非原始变量的方法,如涡量-流函数方法和涡量- 速度方法等. 目前对涡量-流函数方法的研究比较多,如文献[1-6],其中高精度方法也有不 少报道[2-4,6].但对于涡量-速度方法的研究还不深入, 已发展起来的数值格式均是低阶精度 的[7-10] .为此,本文提出了二维非定常不可压涡量-速度变量 N-S 方程组的高精度紧致全隐 格式.其空间为四阶精度，时间为二阶精度，并且是无条件稳定的.对于三个有精确解问 题进行了数值实验，数值实验结果与精确解或文献中的结果吻合的非常的好，从而验证 了本文方法的精确性和可靠性. 1.控制方程 二维非定常不可压 N-S 方程组的涡量-速度形式可以表示为 ξ +uξ +vξ (ξ =+ξ ) / Re （1） t x y xx yy uxx +uyy =−ξy （2 ） v +v =−ξ （3 ） xx yy x 其中, ξ表示涡量, u和v 分别为x 方向和y 方向的速度分量, Re 为雷诺数. 2.数值方法 2.1 涡量方程的四阶紧致差分格式 对二维非定常的对流扩散方程方程 Φ −a(Φ +Φ ) +p(x, y,t )Φ +q(x, y,t)Φ f (x, y,t) t xx yy x y ————————————— 基金项目：国家自然科学基金资助项目10662006） 文献 [13] 给出了空间四阶 , 时间二阶精度的隐式紧致差分格式. 令 a 1/ Re ， φ ξ, p u, q v, f 0 ,即可得到涡量方程(1)的四阶紧致差分格式为： 8 4 4 n+1 n n−1 A ξ B ξ =+ C ξ