文章快速检索  
  高级检索
单孔容腔瞬态充气换热的理论分析方法
丁水汀1,2, 于航1,2, 邱天1,2, 单晓明3, 贺宜红3     
1. 北京航空航天大学 能源与动力工程学院, 北京 100083;
2. 北京航空航天大学 航空发动机气动热力国家重点实验室, 北京 100083;
3. 中国航发湖南动力机械研究所, 株洲 412000
摘要: 换热对于容腔瞬态响应过程有显著影响,而目前缺乏分析容腔瞬态过程换热的通用方法,导致容腔瞬态响应模拟精度较差。针对这一现状,基于自由射流、冲击射流及外掠平板换热理论,提出了一种模拟非绝热单孔容腔瞬态充气过程换热的理论方法。应用该方法模拟了容腔压力和温度的瞬态响应过程,并与试验数据进行了对比。结果表明:该理论方法的模拟结果与试验数据吻合很好,压力最大相对误差不超过3%,温度最大相对误差不超过1%,验证了理论方法的可行性和准确性。而绝热模型的模拟结果与试验数据相比,压力和温度的最大相对误差分别可达12%和14%,等温模型的压力和温度的最大相对误差分别可达6%和7%,说明理论方法显著提高了容腔瞬态响应模拟精度。同时,理论分析方法不仅具有较强的通用性,还能够极大地降低分析容腔瞬态换热的成本,可以有效支撑空气系统非绝热容腔元件建模。
关键词: 容腔    瞬态响应    充气过程    换热    理论分析    
Theoretical analysis method for heat transfer in transient charging of a cavity with single opening
DING Shuiting1,2, YU Hang1,2, QIU Tian1,2, SHAN Xiaoming3, HE Yihong3     
1. School of Energy and Power Engineering, Beihang University, Beijing 100083, China;
2. National Key Laboratory of Science and Technology on Aero-Engine Aero-thermodynamics, Beihang University, Beijing 100083, China;
3. AECC Hunan Aviation Powerplant Research Institute, Zhuzhou 412000, China
Abstract: Heat transfer has a significant effect on the transient response process of the cavity. However, there lacks a general method for analyzing the heat transfer in the cavity transient process, resulting in poor transient response simulation accuracy. In this situation, based on the theory of free jet, impinging jet, and heat transfer for flow over a flat surface, a theoretical method for simulating heat transfer in a non-adiabatic cavity with single opening transient charging process is proposed. The transient response process of the pressure and temperature of the cavity is simulated and compared with the experimental data. The results show that the simulation results based on the theoretical method agree well with the experimental data, the maximum relative error of pressure does not exceed 3%, and the maximum relative error of temperature does not exceed 1%. The feasibility and accuracy of the theoretical method are verified. Compared with the experimental data, the maximum relative error of pressure and temperature based on adiabatic model can reach 12% and 14%. Also, the maximum relative error of pressure and temperature based on the isothermal model can reach 6% and 7%. These findings show that the theoretical method significantly improves the simulation accuracy. Meanwhile, the theoretical analysis method has strong universality, and can also greatly reduce the cost of analyzing the transient heat transfer of the cavity, which can effectively support the modeling of non-adiabatic cavity component of air system.
Keywords: cavity    transient response    charging process    heat transfer    theoretical analysis    

伴随对航空发动机性能需求的不断提高及航空发动机技术的快速发展,发动机强整体、强瞬变、强耦合、强非线性的特征愈发凸显。发动机运行状态复杂多变,对发动机的安全性造成隐患。现阶段,航空发动机尤其是军用飞机发动机在整个飞行剖面内大部分处于多变的运行环境中,发动机经历慢车到最大转速的时间缩短,瞬变特征显著,严苛的运行条件导致发动机在过渡过程中故障频发[1]

航空发动机是一多系统协调工作的实体,其中空气系统是处于发动机主流道内侧、外侧的转子与静子或转静子之间联系不同腔室、间隙、孔缝的空气系统流路总称[2],具有发动机冷却、封严、轴向力控制等重要作用。腔室类元件由于自身质量和能量的存储和释放效应而具有显著的瞬态特征,在发动机过渡过程中起到至关重要的作用。从压气机引气的气流边界发生变化时,各腔室的压力变化存在滞后性,对压力变化的预测直接影响到推算发动机轴向力的动态变化。这一指标在发动机安全性设计中有明确规定:发动机压力平衡系统需提供合适的轴向载荷[3]

