基于非结构动网格的非定常气动力计算

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法耦合求解基于Arbitrary Lagrangian-Eulerian(ALE)描述的欧拉方程,计算了谐和振动NACA0012翼型及矩形机翼的跨声速非定常气动力,计算结果与参考文献提供的结果及实验结果吻合良好。     

多重网格技术在复杂粘性流场计算中的应用及研究

流体力学控制方程的数值求解过程中,当网格加密、粘性效应强或来流马赫数低的时候,流场收敛非常缓慢。本文在多块网格、LUSGS隐式格式基础上采用多重网格技术,计算了NLR7301case#1两段翼型的低速绕流及DLRF4跨声速流场,在复杂二维、三维粘性流场求解中获得了满意的加速收敛效果。文中还比较了V循环、W循环的加速效率,为多重网格的工程应用奠下基础。     

高速磁浮列车会车压力波数值计算研究

根据可压缩粘性流体的N-S方程和k-ε双方程湍流模型,运用移动网格技术,采用有限容积法对磁浮列车以500 km/h速度运行时的会车压力波进行了数值计算研究,计算结果与现车试验结果相差10%以内,表明了本文计算模型和方法的正确性;文中还研究了交会过程中列车周围的流场分布,结果表明:当列车交会时,两车之间的流场相互叠加,而两车之外的流场分布在整个交会过程中变化不大。     

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

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

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

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

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

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

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

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

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

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

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

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

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

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