两种椭圆型方程求源项方法在喷管内流场网格生成中的应用

用椭圆型方程生成二维网格，发展两种Hilgenstock求源项方法，实现对网格与边界间距及与边界正交的控制。     

非交错网格求解传热及流体流动

在非交错网格系统上提出了一种新的以控制容积为基础的有限差分方法，来求解传热及流体流动问题。为了避免压力棋盘问题，压力梯度采用偏迎风压力离散格式进行处理，以保持压力－速度的耦合。利用本算法对数值传热的典型问题进行计算，并将计算结果与采用交错网格的PHOENICS通用程序所得结果相比较，二者吻合甚好。     

用隐式多重网格法计算三维粘性流动

求解高雷诺湍流流动时,边界层法向网格间距较流向和展向相比非常小,因此边界层中存在高度伸缩的网络,这将大大降低多重网格的求解效率。通过谐调的处理多重网格过程的各个细节,既提高了解的精度,又克服了网格展弦比的影响,使得计算效率提高了6～7倍。     

优化程序 服务科研 发动机总公司评出优秀S1流面计算程序

经过近一年的酝酿、筹备工作,49名专家、教授和工作人员于 10月云集四川江油 624所参加 S1 流面计算程序上机考核和专家评审活动。通过10天紧张的调机、上机运行和专家评审工作,评选结果于11月1日正式揭晓。进行如此重大的程序评审活动在我国航空发动机行业尚属首次。评选出的优秀程序是 :涡轮综合第一 :叶栅 S1 流面多重网格欧拉解通用计算程序 (H网格) (中科院工程热物理所刘建军);压气机亚音速第一 :S1流面流函数解计算程序 (中科院北京科能中心华耀南);涡轮跨音速第一 :叶栅 S1 流面多重网格欧拉解通用计算程序 (C网格)。      

转动坐标系中三维跨声速欧拉流的有限体积 TVD格式

在非惯性转动坐标系中,本文采用贴体网格、有限体积法离散和修正数值通量技术,将Harten的一维TVD格式推广到三维。由于转动使方程出现源项,文中通过对源项的巧妙处理,使修改后的格式能用于非齐次双曲守恒律方程组高分辨率的数值计算。为了加速解的收敛,提高显式时间推进的CFL数,本文采用隐式残值光顺技术。三维跨声速带非齐次源项欧拉方程的典型算例表明:捕捉的激波分辨率较高;激波前后没有发现大的数值波动和伪振荡现象;所得的跨声速流场解与实验较接近。     

亚音速复合倾斜透平静叶栅的理论与实验研究

采用三维全势函数多重网格数值方法求解亚音速透平静叶栅流场。以 30 0 MW汽轮机调节级静叶栅的设计为目标 ,采用子午面上外环壁收缩和叶片周向复合倾斜两项措施 ,以求大幅度地减少叶栅二次流损失。实验研究表明新设计叶栅的损失系数比直叶栅降低 42 % ,测出的叶片表面静压分布与理论计算结果十分吻合     

转子分区循环对称接触应力分析

 通过光弹实验表明了端齿连接转子受力的分区循环对称性。针对这~特性建立了各循环对称分区的循环对称接触有限元应力分析方法,并利用叠加原理将各循环对称分区间过渡区化为几个小规模的循环对称问题。     

斜齿轮齿根动应力数值计算与实验研究

运用三维有限元和齿根局部网格密化技术提出了一套计算斜齿轮啮合过程齿根弯曲应力的方法。根据齿轮系统动态特性分析, 求得作用在轮齿上的动载历程, 计算了在动态载荷作用下的齿根应力分布。通过斜齿轮齿根动应力实验, 说明了本文方法的有效性。给出了啮合过程中齿宽不同截面上的动应力分布形态, 讨论了转速和扭矩对齿根动应力的影响。      

多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

二维激波斜壁问题的数值研究

1.前言 在数值模拟激波构造时,通常将流体看作理想流体,即对Euler方程用有限差分法求解。有限差分法的优点是通用性强,但将其用于计算激波问题时,会在不连续波面附近产生数值衰减、过冲或欠冲,难以进行精确计算。FCT(Flux-Corrected Transport Method)法是Boris等开发的一种新的计算方法,它能充分发挥有限差分法的长处,可有效地计算激渡问题。但是,能对二维复杂形状的绕流进行数值研究还是在Thompson等