栅格翼绕流特性的实验研究

研究了栅格翼的绕流特性。采用氢气泡法和丝线法在水洞和风洞中进行实验 ,显示了栅格翼的自由涡系 ,迎角 0°～ 40°范围内各网孔内的流态。还对栅格翼绕流的基本特性以及流动机理作了讨论和分析。     

非结构网格生成技术

描述了非结构网格生成的过程,网格的生成方法基于Delaunay 三角形方法,并对物面附近网格进行拉伸以便获得在边界层内高拉伸比的三角形单元。网格的生成过程分两部分:(1) 按Delaunay 三角形方法生成适于无粘情况下的非结构网格。(2) 对近物面的部分三角形单元进行拉伸,使生成的网格能适用于粘流。使用上述方法本文构造了二维弯管内绕叶栅区域的非结构网格     

一种新的非结构动网格的生成方法

弹簧方法是一种常被采用的非结构动网格生成方法,文中采用一种新的非结构动网格生成方法即弹性体方法(Elastic solid method,ESM)。和弹簧方法相比,虽然弹性体方法的计算量大,但是可以处理大位移的情况,尤其善于处理网格扭转很厉害的情况。另外,采用弹性体方法可以减小网格重新生成的频率,从而可以减小映射误差,同时也可以减小网格重新生成带来的计算量。将两种方法的数值模拟结果进行对比,进一步表明该弹性体方法在处理动网格问题中的优势。     

基于CFD方法的主动襟翼控制旋翼翼型涡特性研究

基于计算流体力学(Computational fluid dynamics,CFD)技术,建立了一个适用于旋翼二维主动襟翼控制(Active flap control,AFC)数值模拟的方法。在满足对翼型参数化分析的前提下,使用Euler方程求解以提高计算速度,并采用嵌套网格方法对AFC旋翼后缘襟翼进行运动控制。应用所建立的方法,首先进行了算例验证计算,然后着重对AFC旋翼翼型进行了数值模拟。在此基础上,进一步开展了AFC旋翼翼型主要参数对后缘涡影响的计算分析。结果表明:提高桨尖马赫数、增加后缘小翼摆动频率能加快涡产生速度;而提高桨尖马赫数、增大后缘小翼摆动幅度和后缘小翼长度能增大涡的强度;但增大后缘小翼与主桨叶缝隙间距仅在一定范围内能够增加涡强度。     

全尺寸栅格舵铰链力矩天平研制

为准确获取带控制系统全尺寸栅格舵的静、动态铰链力矩数据,分析舵控系统性能,需要研制强度、刚度、灵敏度及抗电磁干扰能力均符合试验要求的风洞侧壁支撑天平.栅格舵法向力大,铰链力矩小,且试验平台限制了天平结构的长度,导致天平在载荷匹配与元件布置方面有较大难度.通过有限元分析技术优化天平结构,在适量放大铰链力矩载荷设计指标的基础上,采用法兰盘代替传统的锥连接,天平各元件串联布局等手段,成功研制出了满足试验要求的天平.试验结果显示,模型静、动态试验数据规律良好,天平的动态响应能力及对小量铰链力矩的测量满足试验需求.     

网格加筋壳结构局部受热轴压承载能力分析

弹体结构设计中越来越多地采用网格加筋壳结构,结构局部受热将降低其承载能力。本文首先采用解析方法分析了受热平板的屈曲失效,在此基础上,对网格加筋壳体在轴压和局部热流作用下的承载能力进行数值模拟,表明材料的热膨胀系数是影响结构失稳载荷的主要因素,分析了局部热流对结构局部失稳和整体失稳承载力的影响,其结果对弹体结构局部受热或抗激光设计具有参考意义。     

一种包含运动边界的高精度流场数值计算方法

为了准确快速地模拟运动边界的流场,提出一种反馈力源形式的包含运动边界的非定常流场数值计算方法。该方法采用完全正交的网格,以反馈力源作用点的运动来模拟边界的运动。采用物理量及其各阶导数在边界两侧的突跃修正中心差分格式,使之达到二阶精度,以此离散求解二维不可压Navier-Stokes方程。并且提出了与运动边界相适应的反馈力源构造方法及对边界上速度进行插值的方法。基于此数值计算方法,对低雷诺数的圆柱绕流、静止流体中的振荡圆柱以及昆虫振翅运动的二维非定常流场进行了数值计算,计算结果与以往的数值及实验结果非常吻合,表明本文方法与Peskin的浸入式边界方法在处理运动边界问题时具有同样的高效率,且精度高于浸入式边界方法。     

高阶精度DG/FV混合方法在二维粘性流动模拟中的推广

DG/FV混合方法因其具有紧致性、易于推广至高阶及相比同阶DGM计算量、存储量小等优点,已成功应用于一维/二维标量方程和Euler方程的求解。在此基础上,将该方法推广于二维三角形/矩形混合网格上的NavierStokes方程数值模拟,将格式形式精度提高至4~5阶。物理量的空间重构及离散使用DG/FV混合重构方法;无粘通量计算采用Roe格式;粘性通量计算采用BR2格式;时间方向离散采用高阶显式R-K方法或隐式方法。利用该方法计算了有解析解的Couette流动问题以验证几种格式的数值精度阶,并计算了层流平板流动和定常、非定常圆柱绕流问题等经典算例。计算结果表明DG/FV混合方法达到了设计的精度阶,在较粗的网格上亦能得到高精度的计算结果;定性分析和数值结果表明相比同阶DG方法单步计算量减少约40%。     

DLR-F6/FX2B翼身组合体构型高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6和DLR-F6_FX2B 2种翼身组合体构型的高阶精度数值模拟,计算外形来自AIAA第三届阻力预测研讨会。主要目的是确认WCNS模拟跨声速典型运输机构型和预测局部构型变化引起的气动特性变化量的能力。在固定升力系数条件下,采用粗、中、细3套网格开展了网格收敛性研究,从气动力系数、压力系数分布、表面流态等方面研究了网格规模对DLR-F6和DLR-F6_FX2B翼身组合体数值模拟结果的影响;采用中等网格开展了来流迎角对2种翼身组合体气动特性的影响研究。通过与National Transonic Facility(NTF)的试验结果和CFL3D的计算结果对比,表明采用高阶精度计算方法得到了网格收敛的数值模拟结果,较好地模拟了DLR-F6翼身组合体局部修型引起的微小气动特性变化和翼身结合部流动特性的差异。     

离子推力器栅极系统电子反流阈值的数值分析

阻止束流等离子体中电子反流到加速栅上游区域是离子推力器加速栅负电压的主要作用之一,能够阻止电子反流的加速栅电压最小值称为电子反流阈值。加速栅电压的选择直接影响到离子推力器的工作性能和运行寿命,电子反流阈值电压是确定加速栅电压的重要参考参数。基于PIC方法计算了20cm氙离子推力器加速栅电子反流阈值,并分析了加速栅孔径、栅间距、单孔引出束流电流大小对加速栅电子反流阈值电压的影响,计算结果与试验测量值符合较好。该数值模型为加速栅参数的选择和降低电子反流失效风险方法提供了参考,为下一步电子反流现象对加速栅寿命的预测分析奠定了基础。