30P-30N多段翼型复杂流场数值模拟技术研究

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了30P-30N多段翼型的复杂流场,主要目的是考核湍流模型、转捩位置对多段翼型压力分布和典型站位速度型的影响。本文通过求解任意坐标系下的雷诺平均的N-S方程,采用多块对接结构网格技术,在与相应试验结果对比的基础上,详细研究了SA一方程湍流模型、SST两方程湍流模型、不同的转捩位置对该翼型压力分布和典型站位速度型的影响。本文的研究结果表明,采用全湍流模拟方式可以较好地模拟该多段翼型的压力分布,但对速度型的模拟精度较差;模拟试验的转捩位置可以改善主翼附面层与前缘缝翼边界层尾迹区的模拟精度;采用微吸气技术推迟前缘缝翼的转捩位置,可以进一步提高缝翼尾迹区的数值模拟精度。     

基于需用速度增益曲面的大气层外超远程拦截导引方法

考虑拦截器使用耗尽关机固体推进系统的情况,提出需用速度增益曲面概念,设计了基于该概念的大气层外超远程拦截导引方法。根据Lambert导引,给定拦截时间就有唯一的指令推力方向与之对应,导引过程分为两个阶段:零控拦截到达阶段,选择最优拦截时间,需用速度增益曲面迅速与助推时间-拦截时间平面相切;零控拦截保持阶段,拦截时间不断减小,保持需用速度增益曲面在助推时间-拦截时间平面上滑动,消耗多余的推进剂。利用拟零控拦截概念使两个导引阶段平缓过渡,避免了导引方法切换时推力方向的跳变。导引律计及了J2项摄动对滑行段弹道的影响,导引精度对推进系统参数偏差的鲁棒性强。仿真结果表明本文制导方法用于大气层外超远程目标拦截,脱靶量仅为km量级,推进系统参数偏差对导引精度的影响很小。     

剪切气流驱动液滴运动的受力分析

实验研究了剪切流驱动的液滴在固体表面上起始运动的受力机理。工作中使用一系列液体和固体表面来获得不同的液滴接触角,并在小型风洞中进行实验。实验中对液滴的启动气流速度进行了测量,并综合各种起始时刻的参数信息,建立了一个关于液滴接触线表面张力和剪切气流拖拽力平衡的数学模型,揭示了液滴脱落时刻的受力情况。所建立的模型更适合液滴1变形情况,但对于其它类似情况的剪切气流驱动液滴运动也能够进行合理的描述。     

静压气体轴承气膜力及其与转子耦合动力学特性研究

针对静压气体轴承,进行三维实体建模,采用供气孔区域非结构化网格和非供气孔区域结构化网格相结合的网格划分方法;采用基于有限体积法的计算流体动力学(CFD)商业软件对三维稳态可压缩气体Navier-Stokes(N-S)方程进行求解;根据计算结果,通过数据拟合获得了考虑转子偏心和转速的静压气体轴承气膜支承力模型.基于有限元法建立了气体轴承-转子系统动力学模型,采用Newmark逐步积分法求解了系统的临界转速和不平衡响应.在此基础上进行实验测试,验证了数值仿真结果.研究结果表明:低速、小偏心下,气膜主支承力随偏心呈近似线性变化;高速、大偏心下,气膜主支承力急剧增大,气体轴承的动压效应显著增强;气膜x,y向耦合支承力随转速和偏心呈非线性变化;转子系统一、二阶临界转速对当前结构刚度变化比较敏感,而三阶临界转速对此不敏感.因此,气体轴承气膜支承力的非线性特性及其与转子耦合动力学特性较为复杂,在对气体轴承进行结构设计时,应充分考虑其与转子的耦合,合理设计工作转速范围.      

跨声速压气机转子失速过程的实验

针对一台跨声速压气机转子在不同转速下的失速过程进行了稳态和动态实验研究.通过对原始信号、小波分析结果以及FFT(fast Fourier transform)变换结果的分析得到,该压气机转子在3个不同转速时均为突跃型全叶高失速,流场中存在一个失速团,失速团沿周向的传播方向与转子转动方向相同,其相对静止坐标系的量纲一转速随转子转速的增大而减小,对这一规律产生的原因进行了分析.      

