甲烷/空气预混火焰回火特性数值模拟

为了获得有效的预测回火极限的数值模拟方法,对甲烷/空气预混火焰进行了数值模拟。采用不同的反应机理和燃烧模型,针对一个二维轴对称的本生灯,分别对甲烷/空气预混层流火焰和湍流火焰进行了数值模拟,预测了甲烷/空气预混火焰的回火特性,并与Johnson的实验结果进行了对比。对层流火焰的模拟结果表明,采用Smooke-46反应机理能够较为准确地预测层流回火极限。而单步化学反应所预测的层流回火极限要比实验值低得多,这说明单步化学反应下的火焰更难以发生回火。对湍流火焰的模拟结果表明,涡耗散概念模型比有限速率模型更准确地预测了湍流回火极限。有限速率模型所预测的回火极限略高于实验值。     

碳／碳材料烧蚀对电离边界层的影响

本文给出了碳/碳复合材料热化学烧蚀对层流和湍流电离边界层影响的计算方法。边界层中考虑了20个气体组元,并假设所有化学反应平衡。应用了边界层方程的新解法和简化的输运性质处理。算例表明,烧蚀对边界层的电离特性有很大影响,烧蚀产物中,碱金属的电离并不总占优势。     

湍流燃烧的动态自适应直接数值模拟

为了能以较低的计算成本分辨多尺度的湍流燃烧过程,发展了基于动态自适应网格技术的湍流燃烧直接数值模拟方法。在氢气/空气层流预混火焰传播的模拟中,准确计算出了火焰传播速度,表明该方法具有直接求解火焰以及精确捕捉火焰传播过程的能力。在衰减的各向同性均匀湍流的模拟中,借助FFT分析获得了合理的湍动能谱线随湍动能衰减的变化规律,表明该方法能够直接求解湍流。由于在实际应用中,点火发生的位置与低温机理的速率相关,因此在氢气/空气湍流预混火焰中研究了湍流对低温机理的影响,发现湍流输运以及湍流涡团间强烈的拉伸和剪切对低温机理有明显的加速作用,能显著增强放热。     

有压梯度边界层流中湍流斑形成的机制研究

人们虽然对层流向湍流转捩过程的研究已经付出了许多艰辛的努力,但仍然有一个重要的物理现象还没有弄清楚,即有压梯度边界层转捩过程中湍流斑形成的理论机制以及湍流斑的运动特征是什么?这些问题正有待于人们做进一步的深入探索.本文提出一种以壁面局部脉冲的初始小扰动场来模拟有压梯度边界层流中湍流斑形成的物理模型.采用直接数值模拟的方法研究有压梯度边界层流中湍流斑产生的理论机制和发展规律;数值结果显示,它们在好多方面与湍流斑的基本特征相符.     

湍流射流扩散火焰的层流火焰面模拟

采用双尺度κ-ε湍流模型，标量联合的概率密度函数（PDF）输运方程和层流火争面模型相结合，模拟湍流射流扩散火焰，比较了几种不同κ-ε湍流模型对轴对称湍流射流速度场模拟的结果，建立了混合物分数和湍流频率的联合PDF输运方程，并假定湍流频率和标量 线性关系，得到标量耗散率的概率分布，化学组分的浓度和温度由层流火焰面系综的平均得到，计算结果与直接用化学反应动力学计算的结果及实验数据做了比较。     

涡轮叶片型面换热系数计算

本文用数值方法计算了无气膜冷却涡轮叶片上包括前驻点的整个型面上的换热系数。驻点区采用相似方程计算, 并用圆柱绕流近似表示驻点附近流场。较成功地模拟了边界层转捩和过渡区的换热情况。计算中考虑了来流湍流度的影响。      

细长体大攻角分离流动的数值模拟研究

采用数值模拟方法研究了具有头部微小扰动细长旋成体大攻角分离流动.考察了全层流和全湍流情况下头部扰动位置和大小对背风面流动和截面侧向力的影响.结果显示在给定的扰动形式下,层流和湍流背风面流动存在明显的双稳态特征,侧向力随扰动周向位置呈双周期变化;扰动大小对层流流动的影响明显大于对湍流流动的影响.     

压气机静叶栅层流分离泡转捩与角区分离数值模拟与实验

采用γ-Reθ转捩模型对某可控扩散叶型(CDA)平面叶栅全攻角范围进行了三维数值计算，通过对比数值计算结果与叶栅实验吻合较好。在此基础上，分析了进口来流湍流度和雷诺数变化对叶栅表面层流分离、转捩以及角区分离的影响。结果表明:进口湍流度低于5%时，吸力面存在层流分离，当进口湍流度大于5%后，层流分离移除，但转捩会一直存在；随着进口湍流度或雷诺数增加，吸力面和压力面转捩位置均会前移；随着进口湍流度增加，吸力面角区分离会有所减小，雷诺数增加对角区分离的影响不大。      

