求解代数Riccati方程的一种新算法

本文讨论了Hamiltonian矩阵在辛相似变换下的标准形,由此给出代数Riccati方程存在非负定解的一个充分条件,提出了求解代数Riccati方程的一种新的算法。该方法节省运算量,且精度较高,尤其适合于解阶数不太高、系数矩阵为满阵的Riccati方程。最后给出了一个数值例子,并将该方法与其他方法作了比较.     

李雅普诺夫方程AX-XB=C的误差界

本文对李雅普诺夫方程AX-XB=C及(A+δA)(X+δX)-(X+δX)(B+δB)=C+δC的解X及X+δX给出了相对误差‖δx‖/‖X‖的一个上界。     

双曲型守恒律方程的一个大时间步长二阶TVD差分格式

本文根据Harten,A.的大时间步长差分分格式构造思想,为双曲型守恒律方程弱解计算构造了一个2K 3点大时间步长二阶显式差分格式——LTS-LF格式,得到了其在CFL限制K下为总变差不增差分格式(TVD格式)。文章按照Roe的方法推广格式到方程组情形,并就Burger’s方程和Euler方程组黎曼问题进行数值试验,结果令人满意。     

串列叶栅粘性紊流流场有限分析数值研究

本文将近年发展起来的有限分析方法(Finite Analytic Method)用于曲线坐标系上紊流N-S方程的数值计算。计算了来流雷诺数为2.0×10~5,气流攻角分别为0°,10°,-10°三种情况的二维粘性紊流串列叶栅流场。文中用k—ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。数值计算结果与实验结果相比较,吻合程度令人满意。     

共轴刚性旋翼直升机着舰飞行特性

建立一套能够耦合非定常舰船流场的共轴刚性旋翼直升机飞行动力学模型,以研究该旋翼飞行器在着舰过程中的飞行特性。在舰艉流场模拟方面,采用分离涡方法以获得高精度非定常流场数据;在飞行力学建模方面,引入上下旋翼干扰因子建立共轴刚性旋翼诱导速度模型,并采用等效挥舞运动概念建立其挥舞运动方程;基于“单向耦合”思想构建了计算流体动力学（CFD）向飞行力学模型的数据传递策略,并分别以XH-59A直升机和UH-60A/SFS2组合为算例验证了飞行动力学模型和数据传递策略的正确性。以SFS2舰船模型和XH-59A共轴刚性旋翼直升机组合为研究对象,从直升机操纵余量和非定常载荷水平两方面分析了着舰过程中舰艉流场对共轴刚性旋翼直升机的扰动特征。时均研究结果显示：由于共轴上下旋翼受时均流场的扰动存在差异,飞行员在增大总距以维持高度稳定的同时,还需要减小差动总距以保持机头朝向的稳定。非定常水平研究结果表明：对于共轴构型直升机,舰船流场对拉力及俯仰力矩的非定常扰动,是引起飞行员工作载荷增大的主要因素。     

关节铰柔性多体系统动力学研究

针对关节铰柔性多体系统建立了Ｒ／Ｗ体系的动力学方程，并进一步推广到非树及约束型多体系统，本文采用适合描述弹性变形和进行弹性体离散的刚性连体浮动坐标系，用相对描述方法来分解部件的运动，用部件模态方法对弹性体进行离散，从达朗伯原理出发建立了Ｒ／Ｗ体系的系统动力学方程。方程对于运动的分解，部件的弹性离散以及对系统的结构描述都比较直观和简单，系统的广义坐标数目也小。方程推广到非树系统和约束系统后，采用奇异     

任意回转面叶栅C网格的数值生成

首先尝试采用任意曲面椭圆方程，数值生成任意回转面上叶栅Ｃ网格，经过对边界条件的仔细处理，生成的网格可保证在周期外边界和重合内边界上网格线导数连续，在叶片表面网格线相互垂直，提出了结合初始网格线分布，控制环绕叶片网格线分布的方法。     

二维突扩管流的有限分析解

运用有限分析法求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，获得了二维突扩管流的数值解。对于雷诺数为２５，５０，１００，１５０，２００的来流，预报的流动图谱具有相似性，回流区域及再附距离的计算也获得了比较满意的结果     

翼面吹气对过失速非定常翼面涡的影响

在１ｍ低速风洞中，研究了边条翼在近人速非定常流动控制后的空间流态。实验采用翼面吹气控制翼面非定常流动，通过烟流显示涡轨迹和涡破碎位置，用相锁照相技术记录空间流态，结论表明，在机翼上仰过程中，翼面吹气能延迟前缘涡的破碎，在机翼下俯过程中，吹气有利于前缘涡的生成和发展。     

用 Euler 方程计算某些特殊的二维分离流动

采用非定常Ｅｕｌｅｒ方程、总体时间步长、隐式时间推进，计算了几种特殊类型钝体的二维有分离绕流。一类是光滑的圆柱体，在跨音速时计算出因激波诱导分离而周期性地脱落在尾迹上的非对称涡排。另一类是带尖角的物体，如三角形柱体，计算出从尖角处脱落在尾迹上的分离旋涡，形成卡门涡街。对这两类物体，计算的平均阻力系数，分离涡脱落的Ｓｔｒｏｕｈａｌ数以及流态都与实验接近。结果表明，对于这类绕流问题，用非定常Ｅｕｌｅｒ方程计算，具有一定工程实用意义。