Ni基单晶叶片粘塑性分析

从位错滑移机理出发，建立了一种适用于Ｎｉ基高温单晶结构分析用的晶体滑移粘塑性本构模型。该模型将八面体、十二面体及六面体滑移系作为潜在的开动滑移系，考虑阻应力和背应力两种内应力状态参量。利用某Ｎｉ基单晶叶片材料在７００℃和９５０℃的试验结果对模型进行了考核并标定了模型参数。进而将所建模型编入有限元结构分析软件中，对某型发动机单晶叶片进行了计算分析。     

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

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

准则方程形式的探讨

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

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

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

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

跨音速大迎角Euler方程数值分析

应用有限体积法离散三维Euler方程,选用三步显式格式进行时间推进求解;并通过当地时间步长、残值光顺、焓修正等方法加速收敛。对一维波动方程的von Neumann稳定性分析说明这种三步格式的最大Courant数是2。对三角翼跨音速大迎角流动的数值分析表明文中的三步格式及焓修正方法优于人们经常采用的四步、五步格式以及焓阻尼技术。     

高阶间断有限元法的并行计算研究

根据间断有限元法的数据结构特点,基于METIS网格分区技术,设计并行计算策略,在非结构网格上实现了并行高阶间断有限元法。控制方程的数值通量项使用Local Lax-Friedrichs(LLF)格式计算。设计了并行的牛顿-块高斯赛德尔法(Newton-Block GS)来加速收敛,提高迭代效率。并行性能分析表明,所设计的并行算法能够得到较好的加速比和并行效率,有效地节省计算时间,合理分配内存。这使得采用高阶间断有限元法计算更为复杂的问题成为可能。