几种典型迎风格式的分析与比较

重点分析了AUSM 格式的构造方法,同时将其与Van Leer, Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM 格式,VanLeer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率.获取了对于基础研究和工程应用很有意义的结论.     

Godunov型显式大时间步长格式研究进展

综述了Godunov型显式大时间步长格式的研究进展。首先介绍了显式大时间步长格式的概念、分类和优势。然后重点阐述了Godunov型显式大时间步长格式的构造方法、高阶精度推广方法、多维问题推广方法和收敛特性、分辨率及计算效率等性能，展示了其在典型问题中的应用和验证。最后给出了Godunov型显式大时间步长格式研究进一步可能的发展方向。     

IN 738LC材料Chaboche热粘塑性本构模型的隐式Euler格式

为了提高Chaboche模型的数值积分效率,推导了该模型的Eu ler隐式积分格式,针对IN 738LC材料,将该积分格式编写成FORTRAN程序,利用通用有限元分析软件M arc的用户子程序接口HYPELA 2,将该隐式积分格式与有限单元法结合起来。对实验进行了数值模拟,并将隐式算法和显式算法与实验结果作了对比,结果发现隐式算法是合理的,比显式算法收敛性要好。      

火箭发动机塞式喷管流场的数值研究

用时间相关法和张涵信的无振荡、无自由参数ＮＮＤ显式格式,对火箭发动机塞式喷管二维轴对称定常流场进行了数值模拟。计算得到了有粘和无粘中两种情况下,喷管流场的马赫数压强、温度、流线以及速度矢量分布图,计算网络采用代数方法生成,并通过求解Ｌａｐａｃｅ方程对网格进行修正。计算结果表明：ＮＮＤ格式用于计算塞式喷管流场具有较高的效率和良好的分辨率;无粘和有粘情况得到的流场结构有较大的差别;基础喷管的出口马赫数     

用虚拟黏性法构造Navier-Stokes方程黏性项的隐式格式

提出的虚拟黏性法是一种关于Navier-Stokes方程黏性项(黏性应力项和传热项)的隐式新方法。通过引入虚拟时间和虚拟黏性项,将隐式格式的构造大大简化,从而避免了大型复杂隐式差分方程组的常规求解。在虚拟时间推进过程中所需求解的方程组的系数矩阵是一个三对角矩阵,它具有计算简单且计算量小的优点。用模型方程和Navier-Stokes方程进行了数值仿真,研究各种参数对计算的影响,并在精度和效率上与显式方法进行比较,证实了算法的正确性和优势。最后分析了本文格式的适用范围。     

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

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

三维Euler方程组隐式LU分解的高性能并行计算

许多非定常无粘流体力学问题的数值模拟都需要利用Euler方程组来进行计算,而由于在隐格式下,所选取的时间步长可以比在显格式下时大得多,所以隐格式越来越受到重视,其中隐式LU分解是最常用的方法之一。对三维Euler方程组,采用隐式LU分解进行计算时,网格点所在的各个对角阵面之间存在数据依赖关系,本文分析了采用区域分解且边界上用显格式代替隐格式进行计算的高效性,在长方体建筑物内的爆炸模拟表明,在有112个CPU的某MPP巨型机上,并行计算效率超过60%。本文还分析了计算结果与串行计算时的差异,以及利用区域重叠减小这种差异的方法,同时考虑了对处理器进行合理的逻辑组织,将计算网格映射到处理器网格,以最大限度减少通信开销的方法。文中最后以一个爆炸毁伤的例子实际说明了所述方法的可行性与高效性。     

压气机叶片排内三维紊流流动的隐式高分辨率数值分析

对跨音速压气机叶片排内具有复杂三维激波的紊流流场所进行的数值分析，采用了空间三阶精度无波动格式以保证对紊流特征的正确模拟，利用ＬＵ－ＳＧＳ算法避免了一般隐式算法的繁琐运算，使之与显式方法推进一步的运算工作量基本相当，从而构成了一种既准确又高效的跨音速压气机紊流流动数值求解系统。计算与实验进行了对比，结果是令人满意的。     

点隐式龙格-库塔方法的应用研究

为了提高求解器的效率,在显式龙格-库塔时间推进的欧拉方程求解器之上,发展了点隐式龙格-库塔时间推进格式。给出了其推导过程和非结构网格下中心格式和迎风格式（Roe和Van Leer格式）预处理算子的构造方法。NACA0012翼型和RAE2822翼型的跨音速无粘流动模拟表明：与显式龙格-库塔方法相比,方法能提高求解效率且内存需求相当,具有一定的工程应用价值。     

涡扇发动机加减速特性显式与隐式计算方法

为了探索一种便于进行涡扇发动机过渡态控制规律设计的性能计算模型,提出了基于部件法的涡扇发动机加减速的显式格式和隐式格式计算方法,该方法通过在发动机计算模型中直接给定喘振裕度限制值、燃烧室油气比限制值和涡轮进口总温的限制值,计算出最优的加减速特性,进而获得发动机的最优加减速控制规律。计算模型针对不同的给定值,选择了不同的燃烧室容积效应模型。证明了对一般的涡扇发动机,隐式格式计算模型中,给定压气机喘振裕度算法的解是唯一的。以某涡扇发动机在地面的加减速过程为例,按最优加减速控制规律计算,显式格式算法和隐式格式算法的结果误差小于1.3%.对给定高压转子转速加速率的加速特性也进行了验算,计算结果与最优加速过程的结果误差小于1.7%.本文提出的加减速特性计算方法可为涡扇发动机的过渡态开环和闭环控制规律设计提供便捷的手段。     

