动量增升高升阻比飞行器横航向稳定性研究

考虑高超声速飞行器装填、防热及操稳等多约束条件下的实用化设计需求,提出一种新型下反式高升阻比滑翔飞行器气动布局。借鉴乘波体飞行器的高升力设计方法及动量增升原理,该新型飞行器采用下反式后掠翼构型,上表面为光滑倒圆“Λ”形设计,下表面为内凹的装填空间,为尽可能避免长时间超远距离高超声速飞行带来的防热负担,采用后缘体襟翼及体侧扩张方向舵面设计。采用数值计算方法对所提布局思路进行验证分析,计算结果表明：在飞行高度为40 km,Ma=10的条件下升阻比可以达到4.48左右,在一定迎角范围内均具备很高的气动效率,验证了所提布局的有效性。同时重点针对所提布局共性的横航向稳定性问题基于数值模拟方法探讨了3种不同优化改进方案的效果及可行性,并采用风洞试验对两侧翼梢V尾的控制方案进行横航向稳定性控制效果的试验验证,结果表明：采用两侧翼梢V尾的控制方案是实现横航向稳定性控制的较优方案。     

FlowStar:国家数值风洞(NNW)工程非结构通用CFD软件

计算流体力学（CFD）仿真软件是流体相关的数学物理知识和工程实践经验的数字化表达,是工业数字化转型的重要助推。然而,大型工业CFD软件研发难度极高,需要同时兼顾功能多样、系统稳定、性能优越、交互友好等特征。依托国家数值风洞（NNW）工程,研发出一款通用流场模拟软件NNW-FlowStar,并在航空、航天等工业部门大力推广使用。软件基于非结构有限体积求解方法和大规模并行计算技术开发,结合现代化软件工程思维设计,具备先进的数值方法、高效的计算效率和友好的用户操作界面,可满足各类复杂外形的高效气动模拟。独特的重叠网格技术配合六自由度运动模块,可帮助实现武器分离、舱门定轴转动等各类气动-运动协同仿真需求。多类标模案例和复杂工程应用表明,FlowStar软件算法鲁棒、精度可靠,是一款高精度、高效率、高可靠性的通用CFD仿真软件。通过对软件的架构设计和功能应用进行介绍,使相关从业人员能更好地了解FlowStar软件,最终促进国产自主CFD软件生态的良性发展。     

多场耦合效应对高超声速进气道入口参数影响

为了研究多场耦合效应对高超声速进气道入口参数的影响,采用自主开发的热环境/热响应耦合计算分析平台FL-CAPTER,对吸气式高超声速进气道前体进行了数值仿真研究。介绍了采用的多物理场耦合分析策略及不同物理场求解方法,通过圆管和两级压缩楔外形,初步验证了多场耦合分析方法的可靠性。以此为基础,研究了进气道前体在长时间巡航飞行条件下的结构温升情况和宏观变形量,分析了进气道结构变形对入口参数的影响。结果表明:进气道前体迎风区域和背风区域不均匀的温度分布引起热应力变化,进气道前体压缩面在多场耦合效应作用下上翘约20mm,考虑变形影响后,进气道偏离设计状态,激波边界层干扰效应增强,喉道附近的分离区域有所增大,进气道入口的质量流量增加约4.2%,喉道平均马赫数降低,静压升高,总压恢复系数降低。     

非结构混合网格鲁棒自适应技术

数值格式、湍流模型和计算网格是影响CFD数值模拟精度的3个主要因素。结合流场信息的网格自适应技术具备动态优化计算网格的能力，被NASA列为未来CFD发展的一项关键技术。本文针对非结构混合网格，发展了网格单元分布优化、表面网格几何投影和空间网格协调匹配3项关键技术，建立了高鲁棒性几何保真的网格自适应系统。首先，为了提高自适应方法的鲁棒性和通用性，发展了基于标准面网格的多面体网格单元分布优化方法。其次，发展了仅依赖表面网格信息的局部曲面重构技术，采用参数点映射方法实现了新增表面网格点的几何投影，消除了自适应系统对几何CAD系统的依赖。再次，采用改进的距离函数方法实现了空间网格与投影后表面网格的快速匹配。最后，结合基于流场特征的自适应探测器，采用二阶格式的有限体积方法，开展了30P30N三段翼绕流和三角翼大迎角绕流的网格自适应数值模拟。结果表明，通过网格自适应对网格单元的分布进行优化后，流场求解的收敛性和模拟精度都得到了显著提高。     

分离涡模拟方法在内埋弹舱机弹分离模拟中的应用

