用动态非结构重叠网格法模拟三维多体相对运动绕流

 发展了动态非结构重叠网格技术。本文提出了基于阵面推进的相邻单元搜索算法和适合于高阶空间离散格式的网格间插值关系定义方法，前者提高了非结构重叠网格在模拟多体间具有大幅相对运动的非定常问题时网格处理的效率，而后者保证了网格间边界附近流场求解精度；本文非定常Euler方程的求解采用有限体积法，通量计算采用Osher格式，时间积分采用双时间步长推进。为验证本文提出的方法，分别对定常M6机翼绕流和机翼下外挂物投放过程的非定常绕流进行了模拟,结果表明本文提出的方法提高了非结构重叠网格在模拟多体间具有大幅相对运动问题时的效率，保证了网格边界处的精度，具有较高的工程应用价值。     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

旋翼/机身干扰非定常流场数值模拟(英文)

基于非结构动态嵌套网格方法对直升机旋翼/机身干扰非定常流场进行了三维Euler方程数值模拟。嵌套网格系统由包含旋翼的运动区网格和包含机身的静止区网格组成,两区网格均来自于同一套初始的整体非结构网格,人工边界以及初始的插值贡献单元对应关系在网格分区时同时确定。在非定常计算中两区网格始终存在一定的重叠区域,无需重新建立人工边界。高效的贡献单元搜索算法可以根据初始的插值贡献单元对应关系很快更新每一个真实时间步上的贡献单元对应关系。为了准确计入旋翼桨叶的挥舞和变距变形,在运动区网格中引入线性弹簧方法。采用非定常Euler方程和双时间推进方法模拟了美国佐治亚理工学院旋翼/机身干扰实验外形非定常流场,计算结果与实验值以及文献的计算值相吻合。     

基于非结构网格的气动弹性数值方法研究

基于非结构网格和网格变形技术,研究了耦合求解流体和弹性结构体互相作用问题的数值方法.气动弹性模型采用了耦合求解可压缩雷诺平均的Navier-Stokes方程和线性结构动力学方程的方法.为了能处理在复杂几何域内的问题,采用了基于混合单元的非结构网格有限体积方法.开发的方法在二维翼型简谐振动进行了验证并应用于AGARD 445.6翼型颤振边界的数值模拟.      

基于非结构重叠网格的螺旋桨滑流非定常数值模拟

基于各向异性非结构混合网格及多套网格重叠技术,通过求解非定常雷诺平均Navier-Stokes(URANS)方程,分别研究了螺旋桨动力特性和螺旋桨滑流对涡桨飞机气动特性影响。运用物面相交准则实现多套网格间的挖洞,采用距离权和三线性插值技术传递重叠区变量信息,并针对定轴旋转优化了网格重叠边界。采用双时间步求解控制方程,内迭代计算采用LU-SGS(lower-upper symmetric Gauss-Seidel)方法。模拟了单独螺旋桨旋转的非定常流动,计算的拉力系数和扭矩系数与试验结果一致,表明采用的数值方法和网格技术能有效模拟螺旋桨滑流效应。数值模拟了某涡桨飞机,对比了有无滑流下飞机的气动特性,给出了滑流对飞机影响。结果表明:滑流导致其后部气流加速和旋转,改变飞机气动特性,增大飞机的阻力和俯仰力矩,导致气流下洗,影响机翼及平尾压力分布。      

聚合多重网格法在旋翼前飞流场计算中的应用

基于非结构动态嵌套网格, 将聚合多重网格方法引入到旋翼的前飞流场计算中.从初始的整体非结构网格出发, 快速建立相互重叠的两个子域网格, 同时无需建立背景网格便可以确定初始洞边界和初始插值贡献单元, 且动态嵌套时无需更新洞边界, 并采用了一种无需背景网格便可快速查找新的贡献单元的方法.分区利用聚合法生成多重网格, 利用双时间方法分区运用三层网格V循环求解非定常Euler方程, 两区的流场参数在重叠区实现耦合, 对旋翼的前飞流场进行了计算.结果表明, 多重网格方法能够有效地加快旋翼前飞流场收敛.      

基于动态非结构重叠网格法的直升机前飞非定常流场数值模拟研究

 发展了非结构重叠网格技术，提出适合于直升机旋翼运动的动态非结构重叠网格算法，该算法使得在采用非结构重叠网格模拟直升机流场时实现程序自动处理网格，不需任何人工干预。流动控制方程采用非定常可压缩Euler方程，空间离散采用Jameson有限体积法，时间推进采用双时间步长法，为加速收敛，在双时间步长内循环迭代中采用预处理技术。为验证本发展的预处理解算器以及动态非结构重叠网格的算法的可信性和有效性，开展一系列的数值实验，模拟了不同速度范围的NACA0012翼型绕流,Caradonna旋翼悬停和前飞流场以及Robin直升机机身-旋翼组合的前飞非定常流场，结果表明采用本文算法得到的结果是可信的，网格处理是自动高效的，可以应用于直升机全机流场模拟，为直升机设计提供可信指导数据。     

分层嵌套重叠网格自适应树结构动态组装方法