为预测容腔压力瞬态变化,已开展大量关于容腔瞬态响应规律的研究。Dutton和Coverdill[4]通过试验研究了容腔充放气过程,定性分析了描述瞬态变化的两种理论模型(绝热模型和等温模型)基本假设的适用条件,并使用这两种理论模型预测容腔的压力变化。通过与试验结果进行比较,发现模型与试验数据间存在差异,但是文中没有针对真实响应过程提出更为精确的数值模型。Thorncroft等[5]认为容腔瞬态响应是一多变过程,通过试验修正了多变过程的多变指数,使模拟结果与试验数据吻合。这一结果也表明了响应过程为非确定过程,非绝热也非等温,但只针对特定的试验数据进行模型修正,不具有可拓展性。Gallar和Calcagni[6]描述了静止容腔瞬态响应的数学模型,并借助此模型进行了空气系统网络的建模与分析工作。

在工程上,普遍采用绝热假设,即不考虑瞬态过程中气体与腔壁的换热,基于该方法模拟的压力响应过程与实际存在较大差别[4-7]。另一种较常采用的是等温假设,即在瞬态过程中气体温度始终保持不变[4],这种方法过高地估计了换热作用。两种方法的模拟精度都较差,根本原因在于缺少分析瞬态换热过程的有效方法。为了提高模拟精度,许多学者开展了针对容腔瞬态换热的研究。杨丽红等[8]研究了换热效应对于等温容腔放气过程的影响。郭钟华等[9]则针对真空系统容腔的压力响应建立了考虑换热影响的真空抽吸模型,使用了定努塞尔数换热模型,较好地模拟了真实响应过程。但是瞬态过程流动与换热状态随时间变化,换热分析方法应能适应变工况变化。本文作者团队[10-11]针对非绝热容腔瞬态响应过程进行了建模与试验研究,给出了通用的换热特性准则式,提出了一种研究瞬态换热的可行方法。

目前,开展容腔瞬态换热分析的试验成本和CFD模拟成本都很高,若利用绝热模型和等温模型简化计算,模拟精度又较低。所以,本文提出了一种具有较高精度,且可以大幅降低成本的理论方法。以非绝热单孔立方体容腔为研究对象,从容腔大间距冲击流动和外掠平板换热理论出发,建立了容腔瞬态换热的理论分析方法,并对方法进行了试验验证。

1 理论分析方法

在容腔充气过程中,由于发生质量和能量存储的时间尺度远大于建立稳定流场的时间尺度,所以认为在瞬态过程中的任一时刻,流动与换热仅与当前进气流动参数和腔内气体状态参数有关,与瞬变过程无关,基于这一假设建立瞬态换热的理论分析方法。

首先,需明确容腔充气的物理过程。立方体容腔中心横截面上的充气流动结构如图 1所示,气流自圆形喷嘴射出,射流主要流经3个明显的流动区域,即自由射流区、冲击区和壁面射流区。在自由射流区,流动特性与自由射流相同,在冲击区,气流冲击壁面并在所有方向上均匀地沿径向扩散,气流经历了显著的弯曲,在该区末形成平行于壁面的流动,进入到壁面射流区,随后气流撞击到侧壁面继续流动,最终在容腔内形成循环流动。

图 1 容腔中心截面充气流动结构 Fig. 1 Charging flow structure of cavity central section

对于容腔这类大间距冲击射流结构,且喷嘴截面积远小于容腔截面积,喷嘴出口射流受侧壁及前壁影响较小,可近似为自由射流(除紧邻冲击壁面的位置)。以冲击壁面滞止点为原点建立圆柱坐标系,图 1中纵轴为冲击壁面法向z轴,横轴为径向r轴,自由射流轴向速度满足[12]

$ \frac{{{u_z}}}{{{u_0}}} = 6.2\frac{d}{{L - z}} $ (1)
 

式中:uz为轴向速度;u0为射流速度,在充气过程中,随着容腔内气体压力的逐渐升高而减小;d为圆形喷嘴喉道直径;z为距冲击壁面的法向距离;L为喷嘴出口到冲击壁面的距离。式(1)在自由射流处于充分发展段适用,通常在(Lz)/d>40范围内自由射流处于充分发展,各断面的流速分布具有相似性。同时为了保证射流外部对射流发展没有影响,要求侧壁尽量远离射流,需满足R/d≫1,R为侧壁距射流中心的距离。当z/L≥0.05时,uz/u0d/(L-z)线性变化,可用自由射流轴向速度的变化关系式(1)描述。在冲击滞止点处,速度为0。对于三维轴对称无摩擦流,滞止点附近的轴向速度[13]

