低马赫数流动分区并行计算研究

在由文献[1]提出的低马赫数流动数值模拟方法的基础上,本文开展了该方法的分区并行计算研究,发展了一套通用性较好的低马赫数流动分区并行数值模拟程序.该程序在MPP并行计算机上获得了正确的计算结果和较好的收敛速度,取得了很好的加速比.     

流管实验装置传声器阵列的校准方法

发展了一种应用于流管实验装置的现场传声器阵列校准方法.该方法利用两支标准传声器的测量结果计算管内声场,然后比对待校准传声器的测量信号得到其灵敏度.与传统校准方法相比,该方法具有实时快速、能提供相位校准等优点.在流管装置上进行实验,用模态匹配方法对校准结果进行校核,显示该校准方法得到的声压与计算结果吻合得更好.使用本方法,利用高精度传声器对低精度传声器进行校准并应用,可有效降低测量实验成本.      

工业计算机断层成像系统转台旋转中心的确定

构建一个工业CT(computed tomography)成像系统,转台旋转中心的确定是非常关键的步骤,其定位误差会引起CT图像上的伪影.在详细分析了现存旋转中心确定方法不足的基础上,提出了一种新方法,该方法基于正弦图数据的冗余性,通过计算将正弦图数据总和均分的位置来定位旋转中心.相对现有算法,该方法无须使用模体,亦无需知晓任何几何参数,且基本不受随机噪声影响.实际的实验数据验证了该方法的有效性.      

某涡轴发动机整机逼喘试验研究

根据某型涡轴发动机研制的需要,在整机试验台进行了逼喘试验方法的探索,本文详述了发动机逼喘、试验方法设计、试验装置及试验过程,最后对试验结果进行了分析,研究结果表明利用该发动机上的飞机引气口接入高压气源的方法能够成功地进行该型发动机的整机逼喘试验,从而能够真实获得压气机装在整机上工作时的喘振边界线.      

低耗散高可信度格式及其在气动力热计算中的应用

本文简要介绍了TVD类型格式中近似黎曼解的耗散机理,该耗散机制对于粘性问题的计算有较大的负面影响.通过引入低耗散的格式,所发展的数值方法对粘性流场的气动力热问题的计算具有较高的精度.以该数值方法为基础发展的软件系统,在亚跨超和高超声速领域中的气动力热问题上均得到了验证,计算结果与实验进行了对比验证,表明该计算软件系统具有较高的可信性.     

在大射程导弹上改装捷联惯导系统

主要介绍了捷联惯导系统的基本概念、特点及其优越性,探讨了在大射程反舰导弹上改装该系统的可行性方法,并阐明了影响及意义.     

侵蚀燃速辨识

本文将控制论中的辨识技术用于确定固体推进剂火箭发动机装药的侵蚀燃速关系式,对燃速辨识方法作了较为全面和详细的研究.通过对假想发动机作燃速辨识,对该方法的可行性作了充分的论证.用燃速辨识法对真实型号的发动机作了侵蚀燃速辨识,获得了同一发动机装药在不同初温下侵蚀燃速关系式.由于燃速辨识仅需利用一条发动机试车压力一时间曲线,并且在电子计算机上用最优化方法给出结果,因而该方法具有迅速、经济、准确等优点,能更实际地研究固体火箭发动机工作状态下的装药侵蚀燃烧特性.而且便于观察过去实验研究做得较少的初温对侵蚀燃速的影响.这是一种在工程应用上确定装药侵蚀燃速的切实可行的方法.     

基于NESO的高超声速拦截弹解耦控制研究

针对高超声速拦截弹姿态控制问题,提出了基于非线性扩张状态观测器(NESO)的解耦控制方法.根据Tornambe的分散鲁棒控制理论,把耦合项和不确定性视为广义不确定项,构造基于NESO的估计和补偿信号并加入到闭环控制律中.理论推导证明了该方法可以保证闭环系统跟踪误差一致有界.在高超声速拦截弹模型上进行了仿真验证,并与传统的分通道反馈控制方法进行对比,结果表明所设计解耦控制方法得到了更好的控制效果,在较大程度上消除了通道耦合和不确定性的影响.     

基于DE-RLSSVM算法的航空发动机传感器故障诊断

为了使约简最小二乘支持向量机(RLSSVM)具有更好的稀疏性和泛化能力,利用微分进化(DE)算法选择RLSSVM的支持向量,提出了DE-RLSSVM算法.在benchmark回归数据集上的仿真试验表明该算法具有很好的稀疏性和泛化能力.然后将该算法用于航空发动机传感器故障的诊断,提出了基于DE-RLSSVM算法的航空发动机传感器故障诊断方法.该方法利用DE-RLSSVM算法对传感器故障进行监测,然后进行定位和隔离.数字仿真结果表明该传感器故障诊断系统能够实现对航空发动机传感器硬故障的检测与隔离.      

弹道中段目标质心位置估计

 质心在目标上的位置是弹道中段目标的重要进动特征,也是目标的质量分布特征,可用于识别弹头和诱饵.针对旋转对称目标,提出了基于窄带雷达回波相位特性的质心位置估计方法.首先,选取目标对称轴上的任意点作为质心位置参考点,定义相应的质心位置参数,并构造以该点为相位参考中心的全姿态散射系数模板;然后,对窄带雷达回波采样序列进行相位匹配变换,搜索匹配函数的最大值,从而得到质心位置参数的估计值.仿真结果表明了该方法的有效性.     

用于气动声学计算的非均匀网格紧致差分格式

