高马赫数层流摩阻数值计算精度

飞行器在高空高速飞行时黏性效应显著，摩阻预示精度对飞行器的关键气动性能意义重大。目前，摩阻预示主要依赖数值计算，但高马赫数层流摩擦阻力计算与试验测量结果仍存在差距。以具有高速飞行器典型部件特征的球锥、三角翼为对象，结合风洞试验摩阻测量结果，使用国家数值风洞（NNW）数值计算软件和自研CFD程序研究了数值计算中影响摩阻计算精度的格式数值耗散及壁面温度边界条件等重要因素。研究结果表明：数值耗散越小，表面摩阻的计算精度越高；在速度较低的边界层近壁区内关闭熵修正，将有助于提高表面摩阻的预示精度。此外，在高马赫数流动问题的数值模拟中，壁面温度条件对表面摩阻计算同样具有重要影响。最后，基于分析结果和工程需要提出了对高精度摩阻数值预示的研究需求。     

火箭发动机超声速过膨胀射流气动噪声特性研究

采用大涡模拟LES方法计算了火箭发动机超声速过膨胀射流形态及近场声压分布,研究了入口温度与环境温度的比值(温度比)对声场的影响;将声源分解,基于Ffowcs Williams-Hawkings (FW-H)方程获取了不同位置噪声源的远场噪声,并根据声压级频谱和湍流形态分析了超声速射流噪声的产生机理。研究表明,超声速过膨胀射流气动噪声由湍流混合噪声和宽频激波噪声组成,近场噪声源以马赫波形式向大方位角辐射中高频噪声,下游大尺度湍流向低方位角范围辐射低频噪声,声压级峰值频率随观测角度增大而升高;随温度比升高,马赫波辐射角度增大,噪声指向性发生改变。该研究可为运载火箭发动机地面试车或火箭发射段声学环境设计提供参考。     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。     

液氢贮箱停放过程中的力热分析

为了寻找到运载火箭长时间停放过程中液氢贮箱的最经济液位，用于制定合理的发射流程以及紧急处置方法，采用计算流体力学(CFD)技术，对某型运载火箭停放期间液氢贮箱的力、热情况进行了仿真计算和分析。计算选用了VOF(Volume-of-fluid)两相流模型以及Lee相变模型，为了提高Lee模型在不同压力情况下对相变过程的模拟精度，采用安托因方程修正了该模型。修正后的模型首先由试验数据校验了其精确性，随后开展的液氢贮箱停放过程仿真结果表明：贮箱的竖直方向与径向均存在温度分层的现象，液相内会形成大的漩涡，该漩涡会使得冷热流体不断进行热交换，并导致贮箱内部的液氢出现气化。贮箱停放期间蒸发率最大值超过2 m 3/h，发生在停放4 h左右；而贮箱液位充填至37 m 3以上或17 m 3以下时蒸发率较低，最小值接近1 m 3/h。     

Review and propositions for the sliding/impact wear behavior in a contact interface

A thermal-solid-liquid complex operational environment induces structural interface developing a typical coupling sliding/impact wear behavior. It results in contact damage until systems fail, which may cause significant economic losses and catastrophic consequences. The key point of solving this problem is to reveal the coupling damage mechanism of the sliding/impact behavior in typical systems and life characterization under a complicate evolving environment. This has been a hot topic in the area of mechanical reliability. The main work in this paper can be concluded as follows. Firstly, the main industries in which the “sliding/impact behavior” takes place have been introduced. Then, existing studies on the wear mechanism and degree analysis are presented, which includes surface morphology analysis, wear debris analysis, and wear degree measurement. Meanwhile, existing problems in theoretical modeling and experiments in current research are summarized, so as to point out a bright direction for future research on wear prediction. They include interface contact modeling, mathematic coupling mechanism modeling, wear equation establishment, and wear life characterization, which can provide some new ideas for improving the existing studies on the sliding/impact wear behavior.     

基于成分的燃烧室出口燃气温度计算方法

针对燃烧室不同温度下燃气组成的不同,提出建立不同燃气模型以推算燃气温度的两种计算方法。基于物质守恒、化学平衡和能量守恒原理建立温度计算方程,根据能量守恒方程的不同表达形式给出平均比定压热容法和焓值法,同时给出两种求解非线性方程组的数学解法:三变量迭代法和牛顿法,并采用VB(Visual Basic)编制燃气分析测试程序实现温度算法。列举计算实例比较两种算法的区别,结果显示:随着油气比的增加,两种计算方法得到的燃气温度相差越大,考虑热离解的高温算法结果更小,原因在于燃气中离解组分的摩尔分数随着温度升高而增加,燃烧实际放热量减小,导致燃气温度减小,所以当燃气温度较高时建议采用高温算法更符合实际燃烧情况。      

基于避障伪距离的自由漂浮空间机器人规划跟踪一体化控制

针对自由漂浮空间机器人(Free Floating Space Robot,FFSR)的避障规划与控制问题，提出一种FFSR的避障规划-跟踪一体化控制方法。首先，基于障碍物伪距离技术，采用FFSR逆几何模型求解期望末端位姿下的连杆伪距离估计值，进而通过求解非线性优化问题，获得FFSR避障期望轨迹。其次，将全局轨迹规划与局部在线避障相结合，辅以离散状态黎卡提方程（DSDRE）控制方法实现FFSR的避障规划-跟踪一体化控制。最后，采用6R空间机器人模型验证了所提方法的有效性。仿真结果表明，该方法能够实现FFSR的避障控制，有效克服了传统FFSR控制中末端轨迹规划与控制相分离的问题，提高了FFSR的环境适应性。     

一种部分森林覆盖环境电波传播特性预测方法

针对部分森林覆盖环境中传播媒介的非均匀性,提出一种基于动态网格的抛物方程(PE)预测方法,该方法采用宽角PE法进行大区域电波传播特性预测。为平衡PE求解的准确性与时效性,引入动态网格技术。通过仿真算例,探讨了基于动态网格的PE法在部分森林覆盖环境电波传播特性预测中的有效性。结果表明:在计算精度相似的情况下,相对于均匀细网格,动态网格会大大减少PE所需的计算时间;而相比均匀粗网格,动态网格具有更高的计算精度。     

基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。