共轴对转螺旋桨的非定常气动干扰

为优化共轴对转螺旋桨的气动性能,研究前后桨之间复杂的非定常气动干扰,使用非定常雷诺平均Navier-Stokes （URANS）方程的数值模拟方法,以及滑移网格技术,对6×6构型的共轴对转螺旋桨进行了三维非定常气动干扰计算,获得了共轴对转螺旋桨的气动特性变化规律,并给出了优化气动外形及运行控制的建议。计算结果表明:来流马赫数为0.38的情况下,由于非定常气动干扰导致的整个共轴对转螺旋桨的拉力系数和功率系数的脉动幅值为1.5%左右,推进效率的脉动幅值在0.1%左右。一个旋转周期内,前后桨均出现了12次周期波动,后桨由于能够吸收一部分前桨产生的切向滑流能量,推进效率比前桨更高一些,而且能够提高整个桨的效率。另外通过压力系数以及轴向速度的径向分布情况看,70%左右的桨半径位置上非定常气动干扰强度相对较大。      

Effects of a hot jet on detonation initiation and propagation in supersonic combustible mixtures

The effects of a hot jet on detonation initiation and propagation in supersonic combustible mixtures has been studied with two-dimensional numerical simulations with the open-source program AMROC that uses a block-structured adaptive mesh refinement method. Results indicate that the hot jet could ignite the detonation effectively in supersonic combustible mixtures like a pneumatic ramp. After the realization of the detonation initiation, the hot jet can still play an important role on the detonation propagation during its continuous ejection. For a hot jet with certain diameter, it can result in an overdriven detonation with almost constant overdrive degree. After the shutdown of the hot jet, the stable CJ detonation combustion was realized finally in the supersonic combustible mixtures. With the re-ejection of the hot jet, the failed detonation could be reinitiated quickly. Through the control of the re-ejection of the hot jet, it plays a key role not only in the initiation process, but also in the subsequent continuous detonation combustion period.     

气体动理论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反馈+陷波滤波器的控制算法;进一步分析影响指向精度的因素,设计带有调节滤波器的改进控制器以减小天线视轴指向的周期性误差,提高指向精度与稳定度。最后,给出数值仿真校验的结果并与现有的方法进行比较。结果表明,所提出的控制算法在保证天线指向稳定的同时,可以有效减小天线视轴的指向偏差,并对干扰与振动具有良好的鲁棒性。     

基于切割单元直角坐标网格的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）的大规模纹理、动态顶点数据和灵活的纹理数据操作等特性,有效地降低了系统压力,提高了算法效能。实验证明,该算法满足真三维显示技术的要求。