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

中国科学: 物理学 力学 天文学 2010 年 第40 卷 第12 期: 1542 ~ 1554 SCIENTIA SINICA Phys, Mech Astron SCIENCE CHINA PRESS 论 文 基于预处理的不可压缩N-S 方程拟压缩数值算法 研究 * 钱战森 , 张劲柏, 李椿萱 北京航空航天大学航空科学与工程学院, 北京 100191 *E-mail: qzhsen55130@ 收稿日期: 2010-05-30; 接受日期: 2010-09-01 国家自然科学基金资助项目(批准号: 摘要 通过研究求解不可压缩Navier-Stokes 方程的拟压缩方法与加速刚性双曲型方程时间推进的预处理技 术, 推导了一般曲线坐标系下带预处理的拟压缩 Navier-Stokes 方程特征系统, 并结合有限差分法, 建立了适 用于不可压缩黏性流动计算的拟压缩快速算法. 通过对不可压缩无黏圆柱绕流、平板层流流动、低 Reynolds 数定常圆柱绕流问题的数值模拟研究, 得到了与相关理论与实验测试相吻合的结果, 验证了所建立数值方法 的快速可靠性. 较系统地研究了预处理引入的参数和拟压缩因子的选择对收敛特性的影响. 结果表明, Roe 格 式相对于二阶中心差分格式得到的结果更令人满意; 对拟压缩Navier-Stokes 方程进行预处理能有效提高数值 计算的收敛速度; 自适应的拟压缩因子取值能在很大程度上改善数值解的收敛特性, 且不需要根据具体流动 问题进行人工调节. 最后将本文发展的数值方法用于低Reynolds 数非定常圆柱绕流的数值模拟, 所得结果亦 和实验观测结果及其他文献的计算吻合很好. 关键词 拟压缩方法, 预处理, 迎风格式, 不可压缩流动, Navier-Stokes 方程 PACS: 47.11.BC, 47.10.ab, 02.60.Cb, 02.70.Bf 不可压缩 Navier-Stokes 方程中, 压力仅以空间 数方法, 其基本思想是对 Navier-Stokes 方程进行变 梯度的形式出现, 不存在可压缩流动那样通过热力 换, 写成涡量-流函数的形式, 然后引入涡量方程, 消 学关系将压力和密度、温度相联系的状态方程, 没有 去方程中的压力项, 再引入流函数, 连续方程便可以 可以显式求解压力的方程. 故而不可压缩 Navier- 自动满足. 涡量-流函数方法虽在计算二维不可压缩 Stokes 方程数值求解的难点在于速度场变化和压力 Navier-Stokes 方程中取得了成功, 但将该方法推广到 场变化不耦合, 难以满足散度为零的条件. 近年来已 求解三维问题时却出现了困难. 三维流动的流函数 经发展了许多用于求解不可压缩 Navier-Stokes 方程 无法如二维流动一样直接定义, 需要引入多个流函 的算法, 大体上可以分为两类, 一类是基于非原始变 数, 不但场变量增加, 模型方程也更为复杂. 后人在 量的方法, 另一种是基于原始变量的方法. 此基础上发展起来的如涡量-速度方法[1] 、涡量-矢量 基于非原始变量的方法的典型代表是涡量-流函 势方法[2]等一系列方法, 虽都能用于三维问题, 但在 引用格式: 钱战森, 张劲柏, 李椿萱. 基于预处理的不可压缩N-S 方程拟压缩数值算法研究. 中国科学: 物理学 力学 天文学, 2010, 40: 1542 ~ 1554 中国科学: 物理学 力学 天文学 2010 年 第40 卷 第12 期 选择边界条件时将遇到诸多难题, 目前尚难以普遍 方程, 即对可压缩的 Navier-Stokes 方程进行预