关于一维非定常气—固二相流特征理论的讨论

本文研究了一维非定常气－固二相流的特征，对于目前常用的、包含一定近似的运动方程组，它有两个虚部很小的复根，研究表明这种很小的虚部出现与方程中引人的近似有关。从十多个算例的结果看，这样的虚部不影响人们在工程计算中该方程组求解初值问题。     

二次空气系统一维非定常计算方法

研究了燃气轮机二次空气系统的非定常计算方法.通过将整个二次空气系统划分为两部分,特征"长"元件以及局部损失元件.特征线法用于单个特征"长"元件的求解,通过将连续方程、动量方程、能量方程构成的偏微分方程组转换为常微分方程组并联立求解,可以充分考虑内部流动与换热的相互影响;基于压力修正的流体网络法则用于由多个局部损失元件组成的局部流体网络,首先对动量方程组进行求解,进而对各自的能量方程进行求解.如何在采用这两种方法进行求解的上述两类元件之间进行数据交互将被讨论.此外,所提出的时间步长定义方法以及改进的元件计算模型使得二次空气系统的非定常分析得以实现.计算程序经与某系统定常实验数据及单个元件的非定常计算数据进行对比,在10s的计算时间内,基于特征线法得到结果与CFD数据偏差达到10.2%,结果表明该计算程序可以满足二次空气系统非定常计算的要求.     

推进剂供应管道的集中参数近似模型研究——模态近似模型

提出了用模态近似法建立分布参数湍流管道的集中参数近似模型的新方法,发展了适用于水力和绝热气体湍流管道的两种模态近似模型：Ｉ－Ｒ－Ｃ模态和等效层流模型。模型中考虑了频率相关摩擦损失。利用单根管道阀门关断问题的仿真计算对本文的模型进行了验证,计算结果与考虑频率相关摩擦损失的特征线方法的计算结果十分吻合。     

旋转通道一维非定常计算模型

采用特征线法,对航空发动机二次空气系统的旋转通道结构一维非定常流动及换热的计算模型进行了研究,通过将广义的定截面、有摩擦、考虑对外热交换的旋转通道可压缩非定常流动偏微分方程组转换为常微分方程组,考虑旋转离心力,并采用Colebrook White方程定义其阻力特性,进而获得旋转通道内部沿流动方向压力、质量流量、温度等系统参数的特征线数值解,采用商业软件Fluent对其进行验证.结果表明:特征线法利用压力波传播方式对旋转通道进行求解,能够直观地描述旋转通道一维非定常流动及换热特性,计算误差不超过5%,可以作为二次空气系统旋转通道一维非定常计算研究的有效手段.      

高压段含有颗粒的含灰气体激波管流动数值模拟

本文给出了用类特生法计算的高压室含灰的无限长激波管非平衡流动解。以膜压比Ｐ４１＝５．０、高压室含灰的体积分数α４＝０．０１为例，详细分析了非平衡流场，并同冻结流解和平衡流解作了详细的对比分析，清楚地显示了含灰气体激波管非平衡流动从冻结流向平衡流过渡的各种细节，说明了基于冻结声速的非平衡流类特征线法计算何以能得到无限趋近于平衡流解的原因。本文的结果表明，类特征线方法可以提供非常精细的气固二相流非平衡     

推进剂供应管路内液体瞬变流一维有限元计算

将液路计算中偏微分方程模型化成有限阶数的常微分方程组,由一维液体瞬变管流方程的特征线差分格式,提出了一种计算管腔互联结构形式的液体推进剂供应管路中瞬变流的一维有限元方法,它兼有特征线方法的特点。数值计算表明对单根等截面圆管阀门关断问题的计算结果与特征线方法计算结果吻合。利用本文方法计算了一个模型发动机的脉冲工况动态过程,取得了很好的效果。     

姿控推进系统发动机关机的管路瞬变特性

从一维液体瞬变管路方程出发，采用有限差分格式的特征线方法，针对某姿控推进系统发动机真实氧化剂管路的关机瞬变特性进行了数值研究，对系统不同部位的动态特性进行了计算，喷注前的最大瞬变压力达到储箱压力的3倍，得到的结果与试车数据是一致的，对液体发动机管路工程设计具有参考价值，表明该方法可以方便地对包括许多管路，阀门和推力室的复杂系统进行动态特性分析。     

超声速流场的流线-特征线坐标变换与应用

为了求解二维平面弯曲激波波后流场,讨论并发展了一系列基于流线-特征线坐标系变换的流场代数计算方法。该系列方法依据二维平面流场中气流角、静压沿特征线近似不变的特点,可快速求解平面弯曲激波波后流场。其中,等气流角近似法适合模拟流线轨迹,等静压近似法适合求解波后流场参数。在此基础上,又提出了一种改进的等气流角-等静压混合方法用于计算弯曲激波波后流场。和特征线法对比,等气流角-等静压混合法计算得到的流场基本特征与特征线法得到的结果相同,在来流马赫数Ma=6和Ma=4情况下误差分别仅为0.5%和0.15%,证实了该方法在求解二维平面弯曲激波流场中的适用性。     

二次空气系统盘腔一维非定常计算方法

