烧蚀环境下的圆柱绕流计算模型

基于多学科理论建立了高温高压气流环境下圆柱烧蚀、剥蚀的数理模型。利用离散涡方法计算流场与圆柱表面压力分布;采用三方程烧蚀模型计算热化学控制机制下的烧蚀速率,并将其与扩散控制机制下的结果相比较以确定烧蚀量;引入颗粒轨道模型求解剥蚀颗粒的运动。基于该模型对高温高压燃烧室内圆柱形烧蚀试件的绕流场、烧蚀量以及剥蚀颗粒的运动轨迹进行了编程计算分析。研究表明:低雷诺数条件下烧蚀环境对绕流流场的影响不大,而较高雷诺数条件下烧蚀绕流流场与无烧蚀绕流流场的涡街分布存在显著差别;在低雷诺数条件下流场对烧蚀速率的影响甚微,较高雷诺数条件下烧蚀速率略大于无流场烧蚀情况;涡量的存在决定着剥蚀颗粒的运动轨迹与分布,进而影响尾流热场与烧蚀速率。该模型可为热防护材料的烧蚀实验提供参考,并为烧蚀、剥蚀过程的非线性分析提供数据支持。     

基于离散单元法的二维颗粒阻尼研究

针对自由端带空腔的L型悬臂梁建立了二维颗粒阻尼物理模型,发展并利用离散单元法(DEM)对二维颗粒阻尼进行了数值仿真,重点考查了二维空腔尺寸对颗粒阻尼的影响规律:①高度非线性的颗粒阻尼能够显著地抑制结构振动,阻尼值至少比金属材料阻尼大一个数量级;②当无量纲加速度Γ<1时,颗粒阻尼很小;当Γ＞1时,随Γ增加,颗粒阻尼先增大...     

轮体结构颗粒阻尼器设计方法

利用离散元法及正交设计方法,对典型轮体结构发生伞形振动时的二维等效振动模型进行了空腔尺寸及颗粒填充方案的数值设计,并结合振动试验结果给出了轮体结构颗粒阻尼器的设计流程及设计准则.研究表明:1建立的二维等效振动模型在能够反映轮体结构伞形振动基本特征的同时,也可以较为准确地模拟内部颗粒运动对结构产生的影响;2颗粒阻尼对轮体结构的振动有明显的改善作用,但填充不同材质的颗粒,其减振效果会有较为明显的区别;3空腔的尺寸及颗粒的质量率对颗粒阻尼的影响显著,应首先保证的是较高的颗粒质量率;4由于碰撞间隙的作用,对于固定规格的空腔和颗粒,存在最佳颗粒体积填充率使得颗粒阻尼的减振效果最佳.由试验与数值模拟的最佳方案的一致性可知,发展的设计方法可以较为准确地给出轮体结构颗粒阻尼器的最优空腔尺寸及其对应的最佳填充方法,可用于轮体结构伞形振动减振方案的前期设计.     

协同射流在垂直尾翼流动控制中的应用研究

飞机垂尾的气动特性直接关系到垂尾尺寸和质量,通过流动控制方法提高垂尾的升力可以进一步减小垂尾尺寸,对减小结构质量具有实际意义。采用计算流体力学方法,对垂尾的协同射流主动流动控制进行数值模拟;研究射流动量系数对连续协同射流的影响规律,以及射流动量系数、堵塞比、喷口数量等参数对离散协同射流的影响规律。结果表明：对于连续型协...     

机载220 GHz云雷达衰减订正模拟试验

针对220 GHz云雷达回波信号的衰减问题，结合典型云滴谱和云参数，对水云和冰云模拟取样2 100次，应用离散偶极子近似法计算了220 GHz云雷达探测云滴粒子的衰减系数kc和雷达反射率因子Z的关系。在考虑大气和云对220 GHz云雷达信号衰减的基础上，对水云和冰云的机载220 GHz雷达回波强度进行了数值模拟，根据得到的kc-Z关系，进行衰减订正的模拟试验。结果表明:水滴云比冰晶云的衰减系数大2~3个数量级;相对于小粒子的水云和冰云，大气吸收是造成机载220 GHz云雷达衰减的主要原因。高频段的雷达信号在传播路径上的衰减非常明显，本研究中给出的kc-Z关系适用于220 GHz频率雷达的云衰减订正计算。     

