跨音速叶栅粘流计算的多级广义有限元法

将 Taylor-Galerkin广义有限元法和多级有限元的思想结合起来 ,并在人工粘性的处理上作了改进 ,构成了收敛速度和稳定性均较好的多级广义有限元法。利用该方法 ,并基于 N-S方程研究了二维跨音速叶栅流动。计算结果与实验结果符合较好。通过计算表明 ,该方法计算稳定、收敛较快 ,是叶轮机械内部跨音速流动强有力的计算手段     

矢通量分裂隐式有限体积法解超音速喷流强干扰流场

发展了用于轴对称流场,任意非正交曲线网格下的矢通量分裂隐式有限体积法;统一了二维和轴对称流场计算的表达式。该方法对粘性通量项进行了分解和归类,考虑了包括交叉导数项在内的所有项对隐式增量的贡献。对带中心喷流的流动计算表明该方法收敛速度和计算精度均优于ＭａｃＣｏｒｍａｃｋ显式格式。横向喷流强干扰流场的计算初步揭示了姿控发动机喷流产生间接推力的机理,由此,设计中可降低对主发动机额定推力的要求     

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

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

双轴应力下含裂纹加筋曲板剩余强度的弹塑性有限元计算

 采用考虑了大应变的弹塑性有限元方法,以裂纹尖端张开角为断裂准则,实现了双向受载条件下加筋曲板中裂纹起裂到失稳整个扩展过程的数值模拟计算;还对与加筋曲板有同样尺寸的加筋平板进行了分析和计算。结果表明,所提出的分析方法和所采用的断裂准则适用于对含裂纹加筋结构所进行的剩余强度分析。     

用三维欧拉方程计算机身与机翼绕流

本文用Jameson的三维欧拉方程有限体积法、四步Runge-Kutta时间推进格式,计算机身和大后掠细长机翼的三维可压缩绕流。对钝头机身,用C-O型网格;尖头机身与大后掠细长机翼用H-O型网格。本文介绍钝头旋成体、带座舱前机身和三角翼绕流的计算结果,显示流场等值M线分布,计算压强系数分布与实验比较,以及三角翼大迎角分离涡等。由于用了隐式残值光顺等加速收敛措施,有效地减少推进步数,节省机时。     

固体发动机喷管内两相流的有限元数值模拟

本文给出了固体发动机喷管内两相流的数值模似。从非守恒型两相流甚本方程出发,采用加权积分技术,导出了离散点的流动物理量的有限元数值计算公式。数值模拟分别与JPL喷管的两相流壁面静压分布测量结果以及与差分法计算结果作了比较,比较结果表明:二者是吻合的。本文还以中型喷管为例,计算了粒子的轨迹,粒子的速度滞后和燃气的温度滞后,计算结果表明;入口粒子冲击角对粒子轨迹有明显的影响。     

矩形风洞收缩段流场的计算和分析

本文结合西工大低速风洞设计进行了矩形三元收缩段流场的计算和分析,计算采用差分方法和贴体坐标,以AF1格式进行离散化并利用ADI方法求解。对三种常用的收缩曲线(即Witozinsky曲线,五次方曲线和双三次方曲线)的三元流场进行了具体计算和比较。计算结果正确反映了三元收缩段流动的特点,给出了流动三元性对壁面逆压梯度和出口均匀性影响的数值结果,表明本文方法是三元收缩段设计的有效分析和设计手段。     

有攻角细长旋成体极端地面效应分析

本文用匹配摄动方法研究了细长旋成体贴近地面运动时诱导的不可压缩势流问题,它是经典细长体理论的推广和发展。文中将流场分为内外两个渐近区域,建立了问题的匹配摄动理论,并利用文献[1]中的双圆诱导不可压缩势流的精确解,给出了该问题的首阶解析解。然后又导出了细长体上的吸力和俯仰力矩的计算公式。在此基础上,文中对细长旋成体极端地面效应的运动学规律和动力学规律作了定性和定量两个方面的分析。     

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

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

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

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