基于Euler方程的三维自适应笛卡尔网格应用研究

本文采用非结构自适应笛卡尔网格和有限体积法求解了三维Euler方程。在网格生成过程中，采用了基于模型几何数据、模型表面曲率和流场特征为基础的网格自适应生成及加密方法，采用模型特征恢复技术保证了计算模型与原始外形的一致。在计算中，将以中心差分为基础的Jameson有限体积法在三维笛卡尔网格中进行了推广。本文中对导弹、歼击机等复杂问题进行了数值实验，计算结果与风洞实验结果符合良好，并在工程实际中成功应用。     

非定常旋涡-翼型相互作用问题数值模拟研究

针对直升机旋涡-桨叶相互作用问题,以其二维简化模型旋涡-翼型相互作用问题为研究对象,发展了基于雷诺平均Navier-Stokes方程和自适应混合笛卡尔网格的流场数值模拟方法,通过给定合适的自适应判据,动态网格自适应技术能够准确捕捉旋涡位置和保持旋涡强度。通过两类旋涡-翼型相互作用问题的数值模拟研究,揭示了旋涡与翼型相互作用过程中产生的复杂流动现象和非定常气动力变化。数值模拟结果与相关文献中的试验数据和数值计算结果吻合,且对三维桨叶的气动载荷研究和振动分析具有指导意义。     

共轴式双旋翼悬停流场和气动力的CFD计算

考虑到共轴式双旋翼流场特征高度复杂、桨叶承受非定常气动载荷的特点,为更好地预测共轴式双旋翼的气动特性,把非结构嵌套网格方法和网格的自适应技术相结合,发展了一套适合于共轴式双旋翼流场数值计算的求解器.在该求解器中,采用非结构嵌套网格方案来模拟桨叶之间存在的相对运动,自适应网格技术用来捕捉尾迹对流场和气动特性的影响,求解惯性坐标系下的非定常N-S主控方程来模拟流场的非定常特性.应用该求解器,首先计算了有试验结果可供对比的一副试验旋翼的诱导速度场分布,在此基础上,计算了共轴式双旋翼的桨尖涡轨迹和拉力分布特性,并与单旋翼的计算结果进行了对比和分析,得到了一些有意义的结论.     

基于HyperFLOW平台的客机标模CHN-T1气动性能预测及可信度研究

近年来,针对实际飞行器外形的CFD气动性能预测及可信度研究逐步得到重视,国内也召开了第一届航空CFD可信度研讨会(AeCW-1)。本文首先基于自主研发的CFD软件平台HyperFLOW对NACA0012翼型低速绕流进行了网格收敛性研究,验证了软件对简单湍流问题的模拟能力且具备良好的网格收敛性。其次,针对AeCW-1提供的客机标模CHN-T1,选用其中的两个算例:(1)定升力系数的网格收敛性研究;(2)考虑模型支撑和模型静气动弹性变形的抖振特性研究,研究了计算结果的网格收敛性及模型支撑、静气动弹性变形和湍流模型等对气动特性预测精度的影响。结果表明:观测精度阶和网格收敛性指数显示数值结果具有良好的网格收敛性和可信度;是否考虑模型支撑对力矩的预测精度影响较大,引入尾撑和弹性变形后,数值结果与实验结果吻合较好;对于CHN-T1标模,采用QCR关系式对原始SA模型进行修正对标模力矩特性有一定影响。     

DPWⅢ机翼和翼身组合体构型数值模拟

采用亚跨超CFD软件平台(TRIP)数值模拟了DPW Ⅲ提供的DPW_W1/W2两种机翼构型和DPW-F6/F6_FX2B两种翼身组合体构型,主要目的是通过两种机翼构型和两种翼身组合体构型的数值模拟,研究网格密度对运输机构型气动特性计算结果的影响。数值模拟采用的多块对接网格来自AIAA CFD Drag Prediction Workshop Ⅲ(DPW Ⅲ),采用National Transonic Facility(NTF)的试验结果和CFL3D的计算结果作对比。详细研究了网格密度对两种机翼构型和翼身组合体的总体气动特性和压力分布的影响。采用SST两方程模型计算两种构型均得到了网格收敛结果,网格密度主要影响压差阻力而对摩擦阻力影响较小,计算结果较好地预测了机翼和翼身组合体外形优化前后总体气动特性的变化量。     

