气体动理论BGK格式的网格自适应方法

为了提高气体动理论BGK（Bhatnagr-Gross-Krook）格式在超声速流动问题计算时激波捕捉的准确性与计算效率,提出了一种适用于气体动理论BGK格式的网格自适应加密方法。该方法采用基于四边形的链表技术来描述网格的拓扑结构,在物理量重构过程中,使用了在四边形网格中表现优异的van Leer限制器,以保证粗细网格过渡处物理量重构的精度。用跨声速翼型绕流（马赫数Ma=0.85）、超声速前台阶流（Ma=3）和高超声速圆柱绕流（Ma=8.03）等多个典型算例验证了BGK自适应网格方法。计算结果表明,自适应网格BGK方法在保证数值精度的前提下,可大幅度提高计算效率。这为该方法用于高效地解决复杂问题提供了一种选择。     

Interface flux reconstruction method based on optimized weight essentially non-oscillatory scheme

Aimed at the computational aeroacoustics multi-scale problem of complex configurations discretized with multi-size mesh, the flux reconstruction method based on modified Weight Essentially Non-Oscillatory (WENO) scheme is proposed at the interfaces of multi-block grids. With the idea of Dispersion-Relation-Preserving (DRP) scheme, different weight coefficients are obtained by optimization, so that it is in WENO schemes with various characteristics of dispersion and dissipation. On the basis, hybrid flux vector splitting method is utilized to intelligently judge the amplitude of the gap between grid interfaces. After the simulation and analysis of 1D convection equation with different initial conditions, modified WENO scheme is proved to be able to independently distinguish the gap amplitude and generate corresponding dissipation according to the grid resolution. Using the idea of flux reconstruction at grid interfaces, modified WENO scheme with increasing dissipation is applied at grid points, while DRP scheme with low dispersion and dissipation is applied at the inner part of grids. Moreover, Gauss impulse spread and periodic point sound source flow among three cylinders with multi-scale grids are carried out. The results show that the flux reconstruction method at grid interfaces is capable of dealing with Computational AeroAcoustics (CAA) multi-scale problems.     

大型索网式可展开天线在轨高精度指向控制

以星载大型索网式可展开天线为对象,研究在轨工作过程中天线视轴高精度高稳定度指向控制方法。首先,基于Craig-Bampton法建立表征天线指向的动力学模型;然后针对天线振动对视轴指向的影响,设计PD反馈+陷波滤波器的控制算法;进一步分析影响指向精度的因素,设计带有调节滤波器的改进控制器以减小天线视轴指向的周期性误差,提高指向精度与稳定度。最后,给出数值仿真校验的结果并与现有的方法进行比较。结果表明,所提出的控制算法在保证天线指向稳定的同时,可以有效减小天线视轴的指向偏差,并对干扰与振动具有良好的鲁棒性。     

滚转状态下卷弧翼火箭弹气动特性的数值模拟

采用滑移网格技术对卷弧翼火箭弹的滚转特性进行数值模拟,并对模拟所得的马格努斯力与风洞实验数据进行分析比较。误差分析证明了滑移网格技术模拟火箭弹滚转特性的准确性;模拟结果与风洞实验数据对比表明,采用细长体外形效果显著,升阻力系数受滚转特性影响较小,因而在本次模拟的火箭弹滚速下,可不考虑滚转来模拟火箭弹的升力系数和阻力系数,得出火箭弹滚转运动对气动参数的影响。     

三维凸域中候选卫星点集的确定算法

局部网格生成是无缝有限元并行计算的关键问题,然而,不恰当的局部网格生成算法会产生不一致网格.借助Delaunay三角剖分的唯一性,指出了三维凸域中Delaunay网格局部生成的不一致性,分析其产生的原因,建立避免局部网格不一致性的关键算法--搜索候选卫星点集的优化探索球算法,为网格并行生成做好数据准备.并从理论上证明了任意节点的优化探索球域包含其所有卫星点,从而,保证全域网格没有不一致性.通过算例分析了算法的正确性和局部性.     

基于切割单元直角坐标网格的DSMC方法及优化技术

对传统直角坐标网格进行改进,根据计算机图形学对与物面相交的正方形网格进行切割,使单元与物面贴体,并能根据流场的变化对网格区域进行局部加密。结合直角坐标网格与非结构网格的特点,编写了自适应直角坐标网格的生成程序,并发展了基于切割单元直角坐标网格的直接模拟蒙特卡洛(Direct simulation Monte Carlo,DSMC)方法。同时通过当地模拟分子数、动态时间步长技术优化DSMC方法,提升计算效率。数值结果表明,本方法在计算过程中使搜索效率得到大大提升,同时也保证了物面网格单元的贴体性,提高了物面附近流场的计算精度,还易于实现网格自适应。     

无网格方法关于导数计算的程序验证

为了提高无网格方法对空间导数的数值计算能力,提出并论证了新的适用于直角坐标系、柱坐标系和球坐标系的空间导数无网格算法。基于数值解与理论解的对比分析,对数值误差和收敛性进行了后验评估。评估结果表明:该算法对函数和导数的估计均为2阶精度,支撑域尺度对数值误差的大小有一定影响,但不影响数值计算的精度等级。在离散尺度h为0.01的条件下,对所选函数及其导数的数值计算相对误差不大于0.65%。      

立式捏合机桨叶结构参数对混合釜流场影响的仿真分析

为揭示立式捏合机桨叶几何结构参数与混合釜流场之间的关系,以1 L两桨立式捏合机为研究对象,采用计算流体力学ANSYS Fluent 14.5软件对立式捏合机进行了三维模型仿真研究,分析了桨叶结构参数(桨桨间隙、桨叶螺旋角)的改变对混合物料挤压应力和剪切应力的影响。结果表明,增大桨桨间隙可减小桨叶对混合物料的挤压应力和剪切应力,空心桨螺旋角分别为40°、45°、50°,桨桨间隙由1.0 mm增加到3.0 mm,挤压应力、剪切应力分别平均减小了82.4%、57.23%;增大桨叶螺旋角可增大桨叶对混合物料的挤压应力和剪切应力,空心桨间隙分别为1.0、2.0 mm,桨叶螺旋角由35°增加到55°,挤压应力、剪切应力分别平均增加了92.8%、55.7%。     

变构型航天器参数化前处理方法

针对变构型航天器结构不固定的问题,基于VC 6.0软件对Trail及Session进程文件开发的方式,研究了利用Pro/E及Patran软件对航天器进行参数化建模及网格划分的具体操作方法,通过IGES中间文件的形式,将Pro/E及Patran相连接,利用VC开发工具建立了一套变构型航天器参数化前处理平台。采用两舱构型航天器对参数化前处理方法及平台进行了验证,结果表明:该参数化建模及网格划分方法切实可行,可为变构型航天器的静力学、动力学、热、结构等分析过程的前处理方法提供参考。     

基于硬件加速的真三维显示体素化算法

提出了一种新的体素化算法,利用显示硬件支持二维光栅化算法和新的可编程控制功能,提高了体素化算法的效率。该算法通过逐层处理的方式将三维数据离散成为一系列的切片,利用硬件渲染完成切片的光栅化操作,将最终的光栅化结果组合成为体数据表达。算法利用可编程图形处理单元（GPU）的大规模纹理、动态顶点数据和灵活的纹理数据操作等特性,有效地降低了系统压力,提高了算法效能。实验证明,该算法满足真三维显示技术的要求。