机翼—机身—吊架—短舱的跨声速流计算

本文通过使用多块网格的嵌合全技术，计算了机翼－机身－吊架－短舱的跨声速流场。将机翼－机身－吊架－短舱流场分为机翼机－机身和吊架－短舱两个子区域，根据嵌套网格结构建立机翼－机身－吊架－短舱的组合网格，分别对这两个子区域根据相应的边界条件独立求解，然后采用嵌合体技术传递两个子区域的干扰信息，通过两个流场的耦合迭代得到机翼－机身－吊架－短舱流场的计算结果，数值模拟结果与实验结果进行了比较，两者吻合较好。     

机翼－机身－吊架－短舱的跨声速流计算

本文通过使用多块网格的嵌合体技术，计算了机翼－机身－吊架－短舱的跨声速流场。将机翼－机身－吊架－短舱流场分为机翼－机身和吊架－短舱两个子区域，根据嵌套网格结构建立机翼－机身－吊架－短舱的组合网格，分别对这两个子区域根据相应的边界条件独立求解，然后采用嵌合体技术传递两个子区域的干扰信息，通过两个流场的耦合迭代得到机翼－机身－吊架－短舱流场的计算结果，数值模拟结果与实验结果进行了比较，两者吻合较好。     

飞机/外挂物组合体绕流的数值计算

采用嵌合体埋入方法计算带外挂翼-身组合体复杂外形的流动,将飞机/外挂物外形的流动区域分为机翼-机身和外挂物两个子区域,在机翼-机身子域上求解欧拉方程,在外挂物子域上求解全速势方程,通过两个流场的耦合迭代计算全流场的解。对C-5A飞机及某战斗机的计算表明,数值模拟与实验吻合。实现了不同子域流场使用不同方程的计算方法。     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

翼-身组合体跨声速流场计算

本文把Holst的AF2因式分解格式推广到求解机翼-机身跨声速绕流的全速势方程。利用求解椭圆型偏微分方程组来数值生成贴体的计算网格。算例表明本方法运算可靠,费时较少,可以向工程实用部门推广使用。     

同向旋转螺旋桨对着陆构型气动影响分析

基于分区对接网格技术分别生成翼身区域和螺旋桨区域网格,并通过求解定常N-S方程模拟螺旋桨滑流影响下全机流动。通过对某螺旋桨飞机的着陆构型(机翼+短舱+螺旋桨+机身+平垂尾)进行数值模拟,结果与实验值吻合较好。模拟结果表明:螺旋桨滑流会改变机翼当地绕流迎角、影响压力分布;当螺旋桨同向旋转时,机身两侧襟翼的分离会出现在不同位置并改变升力在机翼展向的分布,进而对全机产生一个滚转力矩。     

基于欧拉方程的机翼跨音速颤振的频域计算

机翼跨音速颤振的频域计算方法是以给定机翼模态分布下机翼上各点的模态值作为运动幅值，以三维非定常Euler方程为控制方程，采用有限体积法和双时间推进，求解三维机翼简谐运动下的非定常气动力。所求得的气动力作为已知值运用于颤振方程，利用v-g法进行求解。对得到的一系列的阻尼、速度和频率进行了线性插值，从而得到颤振速度和颤振频率。     

旋转机翼飞机旋翼模式前飞状态干扰气动特性

和传统直升机相比,旋转机翼（CRW）飞机在旋翼模式前飞时各部件之间存在更为严重的气动干扰。为了获得旋转机翼/机身/鸭翼/平尾之间的非定常气动干扰规律,基于运动嵌套网格技术,通过求解三维非定常雷诺平均Navier-Stokes（URANS）方程,建立了旋翼前飞流场数值模拟方法。首先对传统直升机旋翼/机身干扰模型进行了计算,验证了方法的可靠性,然后对某旋转机翼飞机全机在旋翼模式前飞状态下的非定常流场进行了数值模拟,并对各个气动部件上的非定常气动力和力矩的变化进行了分析。结果表明：飞机在旋翼模式前飞时,机身部件对旋转机翼的干扰较弱,在经过机身上方时拉力峰值仅略有增加;旋转机翼对鸭翼和垂尾干扰较弱,对机身和平尾干扰较强,随着前飞速度增大,旋转机翼对平尾的干扰会产生较大的升力损失和抬头力矩,需要引起重视。计算结果为该类飞行器的总体综合设计提供了参考。     

