间距对凹坑强化传热和流动阻力的影响

对矩形通道内的球面凹坑壁面进行了传热系数和流动阻力测量,研究了通道Re数和凹坑间距的影响.无量纲流向间距分别为1.5,1.2和0.8.实验发现凹坑面上局部Nu数沿流向和横向变化剧烈.与光滑壁面相比,凹坑面的平均换热增强45%左右,流动阻力增大92%左右,平均综合传热性能增强16%.减小凹坑间距使流动阻力增大,平均换热和平均综合传热性能进一步增强.数据显示Re数对凹坑壁面的换热增强和流动阻力增大的影响相对较小.      

采用上游扰动体被动控制减阻技术的实验研究

通过风洞测压和水洞流动显示实验,研究了不同直径的扰动体尾迹下圆柱的流动结构和气动特性。结果表明当扰动体和圆柱串列布置时,流动呈现出两种流动模式:空穴流动和尾流撞击模式。当圆柱处于空穴流动模式时,其阻力可以大幅度降低。当扰动体直径为圆柱直径1/2时,圆柱阻力最大降低98%。综合绕流结构和圆柱表面压力分布的变化,指出阻力减小的原因主要是扰动体对圆柱前体的保护作用以及对后体流动的干扰作用,其中后者占主导地位。随着扰动体直径的增加,扰动体对圆柱前体的减阻效果增加。     

等离子体合成射流对钝头激波的控制与减阻

为了验证等离子体合成射流对超声速流动的流动控制和减阻效果,在考虑热完全气体效应的情况下,工程拟合等离子体热力学属性和输运系数,利用能量源项模型对超声速平板和球头等典型流场结构进行了数值模拟,并取得了与实验较为一致的结果。研究结果表明：对于马赫数为2的超声速平板流动,等离子体合成射流能有效干扰边界层的发展,并诱导产生一系列大尺度结构;对于马赫数为3的球头流场,等离子体合成射流能显著改变激波脱体距离与球头的阻力特性;在放电后第1个周期内,合成射流能使球头平均阻力降低6.3%,而在射流峰值情况下使阻力降低32.0%,激波脱体距离增加2倍,实现了激波控制和流动减阻的预期效果。     

LESS燃烧室非定常旋流流动

结合实验和CFD数值方法,研究了LESS燃烧室旋流器出口附近的非定常流动特征.一维热线测试结果表明:预燃级出口存在着2249Hz周期性速度振荡,而主燃级出口流动未见明显脉动.通过瞬态LES数值结果检验了在旋流器出口附近采用一维热线测量三维流动速度偏差为8%.另外,数值模拟计算得出预燃级出口速度脉动频率为2331Hz,与热线测试结果偏差为3.6%.最后通过频谱分析和相关性分析,得出周期性速度脉动由预燃级内进动涡核主导.      

层板冷却涡轮叶片前缘内部流动与传热特性实验

基于相似理论,对简化的层板冷却涡轮叶片前缘放大模型内部的流动与传热特性进行实验研究,对比了无绕流柱和带菱形扰流柱两种实验模型的流动阻力系数、靶面温度和表面传热系数的分布.实验中采用红外热像技术测量换热面的温度,采用ANSYS软件计算换热面的局部热流密度.结果表明:两种模型的流动阻力随进气雷诺数逐渐增大,带菱形扰流柱模型的流动阻力总体上较大;靶面局部表面传热系数的分布特征基本相同,带菱形扰流柱模型的局部表面传热系数比无扰流柱模型的稍高;靶面平均表面传热系数的差别很小,相同进气雷诺数下带菱形扰流柱模型的平均表面传热系数值最多大7%.      

环形流道流体阻力特性实验研究

