二维任意域内基于节点的局部网格生成算法

凸域内基于节点的局部网格生成算法,克服了基于节点的有限元方法的网格生成可能产生的不一致性。将该基于节点的局部网格生成算法的适用范围拓展到二维任意域。另外,提出了通过使用约束Delaunay路径来划分任意域的区域划分算法,该算法使得在并行实现网格生成的过程中各处理器之间无需通信,从而大大提高了节点给定情形下有限元方法网格生成的并行效率。     

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

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

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

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

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

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

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

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

准则方程形式的探讨

应用积分方程求解换热问题,不仅数学运算较简单而且计算结果具有很好的直观性,容易看出理论解的物理意义,其计算结果使人们一目了然。下面我们利用积分方程,在第一类边界条件求解Pr=0.001～7500的流体绕流平壁层流换热问题。      

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

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

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

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

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

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

二维翼型结冰数值计算

翼型结冰计算是随着结冰时间步长变化的准定常过程,因此在每个时间步长后都要更新各个模块的计算工作,如此反复直到达总结冰计算时间.在各模块计算中,分别采用了泊松方程法生成网格、求解雷诺时均方程计算流场、拉格朗日追踪法求解水滴轨迹、基于Messinger模型计算冰形、采用变距滤点法重构边界.引入动态链接库技术将各模块集成为结冰计算平台,并将三种冰形(楔形冰、槽状冰及混合冰)计算结果与国外软件及试验数据进行了对比.