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

运用全隐式无分裂方法,紊流模型采用修正的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线性解。对圆轨道之间的双冲量转移，给定转移角与转移时间，研究最小变轨速度增量所对应的转移圈数与轨道参数的求解方法，提出满足最小变轨速度增量要求的多圈转移的工程图解法，并从工程应用出发，在飞行时间约束下，按最小速度增量要求，阐述航天器交会两点之间飞行轨道（轨迹）设计方法。这种方法将飞行轨迹划分为初始漂移段、轨道转移段与终端停泊段三部分，应用两点边界值问题的解，选择两次冲量机动时刻，使速度增量之和最小。模拟算例表明，这种方法对航天器交会设计是适用的、有效的。     

二维非定常欧拉方程的摄动解在超声速矩形翼上的应用

本文利用超声速线化理论给出的矩形翼载荷系数解析解作为权函数,对沉浮和俯仰谐振的二维平板翼压强系数进行加权修正,获得矩形翼谐振气动力的近似解,计算结果与线化理论数值结果符合得很好,本方法是超声速和高超声速矩形翼非定常气动力的一种统一和快速的计算方法。     

用椭圆型、抛物型两种方法计算平面叶栅湍流尾迹

本文用连续方程、稳定不可压缩椭圆型和抛物型时均N-S方程及K-ε双方程湍流模型计算了平面叶栅湍流尾迹。计算表明,上述两种方法都能得到较好的数值结果,但抛物方法可节省大量计算时间和内存空间。     

非定常不可压粘性/无粘性耦合方程的一种分步分解方法

给出了数值求解初始变量不可压Navier-Stokes／Euler耦合方程的一种分步块LU分解方法。与传统的时间分裂法不同,该法无需压力中介边条件,从而避免了传统时间分裂法要求的复杂的压力中介边条件逼近。分步块LU分解方法可看做经典的Uzawa算法的改进,后者曾被成功应用于不可压Navier-Stokes／Euler耦合方程的求解。但本文显示分步块LU分解法比经典的Uzawa方法更经济。分析显示该法具有良好的稳定性和高精度,数值结果支持这一理论分析。