基于三维直角网格Euler方程并行计算的气动布局优化研讨.pdfVIP

192 探索创新交流(第二集)——第二届中国航空学会青年科技论坛文集 基于三维直角网格Euler方程 并行计算的气动布局优化研究 王 荣张卫民 (中国航天空气动力技术研究院第一研究部，北京100074) 摘要：本文运用基于三维非结构直角网格Euler方程实现了流场分区并行计算，并以此进行飞行器气动 布局外形的优化研究，通过流场分区并行计算，缩短了流场解算时间，使得以高精度数值模拟流场优化飞 行器外形布局的周期达到了实际工程设计可接受的限度。通过优化无人机外形取得了有意义的结果，优化 后升力系数增加了21％，升阻比提高了6．7％。 关键词：非结构童角网格；Euler方程；并行计算；气动优化 由于计算机性能的飞速提高，与计算机技术紧密相关的计算流体动力学(CFD)在实 际需要的推动和硬件的支持下取得了长足发展，但是伴随着技术的进步，飞行器气动设计 工程应用对计算机性能提出了更高的要求。在CFD应用中，通常的计算资源可以满足单 个飞行状态的预测需求，然而以CFD进行数值优化设计将会上百乃至上千次调用流场分 析软件[1]，这就意味着计算量增加了两到三个数量级，若要在可接受的时间周期内取得 数据，计算机性能就要有相应的增长量，现有的单机水平远不能达到这么高的程度，而且 在短期内也难以实现。采用集群计算机技术与CFD网格分块并行计算是数值优化发展的 必然趋势，也是解决实际问题的一条路径。高效、自动的网格生成技术一直是CFD技术 发展的瓶颈，而本文所使用的非结构直角网格能够快速高质量地离散流场空间，对物面外 形不敏感，适用于任何复杂形状的几何物体，不仅克服了网格生成的难题，而且使利用 CFD数值分析优化任意复杂外形变得简单易行。因此本文运用基于非结构直角网格的Eul— el方程实现了流场分区并行计算，通过流场分区并行计算，缩短了流场解算时间，使得用 高保真度数学模型数值模拟流场的优化设计周期满足了实际工程的时间要求。利用此技术 优化了无人机外形，并取得了有意义的结果。 1非结构直角网格简述 非结构直角网格是一种使用较为方便的网格生成技术，它不必先生成计算的表面网格 再生成空间网格，能够一次性生成计算所需的网格，可以通过直角网格与物体的相交情况 “勾勒”出物体的外形。以三维八叉树非结构直角网格为例，网格生成过程可分为以下几 步‘2叫。 (1)初始网格生成，用一个较粗的长方体网格覆盖包含物体的整个计算域，在 探索创新交流(第二集)——第二届中国航空学会青年科技论坛文集 193 菇、Y、z三个方向划分该长方体区域得到由划分后较小的长方体组成的第一层结构化 网格单元。 (2)判断所有长方体网格单元与物面相交情况。 (3)对与物面相交的或距物面较近的长方体网格单元进行进一步的加细，将其剖分 成八个更小一级的长方体子单元。 (4)对第3步新生成的长方体子单元重复(2)、(3)步直到满足用户设定的要 求为止(比如，达到用户设定的最大网格加密层数或用户限定的最小网格尺寸，网 格剖分结束)。 (5)找出与物面相交的有效网格，这里的有效网格指不再剖分的网格单元，即没有 子单元的所有网格。接下来处理物面与空间网格作用的内边界面，生成贴近物面的一层网 格(物面层网格)。目前较为常用的一种处理内边界的方式是投影法，即通过将与物面相 交的有效网格删除后产生的间隙层的外边界投影至物面形成物面层网格。这一步是整个网 格生成的关键，因为直角网格与给定的物面几何外形作用不仅要生成质量良好的物面层网 格，还得以物面层网格映射在物面上的侧面的集合来模拟物面外形而且应当较好地维持给 定几何外形的保形性，这也是生成网格最复杂的一步。 (6)判断舍弃物面内网格(外流计算)，求出计算域网格(物面外有效的网格与物面 层网格的总和)要求的信息，如体积、侧面面积等。 (7)输出流场解算所需的网格信息，网格生成结束。 从以上非结构直角网格的生成步骤能够看出，相比于结构网格和其他非结构网格，采 用非结构笛卡儿网格具有以下优点日“】：在原始的均匀笛卡儿网格基础上可根据几何外形 特点或流场特点在局部区域内不断进行网格细化，得封精度符合要求、分布又是最理想的 非均匀笛卡儿网格，以此计算复杂几何形状的流场，达到准确模拟外形和捕捉激波等目 的；其