为了研究内埋弹舱非定常流场对武器投放分离特性的影响,基于分离涡模拟（DES）方法和非结构重叠网格技术,建立了内埋武器机弹分离高精度数值模拟方法。在此基础上,首先通过典型空腔流动和机弹分离算例验证了数值方法的有效性,然后采用平均雷诺应力模拟（RANS）和DES两种方法探究了机弹分离过程中弹舱腔内涡结构的演化历程,分析了两种方法在机弹分离流场结构模拟中的差异,研究了内埋弹舱流场结构对导弹分离特性的影响。计算结果表明：机弹分离过程中,高强度涡结构逐渐向弹舱腔体后缘聚集,同时弹舱剪切层流动结构被导弹破坏,导弹出弹舱阶段俯仰方向所受气动力矩急剧变化。相比RANS方法,DES方法捕捉到的高强度小尺度涡结构对导弹压力分布影响较大,使得分离后期导弹俯仰角、俯仰角速度产生明显差异,DES方法更适合内埋武器舱分离特性精细化研究。该研究结果可为内埋武器机弹分离研究的计算方法选择和安全分离策略提供一定的参考依据。     

基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用

为改善高阶间断Galerkin有限元方法(DG)时间推进效率,在三维非结构网格下针对该方法建立了并行广义最小残差(Generalized Minimal Residual,GMRES)隐式时间迭代方法,GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进一步提高GMRES的计算效率,发展了方程组右端项残值雅可比精确计算方法,针对无黏通量Roe格式和黏性通量BR2(Bassi Rebay 2)黏性计算方法,分别解析给出其对守恒变量多项式自由度的雅可比矩阵。基于建立的方法首先采用NACA0012翼型研究了GMRES的重启次数及收敛参数对方法收敛性影响,然后采用无黏及黏性算例对比研究了基于雅可比矩阵不同计算方法的GMRES计算效率,同时对比研究了雅可比矩阵完全近似求解下GMRES和LU-SGS(Lower Upper-Symmetric Gauss-Seidel)的计算效率。结果表明,建立的基于右端项残值雅可比矩阵精确求解的GMRES方法能够大幅提高不同精度DG方法的CFL(CourantFriedrichs-Lewy)数,相比前面提到的其它方法具有更高的计算效率,其收敛速度实现量级以上的提高。     

激波与氦气泡相互作用的实验与数值研究

提出了一种在激波管中实现运动激波与球形气-气界面相互作用的研究方案.以肥皂膜作为气体交界面,高压火花为光源,采用阴影法,在马赫数M=1.2的情况下,在矩形激波管中开展了一系列的实验研究,获得了激波作用下球形氦气泡界面发展变化的相关实验结果.还通过求解二维轴对称的Euler方程,结合网格自适应策略,对实验过程进行了数值模拟.通过与实验结果的对比分析,证实了涡量的产生和分布对于界面的变形和发展有着重要意义.     

火箭子级垂直回收布局气动特性及发动机喷管影响

垂直回收可重复使用运载火箭是运载火箭发展的一个重要方向,大长细比的火箭子级垂直再入过程属于典型的非规则钝头体绕流,与传统低阻力流线体飞行器气动特性差异较大。采用风洞试验辅以数值仿真分析的方法,对基于栅格舵的火箭子级垂直回收构型基本气动特性和非规则钝头体绕流情况进行了研究,获得了发动机外露喷管和栅格舵对火箭子级垂直回收气动特性的影响规律,给出了火箭子级垂直回收布局设计建议。结果表明：火箭子级倒飞状态下肩部区域会在小迎角下产生大分离流动,外露发动机喷管左右两侧诱导出较强的分离涡结构,与火箭尾翼、肩部大分离流动相互作用;垂直回收构型在超声速下阻力会一直处于较高的水平,不同马赫数下压心移动量较大,倒飞时发动机外露喷管会产生较大的干扰静不稳定力矩,其量值与栅格舵提供的静稳定控制力矩基本相当,在火箭子级垂直回收方案设计时需要引起注意。     

流体二次引射推力转向参数影响规律

采用数值模拟手段,对流体二次引射推力转向参数影响规律进行了研究。首先采用二元矩形矢量喷管,结合风洞试验及国外文献计算数据,验证了自主开发的流体推力转向数值模拟软件的可靠性;在此基础上,开展了流体二次引射推力转向的机理研究和各种参数影响规律数值模拟研究,并详细研究了不同主次流压比、引射缝隙位置和缝隙宽度等参数对干扰流场结构及推力转向偏角的影响,获得了各设计变量对喷管性能及内部流态的影响规律,给出了流体二次引射实现推力转向的基本设计原则及较优的参数组合方案,相关结论可为流体二次流引射推力矢量喷管设计提供依据。     

激波与氦气泡相互作用的实验与数值研究

提出了一种在激波管中实现运动激波与球形气一气界面相互作用的研究方案。以肥皂膜作为气体交界面,高压火花为光源,采用阴影法,在马赫数M=1．2的情况下,在矩形激波管中开展了一系列的实验研究,获得了激波作用下球形氦气泡界面发展变化的相关实验结果。还通过求解二维轴对称的Euler方程,结合网格自适应策略,对实验过程进行了数值模拟。通过与实验结果的对比分析,证实了涡量的产生和分布对于界面的变形和发展有着重要意义。