双翼振动的预处理数值方法和性能研究(一).pdfVIP

双翼振动的预处理数值方法和性能研究（一） 双翼振动的预处理数值方法和性能研究 （一） 1 2 3 王涛 ，刘鹏飞 ，杨国伟 1 康迪斯创科技股份有限责任公司，北京，100080 2 加拿大科学院海洋技术研究所，加拿大纽芬兰，120932 3 中国科学院力学研究所，北京 100080 文 摘：本文在多块耦合隐式求解方法的基础上，引入了预处理计算方法，使方程的各个特征值之间互相 平衡， 降低了方程系统的刚性，提高了数值方法的适应能力，通过 5m /s 圆柱绕流的计算结果与试验 数据验证，得到了数值结果误差小于 5% 的工程允许范围。 关键词：NS 方程，预处理方法，数值模拟。 Preconditioning Numerical Method and Performance Investigation of Plunging Dual-Foils (1) 1 2 3 Tao Wang , Liu Pengfei , Guowei Yang 1 CDSTR Tec. Co. Ltd., Beijing, 100080 2 Institute for Ocean Technology, National Research council Canada, NL, A1B 3T5 3 Institute of Mechanics,Chinese academy of science, Beijing 100080 Abstract: A preconditioning method is introduced to the NS solver based on the dual-time iteration for unsteady compressible flows. This method balances the eigenvalues of NS equations, reduces the stiffness of compressible NS solver for low-speed flow, and increases the capabilities of the solver for low speed flow. By comparing the numerical results for the flow around cylinder at 5m/s with the experiment data, a good agreement was obtained with a discrepancy of about 5%. Keywords: NS equations, preconditioning method ，CFD. 引 言 利用可压的非定常多块网格隐式耦合求解程序研究低速流场问题，必须引入预处理方 法，而在时间相关方程中，由于马赫数和雷诺数的影响，系统表现了很强的刚性。对于Euler 方程（当Re =∞ ），刚性的程度取决于条件数（Characteristic Condition Number，定义为最 大和最小特征速度的比值），在M→0 或者M →1 时，条件数由于最小特征速度很小而变成 无穷大，这就使得任何时间步进方法的收敛性大大降低。另外，低马赫数时标准的离散方法 不能保证精度。在NS方程中，粘性耗散尺度迭加到波传播时间尺度，产生了更强的刚性。 预处理方法处理方程的目的在于平衡方程内在的时间尺度，通过改变时间导数，使得用 数值方法来求解方程时能更