基于遗传算法的分段多参气动阻力研究

针对BADA模型未考虑马赫数、仅基于升力系数求得的阻力系数与真实值存在较大误差的问题,考虑到飞机以较大马赫数飞行时气流压缩效应会改变阻力特性,提出了一种基于遗传算法的临界马赫数分段多参数数据拟合方法,建立了更为精确的阻力系数数学模型。以B737-800民机为算例的计算结果表明:该模型计算的阻力系数与真实值的平均相对误差小于5%,与经典BADA模型相比,平均相对误差下降了15.63%,均方根误差下降了0.384%,拟合优度提升了0.344;所提优化方法提高了阻力计算的精确度,可提高ATM飞行仿真的准确性。     

气液两相绕流钝体稳燃器的冷态数值研究

提出并建立了基于可变时间间隔平均方法的多尺度两相湍流模型,通过算例验证,证实了所推导的气液两相多尺度数学模型、所标定的模化参数及描述气液两相作用机理的合理性,及其计算的时间平均阻力系数、脉动升力系数误差、斯特劳哈尔数和回流区大小误差分别为1.45%,0.323%,2.17%,2.33%,计算结果明显优于标准k-ε和重整化群(RNG)k-ε模型.用多尺度湍流模型对气液两相绕流6个不同的钝体稳燃器进行数值计算表明:船形与锥形是综合性能较优的两个钝体稳燃器结构.其中,船形钝体比锥形钝体稳燃器回流区大10.53%,其时间平均阻力系数比锥形稳燃器大4.776%,方均根脉动升力系数比锥形稳燃器小44.73%,通过全方面综合比较,船形是综合性能最优的钝体稳燃器结构.     

人椅系统绕流流场的数值模拟

以某种成熟型号的弹射座椅为研究对象,对其进行合理简化、改进,获得了接近实体的人椅系统几何模型.借助CFD软件FLUENT,采用控制容积法,对N-S方程、壁面函数、标准湍流方程进行了求解.在验证计算数据与实验数据相吻合的基础上,考察了马赫数为0.2,0.6,1.2,攻角为0°60°,偏转角为0°30°条件下可压缩人椅系统绕流流场情况,得到了人椅系统绕流流场的流速、压力、马赫数、表面压力及气动力参数分布情况和变化规律.结果表明:高马赫数、零攻角、大偏转角对人体造成的损害程度最大.      

折叠V形火焰稳定器原型的冷态与稳燃特性研究

为了获得折叠V形钝体这一新型火焰稳定器原型的冷、热态特性,采用风洞实验与混合雷诺平均/大涡模拟结合的方法对10组不同折叠角的折叠V形钝体火焰稳定器模型进行了冷态实验与热态数值模拟研究。其中冷态实验风速10~50m/s,来流雷诺数10000～80000,热态数值模拟工况为进气温度700K、来流速度50~150m/s,来流雷诺数20000～60000,当量比0.2~1。通过实验测得了不同折叠角折叠V形钝体总压恢复系数,通过数值模拟进一步获得了阻力系数、吹熄性能与燃烧效率,并获得了折叠V形钝体下游火焰时均轮廓。研究表明,折叠V形钝体总压恢复系数与阻力系数优于标准钝体,总压恢复系数随钝体折叠角减小而增大,阻力系数随折叠角减小而减小;折叠V形钝体抗吹熄能力与燃烧效率均优于标准V形钝体;折叠V形钝体的下游火焰扩散特性与标准V形钝体有显著差异。上述特性揭示了折叠V形钝体具备改善加力燃烧室性能的较大潜力,在实际应用中,还需对稳定器的空间排布做进一步研究与优化。     

室火轰燃的大涡模拟

