局部网格加密方法计算不可压N—S方程

本文发展了二维非定常不可压N-S方程的局部网格加密方法(LMR):需加密区域预先指定;仅对空间方向加密,分粗细二层网格,细网格覆盖在粗网格上;粗和细网格上分别用显式和隐式差分格式。结合压力修正法类的余量型差分格式,恰当地处理了粗细网格之间的信息传递,使得粗细网格交接面上质量守恒,即满足连续性条件。粗网格通过插值给细网格边界值影响细网格;细网格通过粗细网格压力Poisson方程耦合求解影响粗网格,并且由于压力Poisson方程从动量方程形成,在细网格覆盖下的粗网格上,压力Poisson方程是细网格方程的组合,从而耦合求解时粗细网格压力Poisson方程不需迭代。 本文中计算了二维方腔内的自然对流模型问题,Rayleigh数是10~6。粗网格用显式ULWC格式,细网格用隐式余量型近似因式分解格式。Poisson方程用快速直接算法PO1STG(在FISHPACK中)。粗网格距△=1/16,细网格距△=1/64,加密1/4区域,计算到定常局部网格加密方法所需CPU时间比同等均匀网格(△=1/64)节约一倍以上。     

自力式调压阀变开度流场及特性分析

以研究自力式调压阀性能为目的,利用动态网格方法对阀门内流场进行数值模拟。通过对不同入口流量和弹性系数的多个算例分析得到:稳定压力和入口流量呈直线关系;稳定压力、阀门开度和弹性系数呈二次曲线关系。调压阀在入口质量流量变化时,阀门可以很快稳定,实现对压力和流量进行的调节。     

基于Lattice Bohzmann方法的翼型绕流数值模拟

将通过物理平面的结构化C型贴体网格变换到计算平面的均匀正交网格,引入非均匀网格插值方法,采用基于Lattice-Boltzmann方程的D2Q9模型,对低雷诺数下的NACA0012翼型绕流进行了数值模拟,对五个不同位置的边界层速度型与基于有限体积法的CFL3D程序的结果进行了比较.结果表明,两种方法的计算结果一致.改进后的Lattice-Boltzmann方法适用于曲边边界,非均匀网格,同时对计算速度影响不大.Lattice-Boltzmann方法计算较简单,因而用LB方法来模拟翼型绕流是一种新的选择,并可应用于更复杂的低雷诺数流动中.     

基于非结构网格流场计算的网格重排序

提出基于非结构网格流场计算时的网格重排序,优化了网格信息在内存中的存贮位置,使得计算过程中调用数据的处理速度加快,有效提高了计算效率。以M6机翼及DLR-F4翼身组合体的跨音速无粘流场为算例,分别通过显式四步龙格-库塔推进格式和隐式LU-SGS推进格式来验证网格重排序后对计算效率的提高作用。原始网格均用Delaunay方法生长,并分别使用了两种重排序方法进行比较。计算结果表明:重排序后的计算效率可以提高到未排序时的3到5倍,很具实用价值。     

一类新型的机动目标跟踪算法

针对机动目标跟踪,提出了基于截断正态概率模型的改进自适应目标跟踪算法,该算法具有结构和计算简单,鲁棒性好的特点,通过仿真结果对比,充分说明了文中所提出的跟踪算法能够较好地弥补传统的Kalman滤波方法在跟踪机动目标时的不足。     

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

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

基于滑动网格的螺旋桨复杂绕流数值模拟研究

利用SST湍流模型,求解雷诺平均Navier-Stokes方程数值模拟螺旋桨粘性绕流。为了实现桨叶的相对运动,采用了滑动网格技术,即将计算区域分为静止区域和随桨叶一起运动的旋转区域。以课题组设计的1.2 m螺旋桨为例,分别对不同风速下的桨叶粘性流动进行了非定常数值模拟,计算得到的螺旋桨气动特性与实验值吻合良好,证明了方法的正确性。     

航空弧齿锥齿轮承载传动误差的分析与设计

在计算时变啮合刚度的基础上,推导了承载传动误差简化计算公式,并与有限元计算方法进行了比较,表明两种计算方法得到的承载传动误差的幅值、形状及一阶频谱幅值随重合度的变化趋势基本一致.根据承载传动误差简化公式定性的分析了负载扭矩与实际重合度的关系、承载传动误差的变化规律、齿面印痕和几何传动误差对承载传动误差的影响.最后根据分析结果,提出基于承载传动误差的齿面优化设计的优化策略,优化变量为接触迹线的角度、几何传动误差的幅值,优化目标为承载传动误差在较宽的负载范围内波动最小.      

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

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

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

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