针对环形流道下流体阻力特性较圆管阻力特性存在差异,对不同尺寸间隙的环形流道流体阻力特性进行了实验研究。依据圆管流体流动特性原理,推导出摩擦阻力计算公式,采集实验数据、计算实验结果,并对实验结果进行比较、分析。结果表明,环形通道流动特性与流道间隙大小有关,认为间隙为2.5mm是环形通道流动特性的临界点。研究发现,间隙尺寸大于2.5mm时摩擦阻力系数与理论计算相符,流态转捩雷诺数与常规尺寸基本吻合;间隙小于2.5mm为窄环形通道,窄环形流道流体阻力系数随间隙的减小而减小,流道流体流态转捩点明显提前,并且间隙值越小流态转捩雷诺数越小。     

运输机大上翘机身后体参数影响研究

采用CFD方法研究了外形差异较大的三种大上翘后体机身的阻力与流动特性,分析了不同后体形状及几何参数的影响。压差阻力虽仅占机身总阻力的15~20%,但却决定总阻力的大小。在上翘角和长细比接近的情况下,扁平度是影响流动分离和压差阻力的主要因素。给出了后体设计的参数选择及应采取措施。     

楔形肋片内冷通道传热与流动阻力特性

利用实验和数值模拟方法,分析了楔形肋片几何结构对仿螺旋肋片内冷通道强化传热与流动阻力特性的影响,并将模拟结果与实验数据进行了对比,两者结果吻合良好.研究结果表明,随着楔形肋片导流角的增大,通道的强化传热显著增强,但其流动阻力也明显增大;气流方向的不同对通道总体换热强度影响不是很大,但反吹时通道中的流动阻力却明显高于正吹;从内冷通道的综合换热效果来看,当楔形肋片导流角为7?且气流正吹时是最佳的肋片几何结构参数.      

弹舱对飞翼布局飞机气动特性影响及其控制

以高速风洞气动力测量为研究手段,开展了弹舱开启对飞翼布局飞机气动特性影响及其流动控制试验研究。试验结果表明,对于飞翼布局飞机,弹舱开启主要影响飞机阻力特性,巡航状态下,弹舱开启后使得全机阻力增加60%～110%,Ma=0.8时全机升阻比降低34%。通过在弹舱前缘安装扰流片,对弹舱腔口剪切层施加流动控制,巡航状态下弹舱开启附加阻力最多降低20%,Ma=0.8时全机升阻比提高12.6%。     

扰流柱形状对冲击冷却综合换热效率影响研究

为研究不同形状的扰流柱对冲击冷却系统的影响,对具有菱形、正方形、圆形以及椭圆形扰流柱的冲击冷却进行数值模拟研究,获得了射流Re在1×10~4～3×10~4内冲击靶板的换热特性以及通道内部流场的流动特征,分析冲击射流与不同形状的扰流柱之间的作用机理,并且为了综合考虑换热效果与流动阻力,使用综合换热效率来评价扰流柱在冲击冷却中的作用。研究结果表明:在冷却空气入口条件相同的情况下,菱形扰流柱靶板具有最高的平均努塞尔数Nu以及流动阻力,相比于正方形、圆形和椭圆形扰流柱,Nu分别提高了3.8%,8.9%,10.4%,流动阻力分别提高了9.6%,17.1%,21.3%。圆形扰流柱靶板的综合换热效率相比于菱形、正方形和椭圆形扰流柱分别提高了10.8%,4.9%,1.7%。扰流柱在冲击冷却中所起不同作用的主要原因是其与冲击射流作用后产生涡流的情况及其阻挡横流的能力。当综合考虑换热效果与流动阻力时,在所研究的4种扰流柱形状中,圆形扰流柱是最佳的选择。     

模型超音速燃烧室流场和性能的数值模拟

