轴对称层流到湍流转捩区间速度和温度的振荡统计特性研究

张若凌, 乐嘉陵

张若凌, 乐嘉陵. 轴对称层流到湍流转捩区间速度和温度的振荡统计特性研究[J]. 实验流体力学, 2019, 33(2): 30-42. DOI: 10.11729/syltlx20180117
引用本文: 张若凌, 乐嘉陵. 轴对称层流到湍流转捩区间速度和温度的振荡统计特性研究[J]. 实验流体力学, 2019, 33(2): 30-42. DOI: 10.11729/syltlx20180117
Zhang Ruoling, Le Jialing. The statistics of velocity and temperature fluctuations in axisymmetric laminar-to-turbulent transitions[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(2): 30-42. DOI: 10.11729/syltlx20180117
Citation: Zhang Ruoling, Le Jialing. The statistics of velocity and temperature fluctuations in axisymmetric laminar-to-turbulent transitions[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(2): 30-42. DOI: 10.11729/syltlx20180117

轴对称层流到湍流转捩区间速度和温度的振荡统计特性研究

基金项目: 

国家重点研发计划项目 2017YFB0202400

国家重点研发计划项目 2017YFB0202404

详细信息
    作者简介:

    张若凌(1973-), 男, 湖北钟祥人, 研究员。研究方向:超燃冲压发动机的主动冷却结构设计、计算流体力学和高超声速空气动力学。通信地址:四川省绵阳市二环路南段6号1902信箱(621000)。E-mail:zhangruoling@cardc.cn

    通讯作者:

    张若凌, E-mail:zhangruoling@cardc.cn

  • 中图分类号: V211.1

The statistics of velocity and temperature fluctuations in axisymmetric laminar-to-turbulent transitions

  • 摘要: 热力学平衡系统连续相变的理论方法,被推广用来讨论圆管内轴对称层流到湍流转捩区间的流动和振荡特性。假设在转捩区间径向脉动速度与充分发展区的湍流在数值上完全相同,在每一个截面上转捩流动可以看成是充分发展区的层流和湍流的合成流动。湍流成分的合成比例作为序参数用来定义合成流动。引入合成比例的振荡后,运用最小熵产生准则得到一个可以描述转捩行为的方程。采用相同的处理方法讨论了加热圆管内转捩区间的对流传热特性。在圆管内的流动和对流传热允许相似和独立的转捩过程,在转捩区间宏观振荡同时具有随机性和确定性。最后与实验进行了对比,包括流动和传热实验得到的测量结果。
    Abstract: The procedures for studying continuous phase transitions of thermodynamic equilibrium systems are extended to discuss the laminar-to-turbulent transitional flows in circular tubes. The flow in the transition range is treated as a composition of the laminar and turbulent flows assuming that the radial fluctuating velocity has the same value as that of the fully turbulent flow. The composite ratio of the turbulent flow is used as an order parameter to define the composite flow. The fluctuations of the composite ratios are introduced, and the criterion of minimum entropy production is used to derive an equation which can describe the transition behaviors. The convective heat transfer characteristics in the transition range in a heated circular tube are also discussed adopting the same procedures. Similar and separate processes for the transitions of the flow and convective heat transfer types are allowed in the heated circular tube. The macroscopic fluctuations in the transition range show both probabilistic and deterministic characteristics simultaneously. The agreements with measurements are given including those obtained in flow and heat transfer experiments.
  • 现代机载武器在外形设计、总体布置以及结构形式上除了必须满足机载武器的常规要求外,还要尽可能满足武器内埋以及多弹并列、串列挂载的要求,这些空间和尺寸上的条件约束为其升力面和控制面的设计带来了极大的困难和挑战,同时舵面操纵效能、舵面铰链力矩与小型化舵机的匹配性等矛盾也十分突出。栅格翼[1]作为一种非常规的升力面和控制面,在升力特性、铰链力矩特性和外形尺寸方面都优于传统平板翼,这些都为其在现代机载武器特别是小型化武器上的广泛应用奠定了基础。

    目前,栅格翼虽已成功应用于宇宙飞船救生逃逸系统和一些常规战术导弹上,国内外对其气动特性也开展了大量的实验与数值模拟研究,如Mark S. Miller等通过改变栅格翼边框和栅格壁前后缘形状等开展了减阻特性的试验研究[2];Gregg Abate等通过风洞实验研究了跨声速阶段栅格翼的气流壅塞现象对气动特性的影响[3];而陈少松等开展了栅格翼的减阻特性及格宽翼弦比对其气动特性的影响研究[4-5]等。数值模拟研究方面,Karl S. Orthner等采用非结构网格方法对栅格翼的跨声速气流壅塞现象进行了研究[6];James Despirito等对带栅格翼和平板翼鸭式布局导弹气动特性进行了数值计算,并对栅格尾翼的滚转控制性能进行了研究[7];刘刚等采用结构、非结构混合网格对不同舵偏角下的栅格翼构型进行了数值模拟[8];吴晓军等对多种不同外形单栅格的气动特性开展了数值计算[9]。以上工作广泛地开展了栅格翼相关气动特性的研究,但大多数研究集中于单栅格翼的气动特性或单一类型栅格翼及其尾身组合体的气动特性,对尾身组合体构型下不同类型栅格翼的气动特性研究尚显不足,而且国内有关栅格翼气动特性的试验研究开展得相对较少,这些都阻碍了对栅格翼气动特性的深入认识及其进一步的应用。

    栅格翼因外部框架和内部栅格布置的不同而分为不同的类型,如框架式、蜂窝式等,并且栅格翼的布局形式及外形参数对其气动特性有非常重要的影响[1]。本文针对栅格翼尾身组合体模型,采用常规测力天平,实现了尾身组合体不同类型栅格翼全量气动特性的试验研究,获得了简单框架式、正置蜂窝式以及不同网格密度的斜置蜂窝式4种栅格翼的亚声速和跨声速气动特性,并将其气动特性与平板翼进行了对比,同时研究了弹身迎风侧和背风侧不同位置对栅格翼气动特性的影响。

    图 1给出了试验的不同类型栅格翼尾身组合体模型,由旋成体弹身、×型尾翼等组成。尾翼安装平面与构造水平对称面夹角为45°。试验时,只对1片尾翼进行气动力测量,分别位于1#或3#方位(从尾部向前看分别位于弹体轴yz平面的1和3象限)。尾翼包括1种平板尾翼PF0和4种栅格尾翼(编号分别为GF0+、GF4+、GFX1、GFX2),4种栅格翼的表面积近似相等,翼元格的形状和数量不同,其中GF0+为简单框架式,GF4+为正置蜂窝式,GFX1和GFX2为网格密度不同的斜置蜂窝式栅格翼。试验中尾翼仅作俯仰偏转控制,4片尾翼同时偏转。各类型尾翼的详细结构及尺寸见图 2

    图  1  尾身组合体模型栅格翼测力方案示意图
    Fig.  1  Sketch of model with body and tail and force test of grid fin
    图  2  平板翼与栅格翼外形图
    Fig.  2  Sketch of planar and grid fins

    对于舵面气动特性的测量,试验一般采用专门的铰链力矩天平[10],栅格翼的突出特点是其铰链力矩较小,但同时会带来较大阻力。一般舵面的铰链力矩试验并不测量其阻力,结合本项研究的目的,采用常规的测力天平来测量栅格翼的气动特性,尾身组合体模型采用尾部支撑,待测栅格翼与天平连接,每次测试一个栅格翼的气动力,其余栅格翼起干扰翼的作用。该方法不需要增加额外的翼天平,使用一个天平即可实现不同位置栅格翼全量气动特性的测量,具有模型设计简单的优点,并克服了在较小尺寸模型上安装铰链力矩天平的困难,同时又能测量栅格翼气动力和力矩的全部分量。

    试验在试验段横截面尺寸为0.6 m×0.6 m的直流暂冲式跨超声速风洞中进行,试验马赫数Ma为0.8、1.2,模型迎角范围为-4°~14°,尾翼舵偏角为-10°、-5°、0°。

    本项研究共获得了1种平板翼和4种栅格翼6个分量的气动力特性数据,尾翼气动力和力矩系数分别在相应的舵面固连坐标系中给出。下面对尾身组合体不同类型栅格翼的主要气动特性进行分析。

    图 3给出了1#位平板翼和4种栅格翼在不同马赫数及俯仰舵偏角下轴向力系数随迎角的变化曲线。对于平板翼,轴向力随迎角及舵偏角的变化符合基本规律。δz=0°时,4种栅格翼轴向力在小迎角范围内随迎角变化较小。从不同类型尾翼轴向力的对比曲线可知,与平板翼相比,此时栅格翼表现出了其突出特点,轴向力较大[6],简单框架式栅格翼GF0+的轴向力在4种栅格翼中最小,无舵偏角时,正置蜂窝式GF4+与斜置蜂窝式GFX1的轴向力较为接近,随翼元格数增加,GFX2的轴向力最大。栅格翼的阻力主要取决于其结构和栅格管内的流动状态[7],4种栅格翼具有相同的浸润面积,所以阻力的差别主要来源于管内的复杂流动引起的压差阻力。随着舵偏角的增加,无论是亚声速还是跨声速,4种栅格翼的轴向力有逐渐接近的趋势,Ma=1.2,δz=-10°时,4种栅格翼轴向力的差别明显变小。

    图  3  1#位不同类型尾翼轴向力系数
    Fig.  3  Comparison of axis force coefficients for different fins at 1# location

    图 4是各马赫数和舵偏角下4种栅格翼的法向力系数随迎角的变化曲线。总的来说,对于不同类型栅格翼的法向力系数,在中等迎角以下,法向力随迎角基本呈线性变化,无舵偏角和小舵偏角时,GF0+和GF4+表现出了更好的法向力特性,当舵偏角较大时,2种斜置蜂窝式栅格翼法向力特性较优。4种栅格翼具有近似相等的升力面积,但自较小迎角开始,几种栅格翼的法向力已有较大差别,某些状态下4种栅格翼的法向力在较大迎角时均存在失速现象,且失速迎角均小于平板翼。对于GF0+,Ma=0.8,δz=-10°,在负迎角时,法向力随迎角变化不明显,这是因为此时当地有效迎角增大,栅格翼处于失速的状态,故法向力变化不大,从图中结果可以看出,在相同条件下,GF4+具有相类似的结果,而2种斜置蜂窝式栅格翼则不存在这种现象,可见栅格翼的结构形式对其法向力特性有较大影响。

    图  4  1#位不同类型尾翼法向力系数
    Fig.  4  Comparison of normal force coefficients for different fins at 1# location

    相对于平板翼,栅格翼的优点之一就是具有较好的铰链力矩特性[6],这一点从图 5的试验结果可以看出,栅格翼的铰链力矩特性普遍好于平板翼。与亚声速相比,跨声速时,各栅格翼的铰链力矩随迎角和舵偏角的变化范围较小。总的来说,对于4种栅格翼而言,斜置蜂窝式栅 格翼的铰链力矩特性优于另外2种栅格翼,密网格斜置蜂窝翼GFX2铰链力矩特性最优。

    图  5  1#位不同类型尾翼铰链力矩系数
    Fig.  5  Comparison of hinge moment coefficients for different fins at 1# location

    本文还研究了栅格翼处于尾身组合体迎风侧和背风侧不同位置时的气动特性。为了方便对比,不同位置栅格翼的气动力系数已经转换至统一方向的舵面坐标系中。图 678是1#位和3#位不同类型栅格翼气动特性的对比,亚声速时,对于GF4+,无舵偏角时,在绝大部分迎角下,处于背风侧1#位栅格翼的轴向力、法向力及铰链力矩均较大,迎风侧和背风侧栅格翼的铰链力矩随迎角的变化规律基本相同,随着舵偏角增加,法向力变化规律相似,但轴向力和铰链力矩在部分状态下发生明显变化,如δz=-5°时2个位置栅格翼的铰链力矩,δz=-10°时,大部分迎角下迎风侧3#位的轴向力较大;试验结果表明,对于其它类型的栅格翼,这种变化规律不具有必然性,如GF0+,当δz=-5°时,迎风侧3#位栅格翼的轴向力已经大于背风侧1#位的轴向力,并且不同舵偏角下2个位置栅格翼的铰链力矩相似性更好,量值更为接近;对于GFX2,亚声速时,不同位置时的轴向力变化规律与GF0+相似,但法向力差别较小;与亚声速时相比,跨声速时部分状态下2个位置栅格翼的轴向力和铰链力矩变化规律发生改变。

    图  6  不同位置栅格翼气动特性对比(Ma=0.8,GF4+)
    Fig.  6  Comparison of aerodynamic characteristics for grid fin at different locations
    图  7  不同位置栅格翼气动特性对比(Ma=0.8,GFX2)
    Fig.  7  Comparison of aerodynamic characteristics for grid fin at different locations
    图  8  不同位置栅格翼气动特性对比(Ma=1.2,GFX2)
    Fig.  8  Comparison of aerodynamic characteristics for grid fin at different locations

    对于尾身组合体构型,通过风洞试验研究了平板尾翼及不同类型栅格尾翼的气动特性,讨论了各类型栅格翼气动力的变化规律,对弹身迎风侧和背风侧栅格翼的气动力进行了对比分析,可以得到以下结论:

    (1) 采用了模型设计简单,同时又能测量栅格翼气动力和力矩全部分量的试验方法,该方法克服了在较小尺寸模型上安装铰链力矩天平的困难,不需要增加额外的翼天平,使用一个天平即可实现不同位置栅格翼全量气动特性的测量。

    (2) 对于不同类型的栅格翼,轴向力均大于平板翼,无舵偏角时,轴向力随迎角变化较小,4种栅格翼中,简单框架式GF0+的轴向力最小,随翼元格数增加,斜置蜂窝式栅格翼GFX2轴向力最大。

    (3) 某些状态下,4种栅格翼的法向力在较大迎角时均存在失速现象,且失速迎角均小于平板翼。对于栅格翼GF0+和GF4+,当舵偏角使栅格翼当地有效迎角增大时,栅格翼处于失速状态,法向力随迎角变化不大,而在相同条件下,2种斜置蜂窝式栅格翼则不存在这种现象,可见栅格翼的结构形式对其法向力有较大影响。

    (4) 4种栅格翼的铰链力矩特性普遍好于平板翼,密网格斜置蜂窝翼GFX2铰链力矩特性最优。与亚声速相比,跨声速时,各栅格翼的铰链力矩随迎角和舵偏角的变化范围较小。

    (5) 舵偏角对弹身迎风侧和背风侧栅格翼的轴向力和铰链力矩有较大影响,且不同类型栅格翼的变化规律也不尽相同;位置对密网格斜置蜂窝翼的法向力影响减弱。

  • 图  1   两相自由能的变化

    Fig.  1   Variation of free energy for two phases

    图  2   圆管内的流动示意

    Fig.  2   Schematic diagram of the flow in a circular tube

    图  3   加热圆管内正十烷的密度、粘性系数和热传导系数以及流动的Re

    Fig.  3   Density, viscosity, thermal conductivity of the fluid and Re in the heated tube

    图  4   在转捩区间里r=0处I的大的过冲

    Fig.  4   Large overshoot of I at r=0 in the transition region

    图  5   r=0处u的对比

    Fig.  5   The comparison of u at r=0

    图  6   r=0处u′的对比

    Fig.  6   The comparison of u′ at r=0

    图  7   电加热管外壁温度和Nu的对比[27]

    Fig.  7   Outer wall temperature and Nu comparisons in an electrically heated tube[27]

    表  1   二元混合物相变和圆管内流动转捩的对比

    Table  1   Analogies between the phase transition of a binary mixture and the flow transition

    系统 序参数 热力学函数 处理方法
    二元混合物(平衡系统) 浓度x或其线性函数 自由能,混合自由能 混合自由能取极小值来确定x(Landau理论的精神)
    一个截面上的流体微元(非平衡系统) 合成比例η或其线性函数 平均熵产生,平均合成熵产生 平均合成熵产生取极小值来确定η(利用了Landau理论的精神)
    下载: 导出CSV
  • [1]

    Mullin T. Experimental studies of transition to turbulence in a pipe[J]. Annual Review of Fluid Mechanics, 2011, 43:1-24. DOI: 10.1146/annurev-fluid-122109-160652

    [2]

    Durst F, Ünsal B. Forced laminar-to-turbulent transition of pipe flows[J]. J Fluid Mech, 2006, 560:449-464. DOI: 10.1017/S0022112006000528

    [3]

    Darbyshire A G, Mullin T. Transition to turbulence in constant-mass-flux pipe flow[J]. J Fluid Mech, 1995, 289:83-114. DOI: 10.1017/S0022112095001248

    [4]

    Eckhardt B, Schneider T M, Hof B, et al. Turbulence transition in pipe flow[J]. Annual Review of Fluid Mechanics, 2007, 39:447-468. DOI: 10.1146/annurev.fluid.39.050905.110308

    [5]

    Reichl L E. A modern course in statistical physics[M]. New Jersey:John Wiley & Sons Inc, 1998.

    [6]

    Landau L D, Lifshitz E M. Statistical physics:part 1[M]. 3rd ed. Oxford:Pergamon Press, 1980.

    [7]

    McComb W D. The physics of fluid turbulence[M]. Oxford:Claredon Press, 1992.

    [8]

    Nishi M, Ünsal B, Durst F, et al. Laminar-to-turbulent transition of pipe flows through puffs and slugs[J]. J Fluid Mech, 2008, 614:425-446. DOI: 10.1017/S0022112008003315

    [9]

    Henkel M, Hinrichsen H, Lübeck S. Non-equilibrium phase tran-sitions, Volume Ⅰ:Absorbing phase transitions[M]. Netherlands:Canopus Academic Publishing Limited, 2008.

    [10]

    Cowan B. Topics in statistical mechanics[M]. London:Imperial College Press, 2005.

    [11]

    Huang K. Statistical mechanics[M]. New Jersey:John Wiley & Sons Inc, 1987.

    [12]

    Landau L D, Lifshitz E M. Fluid mechanics[M]. 2nd ed. Oxford:Pergamon Press, 1987.

    [13]

    Eckert E R G, Drake R M Jr. Analysis of heat and mass transfer[M]. Tokyo:McGraw-Hill Kogakusha Ltd, 1972.

    [14]

    Koch H, Tataru D. Well-posedness for the Navier-Stokes equations[J]. Adv Math, 2001, 157(1):22-35. http://d.old.wanfangdata.com.cn/OAPaper/oai_arXiv.org_1302.5785

    [15]

    Durmagambetov A A, Fazilova L S. Navier-Stokes equations-millennium prize problems[J]. Natural Science, 2015, 7(2):54262. http://d.old.wanfangdata.com.cn/NSTLHY/NSTL_HYCC0215037970/

    [16]

    Lifshitz E M, Pitaevskii L P. Statistical physics:part 2[M]. Oxford:Pergamon Press, 1980.

    [17]

    Glansdorff P, Prigogine I. Structure, stability, and fluctuations[M]. New Jersey:Wiley-Interscience, 1971.

    [18]

    Nicolis P, Prigogine I. Self-organization in nonequilibrium systems[M]. New Jersey:Wiley-Interscience, 1977.

    [19]

    Prigogine I. Time, structure, and fluctuations[J]. Science, 1978, 201:777-785. DOI: 10.1126/science.201.4358.777

    [20]

    Bergman T L, Lavine A S, Incropera F P, et al. Fundamentals of heat and mass transfer[M]. 7th ed. New Jersey:John Wiley & Sons Inc, 2011.

    [21]

    Rohsenow W M, Hartnett J P, Cho Y I. Handbook of heat transfer[M]. 3rd ed. New York:McGraw-Hill Book Company, 1998.

    [22]

    Zhang R L, Le J L. Natural laminar-to-turbulent transition inside an electrically heated circular tube[C]//Proc of AIP Conference Proceedings 1770: 030035. 2016.

    [23] 张若凌, 乐嘉陵.电加热圆管内流动的自然转捩过程研究[J].实验流体力学, 2017, 31(2):51-60. http://www.syltlx.com/CN/abstract/abstract11011.shtml

    Zhang R L, Le J L. Natural laminar-to-turbulent transition inside an electrically heated circular tube[J]. Journal of Experiments in Fluid Mechanics, 2017, 31(2):51-60. http://www.syltlx.com/CN/abstract/abstract11011.shtml

    [24]

    Nishi M. Laminar to turbulent transition in pipe flow through puffs and slugs[D]. Der Technischen Fakultät der Friedrich-Alexander-Universität Erlangen-Nürnberg, 2009.

    [25]

    Linne D L, Meyer M L, Edwards T, et al. Evaluation of heat transfer and thermal stability of supercritical JP-7 fuel[R]. AIAA-97-3041, 1997.

    [26]

    Huang H, Sobel D R, Spadaccini L J. Endothermic heat-sink of hydrocarbon fuels for scramjet cooling[R]. AIAA-2002-3871, 2002.

    [27]

    Zhang L, Zhang R L, Xiao S D, et al. Researches on heat transfer correlations of hydrocarbon fuel under supercritical pressure[J]. International Journal of Heat and Mass Transfer, 2013, 64:393-400. DOI: 10.1016/j.ijheatmasstransfer.2013.04.058

  • 期刊类型引用(1)

    1. 田霖,樊文鹏,于建军,杜云鹏,郭光全,唐宏. 基于结构网格的栅格翼子弹翼身组合体气动特性. 科技风. 2023(08): 158-162 . 百度学术

    其他类型引用(1)

图(7)  /  表(1)
计量
  • 文章访问数:  243
  • HTML全文浏览量:  92
  • PDF下载量:  8
  • 被引次数: 2
出版历程
  • 收稿日期:  2018-08-23
  • 修回日期:  2018-09-17
  • 刊出日期:  2019-04-24

目录

/

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