火箭垂直发射燃气流场的数值模拟

数值模拟了火箭垂直发射燃气流场的非定常过程,在计算中采用了二阶精度的ＭＵＳＣＬ格式,在引入单调性限制条件和一定的耗散机制后,该方法适应于Ｅｕｌｅｒ坐标系下的计算。计算结果表明,该方法得到的流场激波结构十分清晰,计算结果和实验规律相一致。     

Gamma过程退化失效可靠性分析

退化型失效产品,可通过退化失效分析对其进行可靠性统计推断。性能退化一般是一个随机过程,并且是单调变化的,为描述产品的退化,给出了一个随机Gamma过程模型,通过求解该过程的首达时分布即可获得产品的失效分布。由于解析方法求解存在很大的难度,给出了一种基于仿真的求解方法。利用所得模型对强激光装置所用的某型金属化膜脉冲电容器进行了可靠性分析,说明了模型的适用性。     

高阶满阵线性代数方程组的快速解法及其应用研究

在非对称、非正定满阵线性方程组的求解过程中,本文以高斯迭代为基础,提出一种残差优化的迭代法.该法在迭代收敛过程中,残差单调下降,大大缩短了线性方程组的求解时间.用该法替代直接解法,增大了求解方程组的阶数,缩短了计算周期,对复杂飞行器气动、隐身一体化优化设计具有重要的工程价值.该方法还可以广泛地推广到其它应用领域.     

固体推进剂动态储能模量主曲线计算应力松弛模量

基于一维线性粘弹性理论中动态储能模量与应力松弛模量的关系,导出了一种利用动态储能模量主曲线计算固体推进剂修正应力松弛模量的方法。通过两类固体推进剂修正应力松弛模量与实测应力松弛模量主曲线对比发现:虽然计算值大大高于实测值,但二者的比例系数随时间单调递减。根据这种关系,利用固体推进剂动态储能模量主曲线可以计算一定时间内的应力松弛模量。最后,对该方法在发动机药柱应力松弛模量监测中的应用进行了初步探讨。     

某型大涵道比涡扇发动机飞行推力确定方法研究

以某型大涵道比涡扇发动机为研究对象,结合动力装置适航取证＂飞行推力确定＂试验科目,对该型发动机飞行推力确定方法进行了分析,并开发了该型发动机飞行推力确定计算分析软件。结合发动机设计单位给出的设计值,对该计算方法进行了验证。结果表明,采用飞行推力确定计算方法得到的计算值与设计值间吻合程度较高。     

奔月转移轨道的快速设计方法研究

 结合重叠圆锥曲线方法,针对着陆型和月球卫星型两种类型奔月转移轨道,给出了一种奔月转移轨道的快速设计方法。该方法是一种无需轨道数值积分的纯代数计算方法,直接从目标轨道参数出发,具有计算速度快、精度高的特点,避免了圆锥曲线拼接法精度差的缺点,可用于这两类奔月转移轨道的初始设计。该方法可为精确轨道计算提供较好的初值,大大缩短了轨道精确计算所需的时间。同时还针对地球扁率( J2 ) 的影响,给出了一种修正方法,提高了本方法的计算精度。     

某型螺旋桨拉力确定试飞参数不确定度分配及验证

对依据某型螺旋桨风洞试验数据所建立的螺旋桨拉力计算模型的准确性进行了验证,利用模型开展了参数的敏感性分析工作。针对以往试飞中存在的问题,提出了一种螺旋桨拉力确定试飞参数不确定度分配方法,该方法能够反映不同参数对计算结果的影响程度,并利用蒙特卡洛法对分配方法进行了验证。结果表明,由该分配方法计算所得的拉力计算的不确定度可满足预定要求,分配方法是可靠的。     

LTRAN2的一种改进方案

本文提出了LTRAN2的一种改进方案BTRAN2用来分析翼型的各种频率的非定常运动。用E-O调转换格式和ADI方法求解了完全的二维非定常跨音速小扰动位流方程,并用单调转换的AF2格式计算了定常跨音速小扰动方程,以此作为非定常计算的初场。本文给出了带后缘正弦振荡襟翼的NACA 64 A006翼型的绕流和做正弦俯仰振荡的NACA 64 A010翼型的绕流计算结果,它们与Euler方程解或实验数据很吻合。     

用柱面坐标法进行三维装药的设计计算

本文介绍了固体火箭发动机翼锥柱药型设计计算的一种方法——柱面坐标法.并针对某发动机给出了该方法的原理、步骤及计算结果与试车结果的比较.该方法以数值解与解析解相结合,计算公式简单明了,处理边界问题方便,计算速度快,提供结果多,是一种比较实用的方法.     

径向基函数在三维Euler方程数值计算中的应用

基于三维Euler方程,对非结构四面体网格给出了一种基于紧支径向基函数重构的ENO型有限体积格式,方法的主要思想是先对每一个四面体单元构造插值径向基函数,而在计算交界面的流通量采用高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法。最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力。     

非结构网格上新型的NND有限元格式

张涵信发展的ＮＮＤ差分格式是由中心差分格式,二阶迎风格式和一阶迎风格式混合组成的杂交型格式。众所周知,和中心差分格式相对应的是Ｇａｌｅｒｋｉｎ有限元格式。通过对中心型有限元格式加修正项的方法本文成功地构造出二阶迎风型有限元格式和一阶迎风型有限元格式。     

