碳基材料氧化烧蚀的快/慢反应动力学问题

国义军, 周述光, 曾磊, 刘骁, 代光月, 邱波

国义军, 周述光, 曾磊, 刘骁, 代光月, 邱波. 碳基材料氧化烧蚀的快/慢反应动力学问题[J]. 实验流体力学, 2019, 33(2): 17-22, 101. DOI: 10.11729/syltlx20180128
引用本文: 国义军, 周述光, 曾磊, 刘骁, 代光月, 邱波. 碳基材料氧化烧蚀的快/慢反应动力学问题[J]. 实验流体力学, 2019, 33(2): 17-22, 101. DOI: 10.11729/syltlx20180128
Guo Yijun, Zhou Shuguang, Zeng Lei, Liu Xiao, Dai Guangyue, Qiu Bo. Study on fast/slow reaction mechanism of carbon-based material oxidation in high speed stream of dissociated air[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(2): 17-22, 101. DOI: 10.11729/syltlx20180128
Citation: Guo Yijun, Zhou Shuguang, Zeng Lei, Liu Xiao, Dai Guangyue, Qiu Bo. Study on fast/slow reaction mechanism of carbon-based material oxidation in high speed stream of dissociated air[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(2): 17-22, 101. DOI: 10.11729/syltlx20180128

碳基材料氧化烧蚀的快/慢反应动力学问题

基金项目: 

国家重点基础研究发展计划项目 2014CB744100

国家重点研究发展计划"大科学装置前沿研究"重点专项 2016YFA0401200

详细信息
    作者简介:

    国义军(1966-), 男, 山东烟台人, 博士, 研究员。研究方向:高超声速飞行器气动热环境, 烧蚀热防护。通信地址:四川省绵阳市二环路南段6号(621000)。E-mail:gyj2236985@sina.com

    通讯作者:

    国义军, E-mail:gyj2236985@sina.com

  • 中图分类号: TP601

Study on fast/slow reaction mechanism of carbon-based material oxidation in high speed stream of dissociated air

  • 摘要: 在前期单/双平台问题研究的基础上[1],就学术上和工程应用中关注的碳在空气中燃烧的"快反应"和"慢反应"问题进行了深入分析,研究发现被广泛使用50余年的"慢反应"并不存在,而被弃置的"快反应"则真实存在,且具有重要的应用价值。采用"快反应"动力学数据,同时考虑CO、CO2两个烧蚀产物,得到的无量纲烧蚀速率随温度的变化曲线存在两个平台,其中,温度稍低情况下出现的第一平台对应的主要烧蚀产物为CO2,温度稍高情况下出现的第二平台对应的主要烧蚀产物为CO,且第一平台对应的无量纲烧蚀速率恰好是第二平台的1/2。过去常被忽略的CO2扮演了重要角色,由它产生的第一平台,将以往文献中看似完全独立、毫无关系的"快反应"和"慢反应"曲线建立了联系。理论分析表明:第一平台之前的快速上升段属于从速率控制区到扩散控制区的过渡区,第一平台及其以后的区域都属于扩散控制区(包括两个平台之间的连接线),它是由反应生成物CO与CO2的分压比δ从0到∞的变化引起的,与表面化学反应动力学条件完全无关。由"双平台"理论得到的从低温到高温、由速率控制区经由过渡区到达扩散控制区的整条烧蚀速率曲线,与实验结果完全吻合。
    Abstract: Based upon our previous research on the single/dual platforms theory for carbon-based material oxidation, it has been found that the generally accepted "slow" reaction model which was presented by Scala in 1962 and has been widely used till now for over fifty years does not exist, while the deserted "fast" reaction model really exists and is of great value in application. Theoretical analysis shows that the "slow" reaction appears only at the boundary-layer-diffusion-rate controlled regime, In this regime, the ablation rate has no business with the chemical reaction rates, so it is a fictitious reaction model and has no physical meanings, although the results from it may be in agreement well with test results. The products of the heterogeneous reactions of carbon with oxygen are CO and CO2, but CO2 is usually neglected when the surface temperature is higher than 1000K in many literatures. However, from theoretical analysis, author find that CO2 plays an important role in spanning the reaction range from the reaction-rate controlled regime to the diffusion controlled oxidation regime, and can't be ignored. It has been shown that for the ablation of carbon-based materials in a high speed stream of dissociated air, the "fast" reaction must be used with both products of CO2 and CO at the same time. There are two different platforms in the diffusion controlled regime. The first platform results from reactions of the predominant CO2 product, and the other is due to the predominant CO product reactions. With the increases of the surface temperature, the ratio of the mass fraction of CO to CO2 at the surface rises rapidly from zero to infinity, which causes the oxidation process to change automatically from the nominated "fast" reaction to the so called "slow" reaction lines. The dual platform theory has been confirmed by several experimental results.
  • 流体壁面剪应力的精确测量是流体力学领域的难题[1]。随着MEMS技术的发展,硅基[2-3]和柔性聚合物基[4]的壁面剪应力微传感器研究发展受到重视,而传感器标定是其测量的基础和前提。根据不同应用需求初步形成了不同类型的动、静态壁面剪应力标定方法和装置[5-9]。本文介绍了平板流动、管道流动和旋转流动等静态标定方法,以及平面波振动、阶梯圆盘旋转运动等动态标定方法,介绍了工作原理和典型结构装置,可为壁面剪应力传感器标定方法的选定提供参考。

    壁面剪应力传感器的静态标定需要施加稳定已知量值的剪应力,读取对应输出。采用流体力学经典模型,通过测量其他物理量从而换算得到对应壁面剪应力值是常用方法。由于之前长期缺乏有效的流体壁面剪应力传感器,相关标定研究不系统,尚未形成统一的标准[10-11]

    目前,常见标定方法主要有平板模型标定法、管道流动标定法[12]和旋转运动标定法[13]

    早在20世纪上叶就有学者和机构对平板模型边界层进行了大量理论和实验研究,并得到平板表面层流和湍流边界层壁面剪应力的理论计算模型。

    在空气中,均匀来流绕过一个厚度很薄的平板时,如果平板迎角为0,可以应用普朗特边界层微分方程解决粘性流动问题,平板边界层如图 1所示[14]

    图  1  平板边界层
    Fig.  1  The boundary layer along a flat plate

    理论上,由于平板绕流压强顺流不变,即,不会发生边界层分离,也没有压差阻力,只有作用于壁面的摩擦阻力。故在这种情况下壁面切应力可由下式计算:

    (1)

    式中:μ表示流体的动力粘度,ν表示流体的运动粘度,U表示平板前端均匀来流流速,x表示距平板前缘的距离。

    随着来流速度增大,雷诺数也将增大。当平板表面雷诺数大于临界雷诺数,绕流场可能由层流状态转变为湍流[15],如图 2所示。

    图  2  平板边界层转捩
    Fig.  2  Laminar-turbulent transition on a flat plate

    通过动量积分方程可以得到湍流边界层中壁面剪应力的计算公式:

    (2)

    式中:xt为湍流边界层前沿到测量点的距离,ρ为空气密度,δ2为动量厚度,可通过下式计算得到:

    (3)

    通过确定测量点到转捩位置的距离xt,就可以得到来流速度U与剪应力τ的函数关系。

    由于平板表面发生转捩的位置不易测得,常常采用强制转捩的方法获得全湍流状态的平板边界层。Tani对转捩带粗糙体高度进行了研究,并得出粗糙体的极限高度计算公式[16]

    (4)

    式中:kcrit为粗糙体的极限高度。

    平板模型标定法能产生数百帕量级的壁面剪应力,标定精度受风洞 (或水道) 来流速度测量误差、测试点距前缘距离测量误差及环境温度测量误差等因素的影响。然而,平板边界层的流场状态不易控制,可结合实验方法 (如升华法、油流法等) 初步判断平板表面流动转捩与分离情况。

    管道流动因其运动受到管壁约束,管道内的流动参数 (如流量、温度和壁面静压等) 较易测量,也是常见的剪应力标定方法。

    当方形管道的宽高比较大时[17],入口处来流稳定的通道内流体运动可近似为二维流动。管道内流体的动量平衡与壁面剪应力作用相关[18]。通过测量管壁沿程压力损失,可以得到壁面剪应力的大小,标定装置如图 3所示。

    图  3  动量损失剪应力传感器标定系统
    Fig.  3  The calibration system of wall shear stress sensors

    定义pinlet/poutlet为管道内入口处压力和出口处压力的比值。对于压力比pinlet/poutlet≤1.1,可忽略管道内流体的压缩性,完全发展的不可压缩层流管道内的压力分布是线性的,剪应力表达式可以简单表示为:

    (5)

    式中:h为管道高度,表示管道内沿流场方向的压力梯度。通过测量Δx距离沿程的压力变化Δp,就可以得到剪应力的值。

    图  4  管道流动剪应力发生装置
    Fig.  4  Wall shear stress calibration device

    马赫数大于0.3的流动问题要考虑流体压缩性。对于压力比pinlet/poutlet≥1.1,管道内压力分布不再是线性变化[19]。通过非线性压力Π=p(x)/poutlet与非线性距离X之间的关系,代入公式 (5) 可以得到大压力比下的管道剪应力计算公式:

    (6)

    式中:Lc为管道的长度,poutlet为管道出口处的压力,X为沿管道方向x位置与管道长度Lc的比值。常数ab能够从Π2=aX+b的曲线计算得到[20]

    1997年,美国麻省理工学院研发了一种基于动量损失原理的空气介质剪应力标定装置[21],量程0~1Pa,如图 5所示。在微型气槽内壁面开凿压力检测孔,由于流体在槽腔运动产生摩擦阻力引起流体动量的损失,导致沿气流运动方向有着不同的壁面压力分布,换算可得压差中点剪应力值。此外,还考虑了空气压缩性,给出了不同压差范围的剪应力计算方法并用此装置对自主研发的浮动元件式剪应力传感器进行了标定。

    图  5  空气介质的剪应力标定装置示意图
    Fig.  5  Air static wall shear stress calibration device

    2002年,美国凯斯西储大学Patel等在研究壁面剪应力传感器时采用了此原理搭建起一套小型的扁平管道标定装置,并对浮动式剪应力传感器进行了标定,标定最大剪应力达到80Pa[18]

    2015年,中国空气动力研究与发展中心李建强等研究建立了管道流动动量损失原理的空气介质剪应力标定装置,增设了精密温控系统,可实现从环境温度到55℃的温度控制[22]

    2014年,中国船舶科学研究中心建立了水介质的动量损失管流标定装置,并利用LDV、PIV等对MEMS热膜剪应力传感器进行了测试标定[7]

    流体在通过极低雷诺数的微型扁薄矩形水道 (具有较大宽高比w/h),可获得可控的槽内纯剪切流[23],如图 6所示。

    图  6  扁形水槽内流体纯剪切流示意图
    Fig.  6  The boundary layer in a thin tunnel

    当管道内的流场形态为均匀层流时,其速度剖面曲线为二次曲线,管道体积流量Q与流动中心处最大流速相关,根据二维定常流动流体壁面剪应力的定义式[24],可以导出管道内壁面剪应力与体积流量之间的关系:

    (7)

    因矩形通道横截面宽度w无法实现远大于截面高度h,在剪应力计算时引入矩形通道的形状修正因子,则修正后剪应力计算公式为[23]

    (8)

    式中:Dh为水力学半径,

    Φ(n) 为矩形通道的修正因子,该修正因子的引入可降低扁薄矩形水道侧壁对中心流速的影响效应[17],

    (9)

    2005年,美国加州理工学院的Xu Y.等人根据该流动原理搭建起适用于MEMS剪应力传感器标定的微型管道测试平台,管道高度仅0.3mm,并实现对传感器0~5Pa范围内的剪应力标定[23]

    2009年,中国西北工业大学依此原理建立了微型扁薄标定水槽。2014年,李雁冰等根据该装置原理,对恒流驱动的热膜剪应力传感器进行了水下标定实验[25],并给出标定模型。王雷涛、严宇超等增加了溢流系统和温度控制系统,并对系统误差进行了分析。

    图  7  微型扁薄水槽剪应力发生装置
    Fig.  7  The micro calibration tunnel of wall shear stress

    平板模型和管道流动标定装置均需要有气源或者水压动力源作为介质流动的动力,尤其搭建专门风洞较为困难,因此有研究者采用圆盘或者圆筒的旋转运动带动周边流场流动从而形成一定的壁面剪切应力。

    2004年,美国国家航空航天局 (NASA) Zuckerwar等采用旋转圆轮建立了标定装置[26],如图 8所示。1个电机速度可控的圆形旋转轮和1个与其同轴心的固定圆弧槽,传感器平齐安装在槽上且靠近轮的边缘。当旋转圆轮与传感器之间的间隙h极小时,圆轮与圆弧槽间的流体形成简单库艾特流动,产生的剪应力为:

    图  8  旋转圆轮标定法示意图
    Fig.  8  Rotating wheel wall shear stress calibration device

    (10)

    式中:r为轮子半径,ω为电机带动的旋转轮的转速,通过调节转速的大小可以获得需要的剪应力。

    Zuckerwar通过对圆轮的精确安装,控制圆轮和同轴圆弧槽之间的缝隙最小达50μm,实现最大发生剪应力达到234.49Pa。

    1992年,美国麻省理工学院的Shajii J[27]等采用锥板式[28-30]的剪应力标定装置对浮动式剪应力传感器标定进行过相关研究,如图 9(a)所示。2006年,韩国浦项科技大学机械工程系的Kim[31]等采用2个平行放置的圆盘建立了壁面剪应力标定装置,如图 9(b)所示,并对热膜剪应力传感器实现0~0.8Pa范围内的标定。另外,美国加州理工学院的Brown G. L.[32]等采用圆筒式粘度计[33-34]的构造,如图 9(c)所示,搭建了剪应力标定装置。2012年,日本名古屋大学的Terashima O.[35]等对圆筒剪应力发生装置进行了改进,将外圆筒改为平板,如图 9(d)所示,更便于传感器的安装。

    图  9  (a) 锥板式标定装置原理图; (b) 静态标定的圆盘示意图; (c) 圆筒式旋转黏度计测量原理; (d) 旋转圆筒-平板标定法示意图
    Fig.  9  (a) Cone-plate wall shear stress calibration device. (b) Disk static calibration device.(c) Coaxial cylinder viscometers measurement device. (d) Rotating cylinder-plat wall shear stress calibration device

    旋转运动剪应力标定装置的搭建较为简单,但要保证产生剪应力的精度,对旋转装置偏心度要求较高,且若要产生较大范围的剪应力,圆盘或圆筒的直径和旋转速度要求很大。通常,此类方法仅用于0~2Pa范围内的剪应力标定。此外,该标定方法流场状态复杂,主要是利用其平均效应。

    动态剪应力标定需要加载一个周期性变化的剪应力。通常动态标定是为了研究传感器在中频或者高频激励下的响应情况,因此需要动态标定装置能产生较高频率的激励信号。

    平波管实现动态测量原理为用扬声器产生声波,声波驱动流体产生动态剪应力信号,利用麦克风测量产生信号的频率实现对剪应力传感器的动态测量。管道尾部加入消声装置以消除回声对测量的影响。

    2000年,美国佛罗里达大学的Chandrasekaran等建立了以平面声波驱动的动态剪应力标定装置,如图 10所示,该装置可产生最大0.1Pa,频率范围100Hz到20kHz的动态剪应力信号,如图 11所示。利用该装置对带有隔热空腔的硅基热敏传感器进行了动态测试[9]

    图  10  佛罗里达大学动态标定装置示意图
    Fig.  10  Dynamic wall shear stress calibration device designed by university of Florida
    图  11  佛罗里达大学动态标定装置在恒定振幅激励下的频率响应
    Fig.  11  Plot of the frequency response of the PWT for constant voltage amplitude of speaker excitation

    2001年,美国麻省理工学院的Sheplak等采用相同原理的平波管动态标定装置,并对浮动式剪应力传感器动态性能进行了标定测试[11]

    平面声波将引起振荡的压力梯度,通过测量压力波动,可以换算得到壁面动态变化的剪应力[11]

    (11)

    式中:p′为压力波动的幅值。

    2006年,韩国浦项科技大学机械工程系的Kim等在研究静态圆盘旋转标定法的同时对转动圆盘进行了改造,采用阶梯式转动圆盘,得到最高激励频率400Hz动态变化的剪应力,其结构如图 12所示。阶梯厚度为δh,2个圆盘的间隙的改变为。在每一个环形库艾特流周期内,圆盘上的m个阶梯提供了m个剪切脉冲。当上部圆盘的转速为n r/min,旋转产生的阶梯脉冲频率为mn/60Hz[31]

    图  12  动态标定圆盘示意图
    Fig.  12  Disk dynamic calibration device

    介绍了国内外流体壁面剪应力标定方法及其原理和装置。通过对比分析,可以得到以下几点启示:

    (1) 采用平板模型标定方法可以依托现成风洞或者水道,可以获得较大发生范围的剪应力 (从0帕到上百帕),实验系统可以适配多种类型的剪应力传感器;

    (2) 微型管道流槽标定是一种紧凑便捷的专用标定装置,截面尺寸及雷诺数很小,水槽高度经常在毫米级以下 (百微米量级)。发生的剪应力范围通常在10Pa以内;

    (3) 动量损失剪应力标定装置也属于专用定制,使用比较方便,剪应力发生范围极大 (可达上百帕);

    (4) 旋转运动标定方法的装置较易搭建,对机械装置及装配等精度要求较高,发生的剪应力范围较小,其流场状态复杂,主要取用平均效应;

    (5) 动态标定方法中研究及使用较多的是平波管法。

    综上,可根据各类流体壁面剪应力传感器的工作要求、尺寸、安装条件 (贴附面为平面或曲面、安装高度等) 和应用环境的不同,选用不同的标定方法或装置。

    目前,不同标定方法流场状态及品质影响,误差和不确定度及精度分析等研究较少,尚需关注。

  • 图  1   实验条件烧蚀组元浓度计算结果

    Fig.  1   Dependence of concentration on surface temperature

    图  2   ATJ石墨烧蚀计算与实验结果比较

    Fig.  2   Comparison of the ATJ graphite oxidation rate between calculated results and experimental data

    图  3   无量纲烧蚀质量流率随温度变化

    Fig.  3   Dependence of un-dimensional ablation rate on surface temperature

    表  1   实验条件下氧化分区所对应的温度范围[9]

    Table  1   The surface temperature corresponds to three oxidation regimes[9]

    Tw/K 氧化速率控制区 过渡区 扩散控制区
    600 700 750 770 800 900 1000 1025 1050 1100 1200 1500 1700 2000
    φch 3.3×105 1721 208 100 33 1.27 0.029 0.01 0.0036 5.3×10-4 1.86×10-5 5.2×10-9 0 0
    下载: 导出CSV
  • [1] 国义军, 代光月, 桂业伟, 等.碳基材料氧化烧蚀的双平台理论和反应控制机理[J].空气动力学学报, 2014, 32(6):755-760. DOI: 10.7638/kqdlxxb-2014.0047

    Guo Y J, Dai G Y, Gui Y W, et al. A dual platform theory for carbon based material oxidation with reaction diffusion rate controlled kinetics[J]. Acta Aerodynamica Sinica, 2014, 32(6):755-760. DOI: 10.7638/kqdlxxb-2014.0047

    [2]

    Scala S M. The ablation of graphite in dissociated air, Part Ⅰ: Theory[R]. AD289298, 1962.

    [3]

    Diaconis N S, Gorsuch P D, Sheridan R A. The ablation of graphite in dissociated air, Part Ⅱ: Experimental investigation[R]. AD290051, 1962.

    [4]

    Maahs H G. Oxidation of carbon at high temperatures: reaction-rate control or transport control[R]. NASA TN D-6310, 1971.

    [5]

    Swann R T, Pittman C M, Smith J C. One-dimensional numerical analysis of the transient response of thermal protection systems[R]. NASA TN D-2976, 1965.

    [6] 黄振中.烧蚀端头的瞬变外形及内部温度分布[J].空气动力学学报, 1981(1):53-65. http://www.cnki.com.cn/Article/CJFDTOTAL-KQDX198101005.htm
    [7] 卞荫贵, 钟家康.高温边界层传热[M].北京:科学出版社, 1986.
    [8] 国义军.炭化材料烧蚀防热的理论分析与工程应用[J].空气动力学学报, 1994, 12(1):94-99. http://www.cnki.com.cn/Article/CJFDTOTAL-KQDX401.014.htm

    Guo Y J. An analysis of a charring ablative thermal protection system with its engineering application[J]. Acta Aerodynamica Sinica, 1994, 12(1):94-99. http://www.cnki.com.cn/Article/CJFDTOTAL-KQDX401.014.htm

    [9]

    Milos F S, Chen Y K. Comprehensive model for multicomponent ablation thermochemistry[R]. AIAA-97-0141, 1997.

    [10] 黄海明, 杜善义, 吴林志, 等. C/C复合材料烧蚀性能分析[J].复合材料学报, 2001, 18(3):76-80. DOI: 10.3321/j.issn:1000-3851.2001.03.018

    Huang H M, Du S Y, Wu L Z, et al. Analysis of the ablation of C/C composites[J]. Acta Materiae Compositae Sinica, 2001, 18(3):76-80. DOI: 10.3321/j.issn:1000-3851.2001.03.018

    [11] 张志成, 潘梅林, 刘初平.高超声速气动热和热防护[M].北京:国防工业出版社, 2003.
    [12]

    Welsh W E, Chung P M. A modified theory for the effect of surface temperature on the combustion rate of carbon surface in air[C]//Proc of Heat Transfer and Fluid Mechanics Institute.: Stanford University Press, 1963: 146-159.

    [13]

    Baron J R, Bernstein H. Heterogeneous rate coupling for graphite oxidation[J]. AIAA Journal, 1971, 9(8):1588-1595. DOI: 10.2514/3.49961

    [14]

    Milos F S. Comparison of ablation predictions for carbonaceous materials using CEA and JANAF-based species thermodynamics[C]. Proc of the 35th Annual Conference on Composites, Materials, and Structures. 2011.

    [15] 益小苏, 杜善义, 张立同.复合材料手册[M].北京:化学工业出版社, 2009.
    [16]

    Clark R K. An analysis of a charring ablator with thermal nonequilibrium, chemical kinetics, and mass transfer[R]. NASA TN D-7180, 1973.

  • 期刊类型引用(5)

    1. 陈纪仲,田海平,丁俊飞,马国祯,陈怡纯. 壁湍流摩阻时空信号无接触测量实验研究. 实验流体力学. 2024(04): 66-74 . 本站查看
    2. 高南,刘玄鹤. 实用化壁面切应力测量技术的综述与展望. 空气动力学学报. 2023(03): 1-24 . 百度学术
    3. 刘丹. 免标定壁面剪应力测量技术研究. 科技创新与应用. 2023(34): 188-192 . 百度学术
    4. 陆连山,李栋. 利用双参数非线性拟合求解湍流边界层壁面摩擦速度. 西北工业大学学报. 2020(02): 253-260 . 百度学术
    5. Xuanhe Liu,Zhuoyue Li,Nan Gao. An improved wall shear stress measurement technique using sandwiched hot-film sensors. Theoretical & Applied Mechanics Letters. 2018(02): 137-141 . 必应学术

    其他类型引用(4)

图(3)  /  表(1)
计量
  • 文章访问数:  287
  • HTML全文浏览量:  131
  • PDF下载量:  12
  • 被引次数: 9
出版历程
  • 收稿日期:  2018-09-13
  • 修回日期:  2018-10-24
  • 刊出日期:  2019-04-24

目录

/

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

重要公告

www.syltlx.com是《实验流体力学》期刊唯一官方网站,其他皆为仿冒。请注意识别。

《实验流体力学》期刊不收取任何费用。如有组织或个人以我刊名义向作者、读者收取费用,皆为假冒。

相关真实信息均印刷于《实验流体力学》纸刊。如有任何疑问,请先行致电编辑部咨询并确认,以避免损失。编辑部电话0816-2463376,2463374,2463373。

请广大读者、作者相互转告,广为宣传!

感谢大家对《实验流体力学》的支持与厚爱,欢迎继续关注我刊!


《实验流体力学》编辑部

2021年8月13日