CFD-FASTRAN理论手册_赵玉新.pdfVIP

CFD-FASTRAN理论手册 国防科技大学航天与材料工程学院赵玉新 序 言 CFD-FASTRAN及其图形用户界面已经发展了很多年，该软件主要用于处理高速流动 问题。熟悉CFD的人都知道，对于可压流动和不可压流动，解决问题的数值方法是完全不 同的，对于不可压流动，方程本身的性质决定了我们主要采用压力速度校正关系实现数值 解法，因此一般被称为基于压力的求解方法。可压流动方程与不可压流动方程的性质却完 全不同，密度的变化在可压流动中十分重要，这就决定了可压流动要有自己的一套基于密 度的解法。CFD-FASTRAN就是基于密度的求解器。其实，随着航空航天的发展，作为流 体力学的一个分支——空气动力学，其发展速度已经远远的超出了流体力学本身。本文将 详细介绍CFD-FASTRAN 的求解方法，从中我们可以看出CFD-FASTRAN所能求解的流动范 围是从低马赫数 （0.1）直至高超声速流动的，因此也主要是针对航空航天问题的求解器。 CFD-FASTRAN最初就是用于求解非移动物体在热完全气体中的高速流动问题的。但 是很多工程问题中是要考虑到运动物体的，而以当时的能力，CFD-FASTRAN无法解决这 类问题。面对困难，CFDRC公司组织人力、物力，将已有的思想——Chimera/Overset及刚 体运动方程集成到CFD-FASTRAN 中。从而实现了Chimera/6DOF/流动求解器之间的结合， 随之对外发布了CFD-FASTRAN V2 。 根据很多内部和外部用户的反馈意见，CFDRC公司认识到CFD-FASTRAN V2在附加的 流动物理模型和运动物体求解方面还有些不足。针对这些意见，CFDRC公司继续扩展流动 求解器和图形用户界面的功能，实现了两大突破：（1）增加了多组分计算能力，从而可以 实现混合和有限速度化学反应的模拟。（2 ）运动模型更加通用化，甚至可以实现规定运 动，增加约束等，从而完善和增强了6DOF算法。增强了这些功能之后CFDRC公司发布了 CFD-FASTRAN V3 。 到CFD-FASTRAN V4 ，主要的提高又包括结构化网格求解器的分布式并行处理及非结 构网格求解器处理通用多面体网格单元的能力。目前，CFD-FASTRAN 已经发展到了V 2003 ，其功能又有了显著的增强，在此不再详述，有兴趣的读者可以浏览CFDRC公司的主 页查阅相关内容。 在阅读本手册以下内容之前需要说明的是：本手册并非简单的介绍，因此需要读者具 有较高的流体力学水平，主要要求读者已经掌握流体力学基础知识、可压缩流动、热传导 以及化学反应流动的基本理论。熟悉差分算法、线性代数、矩阵理论以及计算流体力学的 基础知识。 第一章功能介绍 CFD-FASTRAN理论采用一种模块化的思想，这些模块包括： 1. 流动求解器支持的网格类型：结构和非结构网格 2. Overset/Chimera算法 3. 运动体 （六自由度/规定运动） 4. 流动的物理特征 1.1 网格类型 CFD-FASTRAN有两个求解器，分别为结构网格求解器和非结构网格求解器。结构网格 求解器的功能和支持的物理模型要比非结构网格求解器强大一些。 结构网格求解器采用结构化多区域方法实现复杂流动的模拟。固体阻碍采用障碍物的 概念来实现模拟。目前，CFD-FASTRAN没有自己内嵌的结构化网格生成代码，一般采用 其它网格划分软件生成的网格来实现前处理。主要的网格生成器是CFD-GEOM，这也是 CFDRC公司的产品，专门用来做前处理网格划分的。CFD-FASTRAN支持的格式包括 PLOT3D和DTF两种。 非结构网格求解器采用通常的非结构网格、多区域方法。支持各种类型的非结构网 格：在二维流动中支持三角形、四边形和多变性，在三维流动中支持四面体、六面体、金 字塔、棱柱及多边形。非结构网格求解器在同一区域内可以处理不同类型的单元。该求解 器本身也没有网格生成代码，最常用的是采用CFD-GEOM和CFD-Viscart生成的DTF格式网 格文件。 1.2 Overset/Chimera算法 Overset算法 （重叠网格）只能应用于结构网格求解器。虽然Overset算法可以减少网格 生成的时间和负担，但是要考虑运动体的问题则必须由Chim