$ {u_z} = - 2az $ (2)
 

径向速度为

$ {u_r} = ar $ (3)
 

式中:a为常数;r为径向距离。当z/L≤0.05时,冲击射流轴向速度落在滞止点附近的线性变化区[14],取z/L=0.05位置处的速度u0.05作为线性变化区的特征速度,则根据式(2),常数a

$ a = - \frac{{{u_{0.05}}}}{{2 \times 0.05L}} $ (4)
 

气流冲击到滞止点产生高压区,诱使滞止点临近区域流动为层流[13]。首先定义无量纲温度θ=(T-T)/(Tw-T),其中,T为热边界层外层流体温度,Tw为壁面温度,则层流热边界层内温度型[15]可表示为

$ \theta = 1 - \frac{3}{2} \cdot \frac{z}{{{\delta _{\rm{t}}}}} + \frac{1}{2}{\left( {\frac{z}{{{\delta _{\rm{t}}}}}} \right)^3} $ (5)
 

式中:δt为热边界层厚度。壁面换热热流密度为

$ \begin{array}{l} q = - {\left. {\lambda \frac{{\partial T}}{{\partial z}}} \right|_{z = 0}}{\rm{ = }} - \lambda \frac{{{T_{\rm{w}}} - {T_\infty }}}{{{\delta _{\rm{t}}}}} \cdot \\ \;\;\;\;\;\;\;{\left. {\frac{{\partial \left( {\left( {T - {T_\infty }} \right)/\left( {{T_{\rm{w}}} - {T_\infty }} \right)} \right)}}{{\partial \left( {z/{\delta _{\rm{t}}}} \right)}}} \right|_{z/{\delta _{\rm{t}}} = 0}} \end{array} $ (6)
 

其中:λ为导热系数。利用无量纲温度分布,式(6)可简化为

$ q = \frac{3}{2}\lambda \frac{{{T_{\rm{w}}} - {T_\infty }}}{{{\delta _{\rm{t}}}}} $ (7)
 

对流换热系数为

$ h = \frac{q}{{{T_{\rm{w}}} - {T_\infty }}} = \frac{3}{2} \cdot \frac{\lambda }{{{\delta _{\rm{t}}}}} $ (8)
 

基于上述分析,通过求解热边界层厚度δt即可求得对流换热系数。轴对称、定常圆柱坐标系下的连续性方程及忽略黏性耗散的能量方程分别为

$ \frac{{\partial {u_r}}}{{\partial r}} + \frac{{{u_r}}}{r} + \frac{{\partial {u_z}}}{{\partial z}} = 0 $ (9)
 
$ {u_r}\frac{{\partial T}}{{\partial r}} + {u_z}\frac{{\partial T}}{{\partial z}} = \alpha \left( {\frac{{{\partial ^2}T}}{{\partial {z^2}}} + \frac{{{\partial ^2}T}}{{\partial {r^2}}} + \frac{{\partial T}}{{r\partial r}}} \right) $ (10)
 

式中:α为热扩散系数。

假设热边界层内沿流向传递的热流大小远小于沿壁面法向传递的热流大小,即

$ \frac{{{\partial ^2}T}}{{\partial {z^2}}} \gg \frac{{{\partial ^2}T}}{{\partial {r^2}}}, \;\;\;\;\frac{{{\partial ^2}T}}{{\partial {z^2}}} \gg \frac{{\partial T}}{{r\partial r}} $ (11)
 

则式(10)可简化为

$ {u_r}\frac{{\partial T}}{{\partial r}} + {u_z}\frac{{\partial T}}{{\partial z}} = \alpha \left( {\frac{{{\partial ^2}T}}{{\partial {z^2}}}} \right) $ (12)
 

沿壁面法向0~δt对式(12)进行积分,以无量纲温度表示为

$ \int_0^{{\delta _{\rm{t}}}} {{u_r}\frac{{\partial \theta }}{{\partial r}}{\rm{d}}z} + \int_0^{{\delta _{\rm{t}}}} {{u_z}\frac{{\partial \theta }}{{\partial z}}{\rm{d}}z} = \int_0^{{\delta _{\rm{t}}}} {\alpha \left( {\frac{{{\partial ^2}\theta }}{{\partial {z^2}}}} \right){\rm{d}}z} $ (13)
 

又有

