一种基于数值计算的轴流压气机旋转失速稳定性模型——之二:简化模型及应用

提出了一个基于直接数值计算结果基础上的旋转失速线性稳定性模型,描述了如何从基本的Navier-Stokes方程出发,把轴流压气机流动稳定性问题转化为一个特征值问题的基本理论框架,对压气机不稳定边界的判断通过计算和检验特征值的方式来实现,给出了有限差分法和基于Gauss-Lobatto节点的谱配置方法两种数值离散的方式,并分析了这两种离散方式的收敛特性.     

基于非结构网格有限差分法的扎染算法

应用有限差分法时遵循贴体坐标系下离散等价方程的离散准则,格式中仅用到当地网格点的坐标变换系数,据此可以构建新的非结构网格有限差分法,数值算例表明在空间二维情况下用连接离散点的3条网格线构造出来的一阶迎风格式可以稳定收敛。把这种非结构网格有限差分法推广到三维,在4条网格线基础上进行离散计算,然后提出一种局部区域多重网格、重复计算的新算法。首先在包含物体的整个区域采用直角坐标系下的均匀网格进行计算,然后消除物体内部点对物体外部流场区域的影响,最后在物面和直角网格之间很小的局部区域填补非结构网格进行计算,计算过程类似于中国传统的扎染工艺。处理流场内存在曲面物体边界问题时,相较于常规结构网格差分算法,扎染算法尽管在计算中包含需要特殊处理的无用网格点,也存在重复计算的局部区域,但是以这些计算能力的浪费为代价,换取了编程简便、内存量小、网格生成快捷、易于扩展网格规模等优势。定性分析表明,这种扎染算法非常适合研制大规模并行计算,建立的计算流体力学应用软件可以充分发挥数十万至百万处理器核心的超级计算机的效能。     

再入飞行器烧蚀热防护一体化计算方法

针对再入飞行器烧蚀热防护系统烧蚀与瞬态温度耦合响应预测问题,提出了一体化计算方法,为再入飞行器烧蚀热防护设计提供包括气动热、烧蚀后退、瞬态温度响应在内的动态响应预测依据。该方法采用Sutton-Graves和Tauber-Sutton理论计算驻点的对流热流和辐射热流,通过表面能量平衡整合具有较高精度的烧蚀模型,并通过Landau变换简化烧蚀后退带来的节点删除过程并保证空间离散精度,最后求解瞬态有限差分热传导方程获得烧蚀热防护系统的热环境、烧蚀过程和温度响应。通过对比计算碳-碳材料钝头体地球再入过程和酚醛浸渍基碳烧蚀体（PICA）材料电弧风洞烧蚀模拟,对该方法对于不同材料体系的适用性进行了验证。计算结果表明：对于密度较高的碳-碳材料,本文计算结果与经典的热平衡积分法吻合较好,偏差在7%以内;而对于低密度材料（如烧蚀性能对压力高度敏感的PICA材料）,随着热流和压力的增大,预测偏差逐渐增大。所提出的方法实现了气动热、烧蚀、瞬态温度响应耦合过程的一体化计算,在保证精度的前提下实现快速计算分析,为再入飞行器烧蚀热防护设计提供依据。     

飞行器气动与隐身综合特性数值分析

开展了飞行器气动与隐身综合特性数值研究,编制了相应的计算与分析程序。该程 序的主要 功能包括：飞行器外形参数化和计算网格的生成,气动特性计算（利用有限体积方法求解NS 方程）,隐身特性计算（利用时域有限差分方法求解Maxwell方程）,飞行器气动与隐身综 合特性分析。作为应用的一个实例,对一种带翼钝锥体的气动和隐身综合特性进行了分析, 并设计出两种新的外形,新外形的气动和隐身综合特性与原形相比有明显提高。     

跨流域三维复杂绕流问题的气体运动论并行计算

通过研究求解描述跨流域三维绕流问题的Boltzmann模型方程气体运动论耦合迭代数值格式,分析气体运动论数值计算方法的内在并行性;从变量依赖关系、数据通信与并行可扩展性三方面开展基于离散速度空间区域分解计算的研究,发展求解稀薄流到连续流跨流域三维复杂绕流问题的并行算法。通过对不同Knudsen数、不同马赫数、不同攻角三维球体及返回舱绕流的并行计算,计算结果与实验数据和理论分析吻合较好。研究表明,该并行算法负载平衡和并行可扩展性较好,对不同并行计算机系统具有很好适应性,显示直接求解分子速度分布函数的气体运动论计算方法有良好的并行计算特性。     

振动翼型跨音速绕流的数值模拟

从积分形式的非定常Euler方程出发,在固联于振动冀型的贴体坐标系下,采用有限体积法进行离散,并根据对称型TVD格式构造相应的数值通量,数值模拟了振动翼型的非定常跨音速绕流的变化过程。计算结果表明,该方法能较好地捕捉流场中的运动激波,准确地定出激波位置与强度;并且易于推广到求解任意振动方式的跨音速非定常绕流问题。     

跨音速喷管流动的数值模拟

数值方法采用的差分格式是具有高分辨率和快速稳定收敛性质的全变差衰减(TVD)格式。喷管壁面形状和反压可任意给定,对于典型轴对称与二维喷管(即壁面由收缩角为45°,扩张角为15°的两条直线用圆弧光滑连接而成。圆弧的无量纲半径等于0.625,其参考长度为喉道半宽),进行了有激波与无激波流动的数值模拟。计算结果与实验结果以及其他计算结果符合良好,此法可推广到非定常跨音速喷管流动的计算,并可用于工程中喷管设计。     

二维US-FDTD方法的数值稳定和色散分析

研究了二维US-FDTD方法的数值稳定性和数值色散特性。通过对增长因子的计算，证明了US-FDTD方法的无条件稳定性。利用增长因子的相位，推导出了US．FDTD方法的数值色散关系式。分析了US-FDTD方法的数值色散误差。数值分析表明，与ADI-FDTD方法一样，数值色散误差仍然是决定US-FDTD时间步长选取的关键因素。同时发现，数值色散受时间步长及网格大小的影响。     

热噪声试验中平板温度边界层对近壁面声压分布影响数值分析

为了研究高温平板隔热试件附近温度边界层对近壁面声场的影响,进而对地面模拟试验条件进行适当修正,文章首先构建了热噪声试验中声波穿过温度边界层的简化模型,运用边界层积分法得出常温气体掠过高温壁面时温度边界层内的温度分布;在此基础上,采用有限差分法求解平面波入射非均匀介质时的散射问题,得出高温平板近壁面声压分布。结果表明,温度边界层使近壁面声场发生衰减,在入射声波频率为10 000 Hz时比无温度边界层情况下声压级减小约3 dB,使实际声谱偏离试验条件。     

退火算法在电子元件位置优化上的应用1)

一块电路板上有电子元件，它们不仅有热源，而且相互作用，另外与空气还交换热. 对这类问题建立线性数学模型，从而算出温度稳态分布. 因为许多集成电路系统的可靠性很大程度上取决于温度，所以温度的模拟显得很重要.把温度模型与退火算法相结合可以获得元件最接近优化的布局，并用一个例子加以说明，通过优化能降低板子最高温度，从而提高系统的可靠性.