不同热值生物燃料燃烧特性数值模拟

采用标准k-ε湍流模型、涡耗散湍流燃烧模型、P-1辐射传热模型对二维简化燃烧室燃料甲烷气、沼气、生物质气的燃烧特性进行了数值模拟.通过比较不同截面的燃烧情况,分析燃烧室内的温度分布、NOx质量分数分布、流场结构,研究了不同热值燃料的燃烧特性.结果表明:在入口空气、燃料质量流量相同的条件下,随着燃料热值降低,燃烧温度降低,温度分布更均匀,NOx排放量减少,流动速度降低;适当调节余气系数,能改善燃烧效果,温度分布仍满足高热值燃料燃烧温度高于低热值燃料燃烧温度的规律;在生物燃料混合气中掺入不同热值的可燃成分,可改变混合气的热值,恰当的混合比例能更好地发挥生物燃料的作用.     

弯曲网格上的间断有限元湍流数值解法研究

采用间断有限元方法对雷诺平均Navier-Stokes(RANS)方程进行了数值求解,对Spalart-Allmaras单方程湍流模型进行了部分修正,使得求解器更加鲁棒。构造了分段高次多项式来逼近真实物面,同时物面附近采用多层弯曲网格来避免网格交叉,此外提出了一种快速计算积分点的曲面物面距的方法。采用混合网格对NACA0012翼型以及RAE翼型进行了数值计算,并与实验数据以及前人数据进行了对比。计算结果表明,采用物面弯曲网格结合修正的湍流模型方法在相对稀疏的网格上就能得到比较好的数值解。     

连续尺度变换方程在稳态及非稳态流动模拟中的应用

为了提高涡粘性假设的湍流模型对于非稳态流动的求解精度,同时兼顾其对于稳态流动的求解性能,将雷诺应力项与连续变换方程(CSSE)结合而形成新的应力项,使其根据流场尺度、网格尺度及Kolmogorov尺度来自动调节当地的应力雷诺应力模化水平,避免网格因素在流场模拟中产生不利影响,改正了混合RANS/LES方法的速度型偏离对数率问题;同时,该方法并未引入显式亚格子模型(SGS),因此回避了亚格子系数确定对于流场模拟精度产生的影响,改善了湍流模型对于流动不稳定性的辨识精度。在湍流平板算例中,CSSE方法计算的边界层速度型精度与雷诺平均方法(RANS)相当,而对于圆柱尾迹的模拟则证明了CSSE方法具有混合RANS/LES方法的优点,即能够准确模拟流动的不稳定性特征。     

结构误差对旋转稳定弹丸气动特性影响的数值模拟

基于剪切应力传输(SST)k-ω湍流模型,对考虑结构误差情况下的弹箭模型进行数值模拟计算。首先采用CFD和工程经验公式相结合的方法,得到了不同马赫数、不同转速情况下无结构误差模型的多种气动特性参数。通过与实验数据对比,阻力、法向力、俯仰力矩、压心位置误差在10%以内,验证了该方法的可行性和准确性。然后,建立了考虑结构误差即质量分布不对称、弹体不同轴和无结构误差模型,并分别进行了气动特性模拟计算。结果表明,质量偏心对滚转阻尼力矩系数、马格努斯力矩系数影响很大;弹体不同轴对法向力系数导数、俯仰力矩系数导数、马格努斯力矩系数导数和压心位置都产生很大影响。结构误差使模型的多种气动特性参数产生很大变化,将影响其飞行弹道和稳定性。     

采用钝体式孔板淬熄的富油-淬熄-贫油驻涡燃烧室排放性能试验研究

设计研制了采用钝体式孔板淬熄装置进行淬熄作用的富油-淬熄-贫油(Rich-quenchlean,RQL)驻涡燃烧室(Trapped-vortex Combustor,TVC)模型和试验系统,实现了富油燃烧-快速淬熄-贫油燃烧的分级燃烧。在常压状态下采用RP3航空煤油作为燃料开展了排放试验研究,分析并总结了当量比(包含总当量比和富油区当量比)、进口速度(马赫数)和进口温度等参数对采用钝体式孔板淬熄的RQL-TVC排放性能的影响。研究结果表明:不同进口温度或进口马赫数条件下,随着当量比的增大,氮氧化物(NOx)的排放指数(Emission Index,EI)呈现出先急剧升高,然后急剧下降并趋于平缓的趋势,在富油燃烧区当量比为1.1左右达到峰值,一氧化碳(CO)的EI值先保持平稳而后急剧下降,未燃碳氢化合物(UHC)的EI值呈现连续降低的趋势,燃烧效率持续增大,燃烧效率均高于95%,最大值达99%以上;进口温度较高或进口马赫数较小时,UHC和CO的EI值均较小,燃烧效率较高,而NOx的EI值较大,反之亦然。     