$ \int_0^{{\delta _{\rm{t}}}} {{u_r}} \frac{{\partial \theta }}{{\partial r}}{\rm{d}}z = \int_0^{{\delta _{\rm{t}}}} {\frac{{\partial \left( {\theta {u_r}} \right)}}{{\partial r}}} {\rm{d}}z - \int_0^{{\delta _{\rm{t}}}} \theta \frac{{\partial {u_r}}}{{\partial r}}{\rm{d}}z $ (14)
 
$ \int_0^{{\delta _{\rm{t}}}} {{u_z}} \frac{{\partial \theta }}{{\partial z}}{\rm{d}}z = \int_0^{{\delta _{\rm{t}}}} {\frac{{\partial \left( {\theta {u_z}} \right)}}{{\partial z}}{\rm{d}}z} - \int_0^{{\delta _{\rm{t}}}} \theta \frac{{\partial {u_z}}}{{\partial z}}{\rm{d}}z $ (15)
 

将边界条件代入,式(13)可化为

$ \int_0^{{\delta _{\rm{t}}}} {\frac{{\partial \left( {\theta {u_r}} \right)}}{{\partial r}}{\rm{d}}z} - \int_0^{{\delta _{\rm{t}}}} \theta \left( {\frac{{\partial {u_r}}}{{\partial r}} + \frac{{\partial {u_z}}}{{\partial z}}} \right){\rm{d}}z = \left. {\alpha \left( {\frac{{\partial \theta }}{{\partial z}}} \right)} \right|_0^{{\delta _{\rm{t}}}} $ (16)
 

联立式(9),得

$ \frac{\partial }{{\partial r}}\int_0^{{\delta _{\rm{t}}}} \theta {u_r}{\rm{d}}z + \int_0^{{\delta _{\rm{t}}}} \theta \frac{{{u_r}}}{r}{\rm{d}}z = \left. {\alpha \left( {\frac{{\partial \theta }}{{\partial z}}} \right)} \right|_0^{{\delta _{\rm{t}}}} $ (17)
 

将式(3)和式(5)代入式(17),则有

$ \frac{3}{4}a{\delta _{\rm{t}}} = - \alpha \left( {\frac{{\partial \theta }}{{\partial z}}} \right)\left| {_{z = 0}} \right. $ (18)
 

式(18)的推导是基于滞止点区域的热边界层厚度与径向位置无关的结论[16]。又由式(5),可得

$ {\left. {\left( {\frac{{\partial \theta }}{{\partial z}}} \right)} \right|_{z = 0}} = - \frac{3}{2} \cdot \frac{1}{{{\delta _{\rm{t}}}}} $ (19)
 

联立式(18)和式(19),得

$ {\delta _{\rm{t}}} = {\left( {\frac{{2\alpha }}{a}} \right)^{1/2}} $ (20)
 

代入式(8),则冲击壁面滞止点对流换热系数hs

$ {h_{\rm{s}}} = \frac{q}{{{T_{\rm{w}}} - {T_\infty }}} = \frac{3}{2}\lambda {\left( {\frac{{2\alpha }}{a}} \right)^{ - 1/2}} $ (21)
 

文献[17]在分析中指出冲击区内换热努塞尔数与偏离滞止点的距离无关,故认为整个冲击区域内对流换热系数相同。

气流冲击到壁面后,径向速度从滞止点处零先达到峰值,随后由于能量耗散及壁面阻碍减速,气体由冲击区流入到壁面射流区,流动状态已由层流转捩为湍流[18-19]。文献[14]指出冲击射流壁面射流区范围为r/L>0.22时,壁面射流区速度衰减规律[20]

$ \frac{{{u_r}L}}{{\sqrt K }} = 1.32{\left( {\frac{r}{L}} \right)^{ - 1.1}} $ (22)
 

式中:K为冲击射流的动量通量。

恒定均匀壁面温度条件下流体外掠平板湍流换热局部努塞尔数[15]

$ N{u_r} = 0.032Re_r^{0.8}P{r^{0.43}} $ (23)
 

式中:努塞尔数Nur和雷诺数Rer的特征长度为距起始位置的距离rRer的特征速度为壁面射流区速度; Pr为普朗特数。基于外掠平板换热假设,则壁面射流区局部换热努塞尔数可由局部换热特性式(23)推算。这里重点关注平均换热特性,为便于计算平均对流换热系数,将轴对称流动的环形面积积分区域简化为矩形面积积分区域,则壁面射流区平均对流换热系数$\overline {{h_{\rm{r}}}} $可写为

$ \overline {{h_{\rm{r}}}} = 0.032P{r^{0.43}}\lambda \frac{1}{{0.28L}}\int_{0.22L}^{0.5L} {\frac{1}{r}Re_r^{0.8}{\rm{d}}r} $ (24)
 

