优化程序 服务科研 发动机总公司评出优秀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等     

有限体积法的分析与改进

本文对以Jameson为代表提出的常用的一种有限体积法进行了较深入的分析，指出了它的适用范围，网格对计算精度的影响并提出了改进措施。算例表明本文提出的分析与改进措施是正确、有效的。     

氢气泡网格图像测速法—瞬态二维速度场测量

本文介绍一种二维瞬态多点流场测量方法－网格图像测速法。运用数字图像处理技术由氢气泡时间－染色线流动显示图像自动生成网格图像，网格图像的网格结点为氢气泡块角点的拟合点，由网格结点的位移确定流场速度。     

非结构网格下Euler方程的高分辨率高精度解

提出了一种非结构网格下求解Ｅｕｌｅｒ方程的高分辨率高精度迎风格式。以Ｒｏｅ的矢通量差分分裂为基础,吸收了ＮＮＤ格式的优点,使其具有捕捉激波和滑移线的良好性能;在时间方向上采用Ｊａｍｅｓｏｎ的Ｒｕｎｇｅ－Ｋｕｔａ积分,并结合局部最大时间步长和残差光滑技术加速收敛。最后成功地完成了二维平板激波反射、跨音速Ｌａｖａｌ喷管内的流动和ＧＡＭＭ超音速前台阶绕流等算例,显示了该方法的有效性     

燃烧室三维计算域边界处理自动化

本文叙述了采用正交网格的有限差分法程序进行燃烧室三维数值计算时,其计算域边界处理自动化的必要性,网格划分、阶梯壁产生的约束条件,以及程序设计的一些数学方法和技巧,给出了应用本程序对某型发动机火焰筒计算域边界i,j方向平面网格划分和边界阶梯化图形,以及k向方向分布图形。