研究了简化模型超音速燃烧室流场和性能的数值模拟方法，该燃烧室在支板后缘设置一缝隙式喷嘴，平行于超音速空气流喷射氢气。用椭圆型偏微分方程数学模型，ＭａｃＣｏｒｍａｃｋ差分格式成功地计算了两种不同后绿尺寸的流场及其性能，为了对照比较，还计算了另外两种流场。计算结果表明，支板后的流场存在回流区，因喷氢的压力高于超音速空气流的压力和壁面附面层的影响，燃烧室内将出现斜冲波和膨胀波，压力沿横向变化明显，由此证明流场与“边界层流动”有性质上区别。提出了全流场按性质不同分段进行数值模拟的方法，利用计算得到的流场节点状态参数（如温度、压力等），积分获得燃烧室的性能参数，和文献的实验数据对比，计算结果合理、可靠。     

H2O污染对煤油燃料超燃冲压发动机燃烧室性能影响的数值模拟

在超声速燃烧设备上进行了带凹槽的直连式双模态煤油燃料超燃冲压发动机污染与纯空气来流的对比试验.为了研究来流污染对超燃冲压发动机试验性能的影响,在相同的试验气体参数(马赫数为2,总温约为828K和总压约为800kPa)下,对带凹槽煤油燃料燃烧室进行了一系列对比燃烧试验.采用AHL3D软件对3个试验条件的燃烧室反应流场进行了数值模拟.计算和试验得到的压强分布比较接近且趋势一致.通过数值模拟获得反应流场的细节,探究污染影响的机理.      

Application of Coolant Diffusion Equation in Numerical Simulation of Flow in Air-cooled Turbines

This article deals with time-averaged coolant diffusion equation in arbitrary, body-fitted coordinates. Schmidt number and Lewis number are used to obtain the mass diffusion coefficient. Other than the existing N-S solver, a new solver for the diffusion equation and the influences of cooling air on the thermal parameters, such as gas constant, constant-pressure specific heat, is established. Capable of simulating the flow of multi component fluid, this solver is used to simulate the flow in air-cooled two-stage turbines. The comparison of test results from the modified solver to those from the traditional N-S solver has evidenced that a marked difference present is in thermal parameters and flow field structure.     

飞机热气防冰系统与冰脊预测的数值模拟

基于流固耦合传热的思想建立了一套飞机热气防冰系统的的数值模拟方法,并将其与积冰热力学模型结合起来,实现了热气防冰系统开启时的机翼积冰预测.采用格心格式有限体积法求解N-S方程获得防冰腔与外流场;通过欧拉法在外流场的基础上获得过冷水滴撞击特性;求解三维热传导偏微分方程获得蒙皮的传热特性;采用交接面插值的方法实现防冰腔到外流场的热量传递;建立了考虑三维溢流效应的积冰热力学模型并在此基础上开展了机翼冰脊的数值预测.数值模拟结果表明:热气防冰系统开启时加热机翼表面温度最高可达308K,加热区后的上下机翼表面均有冰脊形成,通过对结果的分析表明该方法是合理可行的.      

燃气轮机双燃料燃烧室流场对比数值研究

针对燃气轮机环管型双燃料燃烧室分别燃用庚烷和裂解气燃料的情况,采用CFD(computa-tional fluid dynamics)技术对其燃烧流场进行对比研究.数值模拟采用了RNGk-ε湍流模型、化学平衡条件下的快速化学反应系统和简单概率密度函数(PDF)燃烧模型、液体燃料的喷雾模型以及SIMPLE算法.模拟对比分析了两种燃料燃烧下的温度分布、燃烧效率、流量分配、壁面冷却效果、空气过量系数等参数,以及它们随工况变化的趋势.所得结论如下:①不同燃料燃烧时,流场内的流量分配基本保持一致;②裂解气燃料燃烧时,其燃烧效率高出燃油状态约1%,但出口温度均匀性较差;③在加入相同化学焓值的燃料进入燃烧室时,裂解气燃料的头部空气过量系数α较大,所得到的出口平均温度低于燃油状态约20～40 K.      

高温富油燃气超燃试验研究