上述壁面射流区平均对流换热系数是基于流体外掠平板局部换热特性推算得到。实际上,在相似的流动状态下冲击射流引起的流体与壁面的换热程度要强于外掠平板。所以,在此引入了换热加强系数j,则实际壁面射流区平均对流换热系数$\overline {{h_{{\rm{r, r}}}}} = j\overline {{h_{\rm{r}}}} $。基于文献[19]的试验研究,随着无量纲冲击距(冲击距离与喷嘴喉道直径之比L/d)的增大,换热加强系数逐渐减小,文献[19]中虽然没有给出适用于容腔大冲击距结构的结果,但是通过对试验结果的合理外推,可以预估适用于更大冲击距的系数大小。文献[19]中指出,当L/d=8时,j=2.3, 而本文容腔冲击距大于30,通过外推近似取j=2。最终,得到冲击壁面平均对流换热系数$\overline {{h_{\rm{b}}}} $

$ \overline {{h_{\rm{b}}}} = \frac{1}{{0.5L}}\left( {{h_{\rm{s}}} \cdot 0.22L + \overline {{h_{{\rm{r}}, {\rm{r}}}}} \cdot 0.28L} \right) $ (25)
 

冲击壁面换热热流密度为

$ q = \overline {{h_{\rm{b}}}} \left( {{T_{\rm{v}}} - {T_{\rm{w}}}} \right) $ (26)
 

式中:Tv为容腔内气体温度。

相较于冲击壁面,侧壁面及前壁面的换热较弱。从流动结构中可以看出,气流在冲击壁面上沿径向流动扩散,以壁面射流区近侧壁速度卷起冲击侧壁面。假设侧壁换热近似为流体外掠平板换热,并且以卷起速度作为来流速度,则根据外掠平板换热特性可以得到其对流换热系数。虽然近壁面流动速度低,但是由于自由射流本身的不稳定性及内部封闭循环流动结构的存在,流场持续受冲击气流影响,流动不稳定,故按湍流特性准则式(23)计算侧壁面换热。实际上,式(23)有其自身的适用范围,即2×105Rer≤5×106,0.6≤Pr≤50。但是对于封闭容腔近壁面的流动换热,不能直接以该适用范围作为适用准则,因为两者在流动细节上并非完全一致。综合考虑两者具有的相似流动特征和容腔近壁面的真实流动状态,选用此特性作为容腔侧壁面及前壁面换热的理论近似。同时,侧壁面流动是由冲击壁面射流区气流撞击引起,撞击点为边界层形成起始点。所以,侧壁面平均换热努塞尔数$\overline {N{u_{\rm{L}}}} $可表示为

$ \overline {N{u_{\rm{L}}}} = \frac{L}{\lambda }\bar h = \frac{{0.032\mathit{P}{\mathit{r}^{0.43}}L}}{\lambda } \cdot \frac{\lambda }{L}\int_0^L {\left( {\frac{1}{l}Re_l^{0.8}} \right){\rm{d}}l} $ (27)
 

式中:h为侧壁面平均对流换热系数;Rel为以侧壁面沿程距离为特征长度的局部雷诺数。

取冲击壁面径向距离r/(L/2)=0.9位置处的速度作为侧壁面来流速度uc,则根据式(27)可以得到侧壁面平均对流换热系数,进而获得壁面热流。由于该径向位置距侧壁面较近,速度量级较小。通过分析,速度选取位置相对偏差10%,壁面换热热流最大相对偏差5%,引起容腔压力及温度瞬态变化最大相对偏差不超过0.1%,说明速度选取位置偏差对容腔瞬态响应过程造成的影响相对较小。前壁面换热以相同的假设进行量化,仍以uc作为来流速度。仿照冲击壁面平均对流换热系数简化处理方法,则前壁面平均换热努塞尔数$\overline {N{u_{\rm{F}}}} $

$ \overline {N{u_{\rm{F}}}} = \frac{L}{\lambda }\bar h = \frac{{0.032\mathit{P}{\mathit{r}^{0.43}}L}}{\lambda } \cdot \frac{\lambda }{{L/2}}\int_0^{L/2} {\left( {\frac{1}{l}Re_l^{0.8}} \right){\rm{d}}l} $ (28)
 

至此,容腔瞬态充气过程各壁面的换热热流均得到确定。上述建立的方法需满足一定的几何条件,具体为自由射流应处于充分发展,满足(Lz)/d>40;侧壁距射流中心足够远,满足R/d≫1。同时,方法中所采用的经验关系式均经过大量的理论或试验验证,在冲击射流领域普遍使用,具有可信度。

