大型连续式高速风洞热交换器设计关键技术研究

赵波, 黄知龙, 陈吉明, 廖达雄, 裴海涛

赵波,黄知龙,陈吉明,等. 大型连续式高速风洞热交换器设计关键技术研究[J]. 实验流体力学,2022,36(5):16-23. DOI: 10.11729/syltlx20210080
引用本文: 赵波,黄知龙,陈吉明,等. 大型连续式高速风洞热交换器设计关键技术研究[J]. 实验流体力学,2022,36(5):16-23. DOI: 10.11729/syltlx20210080
ZHAO B,HUANG Z L,CHEN J M,et al. Research on key technologies of heat exchanger design for large-scale continuous high speed wind tunnel[J]. Journal of Experiments in Fluid Mechanics, 2022,36(5):16-23.. DOI: 10.11729/syltlx20210080
Citation: ZHAO B,HUANG Z L,CHEN J M,et al. Research on key technologies of heat exchanger design for large-scale continuous high speed wind tunnel[J]. Journal of Experiments in Fluid Mechanics, 2022,36(5):16-23.. DOI: 10.11729/syltlx20210080

大型连续式高速风洞热交换器设计关键技术研究

详细信息
    作者简介:

    赵波: (1983—),男,四川绵阳人,硕士研究生,副研究员。研究方向:风洞气动设计及设备传热分析。通信地址:四川省绵阳市涪城区二环路南段6号(621000)。E–mail:zbugo@163.com

    通讯作者:

    廖达雄: E-mail:ldx19631116@163.com

  • 中图分类号: TK172