采用重叠网格可以有效地进行复杂流动的大规模数值模拟,特别是包含运动部件(如旋翼、投弹)的动态模拟。本文将树结构的自适应直角网格用于重叠网格组装过程中的切割和贡献单元的搜索,大大加快重叠网格的组装速度。通过二叉树自适应直角网格对物体外形进行离散,实现切割过程的快速定位;采用八叉树自适应直角网格对流场区域进行离散,高效地搜索贡献单元。使用基于壁面距离准则的重叠区域最小化方法和分层嵌套重叠策略,能提高重叠网格组装的效率和质量。对于具有运动部件的动态重叠网格问题,采用多个二/八叉树减少组装过程中信息更新的冗余计算,从而大幅度减少重叠网格组装的时间消耗。实际算例的重叠网格组装结果说明本文发展的重叠网格组装方法具有很高的计算效率,可以满足运动边界复杂流动问题的动态计算要求。     

多段翼型粘性绕流计算研究

针对多段翼型外形结构复杂,结构网格生成难度大的缺点,运用割补法对重叠网格技术进行了研究,生成贴体及与边界正交的高质量的结构化网格,在识别贡献单元方面,介绍了一种基于kd树的找重方法.流场求解应用雷诺平均Navier-Stokes方程结合Spalart-Allmaras一方程湍流模型在所生成的重叠网格上对多段翼型的粘性绕流进行了模拟,为了兼顾流场计算效率问题,在重叠网格的各个子块上采用多重网格方法.计算结果表明,用该方法可以很好地预计压力分布,并且计算效率得到很大的提高.     

基于非结构动网格的跨声速非定常流场求解

对用于非结构动网格生成的弹簧近似方法进行了研究.通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入防扭转、防挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进.采用改进后的方法并通过耦合求解了基于ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程,模拟了做谐和振动的刚性无后掠机翼及某翼身组合体的跨声速绕流.结果表明,改进后的方法可以大大提高网格变形能力和网格质量,且计算结果与实验结果吻合良好.     

共轴式双旋翼悬停流场和气动力的CFD计算

考虑到共轴式双旋翼流场特征高度复杂、桨叶承受非定常气动载荷的特点,为更好地预测共轴式双旋翼的气动特性,把非结构嵌套网格方法和网格的自适应技术相结合,发展了一套适合于共轴式双旋翼流场数值计算的求解器.在该求解器中,采用非结构嵌套网格方案来模拟桨叶之间存在的相对运动,自适应网格技术用来捕捉尾迹对流场和气动特性的影响,求解惯性坐标系下的非定常N-S主控方程来模拟流场的非定常特性.应用该求解器,首先计算了有试验结果可供对比的一副试验旋翼的诱导速度场分布,在此基础上,计算了共轴式双旋翼的桨尖涡轨迹和拉力分布特性,并与单旋翼的计算结果进行了对比和分析,得到了一些有意义的结论.     

基于非结构动网格的非定常气动力计算

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法耦合求解基于Arbitrary Lagrangian-Eulerian(ALE)描述的欧拉方程,计算了谐和振动NACA0012翼型及矩形机翼的跨声速非定常气动力,计算结果与参考文献提供的结果及实验结果吻合良好。     

基于非结构运动对接网格的旋翼前飞流场数值模拟

发展了一种基于线性弹簧理论的非结构运动对接网格方法,通过求解非定常Euler方程模拟了旋翼的前飞流场。利用运动区网格的转动模拟旋翼的旋转运动,并运用线性弹簧动态网格技术处理旋翼的挥舞和变距运动,方程求解中采用一种精度较高的插值和重构方法。算例表明,本文的方法能够有效地模拟旋翼的前飞流场。     

用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

非结构嵌套网格的直升机旋翼/机身前飞流场数值模拟

建立了一个适用于前飞状态直升机单独旋翼、机身以及旋翼/机身组合体流场的数值计算方法及模型.在该方法中,使用非结构运动嵌套网格描述桨叶之间以及桨叶与机身之间的空间位置变化关系,控制方程采用惯性坐标系下的非定常Navier-Stokes(N-S)主控方程,空间方向上使用二阶迎风格式,时间方向上使用Lower-upper symmetric Gauss-Seidel(LU-SGS)格式.对于求解中遇到的旋翼/机身间距离过近所导致的贡献单元生成困难问题提出了一种新的解决方法.应用所建立的方法,分别对单独旋翼、机身及旋翼机身组合体的前飞绕流场进行了数值模拟与分析,计算结果与试验结果的对比表明本文建立的方法对直升机的前飞绕流场的计算和分析是很有效的.      

基于非结构嵌套网格方法的旋翼地面效应数值模拟

 建立了一个基于非结构嵌套网格的流场求解器,用来精确模拟复杂的旋翼近地流场,为更好地分析地面效应(IGE)对旋翼气动特性的影响提供一套计算方法。在该求解器中,控制方程采用惯性坐标系下的非定常N-S方程,空间方向上采用二阶迎风格式,并用背景网格的一个面模拟地面作用,以方便地面边界条件的处理。应用所建立的模型,首先针对有实验结果可供对比的旋翼无地面效应（OGE）流场进行了数值模拟,以验证计算方法。然后着重计算了IGE下流场中的桨尖涡空间位置和旋翼拉力增益,并对旋翼在小速度前飞状态下的近地流场及地面涡形成过程进行了模拟和分析。在此基础上,得出了一些有意义的结论。