以各壁面热流总和作为容腔零维瞬态模型能量方程中的热流输入,模拟压力和温度的瞬态响应过程。具体的建模方法见文献[11]。

2 方法验证及讨论

通过容腔瞬态充气试验对理论方法进行验证。利用试验获取容腔瞬态充气过程中压力及温度的响应历程,与基于瞬态换热理论分析方法的零维模型模拟结果进行对比,验证方法的可行性。

2.1 试验方法

容腔瞬态充气试验在静止容腔试验台上进行,试验台如图 2所示,其中包括了稳压箱、容腔试验件、数据采集系统等。

图 2 静止容腔试验台示意图 Fig. 2 Schematic diagram of static cavity experimental apparatus

此前,本文作者团队通过瞬态充气换热试验获得了单一进气喷嘴、不同恒定壁面温度条件下的容腔瞬态响应规律,分析了相应条件下的瞬态换热特性[10]。但是本文建立的理论方法的可行性及通用性主要受容腔无量纲冲击距的影响,受壁面温度影响不大。所以,为改变无量纲冲击距,在该试验台的基础上,不改变容腔的几何尺寸,只更换充气喷嘴,进行了多组不同喉道直径喷嘴的容腔瞬态充气试验。

试验装置整体结构及几何尺寸如图 3所示,试验件是立方体容腔,采用金属铝制成,内边长为300 mm,厚度为50 mm,实物如图 4所示。包裹喷嘴管壁外侧的聚氨酯环用于密封及防止腔壁与喷嘴管壁间传热,外层岩棉为隔热材料。

图 3 试验装置 Fig. 3 Experimental apparatus
图 4 容腔实物图 Fig. 4 Photograph of cavity

配备4种不同的渐缩喷嘴,喉道直径分别为3.4、5.0、7.0、9.0 mm,实物如图 5所示。

图 5 渐缩喷嘴实物图 Fig. 5 Photograph of converging nozzle

容腔壁面上布置了3个压力测点,容腔内布置了3个温度测点,压力传感器与容腔内壁面齐平安装,热电偶伸入容腔内部,分别布置在立方体容腔3个顶点附近,该位置与相邻的3个内壁面距离相同,均为35 mm,具体位置如图 6所示。同时,在进气位置上游布置了压力和温度测点用于监测进气边界条件。其中,供气温度是通过布置在气动阀门前30 mm左右位置处的热电偶测试,测点距离喷嘴喉道位置约400 mm。试验开始前,气体处于静止状态,测得的温度为气体总温。对于容腔内气体温度的判定关注的是供气总温,因为对于特定气体,总温就代表了总能量。试验过程中供气温度与环境温度相同,气流流经阀门及喷嘴时,在喷嘴中加速膨胀,此过程中气体温度会略低于喷嘴管壁。但是由于流动过程迅速,在喷嘴中气体与外界基本无换热,可以看作是绝能流动,总温不变,故忽略了气体在喷嘴中的温度变化,认为供气温度与测点位置处温度相同。

图 6 测点布置 Fig. 6 Locations of measurement points
2.2 结果对比与分析

试验进气压力为310 560 Pa,进气温度为298.0 K,容腔壁面热边界为恒温334.0 K(边界的实现见文献[10]),腔内气体初始压力为大气压,初始温度与壁面温度相同。喷嘴喉道直径分别为3.4、5.0、7.0、9.0 mm的试验数据与理论非绝热模型(基于换热理论分析方法)、绝热模型及等温模型的计算结果对比如图 7所示。当容腔内气体压力达到进气压力时,充气试验过程结束。

图 7 不同喷嘴条件下试验、理论非绝热模型、绝热模型及等温模型结果对比 Fig. 7 Result comparisons of experiment, theoretical non-adiabatic model, adiabatic model, and isothermalmodel with different nozzles

图 7中可以看出,压力及温度的试验瞬态响应与理论非绝热模型计算结果吻合很好,压力最大相对误差不超过3%,温度最大相对误差不超过1%。绝热模型与试验数据相比,压力和温度的最大相对误差分别为12%和14%。等温模型与试验数据相比,压力和温度的最大相对误差分别为6%和7%。事实上,绝热模型和等温模型是两种极限模型,分别忽略了换热效应和极大地估计了换热效应。在容腔充气过程中,压力和温度的实际响应过程介于两种模型之间。