CFD-FASTRAN气动热计算模型及网格效应分析

介绍了CFD-FASTRAN软件五种湍流模型。通过B-L湍流模型不同网格气动热计算结果的对比分析,说明近壁面第一层网格间距须保证y＋值处于所选湍流模型要求值范围的1/3左右,才能保证计算结果的准确性。进一步通过直二次圆锥形不同湍流模型计算与试验的对比结果,说明CFD-FASTRAN软件的k-ε及B-L两种模型求解高超声速气动热问题能达到较好的准确性,且B-L模型计算速度快、使用方便,更适宜于整机大型数值计算。     

基于非结构网格求解器的CHN-T1标模气动特性计算研究

高可信度计算流体力学(CFD)技术在大型运输机精细化气动设计中发挥着越来越重要的作用。虽然目前基于RANS方程的数值模拟技术已经应用于大型运输机复杂外形的气动性能计算与分析中,但其可信度仍需进一步确认。针对AeCW-1组委会发布的CHN-T1运输机标模,采用基于RANS方程的非结构网格求解器TNS进行了气动性能计算,并与组委会提供的风洞试验结果进行了对比。此外,还研究了网格量和湍流模型等因素对巡航气动特性计算结果的影响,以及支撑干扰和静气弹变形对跨声速抖振特性的影响。结果表明,TNS预测的巡航升阻比及跨声速抖振特性等气动性能与试验值吻合良好,计算结果完全位于试验值误差带内,验证了所采用的非结构网格求解器在大型运输机复杂外形气动力预测中的可靠性。     

采用TRIP2.0软件计算DLR-F6构型的阻力

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了DLR-F6构型,主要目的是通过计算DLR-F6构型的安装阻力考察TRIP2.0软件的数值模拟精度,并为运输机构型的气动特性计算积累经验.本文数值模拟采用的多块对接网格,测压和测力的试验结果均来自AIAA CFD Drag Prediction Workshop II(DPWII),对比计算结果采用了CFL3D的结果.本文详细研究了网格密度、湍流模型对DLR-F6翼身组合体和翼/身/架/舱复杂组合体两种构型的的总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致.本文采用SST两方程模型计算两种构型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布影响较小.     

采用TRIP3.0模拟CHN-T1模型气动特性

采用TRIP3.0软件平台(TRIsonic Platform version 3.0),对第一届航空CFD(Computational Fluid Dynamics)可信度研讨会组委会提供的运输机标模构型(CHiNa-Transport,CHN-T1)进行了流动数值模拟。本文采用粗、中、细三套网格及百亿极细网格进行了网格收敛性研究;采用相应网格分别从气动特性、压力系数分布曲线、表面流态三个方面分析,研究了支撑装置、机翼静弹性变形以及雷诺数效应对CHN-T1构型气动特性的影响。数值模拟结果与试验结果有良好的一致性。数值模拟得到了网格收敛结果;支撑装置对力矩特性影响较大;机翼静弹性变形对气动特性影响较小;雷诺数效应对气动力特性影响较大。     

梯形翼风洞试验模型数值模拟技术

基于雷诺平均Navier-Stokes (RANS)方程和结构网格技术,采用二阶空间离散精度的MUSCL格式,并结合k-ω剪切应力输运(SST)两方程湍流模型和γ-Re_θ转捩模型,研究了梯形翼风洞试验模型中前缘缝翼、后缘襟翼连接装置对气动特性的影响。简要介绍了本文采用的计算方法;介绍了梯形翼的风洞试验模型及风洞试验结果;在网格收敛性研究的基础上,采用"全湍流"方式和转捩模型研究了梯形翼试验模型连接装置对气动特性的影响。通过与不带连接装置的计算结果的对比,采用"全湍流"模拟方式,计算模型中考虑试验模型的连接装置引起升力系数下降、阻力系数下降、低头力矩减小以及失速迎角提前;通过与试验数据的对比,进一步考虑转捩影响可以提高梯形翼风洞试验模型气动特性的计算结果与试验结果的吻合程度,梯形翼风洞试验模型失速迎角附近的气动特性数值模拟技术还需要进一步的研究。     

