喷管型面对固体火箭发动机性能的影响

用时间相关法数值计算跨音速喷管流动时,其控制方程分别采用守恒型和非守恒型方程.通过对JPL喷管进行计算发现:在相同的精度下,非守恒型方程的收敛速度比守恒型方程的收敛速度快一倍多.采用修正阻尼法(DASHCO法)进行数值计算,其收敛速度均提高六倍以上.最后,采用直接优化法研究了喷管收敛段型面及喉部型面对发动机性能的影响.     

跨声速颤振计算方法研究

采用改进的弹簧近似方法,以ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程为控制方程,耦合结构运动方程,计算了Isogai机翼模型的颤振边界。气动力求解采用双时间推进,外时间为物理时间,与结构运动方程同步;结构运动方程求解采用Rayleigh-Ritz方法,时间推进也采用双时间法,节省了大量计算时间。研究表明,计算结果与参考文献提供的结果吻合良好,证明了所提出的跨声速颤振计算方法的准确性。     

基于LS-DYNA滚动轴承使役可靠度计算方法

利用软件LS-DYNA对滚动轴承使役过程进行了动态有限元计算,基于响应面法（RSM）建立6408滚动轴承输入参数与输出响应的函数关系方程,设定各类随机误差分布,利用矩阵试验法确定抽样点在输入变量抽样空间的位置,进行一系列确定性拟合试验.每次试验借助LS-DYNA进行精确求解,运用求解数据对方程进行最小二乘法回归分析确定方程的组成项及系数.最后利用函数关系方程计算轴承任意时刻的使役可靠度.结果表明:计算得出该轴承在0.020008s时的使役可靠度为97.062%.      

先进飞行器气动隐身一体化设计方法研究现状与趋势

介绍了先进飞行器气动隐身一体化设计的研究现状与发展趋势。首先介绍了设计中的常用的数值模拟技术。针对气动力计算,重点分析对比了位势方程、Euler方程和N—S方程等气动力计算方法;针对隐身性能计算,重点分析对比了高频近似方法、时域法、频域法等雷达散射截面计算方法。在此基础上。讨论了先进飞行器气动隐身一体化设计的策略和步骤,最后对该设计方法的发展趋势进行了分析。     

双欧法与四元数法的应用比较

对解决欧拉方程奇异性的双欧拉和四元数法进行了对比研究，四元数法从理论上讲比较完美，但实际应用中存在较大的累积计算误差，从而影响计算精度；双欧法利用正、反欧拉方程间精华区倒挂关系进行分区交替运算，把精华区扩展到全域，不仅根除了奇异性，而且计算误差小。因此，对于解决欧拉方程奇异性来讲，双欧法要优于四元数法。     

过冷水滴撞击结冰表面的数值模拟

为了计算过冷水滴撞击结冰表面,本文应用雷诺平均Navier-Stokes方程和k-ε两方程紊流模型计算结冰表面外的空气流场。用欧拉法建立了空气中过冷水滴的运动方程,并采用有限控制容积法对方程进行了数值计算。为了提高计算精度,方程中的对流—扩散项采用MUSCL格式进行了离散。计算结果表明过冷水滴的大小和机翼攻角的变化对结冰区的大小和水收集系数的影响很大。      

时间相关法计算喷管跨音速流动

本文用时间相关法完成了定常、无粘跨音速喷管流场的计算,内部点参数用Mac Cormack差分格式计算、壁边界点参数用简化了的特征边界条件计算,轴上点参数用反射原理计算。计算表明:计算是收敛的,计算结果与实验数据符合良好。 定常、无粘、跨音速喷管流动的控制方程是混合型的(即在亚音速区域是椭圆型,而在超音速区是双曲型的),给数值求解带来很大的困难。为了克服这一困难,广泛采用时间相关法,即认为定常流动方程的解是相应的非定常流方程的解在时间趋于无穷时的渐近解。因为非定常流的控制方程不论在亚音速区,还是在超音速区域,都是双曲型的,可以用统一的方法来求解,而且易于求解;另一方面,它可适用于形状比较复杂的喷管。 本文首先阐述了时间相关法的计算方法,然后列出算例喷管的计算结果。通过与实验数据的比较,证明计算是符合实际的。     

基于离散伴随方程的三维雷达散射截面几何敏感度计算

