一种新的用于MHD方程逆风格式的Jacobian矩阵分裂方法

针对理想MHD方程,提出了一种新的基于MacCormack算法的雅可比矩阵分裂方法,克服了原有方法稳定性差的问题,并成功地应用于理想MHD方程的求解.控制方程在非结构混合网格上进行空间离散,其中对流项采用本文发展的逆风向量分裂格式,并引入了双曲型磁场散度清除技术,时间推进为显式5步龙格-库塔方法.对MHD激波管流动和带均匀磁场干扰的二维高超声速钝头体绕流流场进行了数值模拟,得到了与参考文献相吻合的数值结果,表明本文发展的数值分裂方法可以有效地捕捉MHD流场的流动特征,并且具有比MacCormack方法更高的稳定性和计算精度.     

二维爆震波数值模拟

针对爆震发动机的多维效应,利用CFD计算软件对单管脉冲爆震发动机进行了二维的数值模拟。以氢气和氧气混合物为例,仿真了爆震波在管内的传播过程以及传出管外后的流场情况。模拟所得的爆震波参数与文献中的计算结果和实验值基本一致,验证了模拟的正确性。爆震波传出爆震管后退化成一道激波,对该激波的流动情况进行了分析。     

二元高超声速进气道内压缩通道/隔离段曲面构型

对于二元高超声速进气道内压缩通道及隔离段设计,提出了进气道下壁弧形曲面构型方案。在一系列不同收缩比、不同波系配置的平面构型进气道基础上,通过基于N-S方程的数值模拟研究了不同半径的弧形过渡曲面对进气道性能的影响。发现采用弧形曲面过渡可以削弱平面构型方案对气流不必要的膨胀,减小隔离段进口处上侧壁面高压,改善隔离段进口气流均匀性。新构型有助于降低起动马赫数,且弧形过渡半径越大,收缩比越大,降低的程度越明显;还可以大大提高进气道的总压恢复,无须最后一道内压激波打在下壁面肩点上即可获得较高的性能。     

隔离段反压对激波串起始位置的影响

使用商业软件FLUENT对三维矩形截面隔离段流场进行了数值模拟。在来流马赫数2.3条件下,通过改变不同的反压来研究反压对不同长高比等截面隔离段和带扩张角隔离段激波串起始位置的影响。结果表明:隔离段长高比对反压与激波串起始位置之间的关系有非常大的影响,当长高比不是非常大时,激波串起始位置与反压成线性关系,当长高比超过此范围后,线性关系不再成立。     

轴对称喷管内外流场与结构温度场耦合计算

利用有限体积的对流项二次迎风插值格式和重整化群(RNG)k-ε湍流模型,二层增强型壁面函数,同时利用球形谐波法考虑热辐射的影响,以灰气体加权模型(WSGGM)确定气体介质的辐射性质,求解N-S方程、热传导方程、考虑吸收-发射性气体介质的辐射传输方程。采用流固耦合的流动与换热模型,流场与结构温度场互为边界条件交换数据,实现了流场解算与温度场解算的耦合数值分析。采用此计算模型对静止的轴对称收-扩喷管进行了数值模拟,计算得到的喷管壁温与试验数据吻合良好。在此基础上,进一步开展了巡航状态下轴对称收-扩喷管内外流场与结构温度场的耦合数值研究。     

The Expansion of Coronal Plumes in the Fast Solar Wind

Coronal plumes are believed to be essentially magnetic features: they are rooted in magnetic flux concentrations at the photosphere and are observed to extend nearly radially above coronal holes out to at least 15 solar radii, probably tracing the open field lines. The formation of plumes itself seems to be due to the presence of reconnecting magnetic field lines and this is probably the cause of the observed extremely low values of the Ne/Mg abundance ratio. In the inner corona, where the magnetic force is dominant, steady MHD models of coronal plumes deal essentially with quasi-potential magnetic fields but further out, where the gas pressure starts to be important, total pressure balance across the boundary of these dense structures must be considered. In this paper, the expansion of plumes into the fast polar wind is studied by using a thin flux tube model with two interacting components, plume and interplume. Preliminary results are compared with both remote sensing and solar wind in situ observations and the possible connection between coronal plumes with pressure-balance structures (PBS) and microstreams is discussed. This revised version was published online in June 2006 with corrections to the Cover Date.     

钛合金筋板类构件局部加载成形有限元仿真分析中的摩擦及其影响

局部加载等温成形为大型钛合金筋板类构件省力、高性能成形制造提供了有效的途径。然而其成形过程复杂,材料参数、几何参数、工艺参数高度耦合作用,精确预测困难。有限元仿真为此类复杂成形问题的分析和优化提供了有力的工具。而摩擦模型和摩擦条件是数值模拟的重要边界条件之一,影响到预测结果的精度及可靠性。本文着重分析了库伦摩擦模型、剪切摩擦模型、库伦-剪切摩擦模型、库伦-粘着摩擦模型以及摩擦大小对筋板类构件局部加载成形仿真分析结果的影响。研究结果表明:剪切摩擦模型适用于局部加载成形过程的建模仿真;确定精确的摩擦条件可降低载荷预测误差15%以上;主动调控局部区域摩擦条件可改善局部加载成形过程的充填、减少缺陷。     

金属激光增材制造过程数值模拟

金属激光增材制造过程中,热应力导致零件发生形变;气孔与熔合不良等缺陷降低零件的拉伸以及疲劳性能;熔池内的凝固微观组织,尤其是柱状晶等轴晶转变(Columnar to Equiaxed Transition,CET)是影响零件性能的重要因素。针对上述3个方面,回顾了金属激光增材制造数值模拟的发展历史,对其研究现状和存在问题进行了评述,阐述了金属激光增材制造过程中所采用的数值模型和数值方法,包括热应力分析的有限元(Finite Element Method,FEM)方法、模拟熔池金属液流动的计算流体力学方法(Computational Fluid Dynamics,CFD),以及凝固微观组织模拟的相场法(Phase Field,PF)和元胞自动机方法(Cellular Automaton,CA)。在此基础上对金属激光增材制造过程数值模拟的前景及趋势进行了展望。     

齿轮滚轧成形技术及其研究进展

齿轮滚轧成形技术是一种结合齿轮锻造成形优势和轧制技术特点的齿轮塑性成形新工艺.对齿轮滚轧成形技术的基本原理和特点进行了简要介绍,对国内外有关齿轮滚轧成形技术的研究现状进行了综述,重点从工艺基础理论、数值模拟分析以及试验研究等方面给出了当前齿轮滚轧成形技术面临的主要问题,并据此提出了未来齿轮滚轧成形技术研究发展的若干建议.     

大子午扩张涡轮端区的流动传热及端区正弯效果的数值研究

为了研究大子午扩张涡轮端区流动和传热特性,并研究叶片端区正弯技术在大子午扩张涡轮中的气动和传热效果,对某大子午扩张涡轮静叶进行数值模拟。运用SST湍流模型精确捕捉流动结构,并进行了气动和传热预测的有效性实验验证。通过分析结果,对大子午扩张涡轮端区流动和传热特性以及两者相互影响关系进行了深入研究,分析了端区正弯技术在重组大子午扩张涡轮端区流动以及合理分布热负荷的应用效果。结果表明：大子午扩张端壁导致涡轮端壁附面层的强烈分离,通道涡分离点提前约15%,高传热区受马蹄涡和通道涡的强烈影响;端区正弯有效地改善了大子午扩张静叶端壁的附面层分离,减小前缘的热负荷25%,提高涡轮的气热性能。