对流传热场协同原理在高速可压缩边界层流动中的推广

应用理论分析方法对适用于不可压缩层流与湍流流动的对流传热场协同原理进行了可压缩层流与湍流流动上的推广。分析结果表明,与不可压缩流动的对流传热场协同原理不同,可压缩层流与湍流的对流传热取决于流动当地单位体积的动量与总焓梯度的协同。用当地单位体积的动量与总焓梯度的协同研究可压缩流动的壁面传热问题,对层流热流,不但计及了高速流动的密度变化对热流的作用,而且包括了静焓梯度、压力梯度、壁面分子黏性剪切效应对热流的影响;对湍流问题,除了高速流动的密度变化、压力梯度、壁面分子黏性剪切效应对热流的影响外,还计及了雷诺剪切应力对热流的作用。另外,对黏性影响不能忽略的不可压缩流动的对流传热问题,用速度向量与总温（总焓）梯度协同更精确。     

低雷诺数翼型和跨声速翼型边界层计算

本文利用边界层动量积分方程和平均流动能积分方程兼容计算了翼面的层流和湍流边界层流动。文中采用一种e~9型转捩判别公式预测翼面存在层流分离气泡情况的转捩位置。并引入对剪应力系数C_r的滞后方程以体现湍流中雷诺应力的历程效应。为避免在翼面流动分离时边界层方程的奇异,采用反解法计算。算例表明,计算与测量结果吻合良好。     

单面射流冲击淬火参数对7050铝合金厚板冷却的影响

研究了喷嘴直径(D)、相邻喷嘴距离(L)和喷射压力(P)对驻点区和紊流区换热的影响.结果表明:在高压(500kPa)下,紊流区心部换热效果比驻点区心部强;当喷嘴直径和喷射压力一定,相邻喷嘴距离对厚板冷却有影响,距离越小,相邻涡对互相干扰越大,紊流区冷却强度加大;当L/D=8～14时,从470℃冷至300℃的平均冷却速率...     

粗糙壁热交换实验研究

 本文介绍了来流马赫数5的条件下,典型球锥模型的粗糙壁热交换实验结果。模型头部半径R为27.4毫米,底部直径D为60毫米,对五个不同粗糙度的模型进行了实验。模型表面粗糙颗粒直径d分别为0、0.3、0.5、0.7、0.9毫米。风洞前室总压P_t为10～45公斤/厘米。,相应的来流雷诺数R_eD为（O.8～3.6）×10⁶。 实验结果表明:光滑壁模型表面是层流加热,驻点热流与层流理论计算值较一致。粗糙度的影响,在低总压条件下（10公斤/厘米）主要在于促使边界层的转捩和发展。随着风洞总压的提高,物面静压和局部雷诺数的相应增大,粗糙度对热流的影响才明显增强,而严重的区域在端头。在实验最大粗糙度和最大总压条件下（d=O.9毫米、p_t=45公斤/厘米。）,除驻点值外,热流与光滑壁层流驻点值相比(q_i/q_so)的峰值在音点区域且接近4,而在驻点,此模型有别于其它模型,较为特殊,比热流最大值接近6,看来这可能与驻点局部外形变化有关。     

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.      

含湿气流单缝冲击冷却数值模拟

对饱和湿空气夹带液滴的含湿气流单缝冲击冷却进行了数值模拟,比较研究了含湿气流中液滴相变换热/液滴载量比和单缝射流几何参数对强化换热效果的影响,研究发现:气流中液滴的加入对冷却效果的影响十分显著,热流密度为8000W/m2的条件下,空气中注入5%的液滴会比饱和湿空气对驻点区的冷却效果提高90%,当液滴载量比从1%提高到5%,冲击驻点区局部传热系数提高了32%;而合适的缝宽和冲击高度也是影响冷却效果的重要因素,在该条件下冲击高度与缝宽的比值为4.55时冲击冷却效果最佳.      

Experimental investigation on aero-heating of rudder shaft within laminar/turbulent hypersonic boundary layers