采用一套修正的适用于低马赫数浮力流的N-S方程组来描述火灾场,运用Smagorinsky亚格子模型的大涡模拟方法对室内池火的轰燃发展过程进行了数值模拟,把模拟与实验结果进行了对比和分析。结果表明,该方法能较好地预测轰燃。同时模拟结果显示火焰区向房间后部偏转稍大,近壁处温度预测的误差大一些,随着亚格子模型的不断完善和计算能力的提高,用该技术来预测轰燃现象会更加准确。     

激光干涉位移测量系统的环境误差补偿

激光干涉仪在长度测量领域应用十分广泛,然而其测量精度直接受环境的影响,特别是空气折射率变化会引入标尺误差。通过对激光干涉位移测量原理及其误差源的研究分析,建立了环境误差补偿模型,消除了零点误差和空气折射率误差。设计了分别适用于环境条件好、小量程测量的双干涉自补偿系统和适用于复杂环境、大量程测量的多点测量综合修正的补偿系统,对环境误差进行修正。     

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

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

大厚度钝后缘翼型气动性能计算研究

大厚度钝后缘翼型由于其气动和结构的优点,近来被作为大型风力机叶片设计时的内侧翼型成为风力机翼型设计的热点之一EULL由于此类翼型厚度和风洞阻塞度的限制,大雷诺数的风洞实验数据很少。耦合RANS方程和基于线性稳定性分析EULL捩预测方法,进行了大厚度钝后缘翼型气动性能计算研究。使用Gridgen程序生成计算网格,湍流模型为S-A模型。对常用的钝后缘修形方式如直接截断、对称增加厚度等并且有风洞实验结果的几种翼型进行了CFD计算并和实验结果及文献计算结果进行了比较。结果表明升力系数在线性段,计算结果和实验结果吻合很好,但是对失速迎角的捕捉能力差;阻力系数计算结果和实验值吻合较好。不同修形方式的计算结果为钝后缘修形设计提供了参考。     

一个新的可压缩性修正的k-ε模型

考虑结构可压缩性修正的影响,发展了一个同时考虑结构可压缩性修正和膨胀可压缩性修正的k-ε湍流模型,新模型包括Chang可实现性、Heinz湍流动能产生项以及Sarkar可压缩性三部分修正.新模型扩宽了以往发展的可压缩性修正模型的适用范围,适用于高超声速(M>5)复杂湍流流动中.通过对多个复杂超声速横侧射流工况的计算,验证了新模型的预测效果.与实验结果相比表明,几个工况下新模型的预测精度都显著高于标准k-ε模型.流体分离强度越大,新模型的修正效果越显著.与标准k-ε模型相比,新模型计算结果与实验更加接近.     

压缩性修正对γ-Re_θ转捩模型的影响研究

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

Sensitivity analysis of flowfield modeling parameters upon the flow structure and aerodynamics of an opposing jet over a hypersonic blunt body

Implementation of an opposing jet in design of a hypersonic blunt body significantly modifies the external flowfield and yields a considerable reduction in the aerodynamic drag. This study aims to investigate the effects of flowfield modeling parameters of injection and freestream on the flow structure and aerodynamics of a blunt body with an opposing jet in hypersonic flow. Reynolds-Averaged Navier-Stokes (RANS) equations with a Shear Stress Transport (SST) turbulence model are employed to simulate the intricate jet flow interaction. Through utilizing a Non-Intrusive Polynomial Chaos (NIPC) method to construct surrogates, a functional relation is established between input modeling parameters and output flowfield and aerodynamic quantities in concern. Sobol indices in sensitivity analysis are introduced to represent the relative contribution of each parameter. It is found that variations in modeling parameters produce large variations in the flow structure and aerodynamics. The jet-to-freestream total-pressure ratio, jet Mach number, and freestream Mach number are the major contributors to variation in surface pressure, demonstrating an evident location-dependent behavior. The penetration length of injection, reattachment angle of the shear layer, and aerodynamic drag are also most sensitive to the three crucial parameters above. In comparison, the contributions of freestream temperature, freestream density, and jet total temperature are nearly negligible.     