通过对盘腔结构的一维算法展开研究,提出适用于具有非标准结构、包含多个进出口的复杂盘腔的一维建模规则,并采用特征线法,实现盘腔结构的一维非定常流动与换热的计算,通过将广义的变截面、有摩擦、考虑对外热交换的盘腔一维可压缩非定常流动与换热的偏微分方程组转换为常微分方程组,进而获得盘腔内部沿流动方向压力、流量、温度等系统参数的特征线数值解,采用商业计算流体软件CFX对其进行验证.结果表明:所提出的盘腔建模与计算规则能够直观且准确地描述盘腔一维非定常流动及换热特征,仿真结果与CFX计算结果最大偏差约为9.4%,可应用于具有非标准结构的多接口复杂盘腔的模拟,可以作为二次空气系统盘腔结构一维非定常计算研究的有效手段.      

推进剂供应管路充填过程研究

对预充气体和近真空的管道充填过程进行了研究，考虑流体可压缩性和管壁的弹性，以统一的方程形式处理液体和气体的流动；采用带插值的特征线有限差分方法进行数值计算，较好地描述了充填过程中的瞬变特性和水击现象。     

A new accelerating technique for low speed flow: pseudo high speed method

This paper proposes a new accelerating technique for simulating low speed flows, termed as pSeudo High Speed method (SHS), which uses governing equations and numerical methods of compressible flows. SHS method has advantages of simple formula, easy manipulation, and only need to modify flux of Euler equations. It can directly employ the existing well-developed schemes of hyperbolic conservation laws. To verify the technique, several numerical experiments are performed, such as: flow past airfoils and flow past a cylinder. Analysis of SHS method and comparisons with some precondition methods are made numerically. All tests show that SHS method can greatly improve the efficiency of compressible method simulating low speed flow fields, which exhibits in accelerating the convergence rate and increasing the accuracy of the numerical results.     

燃气轮机叶轮非线性瞬态温度场的边界元分析

首先采用克希荷夫变换将包含变物性及混合和辐射边界条件的热扩散方程转换成线性方程，进一步由加权余量法和离散化将方程转变成边界元方程组。用所开发的二维线性—非线性瞬态温度场边界元程序分析计算了燃气轮机叶轮起动工况下的非稳态温度场，与有限元数值计算结果比较表明，边界元法应用于这类非线性工程问题的分析是可靠、初始精度高，数值稳定性好及计算输入数据简单等优点。     

论定常不可压粘流的积分方程解法

本文对定常不可压粘性绕流积分解法中的涡量迭代过程、涡量边界条件的提法以及压强系数的计算等方面都作了改进,使数值解更易于稳定收敛并提高精度。为了验证本文的方法,对低雷诺数下的圆柱绕流进行了计算,结果表明计算效率是高的,迭代过程是稳定的。将本文方法推广到三维流动没有原则上的困难。     

一种求解核函数方程的新方法—辅助核法

本文提出了求解核函数方程的辅助核法。此法克服了偶极子栅格法对网格分布的限制,具有更大的适用范围。数值计算表明,此法对求解绕翼型的跨声速激波流动是很有利的。     

一种基于数值计算的轴流压气机旋转失速稳定性模型——之一:理论模型及算法分析

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

一种高效的多重网格三维N-S方程计算方法及其应用

介绍一种基于 Jameson中心差分和 Runge-Kutta时间推进的三维 N-S方程计算的有限体积方法 ,计算点选取体积单元中心 ,为了提高收敛效率和程序稳定性 ,采用了当地时间步长和多重网格方法 ,特别是针对粘性计算中计算网格形状比高的特点 ,隐式残差光顺采用当地变系数的方法 ,并专门设计了考虑计算网格形状比的人工粘性模型。应用本文方法对收扩喷管的三维流场和NACA0 0 1 2翼型二维流场进行了计算 ,结果表明程序具有良好的收敛能力和较高的稳定性。     

燃气轮机燃烧室中贴体网格的生成与应用

运用贴体坐标系统分别对燃气轮机燃烧室的冷态和热态流场进行了数值模拟。燃烧室的曲线壁面边界用Ｔｈｏｍｐｓｏｎ的非正交贴体坐标网格处理，并把ＳＩＭＰＬＥ算法应用到曲线坐标系下求解各守恒方程。采用双方程ｋ－ε模型来描述紊流特性；采用按简单的单步化学反应假设的Ａｒｒｈｅｎｉｕｓ－ＥＢＵ模型表述紊流燃烧；采用热通量法辐射模型来估算辐射通量和壁面温度。另外，还研究了冷态时不同旋流数和紊流强度对气流结构的影响。     

多级轴流压气机的径向掺混

 随着压气机向着更高的级负荷和更低的展弦比方向发展,径向掺混问题在气动分析和气动设计中越来越重要,当前对决定径向掺混的主要因素众说不一。本文用严格的理论,得到了用统一形式表示含二次流和湍流在内的各个径向掺混因素的多级轴流压气机通流基本方程组,并讨论了二维通流流场湍流掺混与实际流场三维剪切结构的关系。基本方程组的数值结果与试验结果以及其它掺混理论结果比较表呀,该方程组能够满意地模拟多级轴流压气机中的径向掺混特性。     

求解气动弹性系统动态特性的一种准确方法

本文提出一种求解气动弹性系统动态特性的准确方法——正交搜索法,它利用复变函数中解析函数保角映射的特性,沿复频率s的实部和虚部交譬搜索,逐步达到所求复方程的根。该算法能在计算机上实现,直观、简便、实用。     

随机结构孤立特征值的统计特性

研究了随机结构的孤立特征值问题。将材料物理量的随机场扩展为K-L(Karhunen-Loeve)正交展式,采用非正交多项式混沌展式表达孤立特征值,建立了和摄动法类似的一系列确定的递推方程,并通过确定性有限元方法求解了这些递推方程,得到了特征值的均值和方差。在算例中用蒙特卡洛方法验证了本方法的正确性。