用三维欧拉方程计算机身与机翼绕流

本文用Jameson的三维欧拉方程有限体积法、四步Runge-Kutta时间推进格式,计算机身和大后掠细长机翼的三维可压缩绕流。对钝头机身,用C-O型网格;尖头机身与大后掠细长机翼用H-O型网格。本文介绍钝头旋成体、带座舱前机身和三角翼绕流的计算结果,显示流场等值M线分布,计算压强系数分布与实验比较,以及三角翼大迎角分离涡等。由于用了隐式残值光顺等加速收敛措施,有效地减少推进步数,节省机时。     

转动坐标系中三维跨声速欧拉流的有限体积 TVD格式

在非惯性转动坐标系中,本文采用贴体网格、有限体积法离散和修正数值通量技术,将Harten的一维TVD格式推广到三维。由于转动使方程出现源项,文中通过对源项的巧妙处理,使修改后的格式能用于非齐次双曲守恒律方程组高分辨率的数值计算。为了加速解的收敛,提高显式时间推进的CFL数,本文采用隐式残值光顺技术。三维跨声速带非齐次源项欧拉方程的典型算例表明:捕捉的激波分辨率较高;激波前后没有发现大的数值波动和伪振荡现象;所得的跨声速流场解与实验较接近。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

跨音速大迎角Euler方程数值分析

应用有限体积法离散三维Euler方程,选用三步显式格式进行时间推进求解;并通过当地时间步长、残值光顺、焓修正等方法加速收敛。对一维波动方程的von Neumann稳定性分析说明这种三步格式的最大Courant数是2。对三角翼跨音速大迎角流动的数值分析表明文中的三步格式及焓修正方法优于人们经常采用的四步、五步格式以及焓阻尼技术。     

高阶间断有限元法的并行计算研究

根据间断有限元法的数据结构特点,基于METIS网格分区技术,设计并行计算策略,在非结构网格上实现了并行高阶间断有限元法。控制方程的数值通量项使用Local Lax-Friedrichs(LLF)格式计算。设计了并行的牛顿-块高斯赛德尔法(Newton-Block GS)来加速收敛,提高迭代效率。并行性能分析表明,所设计的并行算法能够得到较好的加速比和并行效率,有效地节省计算时间,合理分配内存。这使得采用高阶间断有限元法计算更为复杂的问题成为可能。     

点隐式龙格-库塔方法的应用研究

为了提高求解器的效率,在显式龙格-库塔时间推进的欧拉方程求解器之上,发展了点隐式龙格-库塔时间推进格式。给出了其推导过程和非结构网格下中心格式和迎风格式（Roe和Van Leer格式）预处理算子的构造方法。NACA0012翼型和RAE2822翼型的跨音速无粘流动模拟表明：与显式龙格-库塔方法相比,方法能提高求解效率且内存需求相当,具有一定的工程应用价值。     

在地球极区无奇异的再入飞行器质心运动方程

传统的再入飞行器质心运动方程用经度和纬度描述地理位置，在南北极极点处奇异，在南北极区病态。本文提出利用位置矢量的三个方向余弦，称为“三余弦数”或“三元数”，代替经纬度，推导出无奇异再入质心运动方程。再入飞行器临近极点和到达极点的算例显示，使用三元数的无奇异再入方程完全避免了传统方程有关极点的奇异性。同时，无奇异方程所选状态变量与传统方程可进行很简便的转换。     

超声速平面钝体流动的快速解法

本文介绍了一种平面流动的快速欧拉方程解法。该方法将原参数非定常欧拉方程组重新组合成以广义称曼变量表示的欧拉方程组，再用二点二步迎风格式离散解。针对钝体流动，本文先建立了动网格下的方程，构造了动网格下的算法。提出了一种简单的激波处理方法。计算结果表明，该方法速度快，稳定性好，对初场不敏感。     

工程实践中病态方程组的判定及其对工程实践的作用

本文介绍了病态方程组的判定方法，结合一项具体工程指出了病态方程组对工程实践的影响，分析了出现病态方程组的原因，强调了病态方程组判定对工程实践的作用。     

超静超稳卫星碰振动力学建模

针对分离式卫星载荷模块(PM)受到扰动时可能与服务模块(SM)发生碰撞的问题，综合音圈电机反电动势(back-EMF)和柔性线缆动力学的效应，基于牛顿欧拉法建立了分离式卫星(DFP)载荷模块动力学模型。基于Hertz接触理论，推导了分离式卫星碰撞过程中连续接触力模型，并分析了碰撞过程中产生的接触力对载荷模块指向精度和指向稳定度的影响。数值仿真结果表明，碰撞使得载荷模块指向精度和指向稳定度下降5个数量级，碰撞后载荷模块可再次恢复到超静超稳工作状态，恢复时间超过1400 s。本文建立的碰撞模型对研究分离式卫星碰撞规避和碰撞控制具有重要意义。