钝缘舵高超音速湍流分离特性

给出由半圆柱前缘舵诱导的高超音速湍流分离的实验结果。实验气流Ｍａｃｈ数为７．８,单位长度Ｒｅ数为３．５×１０７ｍ－１。结果表明：钝缘舵诱导的湍流分离极不稳定,分离激波出现大尺度低频振荡,使壁面压力和热流率无量纲标准偏差在主分离线附近达最大值。Ｍａｃｈ数愈高,最大无量纲标准偏差值越大。在前缘区前缘直径是控制分离流场尺度和平均壁面压力、热流率分布的主要参数     

趋于临界马赫数的圆柱跨声速绕流特性分析

通过深化认识趋于临界马赫数Ma_cr的圆柱跨声速绕流特性,明确新型飞行器增升减阻设计的空气动力学理论依据。采用大涡模拟方法数值研究了来流马赫数Ma_∞为0.75和0.85、雷诺数Re为2×10⁵的圆柱跨声速绕流。结果表明:当Ma_∞趋于临界马赫数(Ma_cr≈0.9)时,圆柱的阻力下降且升力系数振荡被抑制;通过力的分解,得知圆柱的阻力减小来自旋涡力的影响,而非可压缩性;圆柱的阻力减小与其背压上升有关;剪切层初始阶段的对流马赫数Ma_c随Ma_∞的增加而增大,而增长率相反,这使得剪切层更为稳定、柱体背压更高。此外,由于Ma_∞=0.85时边界层分离点处的激波和尾迹处的激波向下游推移,使得近尾迹处的湍流脉动减弱,进而导致柱体的表面压力振荡和升力系数振荡被抑制。     

适用变比热的广义气体动力学基本函数

为了计算高马赫数条件下气体动力学参数，引入无量纲的温度平均比热系数m和对数温度平均比热系数n，从气体动力学和热力学基本概念出发，推导建立了能够反映“变比热效应”的广义气体动力学函数和理论动压表达式。广义气体动力学函数在马赫数小于3条件下的简化形式与传统气体动力学函数一致，在马赫数大于3条件下能够有效修正“变比热效应”所导致的传统气体动力学函数计算偏差，通过无量纲系数m、n可以便捷得到总静参数的关系式。     

Detached-eddy simulation of supersonic flow past a spike-tipped blunt nose

Space vehicle in atmosphere travels mostly at supersonic speed and generates a very strong bow shockwave around its blunt nose. Oblique shock and conical separated flow zone generated by a forward disk-tip spike significantly reduce the drag by reducing the high pressure area on the blunt nose. This study employs improved delayed detached eddy simulation to investigate the characteristic flow structures around a spike-tipped blunt nose at Mach number of 3 and Reynolds number(based on the blunt-body diameter) of 2.72x10~8. The calculated time-averaged quantities agree well with experimental data. Characteristic frequencies in different flow regions are extracted using fast Fourier transform. It is found that two distinct instability modes exist: oscillation mode and pulsation mode. The former is related to the foreshock/turbulence interaction with nondimensional frequency at around 0.004. The latter corresponds to the interaction between turbulence and shock structures around the blunt nose, with a typical coherent structure shedding frequency at 0.092.     

Interaction of single-pulse laser energy with bow shock in hypersonic flow

Pressure sensing and schlieren imaging with high resolution and sensitivity are applied to the study of the interaction of single-pulse laser energy with bow shock at Mach 5. An Nd:YAG laser operated at 1.06 lm, 100 mJ pulse energy is used to break down the hypersonic flow in a shock tunnel. Three-dimensional Navier–Stokes equations are solved with an upwind scheme to simulate the interaction. The pressure at the stagnation point on the blunt body is measured and calculated to examine the pressure variation during the interaction. Schlieren imaging is used in conjunction with the calculated density gradients to examine the process of the interaction. The results show that the experimental pressure at the stagnation point on the blunt body and schlieren imaging fit well with the simulation. The pressure at the stagnation point on the blunt body will increase when the transmission shock approaches the blunt body and decrease with the formation of the rarefied wave. Bow shock is deformed during the interaction. Quasi-stationary waves are formed by high rate laser energy deposition to control the bow shock. The pressure and temperature at the stagnation point on the blunt body and the wave drag are reduced to 50%, 75% and 81% respectively according to the simulation. Schlieren imaging has provided important information for the investigation of the mechanism of the interaction.     