天线罩载荷谱编制与验证

天线罩气动载荷为空间力系，我们对气动载荷分别进行三个坐标平面投影，从而使问题等价为三平行力系，分别采用计算机交互分区法进行了载荷分区。计算验证是基于有限元分析技术，对天线罩建立模型并施加气动载荷进行计算，对胶布带的施加区域进行网格局部加密，在胶布带区域施加载荷，从而较准确地模拟了实际加载情况。试验应变监测表明：疲劳试验4倍寿命试验后，通过无损检测。各区域未见异常。     

网格拓扑对DLR-F6构型数值模拟的影响

基于TRIP3.0软件平台和多块结构网格技术,开展了控制方程和网格拓扑等因素对DLR-F6构型数值模拟的影响研究。数值计算采用与试验相同的参数,采用了RANS/TLNS方程,生成了O型和H型2种拓扑的粗、中、细网格进行模拟。分别从气动特性、压力系数分布曲线和表面流态3个方面对结果进行分析。通过与试验数据的对比表明,在模拟小分离流动时,采用H型拓扑网格和RANS方程模型,获得的计算结果会更为准确。     

基于混合笛卡尔网格的翼型绕流数值模拟

在对重叠网格技术研究的基础上,一套以笛卡尔网格为主体、贴体的结构网格为特殊固壁边界的混合笛卡尔网格(Hybrid Cartesian Grid,HCG)生成方法得到了发展,同时基于流场特征的笛卡尔网格自适应技术和有限体积方法的雷诺平均Navier-Stokes(RANS)的数值求解方法也得到了开发。通过对经典跨音速RAE2822翼型、双NACA0012翼型和高升力两段NLR-7301翼型进行绕流数值模拟,获得了较好的结果,验证了方法的可行性与精确性,为研究二维复杂外形流动问题提供了一种解决方法。     

基于笛卡尔网格的非定常气动力计算

在非结构自适应笛卡尔网格基础上构建了一种基于笛卡尔动网格的非定常气动力计算方法。为了适应非定常流场对网格品质要求较高的特点,网格生成采用了基于几何外形、模型表面曲率为基础的自适应生成及加密方法。采用显格式龙格库塔推进方法和中心格式空间离散的有限体积法对NACA0012翼型,三维机翼的跨音速刚性振动进行了数值模拟。对伪时间收敛精度控制,真实时间步长对结果的影响进行了讨论。将计算结果与试验数据进行了比较,得到了良好的效果。     

DLR-F6翼身组合体的高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6翼身组合体的高阶精度数值模拟研究。主要目的是确认WCNS模拟跨声速典型运输机构型的能力。采用粗、中、细3套网格开展了网格收敛性研究,从气动特性、压力分布、表面流态等方面研究了网格密度对DLR-F6翼身组合体气动特性的影响;采用中等网格开展了来流迎角对气动特性的影响研究。通过与试验数据、CFL3D软件和TRIP软件计算结果的对比,表明网格密度主要影响激波位置和压差阻力系数,同时对翼身结合部分离区大小有一定影响;采用高阶精度计算方法显著提高了气动力系数的模拟精度,力矩系数数值模拟结果与试验的差异有待进一步分析。     

基于自适应代理模型的气动优化方法