椭圆的数控加工宏程序编制

根据华中世纪星数控系统及固定循环指令特点,将椭圆型腔的粗加工、精加工程序编制成固定循环格式,用户在加工过程中只需要用指令调用即可。这样在加工过程中大大缩短程序编制时间,提高数控宏程序的集成度,使数控机床高效、高精度的潜力被充分挖掘。     

基于QSS的自适应多步校正算法及其在柔性航天器动力学中的应用

量化状态系统（QSS）算法是一种基于状态变量离散化的数值积分方法,该方法与基于时间离散的传统方法显著不同,QSS通过计算状态变量每次跃迁所需的时间来推进下一步的积分。在求解非刚性常微分方程时,QSS算法比传统算法更具优势,但它不适合求解刚性问题。为此提出一种基于QSS的自适应多步校正算法（AMCQSS）,该算法以QSS为基础、结合隐式多步法思想,在计算过程中可以自适应选择二步法或三步法,以有效提高求解刚性问题的精度和效率。通过对柔性航天器动力学的仿真求解,验证了算法的可行性。将该算法与ODE23tb、ODE15s、ODE45以及QSS等算法进行对比,结果表明AMCQSS算法既能保证求解的效率及精度,又具有较好的收敛性和稳定性。     

基于时间精确自由尾迹方法的旋翼非定常 气动特性分析

 给出了一种显式修正的二阶向后差分格式(PC2B)的时间精确自由尾迹算法,用于分析旋翼非定常气动特性。由于旋翼自由尾迹控制方程具有非线性特征,求解过程中易产生数值不稳定性问题,因此本文推导了修正PC2B算法的修正方程,分析了该算法的求解精度和数值稳定性,并通过数值试验对其稳定性和计算效率加以验证及分析。基于该算法,首先对悬停和总距突增时旋翼非定常气动特性进行了计算,通过与试验结果的对比验证了模型的有效性;其次研究了总距操纵时旋翼非定常气动特性的变化规律。研究结果表明,旋翼自由尾迹几何形状的瞬态变化是造成总距操纵条件下旋翼非定常气动特性波动的主要原因。     

喷管跨声速流场计算

本文用时间相关法求解定常跨声速喷管流场。对守恒型和非守恒型方程采用全场统一时间步长、当地时间步长和修正超阻尼法进行了计算。计算表明:采用当地时间步长可大大提高收敛速度,采用全场统一时间步长时,数值解非守恒型方程收敛较快;采用修正超阻尼法也能提高收敛速度,但与阻尼系数和时间间隔有关。     

基于偏导数的全局灵敏度指标的高效求解方法

在全局灵敏度分析领域,基于偏导数的全局灵敏度指标由于其优良的特性得到了广泛的关注。针对目前求解基于偏导数的全局灵敏度指标计算效率低的问题,提出了一种高效的求解方法。该方法利用乘法降维近似地将响应函数展开为连乘积的形式,从而将求解偏导数全局灵敏度的高维积分问题转化为一维积分的连乘积,在保证计算精度的前提下大大降低了求解基于偏导数的全局灵敏度指标的计算量。在求解特定点的偏导数时,采用了复数步长方法,提高了计算精度。最后,通过算例验证了所提方法的准确性和高效性。     

跨音速大迎角Euler方程数值分析

应用有限体积法离散三维Euler方程,选用三步显式格式进行时间推进求解;并通过当地时间步长、残值光顺、焓修正等方法加速收敛。对一维波动方程的von Neumann稳定性分析说明这种三步格式的最大Courant数是2。对三角翼跨音速大迎角流动的数值分析表明文中的三步格式及焓修正方法优于人们经常采用的四步、五步格式以及焓阻尼技术。     

Numerical simulation of a cabin ventilation subsystem in a space station oriented real-time system

An environment control and life support system(ECLSS) is an important system in a space station. The ECLSS is a typical complex system, and the real-time simulation technology can help to accelerate its research process by using distributed hardware in a loop simulation system. An implicit fixed time step numerical integration method is recommended for a real-time simulation system with time-varying parameters. However, its computational efficiency is too low to satisfy the real-time data interaction, especially for the complex ECLSS system running on a PC cluster. The instability problem of an explicit method strongly limits its application in the ECLSS real-time simulation although it has a high computational efficiency. This paper proposes an improved numerical simulation method to overcome the instability problem based on the explicit Euler method. A temperature and humidity control subsystem(THCS) is firstly established, and its numerical stability is analyzed by using the eigenvalue estimation theory. Furthermore, an adaptive operator is proposed to avoid the potential instability problem. The stability and accuracy of the proposed method are investigated carefully. Simulation results show that this proposed method can provide a good way for some complex time-variant systems to run their real-time simulation on a PC cluster.     

飞行器控制系统的高精度计算方法

对于连续时间线性二次最优控制问题，在动力学方程和价值泛函的基础上，给出其全状态下的微分方程，进而将精细时程积分法引入上述问题，其优点为，放弃了传统的差分算法，使计算过程既简便又稳定；避免了Ｒｉｃｃａｔｉ代数方程的求解；具有非常高的计算精度。通过对某飞行器控制系统的计算，充分说明了上述特色。