喷管结构形式对两级PDE性能的影响分析

为了研究喷管结构形式对两级 PDE 性能的影响及作用规律,以氢气和空气混合物为例,采用 FLUENT 软件中的 k-ε湍流模型,使用非平衡壁面函数、PISO 算法及基于梯度的动态自适应网格加密方法,对收敛、扩张、收扩等不同喷管结构形式的凹面腔内环形向心射流聚心碰撞产生激波会聚起爆爆震波的过程进行了数值模拟,并对两级脉冲爆震发动机产生的推力和冲量进行了计算。通过对不同时间点凹面腔产生的推力和冲量、喷管产生的推力和冲量、两级 PDE 总的推力和冲量三者之间的对比和分析,发现凹面腔产生的推力和冲量是两级 PDE 总的推力和冲量的主要来源,不同喷管结构形式对两级 PDE 推力和冲量性能的影响有较大的差别。以冲量为例,收敛喷管可以提高凹面腔的冲量,但是由于自身同时产生过多负的喷管冲量,因而总冲量较小;扩张喷管虽然对凹面腔冲量的提高作用不大,但是由于其自身可以提供正的喷管冲量,所以总冲量较大;收扩喷管的收敛段也可以较小幅度的提高凹面腔冲量,但同时产生的负喷管冲量会和扩张段产生的正喷管冲量发生抵消,使总的喷管冲量比较低,发动机总冲量也较小。综合考虑各方面因素,扩张喷管性能最优。可见在两级 PDE 喷管结构形式的选择上,要同时考虑凹面腔冲量和喷管冲量,权衡两者,才能选择出合适的喷管。     

倾转旋翼典型飞行状态气动特性的CFD分析

针对倾转旋翼特有的悬停、巡航、过渡3种飞行模式,基于CFD方法分别对倾转旋翼的上述3种典型飞行状态进行了数值模拟.首先,基于运动嵌套网格思想建立了适合倾转旋翼运动模拟的多层运动嵌套网格系统.然后,基于RANS(Reynolds averaged Navier-Stokes)方程建立了适用于非定常流场分析的CFD模拟方法.为了提高整个流场的计算效率,时间推进采用高效的隐式LU-SGS(lower upper symmetric Gauss Seidel)方法,并运用SPMD(single program multiple data)并行计算模式.最后,采用所建立的网格生成方法和流场求解器,进行了不同飞行状态下的倾转旋翼气动特性数值分析,验证了建立CFD方法的有效性并得出了多种飞行状态下的气动特性.结果表明:悬停状态下倾转旋翼相对于直升机旋翼有更加明显的下洗流影响区域.巡航状态两旋翼中心间距为2.5R时,倾转旋翼之间存在明显的气动干扰.过渡状态下,倾转旋翼提供的拉力在倾转角为45°前下降较快,且在倾转角较小时,倾转旋翼桨叶表面拉力分布仍能表现出直升机旋翼的特性.      

NASA共同研究模型翼身组合体数值模拟研究

采用结构网格求解器WiseManPlus软件,对DPW5提供的NASA共同研究模型CRM翼身组合体构型进行数值模拟研究.采用雷诺平均N-S方程,选择SA一方程湍流模型和SST两方程湍流模型,开展了网格收敛性分析及抖振分析计算研究,评估了该软件对此类民机构型的阻力预测能力.按照DPW5的要求,采用会议统一提供的基准结构网格进行计算分析,研究结果表明,计算的巡航设计点的阻力及大攻角下的气动力系数与试验值吻合较好,达到会议统计分析平均水平.计算结果表明,所采用的结构网格求解器WiseManPlus软件计算精度较高,适用于大型民机高速巡航及抖振现象的数值模拟研究.     

一种改进的非定常激波装配算法

基于非结构动网格技术和格心型有限体积方法,提出一种改进的非定常激波装配算法,进一步拓展了其在包含有运动激波的非定常流场的应用范围。首先,针对激波在直/曲壁面传播这类问题,分别建立了壁面间断节点的运动模型;其次,为保证激波在大范围运动时装配阵面不产生失真,基于Bézier曲线拟合方法实现了间断节点分布的自动重构;接着,通过嵌入局部网格自动重构模块,提高了算法的计算效率和自动化程度;最后,对于激波相交点的运动,设计了一种根据位移推算速度的方法进行装配。数值算例表明,所提算法能够有效地处理激波传播问题,相比激波捕捉方法可以提取更多的流场信息,同时可以获得流场间断更加直观清晰的图谱。     