基于黏性伴随方法的跨声速机翼气动优化设计

 进行了基于黏性伴随方法和Navier Stokes方程的跨声速机翼气动优化设计研究。分别推导了适用于三维跨声速机翼气动反设计和减阻设计的黏性伴随方程、边界条件和梯度求解表达式，并研究了伴随方程的数值求解方法。通过将网格生成、流场计算、黏性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合，发展了一种跨声速机翼气动优化设计方法。为了提高计算效率，将多重网格方法应用到方程的数值求解中来加速收敛。跨声速机翼反设计和减阻设计算例验证了本文所发展的方法的正确性。采用本文的方法进行优化设计，一般通过20~30次迭代就能得到满意的结果。     

基于双eN方法的翼身组合体流动转捩自动判断

发展翼身组合体复杂外形流动转捩自动判断方法,对高亚声速民机自然层流（NLF）机翼设计具有重要意义。使用多块结构化网格和三维雷诺平均Navier-Stokes（RANS）方程求解器,耦合边界层方程求解和基于线性稳定性理论（LST）的完全双e^N方法,发展了一套可同时计及Tollmien-Schlichting波和横流不稳定性扰动诱导转捩的翼身组合体流动转捩自动判断方法。对DLR-F4翼身组合体绕流进行了转捩自动判断,将得到的转捩位置与试验结果进行比较,验证了所发展方法的正确性。使用上述方法对配置自然层流机翼的中短程民机翼身组合体外形进行了数值模拟,并将结果与单独机翼的转捩位置进行了对比,结果表明机身三维位移效应增强了自然层流后掠机翼边界层的横流不稳定性强度,导致翼根转捩位置提前至前缘区。     

机翼机身跨声速绕流的计算

本文采用变换与数值网格生成方法，对由椭圆抛物而描述的每一个展向曲面进行网格生成，引入一种用于控制网格疏密分布的函数变换方法，其特点是易于对不同族网格线分别控制，根据物面梯度自适应地调整网格的疏密。用有限体积法计算守恒型全速势方程跨声速绕流，对M-100翼身组合体进行计算。理论计算结果与实验结果表明，本文方法有效。     

基于Delaunay图映射的动弹网格技术及其应用

在多块混合网格的基础上发展了一套基于Delaunay图映射方法的动弹网格技术,实现了针对气动弹性计算的网格变形运动方法,并使用基于多块混合网格的RANS方程的有限体积法耦合结构静平衡方程,针对大展弦比无人机在实际使用中可能出现的显著的静气动弹性问题进行了计算、分析和研究。研究结果表明,所发展的动弹网格方法为粘性混合网格的变形提供了较好的解决方案;RANS方程解算器具有较高的精度和较强的适用性,具备求解粘性流体与结构耦合等复杂航空气动弹性问题的能力。     

一种CFD/CSD耦合计算方法

针对柔性大展弦比机翼气动弹性分析和主动弹性机翼(AAW)设计发展了一种计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法。其主要思想是采用在同一物理时间弱耦合求解CFD/CSD技术。气动力采用非定常N-S方程的双时间有限体积求解技术,结构响应则采用有限元数值求解技术。CFD和CSD耦合计算的边界信息(气动力和网格)由所设计的界面程序传输。网格信息传输采用守恒体积转换(CVT)方法将CSD计算结构响应位移插值到CFD网格点上。变形已有的CFD网格技术用以确定CFD的变形网格。以位移或载荷的迭代误差为判断耦合计算的收敛标准。最后得到了机翼在Ma=0.8395,α=5.06°时CFD/CSD耦合计算的收敛值。针对计算结果分析了机翼受静气动弹性过程中结构响应和气动特性随时间变化的效应。初步研究结果表明:这种弱耦合方法求解非线性气动弹性问题是可行的。     