Research on key technologies of heat exchanger design for large-scale continuous high speed wind tunnel

  • 摘要: 随着大型连续式高速风洞运行功率和精细化测试要求的提高,对风洞热交换器性能提出了越来越高的要求,集中体现在换热压损性能、温度场均匀性和气流扰动特性3个方面。结合近年研究成果,对大型连续式高速风洞热交换器设计中的关键技术及研究成果进行了综述。分析了风洞热交换器的需求特点,总结了影响各种性能的主要因素。介绍了高效低压损设计技术、温度场均匀性控制技术和气流扰动控制技术,给出了热交换器换热效率和压力损失综合权衡的设计原则,阐述了换热芯体排列方式、冷却水流量及进水方式、热交换器段截面形状等对温度场均匀性的影响,以及来流条件和热交换器结构对气流扰动的特性。
    Abstract: With the increase of the operation power and higher requirement on the fine measurement of the large-scale continuous high speed wind tunnel, higher performance of the wind tunnel heat exchanger is required, which is mainly reflected in the heat transfer and pressure loss performance, temperature field uniformity and flow field perturbation characteris-tics. Based on recent years' research, the key technologies and research results in the design of the large-scale continuous high speed wind tunnel heat exchanger are comprehensively discussed. Firstly, the required characteristics of the wind tunnel heat exchanger are analyzed, and the main factors affecting various performances are summarized. Then the high efficiency low-pressure loss design technology is introduced along with the temperature uniformity control technology and the airflow disturbance control technology. The principles of comprehensive trade-off of the heat transfer efficiency and pressure loss of the heat exchanger are presented. And the influences of the heat transfer core arrangement, cooling water flow and its inlet direction, cross section shape of the heat exchanger segment on the temperature field uniformity are discussed. Finally the flow disturbance characteristics downstream of the heat exchanger under different flow conditions and heat exchanger structures are analyzed.
  • 边界层转捩现象对航空领域中的摩阻和热流具有重要影响。对于高超声速吸气式飞行器来说,若流入进气道的边界层为层流,则有可能因逆压梯度和激波边界层干扰而引起流动分离,但湍流边界层在抵抗分离方面具有优势。因此,为了提高飞行稳定性,同时增加燃烧室内燃料和氧气的混合,经常会在进气道前体上安装强制转捩装置,以促进边界层在进入唇口之前转捩为湍流。然而,尽管对边界层转捩现象的研究已经长达一个多世纪,人们仍然不能确切地给出什么情况下边界层会由层流转捩成为湍流,特别是在高超声速复杂构型下,由于受到太多因素(湍流度、雷诺数、马赫数、粗糙度、钝度、壁温、几何构型等)及各种因素之间互相耦合作用的影响,使得转捩预测更加困难。

    现阶段,有多种数值方法可以用来进行边界层转捩研究。DNS和LES被用来研究典型状态下的基本转捩机制和详细的转捩过程[1-2];LST和PSE则被用来研究非稳定波的形成和发展历程;转捩准则则在总结大量实验数据的基础上,用简单的形式给出转捩位置[3-4],如NASP的尖前缘平板转捩准则:Reθ/Mae=305[5]。目前来看,半经验的eN方法在工程上的应用最为广泛,但是eN方法也有几个致命缺陷:由于其基于线性稳定性理论,故难以处理较为复杂的构型,更无法模拟粗糙颗粒诱导的转捩;同时N因子的确定也需要大量的试验数据,而且其在不同马赫数、不同构型下的具体数值也不一样;而且,来流湍流度本是影响转捩的一个重要因素,但在eN方法中却没有任何体现。相反,基于RANS模型发展而来的工程转捩模型,则因其具有较好的鲁棒性且易于与现有的湍流模型框架融合而在最近越来越得到人们的重视,其中,Langtry和Menter[6-8]基于SST湍流模型,提出了一种基于经验耦合的间歇因子转捩模型——γ-Reθ模型。在该模型中,所有的计算都是基于局部变量,没有任何沿流线积分的操作,因此可以轻易地对复杂构型进行模拟,Langtry将该模型应用于大量低速工程算例。由于γ-Reθ转捩模型只是为把基于经验关系式的模型结合进一般CFD方法中提供了一个框架,而具体的物理机制均包含在经验关系式中,因此很多研究者在其框架下针对各种不稳定机制对模型进行修正,从而适应不同情况下的转捩模拟。如Watanabe[9]通过类比Görtler不稳定性和横流不稳定性,引入类似于Görtler数的Kohama参数,Krumbein[10-11]针对三维边界层,通过求解Falkner-Skan-Cooke近似解引入横流转捩动量厚度雷诺数,各自提出了基于该框架下能够模拟横流模态的转捩模型。Langel[12]则通过引入“粗糙幅度”的输运方程发展了能够反映壁面粗糙度影响的修正版本。在高超声速领域,Krause[13]和Zhang[14]分别提出了可以应用于高超声速情况下转捩预测的经验关系式及相应的修正,并都将其应用于Ma8.3的双楔平板模型,得到了与试验吻合较好的壁面压力分布。Cheng[15]则添加了一个针对压力梯度的修正因子,对Ma7.93的7°尖锥进行了模拟,研究了来流湍流度和网格对计算结果的影响。Bensassi[16]模拟了Ma8的7°尖锥模型,得到了与试验一致的壁面热流分布。但是,目前并没有见到该模型在更复杂构型上的应用,如光滑壁面或安装有转捩带的高超声速进气道构型。

    本文基于Langtry和Menter的γ-Reθ模型,进行一些必要的修正,建立一套适用于高超声速情况下的转捩预测方法。通过较为简单的三维构型自然转捩和粗糙颗粒诱导强制转捩的算例对模型进行了验证。将修正的模型应用于X-51A前体边界层自然转捩和强制转捩的研究中,通过与试验结果进行比较,表明本文建立的方法能够较好地进行高超声速复杂构型中的自然转捩和强制转捩预测。

    本文的γ-Reθ模型是基于一套有限体积大规模并行结构网格Navier-Stokes求解器AHL3D建立的。该求解器由中国空气动力研究与发展中心吸气式高超声速技术研究中心开发,能够对理想气体和化学非平衡状态下的混合气体进行求解[17-18]。时间积分采用LU-SGS,无粘通量采用重构-推进方法,重构采用3阶MUSCL插值,推进采用Steger-Warming和AUSMPW+格式,粘性通量采用改进的Gauss定理计算。湍流模型为Menter的SST双方程模型。

    Langtry和Menter[6]提出的γ-Reθ转捩模型主要基于间歇因子γ和转捩起始动量厚度雷诺数Reθt的输运方程。γ方程用来控制转捩程序,其源项由临界动量厚度雷诺数Reθc和最大涡雷诺数Rev之间的关系来控制,而Reθc又与Reθt存在对应关系。Reθt则被当做一个输运标量,主要思想就是使用经验耦合关系式得到边界层外的Reθt,并且允许Reθt通过输运方程扩散到边界层内去。具体形式如下:

    (1)

    (2)

    (3)

    (4)

    (5)

    (6)

    (7)

    (8)

    (9)

    (10)

    (11)

    (12)

    其中,PγEγ分别为γ方程的产生项和破坏项,式中的为应变率张量。Pθt输运方程的产生项, Reθt为边界层外的动量厚度雷诺数, 与当地湍流度和压力梯度有关, λθ在这里表征压力梯度的参数。μμt则分别表示分子粘性系数和湍流粘性系数。Pγ项中有两个关键因素,决定转捩开始的阈函数Fonset和转捩区域长度参数Flength。在上游层流区,Fonset的值为0,当当地涡雷诺数Rev大于2.193倍临界动量厚度雷诺数Reθc时,Fonset被激活,值由0开始转换为1,间歇因子γ逐渐增加,流动由层流逐渐转捩为湍流。而ReθcFlength又都是与转捩起始动量厚度雷诺数有关的参数,Langtry通过大量平板试验数据成功将FlengthReθc这两个参数与关联起来,如式(8)、(9)所示,具体形式可参考文献[7],这两个关联式包含了转捩的内在物理机制,很多基于该转捩模型的修改版本均是对这两个关联式进行了相应的修正。以上两个输运方程通过间歇因子对SST模型中k方程的产生项和破坏项进行修正,如式(10)所示,其中PkDk分别表示SST模型中k方程的产生项和破坏项。γsep则是用来捕获并模拟分离诱导转捩的,可以通过调节式(11)中s1的大小来调整模型对分离诱导转捩的敏感性,Langtry在低速情况下将其标定为2.0。则表示与转捩模型耦合之后k方程的产生项和破坏项。σf, σθt, ca1, ca2, ce1, ce2均为常数,具体数值可参考文献[7]。

    原始的转捩模型是针对低速情况提出的,因此不能直接应用于高超声速情况下的转捩预测,这个模型一经提出,很多研究者均在其框架下针对各种情况对模型的ReθcFlength进行修正,前者用来修正转捩起始位置,后者用来修正转捩区域长度。如Krause和张晓东均提出了该框架下能够应用于高超声速情况下的经验关系式。Krause[13]的修正如下:

    (13)

    (14)

    将来流粘性比RT设置成0.001,因此Tu在到达模型前缘时均衰减到了一个很小的值,通过将Tuinlet加入到经验关系式中来体现不同来流湍流度对转捩的影响规律,经过系列低速平板的标定后,应用于Ma8.3下双楔模型试验的转捩预测中,得到的壁面压力系数与试验吻合较好,显示出高超声速情况下预测转捩的潜力。尚庆[19]也采用该修正关系式对双楔模型进行了更进一步的研究,其结果显示基于该修正关系式的γ-Reθ转捩模型能够在压力系数、热流Stanton数和分离区大小上与试验吻合良好。

    张晓东[14]的修正则如下所示:

    (15)

    (16)

    张晓东还基于Reshotko[20]对高超声速试验中Reθt与边界层外缘马赫数Mae关系的总结,利用当地马赫数对模型进行了马赫数修正,以适应不同Ma数下的转捩预测情况:

    (17)

    张晓东也对双楔模型进行了与试验吻合较好的模拟。

    本文利用文献[21]中的进气道构型,采用二维计算对转捩模型进行考核及修正。进气道外压缩形式为四波系平面顶压,流道全长为550mm,本文模拟其唇口之前压缩面的情况,模型构型及网格如图 1所示。来流条件如表 1所示,表中case3模拟的是在距前缘88mm的第一道压缩面上安装了高度为1mm的钻石型转捩装置的试验结果,转捩带的具体参数可参考文献[22],由于是二维模拟,转捩带构型模化为矩形。计算时为了更好地与试验对比热流值,壁面温度根据实验测量结果分段给定。对比用到的试验数据来自于文献[21-22]。

    图  1  高超声速进气道二维模型及网格示意图
    Fig.  1  Grid of two dimensional hypersonic inlet
    表  1  二维计算时进气道来流条件
    Table  1  Free stream condition for 2-D inlet simulation
    Case Ma 总压/
    MPa
    总温/
    K
    迎角/
    (°)
    Re/
    (107·m-1)
    转捩类型
    1 5.96 2.00 473.2 1.0 1.82 自然
    2 4.97 1.01 357.7 -0.6 2.21 自然
    3 5.96 2.00 473.2 1.0 1.82 强制
    下载: 导出CSV 
    | 显示表格

    试验在FL-31风洞中进行。FL-31属于常规噪声风洞,根据一般的认识,常规噪声风洞的来流湍流度约为1%的量级。在case1情况下,Menter、Krause及张晓东三种版本的经验关系式在不同来流湍流度下的计算结果与试验结果比较如图 2所示。从图中可以看出,Menter版本和Krause版本在来流Tu=0.03时(图中Tu0.03表示来流湍流度为0.03%,RT1表示来流湍流粘性比μt/μ为1),转捩位置仍早于试验值,明显过早地预测了转捩的出现。因此,在本程序框架下,Menter版本和Krause版本将无法正确模拟FL-31风洞湍流度水平下的转捩,更无法模拟静音风洞或飞行试验条件等低湍流度情况下的转捩,而张晓东版本则能在本程序框架下正确反映FL-31风洞的来流湍流度的量级,因此选取张晓东提出的经验关系式作为本文采用的版本。

    图  2  不同计算方法得到的热流分布(case1)
    Fig.  2  Heat flux distribution of different methods (case1)

    将张晓东的经验关系式应用于case2和case3情况,计算结果分别如图 34所示。从图 3中可以看出,当Ma=5.0时,预测的转捩位置要明显晚于试验值。根据文献[20]中的结果,转捩动量厚度雷诺数会随着马赫数的增大而增大,因此转捩位置会随着马赫数的增大而后移,张晓东也是基于其结果添加了马赫数修正。但从图 34的对比结果来看,当Ma6条件下转捩位置与试验结果吻合时,Ma5的修正则过量了,因此模拟得到的转捩位置比试验值要更晚。本文在调试计算过程中发现,直接加入Forsythe[23]的可压缩性修正效果较好,不同马赫数下的修正量有较大差别,与试验结果吻合更好,因此,本文弱化了张晓东的马赫数修正,添加了Forsythe的可压缩性修正。从图 4中则可以看到,其在有粗糙颗粒情况下得到的转捩位置明显迟于试验结果,这是因为在边界层内安装粗糙颗粒后,会给边界层内带来较大扰动并导致提前转捩,但模型中却缺乏体现粗糙颗粒的作用机制。本文通过计算发现,通过增加分离诱导转捩参数s1,提高模型对分离诱导转捩的敏感性,从而可以提高模型对粗糙颗粒的敏感性,增加对强制转捩的模拟能力。此外,张晓东的关系式在Tuinlet=2.5附近时,均会出现奇点,因此本文在经过大量算例计算、修改了可压缩性修正并将分离诱导转捩参数s1由2改为20后,对其原始的关系式也进行了适当的修改,修改之后的计算结果如图 5~7所示。

    图  3  自然转捩热流分布(case2)
    Fig.  3  Heat flux distribution of natural transition (case2)
    图  4  强制转捩热流分布(case3)
    Fig.  4  Heat flux distribution of forced transition (case3)
    图  5  自然转捩热流分布(case1,修正后的模型)
    Fig.  5  Heat flux distribution of natural transition (case1, modified model)
    图  6  自然转捩热流分布(case2,修正后的模型)
    Fig.  6  Heat flux distribution of natural transition (case2, modified model)
    图  7  强制转捩热流分布(case3,修正后的模型)
    Fig.  7  Heat flux distribution of forced transition (case3, modified model)

    修改后的k-w方程为:

    (18)

    其中:

    两个经验关系式为:

    (19)

    (20)

    此外,对Reθt的修正改为:

    (21)

    将修正后的转捩模式在1.3节中的二维构型中进行网格无关性验证。选取1.3节中的case2,并以1.3节中的网格为base网格,网格在第一至第四道压缩面的流向上分别布置了121、61、41、65个网格点,在法向上布置了141个网格点,第一层网格距离壁面0.001mm。总的网格量约为6万。再在此基础上在法向加密1.5倍得到v1.5网格,在法向和流向各加密1.5倍得到v1.5-s1.5网格。3套网格的计算结果如图 2所示。从中可以看出,3套网格计算得到的热流几乎重合,可认为转捩模型此时已经达到网格无关性。

    由于上一节对模型进行的修正都是以二维计算为基础,为了验证其是否能够适应更宽的马赫数、雷诺数范围及三维复杂构型的转捩预测,选取了Ames的全尺寸升力体构型在有迎角情况下的转捩,验证模型对高超声速三维构型自然转捩的预测能力;通过Ma6平板上单个粗糙颗粒诱导转捩的算例,验证模型对粗糙颗粒诱导强制转捩的预测能力。最终,将模型应用于X-51A进气道构型的转捩模拟当中。

    Lockman[24]对Ames研究中心的all-body三维升力体提供了一个完整的高超声速三维实验,模型是一个用于高超声速巡航的通用升力体外形。实验是在Ames的3.5英寸高超声速风洞开展的。模型示意图见图 9。模型全长0.9144m,具有75°后掠角的尖三角翼形状。全体是一个主轴比4的椭圆锥,后体的横截面是具有直边的椭圆形。实验条件如表 2所示。网格拓扑如图 10所示,保持壁面第一层y+ < 1,网格量约为150万。来流湍流度设为0.5%。

    图  8  多套网格计算结果(case2,修正后的模型)
    Fig.  8  Heat flux distribution of three meshes (case2, modified model)
    图  9  全尺寸升力体模型示意图
    Fig.  9  Configuration of the all-body model
    表  2  全尺寸升力体模型风洞试验来流条件
    Table  2  Test conditions for the all-body model
    Ma Re/m-1 总温/K 静温/K 壁面温度/K 迎角/(°)
    7.4 15×106 722 62 300 15
    下载: 导出CSV 
    | 显示表格
    图  10  全尺寸升力体模型网格示意图
    Fig.  10  Grid of the all-body model

    试验测量了迎风面和背风面表征热流的Stanton数,即St,其定义为:

    ρU分别表示来流密度和速度,Hw为壁面上流体的焓值,Ht表示来流总焓。

    图 11为中心线St数分布,对于迎风面,测量到的St数显示转捩大概在距前缘0.25m的位置发生,但是原始的γ-Reθ模型给出的转捩位置却在距前缘0.06m左右,这是由于高马赫数和可压缩性导致了转捩起始动量厚度雷诺数的增加,而这一增加效果在原始的转捩模型中被忽略掉了,因此转捩位置被大大前移。所以原始的转捩模型和全湍流下的SST模型均在x=0.25m的上游区域给出了远高于试验值的热流。而本文修改的模型,不仅准确地捕捉到了转捩的位置,壁面热流也与试验吻合较好。转捩结束位置的热流大约是层流区的3倍,并且要稍高于全湍流计算的结果。对于背风面,由于壁面温度与恢复温度比较接近,热流值都比较小,因此不同模型得到的热流相差不大。图 12所示为迎风面的St数分布和摩擦力线。原始的γ-Reθ和全湍流SST模型给出了相似的St数分布。本文修改的模型计算结果则显示:在上游区域,热流情况基本与层流结果一致, 0.23m < x < 0.28m为热流明显增加的区域,也即转捩区域,之后热流达到全湍流值。摩擦力线显示:层流情况下拐角之后因压力梯度产生的底部横流,在湍流情况下得到了有效的抑制。

    图  11  中心线St数分布
    Fig.  11  Centerline Stanton number distribution
    图  12  迎风面的St数分布和摩擦力线
    Fig.  12  Stanton number distribution and friction force lines of the windward

    Tirtey和Chazot[25-26]在VKI H3高超声速风洞中对单个粗糙颗粒诱导平板转捩情况进行了一系列试验研究。粗糙颗粒构型包括圆柱形、钻石形、斜坡形和半球形,本文选取圆柱形和钻石形粗糙颗粒在高雷诺数下的情况进行模拟。模型如图 13所示,一块长290mm、宽100mm,前缘钝度0.5mm的平板,在距前缘60mm处安装了一个粗糙颗粒,粗糙颗粒细节如图 13(b)所示,其中a=4mm,对于圆柱形粗糙颗粒,k=1.02mm,对于钻石形粗糙颗粒,k=0.83mm。来流马赫数为6,来流单位雷诺数为2.7×107/m,来流总压3.1MPa,总温为500K,壁温在试验开始时均为(292±2)K,试验结束后层流区最大温升为4K,湍流区最大温升为12K,由于这里壁温变化都不大,因此选取300K为壁面温度。粗糙颗粒附近的网格如图 14所示,第一层壁面y+ < 1,网格量为500万左右。来流湍流度设为1.5%。

    图  13  模型示意图
    Fig.  13  Model configuration
    图  14  粗糙颗粒附近网格示意图
    Fig.  14  Mesh of the roughness elements

    此算例中St数定义为:

    Tt表示总温,Tw表示壁温。圆柱形粗糙颗粒尾迹区的红外热试验和计算的St数分布对比如图 1516所示,试验和本文修正的γ-Reθ转捩模型均显示粗糙颗粒诱导的尾迹区的热流明显超出了非尾迹区,而对于原始的γ-Reθ模型,尾迹区和非尾迹区均出现了转捩现象,当来流湍流度被修改为0.5%时,非尾迹区在下游阶段仍出现了转捩现象,而尾迹区的转捩现象却比试验更靠后。这表明原始的γ-Reθ转捩模型不仅对自然转捩的预测比实际情况早,而且对粗糙颗粒对转捩的促进作用模拟不够。而修正后的模型克服了以上两个缺点,在尾迹区和非尾迹区均给出了与试验吻合较好的结果。将修正后的模型计算得到的尾迹区和非尾迹区的St数取出,与试验值以及非尾迹区的全湍流值进行对比,如图 17所示。当流体流过粗糙颗粒后,尾迹区的壁面热流出现了一个快速增长区,之后由于分离涡的作用,其值降低到一个较低值,而后再慢慢增长到全湍流值。St数在尾迹区的值是非尾迹区的2.5~3倍。由于全湍流计算结果中粗糙颗粒后缘的尾迹区也有一个热流值较低的区域,因此可认为转捩在粗糙颗粒后缘就已经完成。试验与计算的对比显示:修正的模型可以模拟粗糙颗粒诱导转捩的整个过程,也能较好地预测壁面热流。但是,在试验中,随着尾迹区的流体向下游流动,非稳定结构在展向发展明显,因此转捩区域越来越宽,而在本文的计算中,转捩区域在展向的发展却比较缓慢。图 16(a)(b)图 17(b)显示的是修改后的模型模拟钻石形粗糙颗粒的情况,结果与圆柱形粗糙颗粒比较类似。

    图  15  红外热试验和计算得到的St数分布(圆柱形粗糙颗粒诱导转捩)
    Fig.  15  Infra-red visualization and calculated Stanton number distribution for cylinder roughness elements
    图  16  红外热试验和计算得到的St数分布(钻石形粗糙颗粒诱导转捩)
    Fig.  16  Infra-red visualization and calculated Stanton number distribution for diamond roughness elements
    图  17  尾迹区和非尾迹区的St数对比
    Fig.  17  Comparison of St number in and out of wake region for cylinder and diamond roughness elements

    图 18则为粗糙颗粒下游流向涡的发展情况,其中展向切片上显示的是流向涡量的大小|ωx|(只显示了|ωx|>2000s-1的情况)。从图中可知,对于两种形式的粗糙颗粒,下游均发展形成了由粗糙颗粒后缘发展起来的主涡和粗糙颗粒两侧形成的次涡,并在继续向下游的发展过程中慢慢耗散掉,流动则转捩为湍流。但本文的计算结果显示其流向涡破碎之后并没有往展向有明显的拓展,因此转捩区域在展向也并未出现明显的展向拓展,考虑到尾迹导致的转捩位置区域很靠前(图 15(a))和很靠后(图 15(c))时,尾迹往展向的发展速度均很缓慢,因此这可能是由于计算方法中耗散过大致使未能捕捉到足够的尾迹流动细节造成的。而过大的耗散则是因雷诺平均引入的湍流粘性和MUSCL插值造成较大的数值耗散共同导致的。

    图  18  粗糙颗粒下游的流向涡发展情况
    Fig.  18  Development of stream wise vortex downstream of the roughness elements

    目前已有的转捩预测方法虽然很多,但是研究的多为平板、圆锥等简单构型,而针对高超声速进气道构型进行完整模拟的还比较少。王亮[27-28]用其提出的三方程转捩模型对X-51A 20%缩比模型在普渡大学Ma6静音风洞中的自然转捩试验结果进行了模拟,结果显示该转捩模型抓住了静音状态下流动一直保持为层流的特征,也较准确地刻画了噪声状态下的转捩起始位置,但是在转捩结束位置上与试验结果还存在较大差别。Orlik[29]则采用了风洞试验与数值模拟相结合的方法,研究了一个高超声速前体的自然转捩情况,其计算采用的是基于线性稳定性分析的eN方法,由于对扰动增长率的积分并不是从中性线上开始的,而是从距前缘1.33mm的固定位置开始,因此积分起始位置的N值是未知的,造成了对N值的刻画出现较大的不确定性。本文希望能够针对进气道模型同时对自然转捩和强制转捩进行模拟,从而考核本文建立的转捩模型对高超声速复杂构型的转捩预测能力。

    Brog[30]在Boeing/AFOSR Mach-6静音风洞中采用20%的缩比模型对X-51A进气道转捩情况进行了一系列试验研究。主要研究了噪声对光滑壁面的自然转捩和粗糙颗粒诱导的强制转捩的影响。模型如图 19所示。模型总长34.42cm。为评估本文的模型对复杂构型转捩的模拟能力,计算了3个算例,并将计算结果与试验数据进行了对比。计算设置如表 3所示。由于没有偏航角,为节省计算量,只计算了一半的模型,在中心面上采用对称边界条件。光滑壁面模型的网格为400万左右,安装了钻石形转捩带的模型网格量为600万左右。壁面第一层网格y+小于1。试验结果显示壁面温度在整个试验过程中变化很小,因此将壁温设为300K。湍流度在噪声条件下为3.0%,静音条件下为0.05%。由于3个算例的马赫数和雷诺数都非常接近,其壁面摩阻和热流差别也很小,因此在计算全层流和全湍流结果时,只计算了case5,并用其来判断3个算例的转捩位置。

    图  19  20%缩比X-51A前体构型示意图
    Fig.  19  Windward side of 20% scale X-51A forebody configuration
    表  3  20%缩比X-51A构型的试验条件
    Table  3  Test conditions for the 20% scaled X-51A forebody configuration
    Case Ma Re/m-1 总压/kPa 总温/K 来流湍流度 转捩类型
    4 6.00 6.59×106 586 418 0.05% 自然
    5 5.80 7.40×106 621 424 3.00% 自然
    6 5.78 7.35×106 614 424 3.00% 强制
    下载: 导出CSV 
    | 显示表格

    图 20为噪声情况下模型迎风面测量得到的壁面温度与计算的摩阻系数。试验与计算结果均在压缩拐角之后开始上升,但测量的壁面温度显示转捩位置在中间最靠后,在两侧区域最靠前;而计算结果显示转捩区域在展向分布较为均匀。从图 21可知,安装了转捩带后,测量的壁面温度在压缩拐角之前就明显大于光滑壁面情况,计算得到的摩阻系数则在转捩带之后立即上升到较高值,表明转捩位置位于压缩拐角和转捩带之间。

    图  20  迎风面测量得到的壁面温度和计算得到的摩阻系数(case5)
    Fig.  20  Measured surface temperature and calculated skin friction for windward surface (case5)
    图  21  迎风面测量得到的壁面温度和计算得到的摩阻系数(case6)
    Fig.  21  Measured surface temperature and calculated skin friction for windward surface (case6)

    图 22为迎风面case4~case6的壁面热流及case5的全层流和全湍流情况。静音条件下,热流分布与层流情况基本一致,表示转捩并没有发生;噪声情况下,拐角之后出现了明显的热流抬升,即转捩现象。而安装了钻石形转捩带之后,转捩位置被推到了转捩带和压缩拐角之间的位置。热流结果与摩阻系数得到的结论吻合。摩擦力线则显示:在层流和静音模式下,由于压力梯度的作用,进气道下表面有明显的横流溢流,会导致进入进气道的流量降低;而对于全湍流或者转捩的状况,这样的横向流动得到了较为明显的抑制,有利于提升进气道捕获流量。在图 23中,噪声条件下,计算的中心线热流显示转捩一经过压缩拐角就开始,而在x=0.28m的位置完成。测量的壁面温度显示转捩起始于压缩拐角,结束于x=0.27m的位置。而静音条件下计算与试验结果都显示没有出现转捩现象。在静音和噪声情况下,试验与计算得到的转捩起始和结束位置均吻合较好。图 24中,在中心线位置,计算的摩阻系数和测量的壁面温度在压缩拐角的位置均已经超过了光滑壁面的情况,进一步表明边界层在压缩拐角之前就已经转捩为湍流。从计算结果可知,转捩带将转捩完成位置前推到了距前缘0.08m的位置。

    图  22  计算得到的迎风面热流对比
    Fig.  22  Calculated heat flux of windward side
    图  23  中心线上测量得到的壁面温度和计算得到的摩阻系数(case4, case5)
    Fig.  23  Measured surface temperature and calculated skin friction of the center line (case4, case5)
    图  24  中心线上测量得到的壁面温度和计算得到的摩阻系数(case6)
    Fig.  24  Measured surface temperature and calculated skin friction of the center line (case6)

    在RANS求解器框架下,通过添加可压缩性修正,修改分离诱导转捩参数,并对经验耦合关系式进行修改,搭建了针对高超声速情况下的γ-Reθ转捩模型。用原始的γ-Reθ模型和修正后的模型对Ma7.4的Ames全尺寸模型和Ma6的单个粗糙颗粒诱导平板转捩进行了计算。结果显示:原始模型在计算自然转捩时,转捩位置过于靠前,而且不能捕捉粗糙颗粒对转捩的促进作用;而修正模型克服了原始模型的这两个缺点,得到了与试验比较一致的结果。最后,将修正后的模型应用于20%缩比X-51A前体进气道模型的自然转捩和强制转捩计算,结果显示模型能够反映高湍流度和粗糙颗粒对转捩的促进作用,表明该模型有望应用于高超声速复杂构型自然转捩和强制转捩的预测研究。

    当然,由于该模型目前仅将来流湍流度作为来流扰动条件,而高超声速来流扰动中的湍流度、扰动频率、扰动模态均可影响转捩位置,因此,如何使模型更有效反映来流扰动,特别是来流压力扰动发展成为第二模态后对转捩的影响,有待进一步研究。

  • 图  1   热交换器性能的主要影响因素

    Fig.  1   Main influencing factors of heat exchanger performance

    图  2   不同换热管排列方式下的温度分布

    Fig.  2   Distribution of temperature in different arrangements of heat transfer tubes

    图  3   某大型连续式高速风洞热交换器模块布局方案

    Fig.  3   Module layout scheme of heat exchanger in a large-scale continuous high speed wind tunnel

    图  4   多排管束气流出口温度分布

    Fig.  4   Outlet gas temperature distribution of multi-row heat transfer tubes

    图  5   某风洞热交换器模块布局方案

    Fig.  5   Module layout scheme of heat exchanger in wind tunnel

    图  6   带包覆外壳体腔体的气流温度分布

    Fig.  6   Gas temperature distribution of heat exchanger with stagnant zone

    图  7   气流温差沿程发展变化曲线

    Fig.  7   Variation curve of gas temperature difference

    图  8   某型板翅式热交换器试验件结构

    Fig.  8   Structure of the plate–fin heat exchanger

    图  9   封条结构对下游流向湍流度的影响[15]

    Fig.  9   Effect of seal on downstream flow direction turbulence[15]

    图  10   拼接缝结构对下游流向湍流度的影响[15]

    Fig.  10   Effect of joint on downstream flow direction turbulence[15]

  • [1] 刘晓波,廖达雄,张巧芸. 大型连续式跨声速风洞热交换器概述[J]. 实验流体力学,2009,23(1):99-104.

    LIU X B,LIAO D X,ZHANG Q Y. A review of heat exchanger study for large-scale continuous transonic wind tunnels[J]. Journal of Experiments in Fluid Mechanics,2009,23(1):99-104.

    [2] 赵波,廖达雄,陈吉明,等. 连续式跨超声速风洞热交换器设计技术初步研究[J]. 空气动力学学报,2010,28(6):696-702. DOI: 10.3969/j.issn.0258-1825.2010.06.014

    ZHAO B,LIAO D X,CHEN J M,et al. Investigation of heat exchanger design for continuous transonic and supersonic wind tunnel[J]. Acta Aerodynamica Sinica,2010,28(6):696-702. doi: 10.3969/j.issn.0258-1825.2010.06.014

    [3] 冯丽丽,杜小泽,杨勇平,等. 椭圆管矩形翅片间空气流动的扰流特征[J]. 工程热物理学报,2011,32(1):119-122.

    FENG L L,DU X Z,YANG Y P,et al. Characteristics of air flow around elliptical tubes with rectangular fins[J]. Journal of Engineering Thermophysics,2011,32(1):119-122.

    [4] 杨涛,袁益超. 管束结构对开缝翅片椭圆管换热器性能的影响[J]. 化工学报,2018,69(4):1365-1373. DOI: 10.11949/j.issn.0438-1157.20170877

    YANG T,YUAN Y C. Effects of tube bundle structure on heat transfer and resistance characteristics of slit finned elliptical tube heat exchangers[J]. CIESC Journal,2018,69(4):1365-1373. doi: 10.11949/j.issn.0438-1157.20170877

    [5] 刘景成,张树有,周智勇. 板翅换热器流道结构改进与流体流动性能分析[J]. 机械工程学报,2014,50(18):167-176. DOI: 10.3901/JME.2014.18.167

    LIU J C,ZHANG S Y,ZHOU Z Y. Analysis of channel structure improvement and its influence on fluid flow in plate-fin heat exchanger[J]. Journal of Mechanical Engineering,2014,50(18):167-176. doi: 10.3901/JME.2014.18.167

    [6]

    KILGORE W A, BALAKRISHNA S, BOBBITT C W, et al. Recent enhancements to the national transonic facility[R]. AIAA-2003-0754, 2003. doi: 10.2514/6.2006-507

    [7]

    KRYNYTZKY A J, HERGERT D W. Boeing Transonic Wind Tunnel upgrade assessment[R]. AIAA-2002-2782, 2002. doi: 10.2514/6.2002-2782

    [8]

    JAUCH C E, LONG D F, HERGERT D W. Performance assessment of a new type of heat exchanger for wind tunnels[R]. AIAA-2002-0740, 2002. doi: 10.2514/6.2002-740

    [9] 李启良,赵兰萍. 矩形翅片椭圆管热交换器流动和换热特性的数值模拟[J]. 流体机械,2006,34(8):67-70. DOI: 10.3969/j.issn.1005-0329.2006.08.019

    LI Q L,ZHAO L P. Numerical simulation of flow behaviors and heat transfer characteristic of rectangular finned ellipti-cal tube heat exchanger[J]. Fluid Machinery,2006,34(8):67-70. doi: 10.3969/j.issn.1005-0329.2006.08.019

    [10] 李启良,杨志刚,赵兰萍. 气动–声学风洞热交换器的气动噪声[J]. 同济大学学报(自然科学版),2012,40(8):1261-1264,1275. DOI: 10.3969/j.issn.0253-374x.2012.08.024

    LI Q L,YANG Z G,ZHAO L P. Aero-noise of heat exchanger in aero-acoustic wind tunnel[J]. Journal of Tongji University(Natural Science),2012,40(8):1261-1264,1275. doi: 10.3969/j.issn.0253-374x.2012.08.024

    [11] 符澄,赵波,徐大川,等. 板翅式及管翅式换热器气流湍流特性研究[J]. 实验流体力学,2019,33(6):22-27. DOI: 10.11729/syltlx20190036

    FU C,ZHAO B,XU D C,et al. Investigation on flow turbu-lent characteristics of plate-fin and tube-fin heat exchanger[J]. Journal of Experiments in Fluid Mechanics,2019,33(6):22-27. doi: 10.11729/syltlx20190036

    [12] 赵波,符澄,裴海涛,等. 椭圆翅片管对风洞气流湍流特性影响研究[J]. 西北工业大学学报,2020,38(2):303-308. DOI: 10.1051/jnwpu/20203820303

    ZHAO B,FU C,PEI H T,et al. Study on the effect of oval tubes on airflow turbulence characteristics in wind tunnel[J]. Journal of Northwestern Polytechnical University,2020,38(2):303-308. doi: 10.1051/jnwpu/20203820303

    [13] 赵波,廖达雄,黄知龙,等. 翅片椭圆管偏斜角对风洞换热器性能的影响[J]. 暖通空调,2020,50(11):120-124,112.

    ZHAO B,LIAO D X,HUANG Z L,et al. Effect of deflection angles on performances of finned oval tube heat exchanger applied to wind tunnel[J]. Journal of Heating Ventilating & Air Conditioning,2020,50(11):120-124,112.

    [14] 廖达雄,陈吉明,郑娟,等. 0.6 m连续式跨声速风洞总体性能[J]. 实验流体力学,2018,32(6):88-93. DOI: 10.11729/syltlx20170086

    LIAO D X,CHEN J M,ZHENG J,et al. General perfor-mance of 0.6 m continuous transonic wind tunnel[J]. Journal of Experiments in Fluid Mechanics,2018,32(6):88-93. doi: 10.11729/syltlx20170086

    [15] 赵波,陈振华,李为民,等. 风洞板翅式热交换器气体流动特性试验[J]. 实验流体力学,2022,36(1):69-76. DOI: 10.11729/syltlx20210081

    ZHAO B,CHEN Z H,LI W M,et al. Flow characteristics experiment of plate-fin heat exchanger in wind tunnel[J]. Journal of Experiments in Fluid Mechanics,2022,36(1):69-76. doi: 10.11729/syltlx20210081

  • 期刊类型引用(6)

    1. 戴梧叶,孙泓朴,吴宁宁,许灵芝. 高超声速边界层转捩研究进展与挑战. 航空动力学报. 2025(01): 180-195 . 百度学术
    2. 易淼荣,张若凌,岳茂雄,李莉,任虎,赵慧勇. 双垂直楔交叉激波与转捩边界层干扰. 实验流体力学. 2024(05): 17-28 . 本站查看
    3. 万兵兵,陈曦,陈坚强,袁先旭,段茂昌,黄章峰,涂国华. 三维边界层转捩预测HyTEN软件在高超声速典型标模中的应用. 空天技术. 2023(01): 150-158 . 百度学术
    4. 黄章峰,张宇琦. 高超声速三维边界层转捩数值研究进展及预测软件. 空气动力学学报. 2023(11): 1-19 . 百度学术
    5. 杜磊,孙波,代春良,卓长飞. 壁面温度对斜爆震发动机进气道流动特性影响的数值研究. 推进技术. 2021(04): 950-960 . 百度学术
    6. 易淼荣,赵慧勇,乐嘉陵,肖保国,郑忠华. 基于IDDES框架的γ-Re_θ转捩模型. 航空学报. 2019(08): 51-65 . 百度学术

    其他类型引用(2)

图(10)
计量
  • 文章访问数:  381
  • HTML全文浏览量:  208
  • PDF下载量:  74
  • 被引次数: 8
出版历程
  • 收稿日期:  2021-08-01
  • 修回日期:  2021-12-01
  • 录用日期:  2021-12-05
  • 刊出日期:  2022-09-30

目录

/

返回文章
返回
x 关闭 永久关闭