一种基于三角形非结构化网格SIMPLE算法的程序设计

为了深入研究解决流动问题的基本算法,提出了一种三角形非结构化网格上SIMPLE算法的计算机程序设计。利用FORTRAN语言的特点,采用自定义变量存储控制单元以及界面信息。单元变量中定义了指向相邻单元和界面的指针,通过对网格文件所提供拓扑信息的读取对这些指针进行赋值,从而建立计算所需的基本数据结构。这样,就避免了采用数组直接存储网格拓扑信息,大大增加了程序的可读性、拓展性和二次开发潜力。通过顶盖驱动流、后台阶流以及圆柱绕流几个经典计算流体力学算例的验证,该程序计算收敛良好,精度较高,所得结果符合物理实际。因而,可在此基础上开发更具实际应用价值的三维程序。

A computer programming of SIMPLE algorithm based on triangle unstructured grid

A computer programming of SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm based on triangle unstructured grid was developed in order to lucubrate the basic algorithm associated with fluid flow.Self-defined variables were employed for storing element and interface information according to the feature of the FORTRAN language.Pointers which point to adjacent elements and interfaces were also defined within these variables.By reading the topological information from the grid files,the essential data structure for programming could be constructed.Thus,storing the topological information by arrays could be avoided.The readability,extendability and the potential of secondary development were largely improved.The program developed in this paper performed well and gave accurate and reasonable results compared to the physical reality when it was tested by several numerical examples such as lid-driven cavity flow,backward facing step flow and laminar flow passing a cylinder.On this basis,three dimensional computer program with higher application value could be developed.