二维H2/O2混合和流动的直接数值模拟

用高阶精度的紧致差分方法,对低速的随时间发展的二维可压缩H-2/O-2混合层流动编制了直接数值模拟的计算机程序,对雷诺数Re为400、来流马赫数Ma分别为0.4(上部的H₂)和1.57(下部的O₂)的H₂/O₂亚跨超声速混合层流动进行了计算,得到了相应的计算结果,并对结果进行了分析和研究.     

三维Euler方程自适应八叉树结构各向异性直角网格算法

采用各向异性直角网格对相应的三维复杂流场进行了模拟 ,提出了物面三角化方法及网格与物面的相交判断过程。用以中心差分为基础的 Jameson的有限体积法对 Euler方程进行求解。计算结果表明 ,各向异性直角网格算法对相应流场进行模拟 ,计算时间较短而效果较好     

基于CFD的静气动弹性优化设计方法

以三维Navier-Stokes方程为控制方程,数值计算弹性机翼跨音速气动力,耦合结构静平衡方程,研究弹性机翼的静气动弹性变形和真实载荷分布,并在此基础上,对机翼的型架外形进行基于静气动弹性的多学科的优化设计(MDO),以及对机翼进行基于总升力不变的飞行姿态确定,以满足弹性机翼在飞行时的品质设计要求.以某后掠机翼为例,设计结果达到预期目标.     

腔内平板对空腔自激励振荡的影响及预估振荡频率方程的改进

运用全隐式无分裂方法,紊流模型采用修正的B-L代数模型,数值模拟了超音速流通过腔内存在平板的空腔流动.研究了腔内有不同长度平板存在的空腔流动的剪切层自持性振荡,并与无平板存在时的现象进行类比.腔内平板将剪切层与空腔隔开,减弱了剪切层振荡与空腔压强脉动的耦合作用;随着腔内平板长度的增加,剪切层振荡的程度有所减小.基于特征长度可变的思想,本文改进了Heller的频率预估方程,改进后的频率方程能更好地预估空腔的振荡频率.与Heller的方程相比,修改后的方程预估所得到的频率值更接近于闭式空腔声学模型方程的预估结果;各个方程预估空腔振荡频率值的比较进一步验证了改进后的方程的正确性及合理性.     

有限元线性方程组改进的全稀疏存贮方法

有限元线性方程组的系数矩阵一般具有稀疏性和对称性的特点,全稀疏存贮方法就是利用这些特点,只存贮对称部分的非零元素,采用链表式管理,即节省存贮空间,又便于动态更改;但在实际求解计算时,查找元素花费的时间较多;这里给出了一种改进的方法,是在链表的基础上生成相应的线性表,能够方便快速查询元素,使存贮结构同时具有链表和线性表的优点.数值算例证明该存贮方案在时间和存贮上都较为占优,可靠高效,能够应用于有限元线性方程组的求解.     

Orbital corrections by a single-input impulsive control applied along the geomagnetic field

This paper examines the possibility of regular orbital corrections for a satellite controlled by an impulsive force applied along the induction vector of the local geomagnetic field. The main purpose of the correcting considered in the paper is to eliminate the drift in secular variation of the relative orbital parameters over the orbital period. The obtained results are applied to the formation flying problem. The developed method is tested using numerical simulation.     

二维欧拉方程的非结构网格自适应方法

给出提高二维欧拉定常解质量的非结构网格自适应方法。网格的生成采用Ｅｅｌａｕｎｅｙ方法,流场的算法为节点型有限体积法,能量的计算采用Ｖａｎｌｅｅｒ格式,自适应采用网格的增加和删减。对流场特性的探测采用密度物质导数的绝对值。最后计算了Ｎａｃａ００１２翼型的跨音速绕流,结果证明自适应效果很好。     

有限元法求解QPNS方程模拟高速旋成体气动力和热

从ＦＮＳ方程中略去流向粘性导数项导出熵变量形式ＱＰＮＳ方程（拟抛物化Ｎ－Ｓ方程）其离散解能自动满足熵增不等式。利用ＧＬＳ（伽辽金／最小二乘法）有限元法构造出ＱＰＮＳ方程的弱解式,所形成的整体方程组用ＧＭＲＥＳ（广义最小余数法）求解。采用空间推进求解ＱＰＮＳ方程,能明显降低计算机内存和运算时间。算例给出典型旋成体的物面压强、剪应力、热流量分布以及流场中密度、温度分布。     

Udimet 720 Li 高温变形特性的粘塑性本构建模研究

采用Chaboche统一粘塑性本构方程,对Udimet 720 Li在700℃时的单调拉伸、循环加载及蠕变特性等复杂的高温变形现象进行了建模研究.探讨了Chaboche本构理论对这些变形现象进行建模的形式,并特别针对快速各向同性软化和非对称循环加载时的平均应力松弛现象,对Chaboche本构理论进行了修正和变化.将经过修正的Chaboche本构理论,与Levenberg-Marquadt非线性优化算法相结合,编制了材料参数优化程序,得到了材料参数值.研究表明,经过修正的Chaboche本构模型可较好地建模镍基高温合金Udimet 720 Li在高温下的各种变形行为.      

航天器交会两点边界值问题

从绝对运动和相对运动两方面讨论近地空间航天器交会中的两点边界值问题。其中，绝对运动涉及多圈Lambert问题，以Lagrange时间方程为研究工具，而相对运动Lambert问题应用C—W线性解。对圆轨道之间的双冲量转移，给定转移角与转移时间，研究最小变轨速度增量所对应的转移圈数与轨道参数的求解方法，提出满足最小变轨速度增量要求的多圈转移的工程图解法，并从工程应用出发，在飞行时间约束下，按最小速度增量要求，阐述航天器交会两点之间飞行轨道（轨迹）设计方法。这种方法将飞行轨迹划分为初始漂移段、轨道转移段与终端停泊段三部分，应用两点边界值问题的解，选择两次冲量机动时刻，使速度增量之和最小。模拟算例表明，这种方法对航天器交会设计是适用的、有效的。