基于离散伴随方程的三维雷达散射截面几何敏感度计算

针对有限差分法计算雷达散射截面（RCS）梯度效率低，采用高精度雷达散射截面评估时计算代价高的问题，提出了一种基于麦克斯韦积分方程离散伴随方程的RCS梯度高效计算方法。基于伴随方程的梯度计算可以通过一次雷达散射截面求解、一次伴随方程求解获得RCS关于所有设计变量的梯度。其中麦克斯韦积分方程离散伴随方程的形式与原方程基本一致，可以采用矩量法（MOM）及多层快速多极子算法（MLFMA）求解。伴随方程求解计算量与直接雷达散射截面评估基本一致，存储量在直接雷达散射截面评估的基础上增加不明显。通过双椎体模型、导弹模型对基于矩量法、多层快速多极子算法的伴随梯度进行验证，证明了基于伴随方法的RCS梯度计算可以实现复杂外形中RCS梯度的高效、高精度求解，为基于梯度的高精度气动/隐身一体化优化提供了基础。     

基于APSE确定高超声速边界层第二模态中性曲线下支

在高超声速边界层流动中,由于存在非平行效应以及激发Mack第二模态时的模态同步机制,因此在进行转捩预测研究时,中性曲线下支的确定需要考虑到这两方面的影响。全文针对高超声速平板边界层,首先采用线性抛物化稳定性方程(LPSE)对考虑感受性的中性曲线下支进行了求解,研究发现该方法存在两大难点:一是计算入口需要给定的离散模态在衰减区求解困难;二是由于存在快、慢模态的同步机制,所以在实际计算时,需要分别计算两次离散模态的演化以确定中性曲线下支,计算效率较低。基于以上问题,本文提出了一种基于伴随抛物化稳定性方程(APSE)来计算非平行边界层稳定性的方法,并以此确定考虑感受性的中性曲线下支位置。该方法从Mack模态的增长区开始进行计算,因此有效地解决了使用LPSE方法计算时所存在的问题。在不同工况下,与线性稳定性理论(LST)以及LPSE给出的预测结果进行了对比。从对比结果中可以看出该方法给出的预测结果与LPSE给定上游主导模态的预测结果一致,同时不必依赖于入口初值的选取且效率更高。该方法有进一步发展的潜力,可以有效地应用于高超声速边界层转捩的预测研究中。     

SMP2仿真引擎的多核并行化

仿真模型可移植性规范2（SMP2）是由欧洲航天局制定并得到广泛应用的仿真模型开发与集     

月球探测器发动机羽流辐射热流计算

火箭发动机羽流的辐射的计算需要求解气体介质的辐射传递方程(RTE).计算模型中基于单线组(SLG)模型考虑气体组分的光谱特性,并采用离散传递法(DTM)求解RTE,对波长积分得到辐射热流密度.通过对圆柱高温CO₂气团辐射热流计算来验证计算程序,然后在流场计算的基础上,计算了月球探测器主发动机和姿控发动机羽流对各部位和仪器的辐射热流密度,根据对不同工况计算结果的分析为探测器的设计提供参考.      

多点喷射模型燃烧室性能CFD分析

采用Fluent的雷诺应力模型、非预混燃烧和非预混平衡化学反应PDF模型、压力-旋流雾化喷嘴模型、P1辐射模型,对多点喷射模型燃烧室的燃烧性能进行了计算。分析了不同余气系数的影响及旋流器下游各截面的温度分布,结果表明,数值计算能够模拟大部分流场结构,余气系数对回流区几乎没有影响,但对主燃区的峰值温度位置有明显作用;对NOx分布规律进行了研究,结果表明,NO的分布与温度分布类似,高温区对应高NO质量分数区;随着流动向下游发展,NO分布趋于均匀。