民机全机气动特性数值计算和验证.pdfVIP

第十四届全国计算流体力学会议论文集 民机全机气动特性数值计算与验证 陈晓东，王光学，王开春，郑 呜，邓有奇，王运涛 (中国空气动力研究与发展中心，四川绵阳621000) 摘 要：研制和验证工程实用的民机全机N-S方程计算软件是减少民机气动方案设计和初步设计阶 段风洞试验次数、缩短设计周期、降低开发成本的需要，是开展民机全机一体化气动分析的需要，也是建 立和完善民机全机一体化气动设计系统的基础。本文详细介绍了近年来我单位在民机数值计算与验证方面 所做的工作，其中包括：计算网格生成研究、全机高速、低速N—S方程计算软件的研制、湍流模型应用研 究以及全机气动特性数值计算和验证等。 关键词：民机；气动特性；N_s方程；数值计算 0 引 言 直到20世纪70年代，在飞机气动力设计工作中，风洞试验是获得可靠的气动数据的主要手段。 对于民机设计而言，由于飞机尺寸大，要求风洞试验的Re数也相应的很高，为获得可靠的气动数据， 必须采用更大尺寸的风洞进行数万次试验，工作量非常惊入，耗资亦十分巨大，设计周期更长。但 是，随着CFD技术和大容量计算机的发展，试验和计算在设计工作中的相对地位不断改变。在方案设 计和初步设计阶段，应用CFD技术已成为主要的设计手段，风洞试验仅仅对最终选定的设计方案进 行的校核；在飞机详细设计阶段，CFD技术又可作为精确校核数值计算与风洞试验数据的修正及风洞 与飞行的相关性研究的重要手段，指导和补充风洞试验。由此可见，CFD技术的使用，大大减轻风洞 的负担，使风洞试验用于解决更为重要方面的问题及进行必要的验证，大大缩短了设计周期，节约了 设计成本，确保了飞机的高性能设计指标，对飞行器气动力设计方式带来了深刻的变化，产生了巨大 的社会经济效益。因此，研制和验证工程实用的民机全机卜S方程计算软件是减少民机气动方案设 计和初步设计阶段风洞试验次数、缩短设计周期、降低开发成本的需要，是开展民机全机一体化气 动分析的需要，也是建立和完善民机全机一体化气动设计系统的基础。 1计算方法简介 对于计算全机低速气动力问题(计算起飞着陆构形)，我们采用有限体积法求解雷诺平均的不可 压N_S方程模拟低速粘性流动，湍流模型采用低雷诺数的k一￡两方程模型。不可压N—S方程求解采 用目前国际上比较成熟和先进的SIMPLE系列算法，在保留原有多盲区网格系统计算软件优点的基础 上，从软件的通用性及提高处理复杂外形的能力角度出发，采用多块网格技术和边界多窗口技术研 制SIMPLE算法求解不可压流动的计算软件。 对于计算高速气动力问题，我们采用三维任意曲线坐标系下的雷诺平均N-S方程，所涉及的具 体边界条件有两种，即远场边界条件和物面边界条件。亚声速远场边界利用一维流的特征关系来处 理；超声速远场边界中，入流边界用来流值，出流边界用内场外推得到。物面采用无滑移条件。在 SST二方 方程的求解中采用隐式LU分解的方法求解，湍流模型采用两方程k-‘_)SST湍流模型。k-co 程模型通过混合函数将k一￡和k-o)二方程模型结合起来，利用混合函数E改进涡团粘性系数pt在壁 面逆压流动区域的结果，充分发挥了k一￡模型处理自由流，k-m模型处理壁面约束流动的特长，是目 前广泛采用的两方程湍流模型。 另外，我们还开发了结构／非结构混合网格Euler／№S方程流场解算器，解决了结构／非结构网 格边界的匹配问题和界面流场传值问题。 第十四届全国计算流体力学会议论文集 有关计算方法方面的详细情况参见文献[1～4]。 2网格生成 由于全机外形复杂，如果采用常规的对接网格进行处理，将会遇到网格规模庞大和网格生成困 难等多个问题，为此我们采用了一种比较灵活的处理复杂外形的网格结构形式一对接／拼接网格。拼 接网格的特点是在网格块的公共面上，网格线不需要连续对接，因此各块网格之间的拓扑关系在交 接面上被切断，各块网格之间是基本独立的。在网格生成时，对利用拼接处理的网格块可以单独生 成网格，不需要顾及相邻网格块的拓扑关系，从而使得网格生成的难度大大地降低。另外，由于切 断了各块网格之间的网格线，因此对于局部进行加密的网格点就不会传递到不需要加密的空间或其 它网格块中，从而节省了大量网格单元。另