目标函数误差估算及网格自适应处理

发展了一种基于目标函数误差估算的网格自适应准则,进而通过网格自适应处理提高目标函数计算的效率和准确性。首先描述了目标函数的误差估算及修正的方法,该方法通过伴随方程将原方程的残值误差与目标函数联系起来。然后,基于误差估算建立网格自适应准则,以减少目标函数修正后的剩余误差,提高目标函数计算结果的准确性,并将方法进一步发展至多目标问题。最后将该准则应用于Euler方程控制的NACA0012翼型无粘可压流动的网格自适应模拟。数值计算成功地捕获了与升力、阻力和力矩等目标函数相关的特征流动区域,计算出符合指定精度要求的目标函数值,验证了本文所发展的方法。     

悬停状态倾转旋翼/机翼干扰流场及气动力的CFD计算

基于一套高效通用的多层运动嵌套网格技术,建立了适合悬停状态倾转旋翼机旋翼/机翼气动干扰特性分析的高效混合计算流体力学（CFD）方法。在倾转旋翼/机翼贴体网格区采用可压雷诺平均Navier-Stokes（RANS）方程作为主控方程,过渡/背景网格区选用Euler方程,湍流模型选用Spalart-Allmaras模型。时间推进上采用双时间推进格式进行非定常求解,并在方法中运用了SPMD（Single Program Multiple Data）模式的并行加速技术。在此基础上,首先分别采用UH-60A直升机旋翼及XV-15倾转旋翼机旋翼作为数值算例,验证了CFD方法的有效性。然后着重对倾转旋翼/机翼的非定常干扰流场及气动力分布特征进行了数值研究,模拟得到与实际情况相符的“喷泉效应”干扰现象。计算结果表明,干扰作用使得倾转旋翼相对于孤立旋翼拉力数值减小了3%,但总的拉力系数损失达到了17%,证明悬停状态下气动干扰对飞行器气动性能有重要影响。     

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。     

高超声速飞行器跨流域气动力／热预测技术研究

分析了跨流域流动特点、工程背景需求、国内外研究现状,针对高超声速飞行器跨流域气动力／热预测存在的主要问题,从粘性干扰参数与克努森数搭接的跨流域模拟准则、微量天平结构优化及测力技术、红外大面积中低量值热流测量技术、流场显示与测量技术和 N-S／DSMC 自适应紧耦合数值模拟方法等方面,总结了近年来取得的研究进展,并讨论了下一步研究方向。     

基于改进粒子群算法的串列叶型优化设计

为了提高串列叶型优化设计的质量,设计了一套基于改进粒子群算法的串列叶型自动优化系统。研究了原始粒子群算法,提出了一种粒子群算法的改进方法。结果发现,改进粒子群算法的收敛速度和收敛精度明显优于原始粒子群算法和遗传算法。以50°大弯角串列叶型为研究对象,使用程序对串列叶型参数化。以叶型参数和串列叶型相对位置的参数作为优化变量,结合BP神经网络和改进粒子群算法对串列叶型进行优化设计。优化结果表明,优化后的叶型在设计攻角时的总压损失系数降低了22%,静压比升高了0.6%,在负攻角时,优化后的叶型的流动性能得到了明显改善。适当减小串列叶型前后缘的半径可以减小叶型的损失,合理的缝隙结构可以有效减小前排叶型压力面和后排叶型吸力面附面层的分离损失。     

预设堵块法检测进气道自起动能力的数值研究

在进气道内预先设置轻质堵块迫使进气道不起动,堵块被气流吹出后流道恢复畅通,为激波风洞提供了一种检测进气道自起动能力的方法。为了深入认识预设堵块方法的检测过程,将堵块简化为一个自由度的刚体运动,采用k-ωSST湍流模型,结合铺层动网格技术,对该检测方法进行了二维非定常数值模拟。在Ma_∞=5.9条件下,采用改变堵块质量的方式,获得了三种典型的检测过程。通过分析堵块在气动力作用下的运动与进气道非定常波系演化相互耦合的过程,揭示了堵塞作用产生和消除的流动机理。结果表明,堵塞时间没有显著地改变进气道的自起动能力,预设堵块方法成功检测到进气道在Ma_∞=4.9自起动。此外,建立了预估堵块运动的匀加速模型,预估的堵塞时间与数值模拟结果较为符合。在Ma_∞=5.9条件下,堵块的临界质量约为3.8g。考虑到激波风洞实验时间短暂,应选择合适的堵块。     

基于代理模型的二元收扩喷管流道型面优化设计

在二元收扩（2D-CD）喷管设计参数对喷管气动性能影响研究的基础上,以获得尽可能高的喷管推力系数为目标,以喉道宽高比、喉道型面半径比、收敛半角和扩张半角为设计变量,对二元收扩喷管的流道型面进行了优化设计。设计过程中,利用正交试验设计方法确定初始样本点,建立喷管推力系数与设计参数间的Kriging代理模型,采用自适应模拟退火算法（ASA）对代理模型进行分析求解。结果表明,二元收扩喷管的优化型面参数为:喉道宽高比为6,喉道型面半径比为0.3,收敛半角为15°,扩张半角为5.64°,此时最大推力系数为0.97847,流量系数为0.98778.      

基于阶跃响应的飞行器单独静、动导数数值求解方法

基于阶跃响应方法,结合刚性动网格技术,对飞行器的单独静、动导数的精细化数值计算进行了研究.以纵向为例,通过给物面施加恒定附加攻角,求解得到阶跃响应运动过程的非定常气动力,求导得到静导数.同样给物面施加恒定的俯仰角速度,并同时强迫物面平动以抵消俯仰转动产生的附加攻角影响,可由非定常气动力求导得到动导数值.分别利用NACA0012翼型和三维SACCON飞翼无人机进行了计算验证,各攻角下的静、动导数值与文献、试验结果吻合得很好,最大误差不超过5%.结论表明:基于阶跃响应的单独静、动导数直接模拟方法计算耗时仅为传统强迫振动方法的21%,效率相对较高,且可推广到横航向的动导数计算,为飞行器的稳定性研究提供参考.      

曲线物面边界条件在DGM中的应用研究

基于二维非结构网格,采用高精度间断Galerkin有限元方法数值求解定常Euler方程。为了有效克服传统的反射物面边界条件在物面处易于生成伪熵降低计算精度的缺陷,提出一种曲线物面边界条件,从而实现对物理模型而不是数值模型进行数值模拟。对经典圆柱亚音速无粘绕流进行数值模拟,结果表明:采用曲线物面边界条件之后,流场解的精度得到很好的提高;此外,在非常稀疏的网格上,通过提高基函数的阶次仍可以得到高精度的数值解。     

基于遗传模拟退火算法的无人机航迹规划

航迹规划技术是无人机任务规划系统中重要的核心技术之一,无人机飞行空间广阔,需要一种快速搜索最佳路径的方法.首先在飞行区域中建立数字地图模型和防空威胁区模型,在满足无人机飞行约束条件的情况下,为无人机航迹规划提供一种遗传模拟退火算法,充分利用模拟退化算法的概率突跳特性和遗传算法强大的快速搜索能力.仿真结果表明,使用该算法无人机能够自动避开模拟数字地图的威胁区,搜索出一条安全有效航迹,并保证航线的完整性和最优性.