基于预处理方法的流动数值模拟分析-力学、流体力学专业论文.docxVIP

南京航空航天大学硕士学位论文 南京航空航天大学硕士学位论文 基于预处理方法的流动数值模拟研究 基于预处理方法的流动数值模拟研究 摘 要 本文研究预处理方法用于数值模拟无粘和有粘低速流动问题。首先，在二维结构网格上， 研制了求解 Euler 方程的计算程序，空间采用有限体积中心格式离散，时间采用四步 Runge-Kutta 格式推进求解，并通过引入当地时间步长和隐式残值光顺等加速措施，成功地用于翼型等典型 可压缩流动算例。为了使发展的程序可用于低速流动，在此基础上，开展了预处理方法的应用 研究，文中对通过预处理降低迭代矩阵刚性，达到加快收敛等基本特点进行了阐述，并选用不 同的预处理方法对发展的可压缩程序进行了预处理，成功地模拟出低速无粘翼型绕流和凸包管 道流动。接着，本文把发展的方法推广到求解 Navier-Stokes 方程，模拟出圆柱和空腔等低速粘 性流动，计算结果都与文献结果进行了比较，展示出不同预处理方法的收敛特性和计算效果。 关键词：预处理，低速，有限体积，Euler 方程，Navier-Stokes 方程 I II II PAGE PAGE VII Abstract Preconditioning method was studied in this paper, for simulating inviscid and viscous low speed flows. Euler solver was developed based 2d structural grid, in which finite volume center method and four-stage Runge-Kutta time-stepping scheme were applied. Convergence rate was accelerated by means of local time-stepping and implicit residual smoothing. The solver simulates compressible flows such as flow past an airfoil successfully. To make the solver can be used in low speed flows, preconditioning method was investigated, the basic characteristics of which were decribed in this paper. The solver was preconditioned by different preconditioning methods. The developed codes were used to simulate the steady inviscid low speed flows past an airfoil in free stream and a bump through a channel. Then, the method was extend to the Navier-Stokes equations to simulate viscous flows past a cylinder and cavity. Calculated results were compared with the literature, demonstrating the convergence of different preconditioning methods and calculation results. Keywords: preconditioning method, low speed, finite method, Euler equation, Navier-Stokes equation 图表清单 图 3.1 相邻网格单元…………………………………………………………………………………13 图 3.2 收敛历程（M=0.8，α = 1.25 o ）………………………………………………………………20 图 3.3 等马赫线图（M=0.8，α = 1.25 o ）……………………………………………………………20 图 3.4 等压线（M=0.8，α = 1.25 o ）…………………………………………………………………20 图 3.5 翼型表面压力系数（M=0.8，α = 1.25 o ）……………………………………………………21 图 3.6 NACA0012 翼型绕流收敛历程（M=0.001）………...…………………………………………