针对有限差分法计算雷达散射截面（RCS）梯度效率低，采用高精度雷达散射截面评估时计算代价高的问题，提出了一种基于麦克斯韦积分方程离散伴随方程的RCS梯度高效计算方法。基于伴随方程的梯度计算可以通过一次雷达散射截面求解、一次伴随方程求解获得RCS关于所有设计变量的梯度。其中麦克斯韦积分方程离散伴随方程的形式与原方程基本一致，可以采用矩量法（MOM）及多层快速多极子算法（MLFMA）求解。伴随方程求解计算量与直接雷达散射截面评估基本一致，存储量在直接雷达散射截面评估的基础上增加不明显。通过双椎体模型、导弹模型对基于矩量法、多层快速多极子算法的伴随梯度进行验证，证明了基于伴随方法的RCS梯度计算可以实现复杂外形中RCS梯度的高效、高精度求解，为基于梯度的高精度气动/隐身一体化优化提供了基础。     

紧致格式数值模拟超音速粘性绕流问题

通过数值求解可压Navier-Stokes方程的方法计算了球锥粘性绕流问题。方程中粘性项按通常办法逼近,为改善收敛速度利用了文[1,2]中之算子附加修正方法。方程中之流向导数利用了四阶精度的紧致差分。这是由于该方向上流动参量光滑和边界条件易于处理的缘故。在法向方向上利用了中心差分。对差分方程之隐式部分利用了近似因式分解法。由于选用了文[3]中之特殊Jacobiall系数矩阵分裂法而使计算工作量大为减少。用这种方法计算了超音速粘性绕流问题。计算结果与实验结果相符很好。     

改进节点积分的无单元Galerkin法及其在流动问题中的应用

无单元Galerkin法需要在背景网格上积分,计算量大,且在求解对流占优问题时会出现非物理的数值伪振荡现象。为此,基于局部Taylor展开思想,采用节点处的局部Taylor展开计算积分,建立了局部Taylor展开积分无单元Galerkin法。该方法同时解决了标准的无单元Galerkin法计算量大和对流占优时会出现数值伪振荡的问题。一维定常对流扩散方程和二维Burgers方程的求解说明了该方法的有效性,且计算效率远高于无单元Galerkin法。     

Finite element analysis of dynamic stability of bearingless rotor

Dynamic stability equations of bearingless rotor blades were investigated using a simplified model.The aerodynamic loads of blades were evaluated using two-dimensional airfoil theory.Perturbation equations were obtained by linearization of the perturbation.A normal-mode approach was used to transform the equations expressed by nodal degrees of freedom into equations expressed by modal degrees of freedom,which can reduce the dimension of the equations.The stability results of rotor blades were presented using eigenvalue analysis.The shape function matrix was obtained using spline interpolation,which simplified the analysis and made assembly of the inertial matrix,damping matrix,and stiffness matrix a simple mathematical summation.The results indicate that the method is efficient and greatly simplifies the analysis.     

低马赫数流动中的预处理Euler方程的收敛特性

对预处理后Euler方程在低马赫数流动数值模拟中的收敛特性进行了细致研究。首先采用量级分析法研究了Weiss-Smith预处理方法对低速情形下Euler方程收敛性的影响;然后针对低速情形下预处理后的Euler方程组中各方程的收敛性,提出了一个改进的预处理矩阵。数值模拟结果表明：Weiss-Smith预处理方法可以很好地改善Euler方程在低速流动时的收敛特性,并且改进后的预处理矩阵在连续方程的收敛特性方面效果明显。     

统一粘塑性本构方程一致切线刚度矩阵的研究

统一粘塑性本构方程是微分或增量形式的方程组,内变量是隐式的,一致切线刚度矩阵很难给出明确的表达式.针对包含热恢复项的统一的粘塑性本构方程,在向后欧拉积分的基础上,为本构模型推导和确定了一个新的考虑率相关塑性的一致切线刚度矩阵(Consistent tangent modulus)的表达式,并在通用的有限单元分析程序MARC中实现.通过一些算例的有限元分析,检验了一致切线刚度矩阵的可行性.      

火箭发动机两相喷雾燃烧的并行虚拟机仿真