气氢/气氧同轴剪切单喷嘴燃烧室壁面热流计算与分析

以气氢/气氧为推进剂,采用数值模拟方法,研究了同轴剪切喷嘴设计参数——氢氧速度比和氧压降比对单喷嘴燃烧室内燃烧过程和壁面热载的影响,并将绝热壁面条件、等温壁面条件的计算结果与试验结果进行了对比分析,结果表明:氢氧速度比增大,燃烧性能提高,壁面热载增加;氧压降比增大,燃烧性能下降,壁面热载减小;相比采用壁面绝热燃气温度,采用热流预示燃烧室壁面热载与真实情况更为接近.      

轴流风机中雷诺应力与应变滞后涡黏模型 的PIV实验

采用粒子成像测速(PIV)技术,精细测量了低速轴流风机内部叶片60%跨度处压力面与吸力面的切向和轴向速度,以及风机出口径向平面的切向和径向速度,基于叶片的周缘速度和叶片弦长的雷诺数为370000.应用周期相位平均技术分析雷诺剪切应力与平均速度剪切应变的时空演化过程,并应用互相关技术研究了两者之间的相位关系.结果表明雷诺剪切应力与平均速度剪切应变的极值发生在不同时空相位,存在相位滞后.这一现象广泛存在于三维非平衡复杂湍流中,显然,工业领域广泛应用的线性k-ε模型具有一定的缺陷,应考虑雷诺应力张量与平均速度应变张量之间的时空滞后现象,从而得到符合物理规律的计算结果.因此,雷诺应力张量与平均速度应变张量滞后的涡黏模型将成为一个很有发展前景的封闭模型,来更加精准地预测工业领域中广泛存在的非平衡复杂湍流.      

气体分配方式对多隔仓燃油箱地面惰化的影响

在各仓温度和压力为定常边界条件下,采用数值积分法,建立多仓燃油箱冲洗的数值模型,应用波音747中央翼燃油箱,通过将计算结果与国外文献中公布的实验数据进行比较,验证了模型具有较高的计算精度;然后,选择波音737中央翼燃油箱作为研究对象,采用惰化气体均匀和非均匀两种分配方式,在惰化气体流量按体积和数量两种分配方式下,计算出各隔仓氧气体积分数随时间及换气次数的变化关系,结果表明:要使各隔仓均惰化至目标氧气体积分数,按体积分配富氮气体的均匀进气方式所需的惰化时间最短,按数量分配的均匀进气方式所需的时间最长;由左舷、中弦、右弦整体平均氧气体积分数随换气次数的变化,可知非均匀进气方式优于均匀进气方式.      

高前进比旋翼气动特性分析方法研究

提出一种前进比高达0.8的旋翼气动特性分析方法,该方法针对高前进比旋翼前行桨叶压缩性、后行桨叶失速效应严重以及桨叶偏流作用和反流区大的特点,建立了旋翼气动力模型以及与之相适应的旋翼诱导速度时变非均布模型与桨叶非定常挥舞运动模型,然后根据高前进比旋翼气动力、旋翼诱导速度和桨叶挥舞运动三者之间的内在耦合关系提出了高前进比旋翼气动特性的动态响应计算方法,最后以H-34旋翼为例计算了该旋翼高前进比状态的气动特性,并将计算结果与风洞试验数据进行对比验证,结果合理。     

基于黏性涡模型的旋翼流场数值方法

 建立了一种适用于旋翼非定常流场特性分析的黏性涡数值方法。在该方法中:流场中的大尺度涡被离散为若干微小的涡元,通过求解涡量-速度形式的Navier-Stokes方程模拟涡元的输运等过程;黏性扩散效应采用高精度的粒子强度交换法进行计算,而桨叶附着涡以及新生涡环量采用了Weissinger-L升力面理论进行求解;为显著提高计算效率,在诱导速度及其梯度的计算中还引入了快速多极子算法(FMM)。应用上述方法,对悬停和前飞状态下的多个旋翼流场算例进行了计算,通过对比旋翼尾迹涡量特征和诱导速度分布等,验证了该方法的有效性。此外,还将本方法与旋翼计算流体力学(CFD)方法及传统的自由尾迹方法进行了比较,结果表明黏性涡方法在兼顾效率的同时,还能够更好地捕捉旋翼尾迹运动。