气氧/煤油富燃燃气发生器燃烧试验研究

为了研究航天煤油的富燃燃烧特性,设计了一个气氧/煤油富燃燃气发生器,并进行了室压范围1~4MPa、混合比范围0.25~0.52的燃烧试验,获得了混合比和室压对特征速度、燃气温度及燃烧效率的影响规律。试验结果表明,采用火花塞点火方式可以实现该燃气发生器在高度富燃条件下的可靠点火;在本文研究范围内,试验测得的特征速度和燃气温度明显低于化学平衡理论预测值,燃烧效率介于63.9%~88.3%之间。混合比是特征速度、燃气温度及燃烧效率的最大影响因素,室压的影响远小于混合比。燃烧效率与实际燃气温度与理论燃气温度比值的平方根成正比。     

转捩模型在翼型气动性能预测 中的应用研究

边界层转捩位置的准确预测对于提高飞行器气动性能的预测精度具有重要意义.选取与k-ω SST湍流模型相耦合的γ-Reθt模型,以零压力梯度平板为研究对象,通过求解基于有限体积法的雷诺平均N-S方程验证该模型自动捕捉流动转捩的准确性;将该模型应用于传统有压力梯度的NACA 0012翼型的流场特性和气动性能的研究中,并与原始k-ω SST模型的计算结果及全湍流试验数据进行比较.结果表明:远场边界距离对翼型阻力系数有较大的影响;与无转捩模型相比,γ-Reθt转捩模型对翼型阻力系数的预测精度有一定程度的提高;对于二维模型,γ-Reθt转捩模型难以捕捉翼型表面的三维效应和非定常分离特性.     

带迎角钨杆斜侵彻铝板数值仿真研究

为研究迎角对钨杆斜侵过程的影响,开展了带迎角钨杆斜侵彻铝板数值仿真研究。仿真中采用SPH方法,Shock状态方程和Steinberg本构模型。钨杆共有两种尺寸,分别为Φ10mm×23mm和Φ10mm×46mm。铝板厚10mm,撞击速度为2km/s,撞击角度为60°。为简化分析,将迎角分解为俯仰角和偏航角,其中俯仰角范围为-90°～90°,偏航角范围为0°～90°。对钨杆侵彻过程的分析结果表明：斜侵彻存在临界迎角,在临界迎角范围内钨杆所受力矩可忽略不计;在临界迎角范围外,钨杆将发生弯曲甚至折断。对钨杆的剩余质量分析结果表明：当偏航角为0°且钨杆未折断时,剩余质量只在较小的范围内变化;当俯仰角为0°时,剩余质量随迎角增加线性减小。对钨杆的剩余速度分析结果表明：偏航角对剩余速度的影响与俯仰角对剩余速度的影响基本一致,剩余速度曲线关于0°迎角基本对称。根据仿真结果,采用曲线拟合的方法给出了钨杆斜侵彻铝板的剩余速度经验公式。     

有阀式脉动喷气发动机出口瞬态速度场测量

通过使用TR-PIV(time-resolved-patricle image velocimetry)系统对某亥姆霍兹型有阀式脉动喷气发动机出口瞬态速度场进行测量.实验结果表明:有阀式脉动喷气发动机出口压力及速度场呈明显的周期性变化过程,时均速度最高的区域并不在尾喷管出口截面,排气阶段的旋度及持续时间均大于吸气阶段.有阀式脉动喷气发动机出口未完全膨胀的高温燃气在排出尾喷管后仍会继续膨胀,对于提高性能不利.      

可压缩粘性流动笛卡尔网格虚拟单元方法研究

以二维高雷诺数可压缩粘性流动问题为背景,提出了一种全新的笛卡尔网格虚拟单元方法。基于壁面函数基本假设,构造了壁面函数-虚拟单元方法(WF-GCM),用于定义湍流壁面边界条件。引入参考点的概念计算虚拟单元上的基本变量与湍流变量值,定义了"非贴体"笛卡尔网格下的湍流壁面边界条件,并通过壁面函数模型修正近壁面单元与界面单元。基于自适应笛卡尔网格体系,采用发展的具有二阶精度的格心格式有限体积求解器,数值模拟了跨音速RAE2822翼型绕流问题与超音速圆柱绕流问题,计算结果与实验值吻合良好,显示了WF-GCM对高雷诺数可压缩粘性问题是有效的。