为克服传统紧致差分格式在数值求解非均匀网格问题时产生的寄生波,构造了一种新的高精度紧致差分格式.通过泰勒展开分析方法,详细给出了格式系数的通用形式;利用傅里叶分析方法,分析了数值耗散、色散误差.以3对角6阶精度紧致差分格式求解均匀扰动网格问题为例,计算表明:色散值和耗散值随扰动因子的增加而更加趋近于精确值;当扰动因子大于0.213时,格式不稳定,当扰动因子小于等于0.213时,格式渐近稳定;对一维对流波和二维波传播的模拟计算所得数值解与精确解吻合,体现了该格式在求解非均匀网格问题时的优越性.      

气体动理论BGK格式的网格自适应方法

为了提高气体动理论BGK（Bhatnagr-Gross-Krook）格式在超声速流动问题计算时激波捕捉的准确性与计算效率,提出了一种适用于气体动理论BGK格式的网格自适应加密方法。该方法采用基于四边形的链表技术来描述网格的拓扑结构,在物理量重构过程中,使用了在四边形网格中表现优异的van Leer限制器,以保证粗细网格过渡处物理量重构的精度。用跨声速翼型绕流（马赫数Ma=0.85）、超声速前台阶流（Ma=3）和高超声速圆柱绕流（Ma=8.03）等多个典型算例验证了BGK自适应网格方法。计算结果表明,自适应网格BGK方法在保证数值精度的前提下,可大幅度提高计算效率。这为该方法用于高效地解决复杂问题提供了一种选择。     

浸入式贴体网格边界方法

提出一种基于重叠网格的浸入式边界方法,用以模拟复杂的三维黏性跨声速流动。该方法采用一套固定的笛卡儿正交网格,用以主流的求解;采用一套可以移动的贴体等距面网格,用以拟合或者离散物面的作用力,通过空间插值,实现两套网格重叠部分的信息传递。分析了浸入式贴体网格边界方法的优势,介绍了贡献单元的寻找策略和物理通量的插值方法。流场求解采用Spalart-Allmaras带湍流模型的Navier-Stokes方程组,其中对流项采用流通矢量分裂和5阶WENO(weighted essentially non-oscillatory)-Z格式离散,黏性项采用6阶中心差分格式离散,时间项采用龙格库塔显式格式离散。数值验证算例表明:该方法具备高于4阶的空间求解精度,并适用于刚体动网格非定常流场模拟,且无需更新网格形状。等距面贴体网格生成过程简单,避免了繁琐的人工调节过程,与笛卡儿网格结合,可提供足够的壁面附近网格密度,同时有效减少了网格总量需求。     

一种新的用于MHD方程逆风格式的Jacobian矩阵分裂方法

针对理想MHD方程,提出了一种新的基于MacCormack算法的雅可比矩阵分裂方法,克服了原有方法稳定性差的问题,并成功地应用于理想MHD方程的求解.控制方程在非结构混合网格上进行空间离散,其中对流项采用本文发展的逆风向量分裂格式,并引入了双曲型磁场散度清除技术,时间推进为显式5步龙格-库塔方法.对MHD激波管流动和带均匀磁场干扰的二维高超声速钝头体绕流流场进行了数值模拟,得到了与参考文献相吻合的数值结果,表明本文发展的数值分裂方法可以有效地捕捉MHD流场的流动特征,并且具有比MacCormack方法更高的稳定性和计算精度.     

振动翼型跨音速非定常Euler方程数值解方法

通过坐标变换使变化的物面边界条件及无穷远边界条件变为固定的有限范围内的边界条件。因此可在不增加计算机时的情况下,边界条件可以精确地满足。由于主控方程是Euler方程,因此,凡适合于Euler方程解的问题（跨音速、大扰动、非定常、高简化频率）均可应用此法。同时,在计算过程中采用贴体、正交网格与可控网格技术获得最佳网格分布。算例计算结果是令人满意的。     

关于 CFD 高精度保真的数值模拟研究

本文由四部分组成。第一部分研究了利用高阶格式、网格、物理模型和高性能计算机求解 NS 方程时,计算数据的高精度保真性的要求以及高精度保真的指标。第二部分研究了利用大规模计算开展大涡模拟和直接数值模拟存在的问题。第三部分分析了在高雷诺数下 NS 方程计算方法和网格的关系。最后提出了建立计算数值验证、确认的新方法。     

基于激波拟合与柔性壁的反设计数值方法研究

将激波和待求激波压缩面作为运动间断边界拟合,研究了超声速来流下的型面反设计数值计算方法.算法秉承拟合法思想,采用二阶格式时间推进欧拉方程组,在运动边界和动网格条件下,按预估的波系结构划分网格块,反求可产生指定压力分布的三维型面,并获得精确的激波型面.通过包括凸包进气道在内的三个反设计算例验证了本文算法的合理性及其对网格质量的适应性,证明了在型面反设计中,采用激波拟合法可大大提高型面精度.本算法可用于超声速进气道压缩面和唇口完全三维反设计.     

非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

A SELF－ADAPTIVE FINITE ELEMENT METHOD FOR SOLVING 2－D EULER EQUATIONS ON THE UNSTRUCTURED MESHES

ＡＳＥＬＦ－ＡＤＡＰＴＩＶＥＦＩＮＩＴＥＥＬＥＭＥＮＴＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧ２－ＤＥＵＬＥＲＥＱＵＡＴＩＯＮＳＯＮＴＨＥＵＮＳＴＲＵＣＴＵＲＥＤＭＥＳＨＥＳＺｈｏｕＣｈｕｎｈｕａ；ＹａｎｇＺｕｏｓｈｅｎｇ（ＮａｎｊｉｎｇＵｎｉｖｅｒｓｉｔｙ...     

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。