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

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

尖端具有非线性吸附接触的界面裂纹的动态特性研究

 研究一类尖端具有非线性吸附接触型裂纹模型的弹性波散射问题。利用积分变换和积分方程方法推导了确定这类问题基本变量的奇异积分方程组,采用围道积分技术和切比雪夫多项式展开技术得到了待定系数的非线性代数方程组。最后给出了裂纹尖端接触区的大小和接触应力的数值结果,揭示了这种接触型裂纹模型的动力学特征及物理上的合理性。     

Study on the breakup lengths of free round liquid jets

IntroductionBreakup and atomization of roundliquidjetshad been widely used in internal combustion en-gines,gas turbines,and liquid rocket engines.Recently,with the development of scramjet en-gines and the deepening of fuel atomization re-search,it attract…     

求解非线性方程及方程组的模拟退火算法

探讨了用模拟退火算法求解非线性方程及方程组的方法,该方法把方程组求解问题转化为函数优化问题.计算中不需要使用目标函数的导数信息,在Matlab环境下实现了该算法,数值实验结果表明了该算法的有效性.     

叶轮机械三维不可压Euler方程的直接解法

本文用近似因式分解交替方向隐式格式直接求解叶轮机械中原参数三维不可压Euler方程组。在连续方程中引入“拟压缩性”之后,将原来类型不确定的控制方程变成双曲型方程,给定初边值之后用时间推进法求得稳定解,对一个单级压气机转子内部流场进行了计算,结果与实验值比较符合。     

边条翼布局双垂尾抖振的数值模拟

对边条翼布局的双垂尾抖振进行了较为深入的数值模拟研究。模拟来流马赫数为0.2,迎角为10°~40°。通过非定常Euler方程计算各迎角下的非定常流场及垂尾根部弯矩系数。并将随时间脉动的根部弯矩系数进行计算得到根部弯矩系数均方根值,从而得出根部弯矩响应大小随迎角的变化曲线。结合流场特性对该布局双垂尾抖振的发生机理及抖振响应随迎角的变化规律作了深入分析。结果表明:该边条翼布局双垂尾抖振主要是由于边条破裂涡作用在垂尾上的脉动载荷引起的。最后,将垂尾根部弯矩响应的计算结果与该模型双垂尾抖振的风洞实验结果作了比较,结果符合得较好。     

用于棘轮变形预测的棘轮演化统一模型研究

 在大量单轴棘轮实验基础上,研究了均值、幅值、峰值和谷值应力对 304不锈钢的饱和棘轮应变的影响规律,首次提出了棘轮应力σ_r和棘轮门槛值σ_rth 的概念,建立了基于单参数控制的、用于饱和棘轮应变预测的 SRM饱和棘轮本构模型和用于独立循环应力工况下棘轮应变演化预测的 REM棘轮演化模型,并由此发展了全面描述任意循环应力工况下棘轮应变演化规律的 URM棘轮演化统一模型。抛物律模型 SRM和幂律模型 REM对实验数据拟合的相关系数均超过 0.98。URM建模容易,只需 4～ 6个试样的单轴棘轮实验数据。此外还讨论了获得 SRM本构模型的单试样实验法。     

Ni基单晶叶片粘塑性分析

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

高性能涡轮发动机润滑系统设计

简述了润滑剂对发动机设计的影响。阐明了使用酯类润滑油的设计问题及润滑系统通用的设计要求。     

高雷诺数时串列双圆柱脉动压力的实验研究

本文通过风洞实验给出了二维串列双圆柱在高雷诺数时(分别为Re=3.25×10~5和Re=6.5×10~5)的脉动压力分布。实验在低湍流度的均匀流中进行,模型表面光滑。实验结果表明:在这二种分别属于高亚临界和低超临界雷诺数下,串列双圆柱不同间距比时的脉动压力分布差别很大。Re=3.25×10~5时前后圆柱间的相互干扰对圆柱表面脉动压力分布的影响要比Re=6.5×10~5时大得多,而且二者干扰的规律和特性也截然不同。