从压力瞬态响应规律可以看出:①随着充气过程的进行,气体压力逐渐升高,最后稳定,稳定腔压等于进气压力;②不同喷嘴喉道直径条件下的瞬态响应快慢不同,喉道直径越大,充气流量越大,响应时间越短;③随着喷嘴喉道直径的增大,压力瞬态响应更接近绝热模型,偏离等温模型。喉道直径大意味着充气流量大,即进气焓流量大,转化为储存的气体内能多,引起气体温升大,等壁温条件下换热热流密度大。虽然换热热流大,但是瞬态响应过程更接近绝热模型,偏离等温模型。判断响应过程更符合哪种假设,不是以换热热流的大小为判断依据,而是以进气焓流与换热热流之比作为依据。在喷嘴喉道直径大的条件下,进气焓流与换热热流之比更大,所以能量较多地储存在容腔中,响应过程也更接近绝热模型。

从温度瞬态响应规律可以看出:①随着充气过程的进行,气体温度先升高,后缓慢下降。初始时刻,进气焓流要大于壁面换热热流,进气焓更多地转化为气体内能,导致温度升高。随着壁面换热增强,当进气焓流与换热热流之比小于1时,温度开始下降。②在绝热模型中,进气焓全部转化为气体内能储存在容腔中,模拟的温度持续上升,背离真实变化趋势。③随着喷嘴喉道直径的增大,温度升高过程中达到的峰值越大。随着喷嘴喉道直径的减小,温度下降的趋势则更明显。从前面的分析知道,喉道直径小响应规律偏离绝热模型,接近等温模型,气体温升小,换热热流小。虽然小喉道直径下换热热流小,但是在温升阶段,小温升导致的气体内能增量也偏少,两者综合作用导致温度下降的趋势更明显。

3 结论

本文基于自由射流、冲击射流及流体外掠平板换热等理论,建立了分析容腔瞬态充气换热的理论方法,模拟了非绝热单孔容腔瞬态充气过程的压力和温度变化,并通过试验进行了验证。主要结论为

1) 对容腔瞬态充气过程换热的理论分析方法进行了研究,并基于理论方法模拟了容腔压力及温度的瞬态响应过程。模拟结果与试验数据吻合很好,其中压力最大相对误差不超过3%,温度最大相对误差不超过1%,验证了该方法的可行性及准确性。

2) 在相同的物理边界条件下,喷嘴喉道直径越大,充气过程越接近绝热模型,偏离等温模型,响应速度越快。

3) 相较于绝热模型或等温模型,利用该理论方法不仅提高了模拟精度,而且分析方法更具通用性。将封闭容腔内部流动抽象为基本的自由射流、冲击射流、壁面射流及外掠平板流动单元的组合,该方法适用于具有这类相似流动结构的换热分析。虽然抽象结构与真实流动存在差异,牺牲了一定精度,但提高了方法的可拓展性,使其具备了解决相似问题的能力。同时,理论方法对于研究空气系统容腔类元件的瞬态换热具有一定的指导意义。

