基于预处理的不可压缩N-S方程拟压缩数值算法研究  

作　　者：钱战森[1] 张劲柏[1] 李椿萱[1] 

机构地区：[1]北京航空航天大学航空科学与工程学院,北京100191

出　　处：《中国科学：物理学、力学、天文学》2010年第12期1542-1554,共13页Scientia Sinica Physica,Mechanica & Astronomica

基　　金：国家自然科学基金资助项目(批准号:50836006)

摘　　要：通过研究求解不可压缩Navier-Stokes方程的拟压缩方法与加速刚性双曲型方程时间推进的预处理技术,推导了一般曲线坐标系下带预处理的拟压缩Navier-Stokes方程特征系统,并结合有限差分法,建立了适用于不可压缩黏性流动计算的拟压缩快速算法.通过对不可压缩无黏圆柱绕流、平板层流流动、低Reynolds数定常圆柱绕流问题的数值模拟研究,得到了与相关理论与实验测试相吻合的结果,验证了所建立数值方法的快速可靠性.较系统地研究了预处理引入的参数和拟压缩因子的选择对收敛特性的影响.结果表明,Roe格式相对于二阶中心差分格式得到的结果更令人满意;对拟压缩Navier-Stokes方程进行预处理能有效提高数值计算的收敛速度;自适应的拟压缩因子取值能在很大程度上改善数值解的收敛特性,且不需要根据具体流动问题进行人工调节.最后将本文发展的数值方法用于低Reynolds数非定常圆柱绕流的数值模拟,所得结果亦和实验观测结果及其他文献的计算吻合很好.

关 键 词：拟压缩方法 预处理 迎风格式 不可压缩流动 Navier—Stokes方程 

分 类 号：O35[理学—流体力学]