高精度差分求解气动方程的几个问题

本文探讨了发展高精度格式的必要性，研究了高精度格式与网格以及熵增条件的关系，并发展一种半离散化的空间为三阶精度的格式。从此半离散化的格式出发，可建立多步显式格式及隐式格式。模型问题的计算表明，该三阶精度的格式具有好的精度，且激波附近基本上没有波动。     

求解双曲型守恒律的高精度,无波动样条逼近有限体积方法

本文从三次样条逼近出发，提出了求解流体力学双曲型守恒律的一种高精度、无波动的数值方法。针对样条函数的特性，本文提出了一种新的通量限制技术，使该方法在光滑区可以达到四阶精度，在流动参数的空间分布出现拐点或极值点时分别退化为二阶或一阶精度的格式。数值实验表明，该方法对流场中的激波和接触间断有很高的分辨率，优于二阶精度的ＴＶＤ格式。     

变结构系统的递阶切换模式

 本文研究关于线性多变量变结构系统的递阶切换模式及其切换函数c与变结构控制u（x,s）的设计问题。文中提出了两种递阶一切换模式，并按模式的到达条件及对滑动模态的极点配置的原则给出了各种情况下切换函数及变结构控制的设计方法。此种模式的最大优点是系统的控制容许失误。     

一种混合型四阶格式、基于特征的边界条件及其应用

本文描述并应用差分格式的构造原则,采用了时空转换的格式构造方法构造了一个四阶精度的差分格式,文中还进一步提出将四阶格式和二阶ＮＮＤ格式混合的方法,即在光滑区格式保持四阶精度、在激波附近降为二阶ＮＮＤ格式,来增强格式捕捉激波的能力。同时,本文基于特征传播的思想,建立了一种开边界的数值计算方法,并采用将边界条件嵌入控制方程的方法以及新的简化形式的边界处理方法,完整给出了无粘、粘性物面和开边界的的边界计     

高精度两步TVD格式的构造及数值检验

基于通量分裂的思想,利用Taylor级数理论和第二步限制器函数的控制,构造了在光滑区域空间为三阶、时间为一阶精度的高分辨率全变差递减(TVD)格式.该格式在激波过渡点降为一阶迎风格式.并从单个线形方程推广到非线性方程及方程组情形.通过几个典型算例的计算,并与二阶TVD格式作了比较,表明该方法对流场中的激波有较高的分辨率,且是无波动的.     

新型高低速流统一算法及其叶轮机应用

 对于可压缩/不可压缩混合流动或高低速混合流动问题，目前没有较好的解决方法。基于Roe格式及预处理方法，提出了全速度Roe格式。全速度Roe格式是将Roe格式与低速Roe格式通过一个与本地马赫数相关的函数联合而得到的。与标准Roe格式相比，对亚音速流动，全速度Roe格式中声速对数值耗散的影响随着马赫数变化而变化。与预处理Roe格式相比，全速度Roe格式简单易于实现，而且没有计算稳定性的限制，不需要采用全局截断处理，从而在理论上对高低速流动具有一致精度。事实上，全速度Roe格式意味着可以通过直接改变Roe格式的非线性特征值，解决可压不可压流动统一计算精度的问题。这是对统一算法的新观点。数值模拟表明,全速度Roe格式具有与预处理Roe格式一样的激波捕获能力、极低速流动准确模拟能力。     

DCS5在计算气动声学中的应用研究

对流项采用五阶线性耗散紧致格式(DCS5)离散,粘性项采用6阶中心紧致格式离散,时间项采用5级5阶龙格-库塔时间推进格式求解流体动力学方程组,并由此建立近声场的直接模拟方法。为了应用结构网格在复杂几何形状条件下对流场进行计算,发展了高精度的多块特征对接边界技术。完成了低速无粘振荡圆柱和粘性圆柱绕流等标准计算声学算例,获得了令人满意的结果。     

一种求解透平叶栅三维流场的高精度TVD格式

本文提出了一种应用三阶精度Chakravarthy-OsherTVD格式求解环形透平叶栅三维定常流场的隐式近似因子分解方法。控制方程为柱座标速度分量表示的任意曲线座标系下的三维欧拉方程。采用Beam-Warming近似分解, 对隐式项进行迎风差分, 得到了时间精度为一阶、空间精度为5点三阶的隐式格式。算例对5点三阶和5点二阶精度的格式进行了比较。表明空间三阶精度隐格式对计算结果没有明显改进, 但收敛残差有较大降低。      

基于5阶精度格式WCNS-E-5的p-multigrid方法研究

p-multigrid方法的基本思想是:在保证收敛结果为高阶精度的同时,利用低阶精度格式耗散大的特点,来改善高精度有限差分格式在迭代计算时收敛速度慢的弱点.本文基于5阶精度WCNS-E-5差分格式,引入1阶精度迎风格式和3阶精度加权格式,构造了p-multigrid方法,在迭代过程中采用了V循环、W循环、S循环、PreV和FMG循环等不同方式来应用这三种格式,并通过典型算例考察了这些循环方式对收敛速度的影响,初步数值试验表明,采用恰当的循环方式,本文所设计的p-multigrid方法能够加快收敛速度,并保证了最终收敛结果与5阶精度WCNS-E-5差分格式的一致性.