仿生正弦前缘对翼面动态失速的影响

动态失速导致叶片气动载荷急剧变化，造成振动载荷激增，桨叶寿命大幅衰减。针对动态失速问题，从座头鲸胸鳍在动态倾转下取得良好的流动特性获得启示，据此模化出仿生正弦前缘翼面（包含3种波峰和2种波长），旨在实现动态失速控制。借助三维非定常数值模拟方法，采用运动网格技术，基于SC1095旋翼翼型，研究了仿生前缘动态失速流动控制机理及运动参数和来流速度的影响。结果表明：正弦前缘大幅度降低俯仰力矩系数峰值和阻力系数峰值；前缘波峰越大、波长越小，阻力系数峰值与俯仰力矩系数峰值的抑制效果越明显，虽然升力系数峰值减小，但其减小量远小于前两者，例如其中一种仿生翼使俯仰力矩系数峰值减小了47.7%，阻力系数峰值减小了36.4%，升力系数峰值减小14.1%；在最大迎角附近，正弦前缘能够缓和失速特性，使载荷变化更为平缓；在高平均迎角、低俯仰频率、低马赫数下，仿生翼动态失速控制效果更强，相比较而言迎角振幅的影响较小。     

二元收-扩喷管与球形收敛调节片喷管的性能对比研究

为了对比二元收-扩喷管(2DCD)与球形收敛调节片喷管(SCFN)的性能,基于CFD数值计算软件对2种喷管在相同工作点下进行了内外流场计算与分析。数值模拟结果表明:在相同高度及飞行马赫数条件下,SCFN的流量系数相较于2DCD呈现出较大不同,推力系数则略差于2DCD;在低落压比条件下,2种喷管内流对后体阻力特性影响不大;在高落压比条件下,2种喷管内流对后体阻力特性影响较大。     

基于边界层转捩的高超声速进气道特性研究

为了探索边界层非强迫转捩对进气道性能的影响,采用数值计算的方法开展了边界层转捩对轴对称混压式高超声速进气道流场特性的研究。研究表明:随着进气道中心锥锥尖钝化半径增大,边界层转捩先推迟。当锥尖钝度大到一定程度时,边界层转捩位置前移。随着钝化半径进一步增大,边界层转捩再次推迟,转捩位置逐渐后移。来流湍流度越大,边界层越不稳定,边界层转捩越易发生。与湍流边界层相比,考虑边界层转捩时进气道的总压恢复系数及流量系数较高、热载荷及阻力系数较小,Ma=6.5时喉道处总压恢复系数最高上升17.3%,进气道阻力最大下降17.4%。边界层转捩对壁面热流密度分布影响较大,但对壁面压力分布影响较小。钝化影响进气道的自起动性能,随着钝化半径增大,自起动马赫数升高,而边界层转捩对进气道自起动性能影响较小。     

结冰对翼型流场影响的研究

本文以NACA23012翼型为研究对象,采用CFD软件模拟了结冰对翼型流场的影响。结果表明：结冰对翼型的气动性能有着明显的影响,尤其是在大攻角时,升力系数显著下降,而阻力系数则明显增加;对于光滑翼型和结冰翼型．计算值与试验值在‘定攻角范围内吻合良好：马赫数对两种翼型的升力系数和阻力系数随攻角的变化影响均比较小,雷诺数也一样。