气动优化设计中,为了减少优化系统的计算周期,提高搜索效率,引入结构简单、计算量较小的代理模型,而运用有效的插值和选样方法(自适应选样)可以大大减少建立代理模型的时间。因此本文提出了一种基于自适应代理模型的气动优化方法。首先对自适应代理模型进行研究,建立了 Kriging 自适应代理模型和支持向量回归自适应代理模型,这两种自适应代理模型在相同样本点情况下比一般代理模型拥有更高的预测能力,然后将这其应用到翼型优化设计中,取得了良好的优化效果,从而表明这两种自适应代理模型不仅简单实用,而且明显提高了气动分析的计算效率。     

CRM-WBPN风洞试验模型数值模拟

基于雷诺平均Navier-Stokes（RANS）方程和多块对接结构网格技术,采用亚跨超CFD软件平台（TRIP 3.0）,研究了共同研究模型（CRM）风洞模型支撑装置对CRM翼/身/架/舱（CRM-WBPN）模型压力分布及气动特性的影响。构造了一个包含小、粗、中、细网格的4套网格,并进行了网格收敛性研究,在升力系数为0.50条件下,开展了固定升力系数下模型支撑对CRM-WBPN模型和CRM-WB模型压力分布及挂架短舱阻力增量的影响,以及不同来流迎角下模型支撑对CRM-WBPN模型气动特性的影响。通过与NASA国家跨声速风洞设备（NTF）试验结果和CRM-WB模型数值模拟结果的对比分析,表明模型支撑对固定升力系数下挂架短舱阻力增量影响较小;在计算构型中考虑模型支撑的影响,引起升力系数和阻力系数下降、俯仰力矩系数增加,数值模拟结果更加接近相应的试验结果。     

基于共模型技术的大展弦比飞机气动/结构耦合性设计

在考虑弹性变形情况下,以大展弦比飞翼飞机为例,对气动/结构耦合设计方法进行研究提出了几何/结构/气动的共模型技术.首先是基于共模型技术生成表面气动网格和空间气动网格以及结构有限元网格,并实现气动载荷在结构模型和气动模型间的数据传递,最后进行气动/结构耦合计算分析,在进行气动/结构耦合计算分析时要进行多次迭代,直到结构变...     

基于非结构网格CFD技术的旋翼气动噪声计算方法研究

将基于非结构网格技术的旋翼流场CFD计算方法与基于FW-H和Kirchhoff方程的声学方法相结合,建立了一套既适合于直升机旋翼厚度、载荷和桨-涡干扰噪声,又适合于跨声速高速脉冲噪声的综合计算模型。为提高旋翼流场及桨叶表面气动载荷计算的精度,主控方程的求解采用了三维可压非定常的N-S方程,网格划分则使用非结构运动嵌套网格方法。在噪声计算中,通过FW-H方法计算旋翼的厚度噪声和载荷噪声,并选取能够包含流场非线性区域的旋转面作为Kirchhoff积分面,由Kirchhoff方法计算包含四极子项的高速脉冲噪声。应用该模型,以AH-1旋翼为算例,计算了不同飞行状态下的旋翼气动噪声,并与可得到的试验结果进行比较,验证了方法的有效性。然后,着重对两种声学方法对计算结果的影响进行了对比研究,并分析了旋翼厚度噪声、载荷噪声和四极子噪声的特性。     

HiLiftPW-3高升力构型数值模拟

基于雷诺平均Navier-Stokes方程和拼接网格技术,采用MUSCL-Roe格式和Spalart-Almaras一方程湍流模型,对第3届高升力构型性能预测会议提供的两组高升力标模进行了数值模拟,主要目的是确认本文计算方法模拟复杂高升力构型的能力。研究内容主要包括高升力构型网格生成技术、网格收敛性研究及气动特性数值模拟。通过与JAXA （Japan Aerospace eXploration Agency）提供的测压、测力试验结果的对比分析,表明,在失速迎角之前,数值模拟得到的气动力系数和压力分布均与试验结果吻合;较好地模拟了局部外形变化引起的气动特性差量。本文建立的数值模拟方法对典型运输机三段翼布局的低速问题具有良好的适用性,可以为大飞机低速构型的气动设计及评估提供技术支撑。