S_1流面H型与C型双重网格计算

本文发展了一种S_1流面的H型和C型双重网格计算方法,即首先采用H型网格进行全流场计算,再在叶片前缘附近采用局部C型网格进行加密计算,在生成C型网格时采用了一种组合二次曲线拟合方法以便适用于各种压气机及透平叶型.文中给出的算例表明,本方法是提高S_1流面亚声速绕流计算精度和预测叶栅变工况性能的有效手段。     

高架污染绕建筑对流扩散风洞实验与数值仿真

采用激光衍射瞬时浓度场测量系统和图像处理系统,对不同高架与周围建筑高度比下污染点源绕建筑对流扩散进行了风洞实验.Fluent双相耦合模型与实验进行了相互校核,得到了建筑的各个高度层面污染分布、不同高度比下的污染扩散及污染分布随时间变化的规律.不同高度比对污染对流扩散具有明显的影响,建筑各高度层面污染浓度分布也存在明显差别.当高度比较小时,污染将堆积在建筑物底层.污染分布随时间变化有助于了解污染对流扩散的时间效应.结果可以为高架与周围建筑规划提供参考与依据.     

随机过程数值仿真的精度实验与分析

为了考核随机过程数值仿真中采用的混沌多项式展开,小波动态自适应网格等的可靠性和近似程度,本文构造了一个具有解析解的随机过程,随后进行数值仿真。比较了数值仿真结果与解析解,证明两者是相互一致的。说明所采用的数值仿真方法是成立和近似度是可以接受的。所进行比较的随机过程解析解是相对比较简单的,当对更复杂的情况需要进一步数值实验验证。在文中并讨论了影响数值仿真精度和增加数值仿真工作难度的因素。     

一些完全非线性函数的原像分布值

确定完全非线性函数的原像分布值,是决定和分析完全非线性函数以及构造相应线性码的公开问题和重要课题之一.本文讨论了完全非线性函数的原像分布所满足的基本方程的求解问题,完全解决了该方程当m=5以及m=6时的情形.     

使用任意非正交曲线座标系的任意翘曲S_1流面正问题流场矩阵解

本文包括了两个部分内容(1)根据给定的几根流线上的离散点拟合一个光滑的、任意翘曲的S_1流面,并在此曲面上进行计算网格加密.(2)在此任意翘曲的S_1流面上进行亚声速叶栅绕流计算,并给出了计算实例.本文将用此方法拟合的曲面同解析解进行了比较,其误差为0.004％,可见是相当满意的,而且拟合并不复杂,适合于工程上使用.此曲面拟合的方法也适用于其它工程计算。     

给定激波模型的叶栅跨声流的计算

本文提出了一个已知激波形状的情况下计算跨声速叶栅流场的方法,通道中只有一道激波并且激波后的气流速度低于声速。对超声区和亚声区分开进行计算,两区连接处满足激波关系.典型算例表明,本方法能迅速而精确地决定整个叶栅流场,得到的压力系数和唯一进气角都较好地符合实验值。计算结果不只提供了整个流场的流线分布和马氏数分布,并且还提供各流线上气流各参量变化的情况,特别是通过激波时的突跃变化。     

任意迴转流面叶栅出气角的确定

积分区域ABGD参见图1。 (1a),(2a),(3a)构成封闭方程组。 用其进行迴转面叶栅绕流计算,实际是求解椭圓型微分方程的边值问题。 类似于文献[2]迭代求解,为在计算过程中求解出气角因而增加了出气角渐近循环过     

使用非正交曲线座标与速度分量S_1流面正问题流场矩阵解

本文基于吴仲华提出的使用非正交曲线座标与相应的非正交速度分量的叶轮机械三元流动气动热力学基本方程组,引入流函数,得出求解的主要方程:流函数的二阶拟线性偏微分方程.除了与密度有关的项以外,流函数的各阶导数都置在方程的左端.这样加快了收敛的速度.用中心九点差分格式,将微分方程离散化后,所得的线性代数方程组用矩阵[L][u]分解直接求解.这种解法收敛速度较快.系数矩阵为对角线带状稀疏矩阵.采用了:(1)非零元素按对角线编号;(2)增设虚点两项办法.大量减少了计算机内存量.由流函数求密度时采用了内存密度函数表插值方法.简单地讨论了松弛因子的选取.用此程序对一些压气机和透平的叶型进行了计算,同实验结果及理论解析解进行了比较,相互是一致的.     

二维对流扩散方程的高精度全隐式多重网格方法

提出了数值求解二维非定常变系数对流扩散方程的一种时间二阶、空间四阶精度的三层全隐紧致差分格式。为了加快迭代求解隐格式时在每一个时间步上的收敛速度，采用多重网格加速技术，建立了适用于本文高精度金隐紧致格式的多重网格算法。数值实验结果验证了本文方法的精确性、稳定性和对高网格雷诺数问题的强适应性。     

跨声速叶栅流动中激波与层流边界层干扰的分析与计算

基于实验结果,对流动模型进行了分析。提出了叶栅通道内激波系中“主要分支”的概念。此主要分支位于激波系的前沿。边界层外部流动M数在此发生最剧烈的下降。由于激波系的不稳定,叶栅实验中的非二维性,光学玻璃窗的折射等原因,使纹影照片中激波阴影加宽,而阴影的宽度并不体现激波强度。对包括激波与层流边界层干扰的跨声速边界层的特点进行了讨论,并编制了一个计算程序。用该程序计算了一个典型的超声速入口,M_i=1.34,跨声速叶栅边界层的发展。其结果与实验数据相比,有很好的一致性。进一步证实了所提出的主要分支概念,该计算简便,在UNIVAC-1100机上仅需2—3分钟CPU时间。