用于复杂流场数值模拟的重叠／嵌套网格方法

介绍了一类用于复杂流场数值模拟的重叠／嵌套网格方法。它是一类多块结构化网格,可用于有或无多重网格加速收敛技术的有限体积、有限差分等计算。其基本思想是把复杂流场的计算区域化分成相对简单的子区域,每个子区域上的网格可以独立生成;网格界面之间的流场信息交换通过三线性插值来实现,通过修改流场解算方法来排除洞中的点。     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

基于超椭圆方程和序列响应面法的回转壳开孔形状优化

含有多种功能开孔的轻质回转壳结构是航空航天结构系统内的常用支撑构件,其开孔形状直接影响结构的静动态性能。以回转壳结构为对象,基于超椭圆方程和坐标映射变换推导了回转壳开孔边界的参数化表达,开展了开孔形状动力学优化研究。为提高结构优化计算的精度、效率和收敛性,提出了准等弧长方法和基于均匀设计的序列响应面近似建模方法(SRSM),以分别实现空间超椭圆曲线的精确逼近、减少结构有限元分析成本和加快迭代收敛。以非支配排序遗传算法Ⅱ(NSGA-Ⅱ)作为响应面模型求解算法,结合有限元分析构建了回转壳开孔形状优化设计流程,开展了最大化结构一二阶频率带隙的典型回转壳结构开孔形状优化设计。结果表明,基于超椭圆方程和序列响应面法的开孔优化方法获得了有效改进结构动态特性的回转壳开孔形状,对开展计算耗时工程结构形状优化设计具有一定应用价值。     

用于粘性可压缩流动数值计算的SIMPLE方法

对SIMPLE方法加以推广,使之适用于可压缩粘性流动的数值计算。采用非交错网格技术,推出三维任意曲线坐标系下可压缩形式SIMPLE方法的计算方程。对亚音速、跨音速及超音速等5个流场进行了数值计算并与有关文献及实验数据进行比较。      

翼-身-尾组合体的网格生成及 Euler 方程数值模拟

采用代数方法快速生成翼-身-尾组合体三维O-C型贴体网格。该代数网格生成方法使用双边界法,可保证物面附近近似正交,真实地模拟了机身头部及与机翼、尾翼交界处的形状,提高了代数方法生成网格的质量。并用中心格式有限体积法求解跨声速Euler方程,以某翼-身-尾组合体模型为例,计算结果与实验符合良好。     

Aerodynamics of smart flap under ground effect

Aerodynamic characteristics of two-dimensional smart flap under the ground effect have been assessed by a numerical simulation. In this process, a pressure-based implicit procedure to solve Navier–Stokes equations on a nonorthogonal mesh with collocated finite volume formulation is used. The boundedness criteria for this procedure are determined from the Normalized Variable Diagram (NVD) scheme. The procedure incorporates the k–ε eddy–viscosity turbulence model. Cantilever beam with uniformly varying load with roller support at the free end is considered for the configuration of the smart flap. The method is first validated against experimental data. Then, the algorithm is applied for turbulent aerodynamic flows around airfoil with smart and conventional flaps for different attack angle, flap angle and ground clearance where the results of two flaps are compared. The comparisons show that the quality of the solution is considerable.     

计算机翼跨声速压力分布的格林函数方法

本文给出一种用格林定理将跨声速小扰动速势方程化成积分方程式,然后进行数值求解的方法,文中对场元强度进行松弛迭代,一般只需20次左右即达收敛。算例表明,该方法结果与差分法结果一致。