基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

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

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

液体火箭燃烧尾焰冲击干扰效应数值研究

火箭发射时其燃烧尾焰的冲击干扰效应对发射稳定性和发射架、导流槽等地面设施有重要影响。文章采用压力隐式算子分裂算法,通过求解Navier-Stokes方程,对氢氧液体火箭发动机燃烧室内燃烧过程与尾焰流场进行了一体化数值计算,得到了火箭发射后尾焰与地面撞击产生的冲击流场。结果表明:尾焰流场计算模型、方法与结果合理;尾焰冲击干扰效应会大幅度提高地面附近的压力和温度;火箭尾焰撞击地面后,高温区出现在离地面一定距离的高温层内,此时地面附近为低速区;尾焰对其正下面的地面区域产生冲击最大,主要干扰区域集中于半径为15 m的圆形区域。     

C/C-SiC复合材料致密度影响因素

利用先驱体转化法制备C/C-SiC复合材料,对试样进行微观结构分析和性能测试,研究渗硅温度、保温时间、真空度和裂解周期对C/C-SiC复合材料致密度的影响。结果表明:随着渗硅温度的升高,材料的致密度呈先加速升高后快速下降趋势;随着保温时间的延长,材料的致密度先快速升高,保持一段时间稳定后再缓慢降低;随着烧结真空度的提高,材料的致密度加速升高;随着裂解周期的增加,材料的致密度不断增大,但增速逐步降低。经过11周期的“浸渍-固化-裂解”过程后,所制备的C/C-SiC复合材料获得最大密度2.09 g/cm3、最小孔隙率7.6%,其综合力学性能最为优异:弯曲强度468 MPa、拉伸强度242 MPa、断裂韧度19.6 MPa?m1/2、维氏硬度17.2 GPa。     

某型飞机前起落架收放载荷当量化研究

飞机起落架的收放大部分是在飞机飞行时进行的,起落架主要承受着飞行时的气动载荷、质量力和惯性载荷,这些载荷的大小或方向随着飞机的飞行速度和起落架的收放不断发生变化。在地面进行起落架收放系统可靠性试验时,为真实反映起落架收放时的收放载荷,施加多大的收放载荷以及如何施加收放载荷成为起落架收放系统可靠性试验中的关键技术。本文对某型飞机起落架收放载荷进行了研究,提出了起落架收放载荷当量化处理的一种方法,并采用动力学软件对当量化结果进行了模拟分析,分析结果与飞行实测结果十分吻合,而且该当量化方法简易可行,便于在起落架收放系统可靠性试验中施加载荷。     

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。     

飞行器外挂投放数值模拟

采用非结构重叠网格Euler方程的CFD方法数值模拟了外挂投放。该方法在重叠网格"挖洞"时,以"离壁面距离"为挖洞判据;在建立重叠区域插值关系时,对Neighbor to Neighbor(N2N)和Alternating Digital Tree(ADT)两种方法的效率进行了比较。采用机翼/挂架/带舵外挂物模型开展了验证计算。在此基础上,完成了某飞行器外挂投放的初步计算。     

燃烧室壁冲击-逆向对流-气膜冷却特性的数值研究

采用FLUENT软件对冲击-逆向对流-气膜冷却结构进行了流固耦合计算, 得到了对流通道和气膜形成区流场、流体温度场以及固壁温度场分布, 分析了该冷却结构内部复杂的流动和换热情况.计算研究目的是摸索一种火焰筒壁温的计算方法, 以较好地获得火焰筒壁温分布.计算结果和试验结果进行了对比, 两者有较好的吻合性.分析表明采用流固耦合计算获取冲击-逆向对流-气膜冷却结构的壁温分布是可行的.      

多斜孔壁气膜冷却绝热温比传热传质类比实验研究

采用传热传质类比的方法,对多斜孔壁气膜冷却的绝热温比进行了实验研究。多斜孔为规则的叉排形式,小孔与平板表面在主流方向上夹角α均为３０°。在吹风比为１．０～４．０的范围内,研究了具有不同展向间距比ｐ／ｄ,流向排距比ｓ／ｄ的实验板绝热温比分布;对实验数据进行分析比较,得出了孔距与排距的最佳比值,以及吹风比对绝热温比的影响。     

非结构混合网格自适应并行技术

计算流体力学（CFD）模拟实际工程问题所采用的网格规模可达千万量级，并行技术是减少计算时间的有效方法。耦合流场信息的网格自适应技术能有效动态优化计算网格，被NASA视为一项亟待发展的CFD关键技术。混合网格自适应系统包含网格分布优化、表面网格投影和空间网格匹配等关键技术。针对以上3项关键技术分别建立了高效的并行算法。首先，提出了"先唯一后同一"的两步法策略实现了网格单元分布优化过程的并行相容性；其次，基于局部曲面拟合思想，实现了曲面重构和新增物理网格点投影的完全并行；再次，提出了空间网格匹配技术的半并行算法，快速解决了网格单元交错问题。为了提高后续流场计算的并行效率，发展了基于并行重分区-网格数据迁移方法的动态负载平衡技术，并采用圆柱激波流场自适应模拟对动态负载平衡技术进行初步验证。最后，采用三角翼自适应加密测试了自适应系统的并行效率。结果表明，建立的混合网格自适应系统并行效率较高，且相比流场求解耗费总时间的比例低于1%。