涡轮机械叶片的流固耦合数值计算方法

发展了一种流固弱耦合数值方法来判断涡轮机械叶片气弹稳定性.通过一种数据交换方法将计算结构动力学(CSD)的节点振动位移施加到计算流体力学(CFD)的叶片表面网格点上,CFD分析采用多层动网格技术实时更新可动域的网格点坐标,并通过有限体积法求解了用k-ε湍流模型封闭的雷诺平均Navier-Stokes(RANS)方程.以NASA 67转子叶片为例,在设计状态下,对比了叶片前4阶模态下一个振动周期内的气动功及模态气动阻尼比,初步分析了叶片的气弹稳定性.计算结果表明:所发展的流固弱耦合数值方法用于判断涡轮机械叶片气弹稳定性是可行的.      

电荷交换离子对栅极系统束流影响的数值研究

采用二维网格质点法(PIC)计算离子在离子发动机栅极系统中的运动, 通过在模型中添加离子和中性粒子电荷交换的Monte Carlo碰撞模块, 得到了电荷交换离子在栅极周围的分布及电荷交换离子的运动规律.计算结果表明:考虑电荷交换离子后, 屏栅极电流较不考虑电荷交换离子情况时增大了1.42%, 所受影响不大, 加速栅极电流由0增大到主束流电流的1.41%.模拟结果表明:加速栅极下游较远处产生的电荷交换离子, 是造成加速栅极下游面腐蚀及加速栅极电流的主要原因.      

高效确定重叠网格对应关系的距离减缩法及其应用

提出了“距离减缩法”,可以通过搜索有限几个数目很少的点来寻找到目标单元,从而可以快速确定重叠网格之间的对应关系,应用于重叠网格洞点识别与插值单元寻址两个方面,并且还可以结合Inverse Map方法高效建立网格之间的映射关系。该方法将三维搜索化为准一维搜索,大幅度提高搜索效率,且算法通用性强、易于实现。应用表明：利用该方法可以迅速识别洞点,与等参变换结合可以有方向地快速锁定插值单元,找点效率高,适用于计算流体力学(CFD)三维非定常计算中有相对运动的重叠网格系统在每一个时间站位都需要重新进行洞点识别与插值单元寻址的计算,可以有效节省计算时间,具有工程实用价值。     

涡轮叶片叶面B-spline曲面层次化拟合方法

航空发动机涡轮叶片表面形状复杂,其三维重构是叶片成型精度分析以及质量控制的重要手段。涡轮叶片叶面是一个复杂的自由曲面,其拟合是重构涡轮叶片的关键步骤之一。本文以涡轮叶片叶面拟合为研究对象,提出一种B-spline曲面层次化拟合方法,该方法首先构建网格曲面的多分辨率模型,然后完成该多分辨率三角网格模型的参数化,最后在多分辨率三角网格的基础上通过粗略拟合、精确拟合以及最终拟合3个主要步骤完成曲面拟合。该方法可以完成高阶连续自由曲面的拟合,适用于海量测量数据的快速拟合,具有较高效率并且拟合结果精度高。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

网格计算技术在大规模CFD计算上的应用

针对大规模数值计算规模化管理的难题,采用网格计算技术将网络环境下分散的计算资源组建了CFD网格计算环境.采用主动接入方式,以中心服务器为核心的网格计算管理和服务模式,构建CFD网格计算环境,实现了CFD计算在网格环境内的自动化计算管理,提高了计算资源的利用效率,降低了设计人员的劳动强度.     

基于动态边界控制的优化方法在高升阻比气动外形优化中的初步应用

本文是利用动网格技术和有限体积法来研究基于动态边界控制的气动优化的方法.该优化方法是通过在求解非定常流动的过程中计算优化的目标函数随设计变量变化的梯度.而利用参数化的结构动网格技术能有效、快速地生成与优化过程出现的外形对应的网格并进行计算,具有较高的效率.本文以典型的高升阻比外形优化为例,用非定常NS方程为流动控制方程进行优化,取得了初步的结果.     

基于准定常的飞行器动导数的高效计算方法

结合工程设计经验和工程计算方法,针对飞行器动导数计算方法展开研究。基于求解准定常欧拉方程发展了一种适合于亚跨超声速状态下,任意气动外形气动导数的高效计算方法,并将计算结果与一般非定常方法的求解结果以及国际标准实验数据进行了对比,研究结果表明该方法结果可靠、效率较高。     

利用CFD计算技术进行机翼模型气动特性分析

计算流体动力学方法（CFD）已经成为空气动力学领域重要的计算方法。利用该方法进行机翼模型Naca0012的数值模拟,包括计算模型的气动力系数以及气动导数。采用了效果很好的新型网格控制算法,通过自有程序处理流固耦合边界问题以及商用软件的求解器模块进行整体区域的数值计算。在证明了方法的有效性之后,又通过多种振幅奈件下的数值模拟结果,对以往广泛采用的理论方法的有效性进行了探讨,研究了机翼模型气动导数随振幅变化而变化的情况。     

高雷诺数下的翼型绕流LBM数值模拟

采用格子Boltzmann方法(LBM)的二维9速度(D2Q9)模型和贴体网格,通过引入非均匀网格插值方法和非平衡态外推边界处理,分别结合Baldwin-Lomax(B-L)湍流模型和Spalart-Allmaras(S-A)湍流模型,对高雷诺数Re≥5×105下的NACA0012翼型绕流进行了数值模拟和对比研究,两者的结果与CFL3D的结果和实验结果均吻合的很好,相比之下,采用S-A模型能更好地预测失速迎角,其处理分离流动的能力要强于B-L模型。改进后的LBM适用于非均匀贴体网格,曲边边界,计算简单,并可应用于更复杂的高雷诺数流动中。