直升机前飞旋翼桨尖涡定常质量射流控制计算

针对直升机前飞状态开展了旋翼非定常桨尖涡的模拟,采用定常桨尖空气质量射流(TAMI)控制方式对旋翼桨尖涡进行了控制。建立了一个适用于前飞旋翼桨尖涡高精度捕捉和质量射流控制模拟的数值方法,在该方法中,采用有限体积法进行空间离散,应用5阶Roe-WENO(weighted essentially non-oscillatory)格式进行流场重构及控制面对流通量计算;采用双时间方法进行时间推进,伪时间步上应用隐式LU-SGS(lower upper symmetric Gauss-Seidal)格式;引入射流边界条件对质量射流进行模拟;采用运动嵌套网格方法对前飞旋翼桨尖的挥舞运动进行模拟,并对桨叶网格和背景网格进行针对性加密。基于所建立的方法对前飞状态旋翼非定常桨尖涡及其质量射流控制进行了模拟,计算结果表明:前飞状态下旋翼桨尖涡存在较大的前后差异,桨盘前侧的桨尖涡涡核强度远小于桨盘后侧;桨盘前侧旋翼桨尖涡的涡核强度很难由定常质量射流控制来减弱,而桨盘后侧的旋翼桨尖涡则可以通过定常质量射流得到有效控制。     

结合等值面绘制与体绘制的电磁环境可视化方法

为提高现有的电磁环境三维可视化效果,提出一种将等值面绘制融合到光线投射体绘制中的网格投影算法。将等值面网格化,利用空间网格表现等值面并融合进体绘制生成的电磁环境图像中。将平面均匀网格平行投影到电磁环境体数据中具有同一数据值的数据采样点上,由此生成网格化的等值面。投影到体数据中的网格沿投射光线的反方向确定其对应于光线投射生成图像上的像素点,并将网格颜色值融合到此像素点中,实现融合绘制电磁环境。在统一计算设备架构(CUDA)下并行实现本文算法后,不仅电磁环境的绘制效果比融合前有提升,而且更新绘制速度能够达到实时,支持用户实时交互控制。     

数据驱动的曲面构件形状-拓扑协同优化方法

针对紧凑设计空间下曲面构件极致轻量化需求,提出一种数据驱动的曲面构件形状-拓扑协同优化方法,主要包括离线、在线和更新3个阶段。离线阶段中,首先通过拉丁超立方采样方法在设计空间内采样,并利用网格变形技术进行参数化建模,获得样本点对应的网格模型。然后对上述网格模型分别进行拓扑优化,获得优化后的应变能。基于上述步骤获得的样本数据,训练径向基函数代理模型,其中形状设计变量为输入,拓扑优化获得的应变能为输出。在线阶段中,将基于离线阶段获得的代理模型开展优化设计,采用协方差矩阵自适应进化策略来提高优化效率。更新阶段中,计算代理模型优化结果的真实响应,并将其加入样本数据集进行代理模型更新。最后,通过简支梁和航天器舱门算例开展算法验证。结果表明相比针对固定形状获得的拓扑优化结果,提出方法优化获得的应变能结果分别降低了20.08%和37.93%,表明提出方法具有更优的设计能力。     

基于网格的分布式仿真技术

介绍了网格技术及其国内外研究情况、基于网格的分布式仿真的关键技术 ,包括 :体系结构、自适应网格动态变化技术、负载平衡技术、通信质量控制技术。     

基于混合网格上的高温非平衡流的数值模拟

新一代再入飞行器及空间传输系统需要了解其飞行时的气动问题及热力学问题。数值模拟高空、低密度、高熵及非平衡流动是一具有挑战性的问题。本文针对带座舱飞船高超声速再入大气过程中存在的严重气动加热现象,利用混合网格及Osher逆风格式,数值求解了三维化学非平衡Navier-Stokes方程,其中化学模型是5组分1 7个化学反应的空气化学模型,对带座舱飞船再入高度为40 km和马赫数为2 0 ,1 0的化学非平衡流场进行了数值模拟,给出了迎角为0°和2 0°情况下的各个组分的密度分布、压力等参数,并与量完全气体的计算结果进行了比较。