采用计算流体动力学方法对液体火箭发动机内部多维喷雾两相燃烧过程进行了数值模拟。气相控制方程组用欧拉坐标系下的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组描述,液相控制方程组在Ｌａｇｒａｎｇｉａｎ坐标系下进行描述。气、液两相作用通过方程组的源项互相耦合,编制了串行和并行程序,并在并行虚拟机环境下进行了测试。计算结果显示并行计算的效率较高。     

改善航空发动机特性计算收敛性的方法

详细分析了航空涡轮发动机部件共同工作的非线性方程组的求解收敛性问题,统计出不收敛的类型并分析了不收敛的机理.针对常见的不收敛问题,以牛顿迭代法为基础,分别提出了以发动机部件特性图为基础的独立变量值限制法、独立变量初值拟合法、变步长牛顿法以及部件特性扩展法等解决措施.特性计算程序中采取这些改进的方法后,对某型定几何单轴涡喷发动机和某型尾喷管喉部面积可调的变几何双轴涡扇发动机特性计算结果表明,发动机特性计算的收敛性和收敛速度得到大幅度提高.      

用广义乘子法求解航天器最优平面再入轨迹

从纵向平面无量纲航天器质点运动方程出发，引入了能量参数对运动方程进行推导，使得运动方程和优化问题易于处理。通过将优化变量转化为分段连续的线性函数和Runge—Kutta数值计算方法，将轨迹优化问题转化为非线性规划问题，应用广义乘子法对其进行数值分析，给出了不同航程要求下．总加热量最小，满足终端要求以及过载、动压、热流约束的航天器最优轨迹，并分析了其特点。     

多转子系统振动分析的子结构耦合矩阵法

提出了一种用于多转子耦合系统振动分析的建模方法.将耦合的多转子系统在耦合处分开,把耦合力作为单转子系统上的外力,利用集总参数法和节点的自由振动微分方程通式建立单转子系统的振动方程.通过分析耦合部件作用在耦合节点上的力,得到了多转子系统耦合部件的耦合矩阵;利用耦合矩阵将单转子系统振动微分方程耦合在一起,建立了整个多转子系统的振动方程.通过分析多转子系统的振动方程,解出了多转子系统的固有频率、临界转速、幅频响应曲线、涡动轨迹等.通过两个算例,证实了多转子系统振动分析的子结构耦合矩阵法是可行的,计算结果是可靠的.      

抛物化Navier-Stokes方程在求解高超声速流动稳定性分析 的基本流中的应用

通过在边界层内部使用边界层上缘的流向压力梯度代替该处的流向压力梯度,对传统抛物化Navier-Stokes(PNS)方程求解中的流向压力梯度的处理方法进行改进,从而使它可以求解用于稳定性分析的高超声速流动的基本流,并分别计算了平板、零攻角钝锥和小攻角钝锥三种典型算例.对于流向没有分离,且横向流动不强的流动,使用PNS方程计算高超声速流动稳定性分析用的基本流是可靠的.特别是改进后的方法求得的基本流的稳定性分析结果与直接求解Navier-Stokes(N-S)方程的结果非常吻合,但计算的时间消耗和空间占用都减少了一个数量级.     

气体涡轮燃烧室流动与燃烧三维数值模拟

用三维湍流N -S方程和单步快速不可逆化学反应描述涡轮燃烧室内的流动与燃烧过程。采用非交错网格中的SIMPLE算法求解控制方程并对涡轮燃烧室内的流动与燃烧过程进行三维数值模拟 ,得到了有代表性的圆筒形燃烧室内燃气参数的详细分布情况。模拟结果表明 :对圆筒形涡轮燃烧室内的流动和燃烧而言 ,预计结果与测量结果趋势相同。本文是研究涡轮燃烧室的初步尝试     

大涡模拟超声速两相流混合层的非线性分析

为研究超声速两相流混合机理,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构.用二维盒式滤波法处理可压缩气相N-S方程,得到亚格子形式的LES控制方程;用轨道模型仿真液滴运动,并用单液滴蒸发模型模拟液滴蒸发过程.利用伪邻近法和互信息量法确定了相空间重构的嵌入维数和时间延迟,对混合层压力时间序列进行了相空间重构.分析了仿真结果中混合层的涡结构特性和液滴对流场的反作用效果,并讨论了混合层失稳过程的非线性特征.发现混合层发展的不同阶段与其系统的非线性特征值密切相关.