在空气流量１．２ ｋｇ／ｓ 左右的地面连管试验台上, 进行了模拟飞行Ｍａ＝ ４, ５, ６的三个气流总温状态的碳氢燃料（煤油） 超燃试验。试验用双燃烧室方案, 由突扩型亚燃燃烧室燃烧产生的高温可燃气以马赫数１．２５喷入超燃室, 超燃室空气流马赫数为２．１５ （或２．１３）。不同空气流总温状态下燃料当量比对亚燃燃烧室和超燃燃烧室的试验结果表明, 双燃烧室方案实施煤油的超声速燃烧是可行的。若进一步采取混合增强和合理控制油量分配等措施, 则可提高双燃燃烧室超燃效率。     

气冷涡轮叶片三维换热问题计算

建立了三维气冷涡轮叶片换热计算的热-流耦合计算模型, 将辐射热量作为源项加入到方程中, 采用热-流耦合方法完成了某气冷涡轮叶片带辐射的三维换热计算.辐射热量的计算采用了P-1辐射模型, 比较了计算辐射和不计辐射时叶片表面的温度值及温度分布的差异.计算表明对于1600 K这样的高温燃气, 辐射热量不可忽略, 考虑辐射热量后, 叶片温度明显升高, 温度分布也产生了显著差异.      

二次空气系统一维非定常计算方法

研究了燃气轮机二次空气系统的非定常计算方法.通过将整个二次空气系统划分为两部分,特征"长"元件以及局部损失元件.特征线法用于单个特征"长"元件的求解,通过将连续方程、动量方程、能量方程构成的偏微分方程组转换为常微分方程组并联立求解,可以充分考虑内部流动与换热的相互影响;基于压力修正的流体网络法则用于由多个局部损失元件组成的局部流体网络,首先对动量方程组进行求解,进而对各自的能量方程进行求解.如何在采用这两种方法进行求解的上述两类元件之间进行数据交互将被讨论.此外,所提出的时间步长定义方法以及改进的元件计算模型使得二次空气系统的非定常分析得以实现.计算程序经与某系统定常实验数据及单个元件的非定常计算数据进行对比,在10s的计算时间内,基于特征线法得到结果与CFD数据偏差达到10.2%,结果表明该计算程序可以满足二次空气系统非定常计算的要求.     

气体静压止推轴承静动态性能及振动抑制

为探讨抑制气体轴承自激振动的新手段,研究了圆盘形单个小孔节流的气体静压止推轴承振动问题。通过两平行圆盘间的气体流动方程,推导建立了气体轴承静动态特性分析模型。通过数值分析的方法,研究并总结了供气压力和气膜厚度对气体轴承静态承载力、静态气体质量流量、静动态刚度以及动态阻尼的影响,计算得到了一系列气体轴承静动态特性变化曲线。然后针对系统出现负阻尼情况,通过引入非线性能量阱(NES)来抑制自激振动的问题。研究结果表明,提高供气压力可以有效地提高气体轴承的静态性能;气体轴承的自激振动主要源于气膜挤压的负阻尼特性;同时,当NES的阻尼超过临界阻尼时,可以使系统的振动得到有效的抑制。     

Density and Temperature in Wind Tunnels Using Electron Beams

In a paper presented at the last meeting of the International Congress on Instrumentation in Aerospace Simulation Facilities, a method of measuring gas density by radiation scattering was presented. Since that time, the method has been further developed using scattering of 50-keV electrons and has been applied in a wind tunnel using air as a test gas. Test section temperature and density were in the neighborhood of 150°K and 6×10-8 g/cm3. In this application, accuracies better than 1.0 percent appear to be attainable, together with excellent spatial and temporal resolution. The application is described, and comprehensive design data are presented. In addition, the density measurements have been supplemented by measurements of the local rotational temperature, utilizing nitrogen molecular band radiation excited by a second electron beam. Comparisons between the density and temperature measurements and flow properties obtained by standard gas-dynamic calculations are made.