高超声速混合流动的粒子模拟耦合算法

通过采用基于三维非结构网格的子网格技术,对适用于高超声速混合流动的自适应时间步长粒子模拟耦合算法(Improved Hybrid Particle Simulation Method,IHPSM)进行了改进,在保证算法计算效率的同时降低数值误差。通过对三维双曲钝锥外形的数值仿真及与DSMC(Direct Simulation Monte Carlo)计算结果的对比分析,证明改进的IHPSM算法具有较高的计算精度,且能较大幅度提高计算效率。基于改进的IHPSM算法,文中针对双曲钝锥外形进行了稀薄气体效应和飞行马赫数对高超声速流动影响规律的研究。结果表明,气体稀薄程度的增加会减缓双曲钝锥前端流场宏观物理量的变化梯度,使流场激波结构变厚;来流马赫数的增大会使激波明显增强,但对激波厚度与结构的影响较小。     

压缩性修正对γ-Re_θ转捩模型的影响研究北大核心CSCD

将基于流场当地变量的γ-Reθ转捩模型加入可压缩RANS方程计算程序,对平均流场控制方程、湍流和转捩控制方程进行基于LU-SGS的隐式紧耦合求解。结合三种压缩性修正方法对高超声速平板、双楔和圆锥流动进行了数值模拟,给出了壁面斯坦顿数、热流值和实验值的比较,以及湍动能和间歇因子分布云图。数值计算结果表明,相同来流湍流度下不同压缩性修正方法得到的转捩位置差别较大。相比于未修正的模型,基于当地马赫数压缩性修正后的转捩位置最靠后,其对带有分离的双楔流动预测较为准确;基于来流马赫数的压缩性修正仅使得转捩位置稍有延迟;对湍动能源项的压缩性修正可提高模型在转捩后湍流段的热流预测精度。     

考虑多种物理效应的钝锥俯仰稳定性参数影响分析

首先采用二阶ROE格式计算了非定常NS方程下高超声速钝锥体的俯仰阻尼导数,通过与文献中内伏牛顿法和实验结果对比验证了本文所采用非定常计算方法的正确性.同时采用二阶ROE格式及双时间步长、刚性旋转动网格技术等,数值分析了NS方程计算条件下高空高马赫数钝锥体模型的强迫振动非定常流场特征,并运用最小二乘法辨识出其俯仰静稳定性导数与俯仰阻尼导数.通过计算结果对比分析了不同高度(含稀薄气体效应)、Maxwell滑移边界条件及五组元化学非平衡模型对钝锥体模型非定常流动特征及静、动导数的影响.结果分析表明,本文所计算的钝锥模型在所假设强迫振动下的俯仰稳定性参数随高度变化剧烈,且考虑滑移边界条件与化学反应会使得俯仰力矩迟滞曲线形态发生较大程度改变.准确计算高空、高马赫数条件下飞行器俯仰阻尼导数应考虑化学非平衡效应和稀薄气体效应的影响,是否考虑这些物理效应会关系到飞行器静、动稳定的程度.     

9310钢螺旋锥齿轮喷丸强化残余应力场计算仿真研究

喷丸加工诱导零件表面产生残余压应力场以提高零件疲劳寿命,是螺旋锥齿轮的关键强化工艺。为准确计算喷丸后齿面残余应力场,基于离散元与有限元耦合的方法提出一种螺旋锥齿轮喷丸工艺计算仿真模型。模型计算结果与试验结果误差在10%以内,表明模型可准确预测齿面残余应力分布。基于该模型,以AISI 9310材料螺旋锥齿轮为研究对象,探讨了喷丸工艺参数与残余应力场特征参数的关联规律。研究发现,本文所用工艺参数加工：（1）喷丸工艺主要影响轮齿表层50μm内的残余应力场;（2）喷丸覆盖率为200%时,弹丸直径和速度的改变,对表面残余压应力影响较小;（3）当弹丸速度和直径提高时,靶板表层最大残余压应力值和最大残余压应力深度都得到明显提高,其中最大残余压应力值可提高到–1251.5 MPa,最大残余压应力深度可提高到40μm。本文建立的残余应力计算模型为螺旋锥齿轮的喷丸工艺参数优选提供了计算工具与方法,把依赖试错迭代的工艺方法上升到可计算、可预测的层面。