The aero-heating of the rudder shaft region of a hypersonic vehicle is very harsh, as the peak heat flux in this region can be even higher than that at the stagnation point. Therefore, studying the aero-heating of the rudder shaft is of great significance for designing the thermal protection system of the hypersonic vehicle. In the wind tunnel test of the aero-heating effect, we find that with the increase of the angle of attack of the lifting body model, the increasement of the heat flux of the rudder shaft is larger under laminar flow conditions than that under turbulent flow conditions. To understand this, we design a wind tunnel experiment to study the effect of laminar/turbulent hypersonic boundary layers on the heat flux of the rudder shaft under the same wind tunnel freestream conditions. The experiment is carried out in the ?2 m shock tunnel(FD-14 A) affiliated to the China Aerodynamics Research and Development Center(CARDC). The laminar boundary layer on the model is triggered to a turbulent one by using vortex generators, which are 2 mm-high diamonds. The aero-heating of the rudder shaft(with the rudder) and the protuberance(without the rudder) are studied in both hypersonic laminar and turbulent boundary layers under the same freestream condition. The nominal Mach numbers are 10 and 12, and the unit Reynolds numbers are2.4 × 10~6 m~(-1) and 2.1 × 10~6 m-1. The angle of attack of the model is 20°, and the deflection angle of the rudder and the protuberance is 10°. The heat flux on the model surface is measured by thin film heat flux sensors, and the heat flux distribution along the center line of the lifting body model suggests that forced transition is achieved in the upstream of the rudder. The test results of the rudder shaft and the protuberance show that the heat flux of the rudder shaft is lower in the turbulent flow than that in the laminar flow, but the heat flux of the protuberance is the other way around,i.e., lower in the laminar flow than in the turbulent flow. The wind tunnel test results is also validated by numerical simulations. Our analysis suggests that this phenomenon is due to the difference of boundary layer velocities caused by different thickness of boundary layer between laminar and turbulent flows, as well as the restricted flow within the rudder gap. When the turbulent boundary layer is more than three times thicker than that of the laminar boundary layer, the heat flux of the rudder shaft under the laminar flow condition is higher than that under the turbulent flow condition. Discovery of this phenomenon has great importance for guiding the design of the thermal protection system for the rudder shaft of hypersonic vehicles.     

转动园盘上可压缩流动换热的计算研究

本文将边界层方程无因次化,用“BOX”算法求解。研究计算了转动园盘上层流区、过渡区和紊流区的扭矩系数和局部换热系数以及变物性的影响,并研究了适合计算转盘问题的双层紊流模型。     

Effects of temperature distribution and level on heat transfer on a rotating free disk

In gas turbine engines, with the existence of the intense forced convection and significant buoyancy effects, temperature distribution and level on turbine or compressor disks affect the heat transfer characteristics strongly. In this paper, numerical simulations were performed to analyze these influences for a free disk, with the laminar and turbulent flow respectively. The influences of temperature distribution on the heat transfer were observed by using incompressible cooling air, and temperature profiles of nth order monomial and polynomial were assumed on the disk. The analysis revealed that the heat transfer for two flow states on the free disk is determined by the exponent n of the monomial profile when specifying the rotating Reynolds number; for an arbitrary polynomial profile, the local Nusselt number can be deduced from results of monomial profiles. To study the effects of temperature level on heat transfer singly, monomial profiles were used and the local Nusselt number of compressible and incompressible cooling air were compared.And both for two flow states, the following conclusions could be drawn: the relative difference of local Nusselt number is mainly controlled by nondimensional local temperature difference, and almost independent of the monomial's coefficient C, exponent n and the rotating Reynolds number. Subsequently, a correction method for heat transfer of the free disk is presented and verified computationally, with which the local Nusselt number, obtained with a uniform and low temperature profile, can be revised by arbitrary distribution and high temperature magnitude.      

MHD控制激波诱导湍流边界层分离的机理分析

为了研究磁流体动力学(Magnetohydrodynamics:MHD)加速边界层对激波-湍流边界层相互作用的影响,用高阶有限差分法求解了小磁雷诺数近似的MHD湍流方程。其中,无粘通量采用WENN格式离散、粘性通量采用Roe平均中心差分离散,时间采用半隐式推进,并采取追赶法求解。计算给出了湍流、电场、磁场和电导率等参数对边界层分离的影响,数值结果显示:在同样的逆压梯度下,湍流边界层分离能更快地趋于稳态流场,且分离区比层流小;通过施加洛仑兹力加速,边界层速度型面变得更加饱满、位移厚度减小、分离点和再附点向激波与固壁的交点靠近,分离区尺寸减小甚至最终被消除。     

客机APU舱温度场的数值计算

为研究大型客机APU舱传热过程和进行结构热分析,采用数值方法计算了APU舱温度场。采用商业CFD软件Fluent,选用Realizable k-ε湍流模型和S2S热辐射模型建立了APU舱内流动和传热数值计算模型,得到了流场、温度场和热流信息。APU舱中主要传热过程是辐射传热;APU防火罩和排气管隔热罩起到了热辐射遮热罩的作用;防火罩内的APU冷却空气对于降低防火罩温度有重要作用;排气管内的引射冷却空气有效降低了壁面温度。