参考文献
[1] SHIN J. The NASA aviation safety program: Overview: NASA/TM-2000-209810[R]. Washington, D.C.: NASA, 2000.
[2] 王华阁. 航空发动机设计手册——空气系统及传热分析[M]. 北京: 航空工业出版社, 2001: 1-5.
WANG H G. Design manual of aero-engine-Analysis of air system and heat transfer[M]. Beijing: Aviation Industry Press, 2001: 1-5. (in Chinese)
[3] 中国人民解放军总装备部.航空涡轮喷气和涡轮风扇发动机通用规范: GJB241A-2010[S].北京: 总装备部军标出版发行部, 2010: 25-26, 75-76.
General Reserve Department of PLA. General specification for aircraft turbojet and turbofan engine: GJB241A-2010[S]. Beijing: General Reserve Department Military Standards Publishing Department, 2010: 25-26, 75-76(in Chinese).
[4] DUTTON J C, COVERDILL R E. Experiments to study the gaseous discharge and filling of vessels[J]. International Journal of Engineering Education, 1997, 13(2): 123-134.
[5] THORNCROFT G, PATTON J S, GORDON R. Modeling compressible air flow in a charging or discharging vessel and assessment of polytropic exponent[C]//ASEE Annual Conference, 2007.
[6] GALLAR L, CALCAGNI C. Time accurate modelling of the secondary air system response to rapid transients[J]. Journal of Aerospace Engineering, 2011, 225(8): 946-958.
Click to display the text
[7] 高文君, 刘振侠, 朱鹏飞, 等. 航空发动机静止盘腔瞬态特性数值与实验研究[J]. 推进技术, 2019, 40(3): 496-503.
GAO W J, LIU Z X, ZHU P F, et al. Numerical and experimental investigation of transient response of static disc cavity in aero-engine[J]. Journal of Propulsion Technology, 2019, 40(3): 496-503. (in Chinese)
Cited By in Cnki
[8] 杨丽红, 沈航明, 宋元明. 等温容器放气过程中对流换热模型的研究[J]. 中国机械工程, 2014, 25(18): 2489-2495.
YANG L H, SHEN H M, SONG Y M. Study on convection heat transfer model of isothermal chamber during discharge[J]. China Mechanical Engineering, 2014, 25(18): 2489-2495. (in Chinese)
Cited By in Cnki
[9] 郭钟华, 李小宁, 香川利春. 考虑热传递的真空容腔压力响应研究[J]. 真空科学与技术学报, 2015, 35(1): 74-78.
GUO Z H, LI X N, TOSHIHARU K. Impact of heat transfer on pressure response in vacuum chamber of pneumatic vacuum system[J]. Chinese Journal of Vacuum Science and Technology, 2015, 35(1): 74-78. (in Chinese)
Cited By in Cnki (1)
[10] DING S T, YU H, QIU T. Modeling of the cavity response to rapid transient considering the effect of heat transfer: GT2018-75264[R]. New York: ASME, 2018.
[11] 丁水汀, 于航, 邱天. 非绝热单孔容腔瞬态响应的零维建模[J]. 北京航空航天大学学报, 2018, 44(2): 215-222.
DING S T, YU H, QIU T. Zero-dimensional modeling for transient response of non-adiabatic cavity with single opening[J]. Journal of Beijing University of Aeronautics and Astronautics, 2018, 44(2): 215-222. (in Chinese)
Cited By in Cnki | Click to display the text
[12] ALBERTSON M L, DAI Y B, JENSEN R A, et al. Diffusion of submerged jets[J]. Transactions of the American Society of Civil Engineers, 1950, 115(1): 639-664.
[13] SCHLICHTING H. Boundary layer theory[M]. 9th ed. Heidelberg: Springer, 2017: 118-119.
[14] BELTAOS S, RAJARATNAM N. Impinging circular turbulent jets[J]. Journal of the Hydraulics Division, 1974, 100(10): 1313-1328.
Click to display the text
[15] LIENHARD J H IV, LIENHARD J H V. A heat transfer textbook[M]. 3rd ed. Cambridge, MA: Phlogiston Press, 2006: 301-302.
[16] HARTNETT J P, IRVINE T F. Advances in heat transfer:Volume 13[M]. New York: Academic Press, 1977: 4-5.
[17] BEITELMAL A H, SHAH A J, SAAD M A. Analysis of an impinging two-dimensional jet[J]. Journal of Heat Transfer, 2006, 128(3): 307-310.
Click to display the text
[18] GARDON R, AKFIRAT J C. The role of turbulence in determining the heat-transfer characteristics of impinging jets[J]. International Journal of Heat and Mass Transfer, 1965, 8(10): 1261-1271.
Click to display the text
[19] KATTI V, PRABHU S V. Experimental study and theoretical analysis of local heat transfer distribution between smooth flat surface and impinging air jet from a circular straight pipe nozzle[J]. International Journal of Heat and Mass Transfer, 2008, 51(17-18): 4480-4495.
Click to display the text
[20] POREH M, TSUEI Y G, CERMAK E. Investigation of a turbulent radial wall jet[J]. Journal of Applied Mechanics, 1967, 34(2): 457-463.
Click to display the text
http://dx.doi.org/10.7527/S1000-6893.2019.23221
中国航空学会和北京航空航天大学主办。
0

文章信息

丁水汀, 于航, 邱天, 单晓明, 贺宜红
DING Shuiting, YU Hang, QIU Tian, SHAN Xiaoming, HE Yihong
单孔容腔瞬态充气换热的理论分析方法
Theoretical analysis method for heat transfer in transient charging of a cavity with single opening
航空学报, 2020, 41(1): 123221.
Acta Aeronautica et Astronautica Sinica, 2020, 41(1): 123221.
http://dx.doi.org/10.7527/S1000-6893.2019.23221

文章历史

收稿日期: 2019-06-14
退修日期: 2019-07-09
录用日期: 2019-07-29
网络出版时间: 2019-08-12 16:51

相关文章

工作空间