旋转直管流场的数值模拟

根据旋转坐标系下的 Favre平均 NS方程以及旋转条件下的流动特征 ,修改了 B/ L湍流模型 ,并提出了一种旋转条件下出口边界压力分布的给定方法 ,对已有的 NAPA通用湍流计算软件进行了改进和发展 ,并计算了旋转直管中的流动。通过与实验数据的比较 ,文中的计算结果与实验结果吻合较好 ,正确反映了流动特征和物理现象 ,验证了修改后的 NAPA软件的可靠性。     

低速高湍流度90°弯管流动数值模拟

分析了曲率对弯曲管道流动的影响 ,给出了曲率修正的双层 k-ε湍流模型 ,并数值模拟了非均匀来流低速高湍流度 90°弯管内流动。文中比较了曲率修正双层 k- ε湍流模型的计算结果与实验结果 ,以及双层 k- ε湍流模型和逆压力梯度修正双层 k- ε湍流模型的计算结果。通过比较 ,发现经过曲率修正后的双层 k- ε湍流模型更好地模拟了非均匀来流低速高湍流度 90°弯管内流动。同时比较了壁面参数的插值方法 ,发现在弯曲管道流动模拟中对壁面参数进行一阶插值比零阶插值能更好地模拟真实流场     

一类非线性比例积分控制规律及其应用的研究

提出一类非线性比例积分控制规律，理论分析和应用结果表明，它能显著改善控制系统的静态和动态性能。     

S弯扩压器内旋流的自动抑制

旋流的自动抑制意在建立一控制系统对飞行中进气道内旋流进行实时控制,以求在各种飞行状态下,进气道内旋流均较小.文中首先对旋流的监控进行了广泛、深入的探讨和分析,给出了监控参数与旋流之间的关系及控制措施对其关系的影响.接着建立了一套简单的旋流自动抑制系统,其中以进口段可调导流叶片作为控制措施,用步进电机带动叶片转动,通过对旋流状态的监测,由一台386微机计算、比较并调整叶片的角度.实验证实,该系统工作正常,当气流攻角不断变化时,自动抑制系统不断调整叶片角度.对抑制前后旋流的测量表明,该系统可在各种攻角情况下对旋流进行有效抑制,抑制后的旋流不再有明显的单涡旋流特征,当攻角|a|≥40°时旋流可减少70%左右.     

一种改进的MUSCL格式

给出了一种改进的二阶精度的ＭＵＳＣＬ格式，并用来求解Ｅｕｌｅｒ方程，改进后的ＭＵＳＣＬ格式在线化插值步中，用Ｖａｎ Ｌｅｅｒ的插值斜率对特片变量进行插值，并对与熵波、剪切波对应的特片变量使用了小耗散的插值斜率，在迎风步采用Ｒｏｅ的近似黎曼解求得迎风通量。比较一维和二维的计算结果均可以看出，改进后的ＭＵＳＣＬ格式在消除数值振荡的同时，明显提高了流场间断的分辨率。     

S变扩压器内旋流的自动抑制

旋流的自动抑制意在建立－控制系统对飞行中进气道内旋流进行实时控制，以求在各种飞行状态下，进气道内旋流均较小。文中首先对旋流的监控进行了广泛、深入的探讨和分析，给出了监控参数与旋流之间的关系及控制措施对其关系的影响。接着建立了一套简单的旋流自动抑制系统，其中以进口段可调导流叶片作为控制措施，用卞进电机带动叶片转动，通过对旋流状态的监测，由一台３８６微机计算、比较并调整叶片的角度。实验证实，该系统工作     

DSMC方法模拟微尺度流动时的特征参数选取

利用DSM C算法对微尺度平板附面层流动进行了数值模拟研究,分析了算法中3个主要特征参数:采样循环数、单位网格粒子数和子网格数对模拟结果的影响。分析结果表明,当采样循环数增加到一定程度后,计算结果趋于一致,考虑到计算机时,采样数取8左右为宜;单位网格粒子数过多或过少都会使计算机时增加,太少了还会影响计算精度,建议4~5个左右合适;子网格数对运算时间影响很大,分析认为每个网格内的子网格数为4时最为节约计算时间。最后通过微通道内流动计算进一步验证了本文参数选择方法的可行性。     

对称带状矩阵的并行Cholesky分解及相应线性方程组的并行计算

大型结构问题所导出的方程组系数矩阵阶数往往非常浩大,传统的串行计算机受存储容量与计算速度限制往往难以处理。本文给出适合寄存器—寄存器加工方式流水线向量机上对称带状矩阵三角分解的并行算法MPLDLT和对称带状线性方程组求解的并行算法MCSA。在YH—1机上通过对实例的计算表明,算法是高效的。当矩阵的阶数仅力1666阶时,算法MPLDLT比相应串行算法计算速度快25倍,算法MCSA比相应串行算法计算速度快47倍。若结合YH—1机的特点,使用向量“链接”技巧,则算法MPLDLT比相应串行算法的计算速度快74倍。