CN108595746A - 一种空心砖蓄热式加热器的总体设计评估方法 - Google Patents
一种空心砖蓄热式加热器的总体设计评估方法 Download PDFInfo
- Publication number
- CN108595746A CN108595746A CN201810189272.9A CN201810189272A CN108595746A CN 108595746 A CN108595746 A CN 108595746A CN 201810189272 A CN201810189272 A CN 201810189272A CN 108595746 A CN108595746 A CN 108595746A
- Authority
- CN
- China
- Prior art keywords
- heater
- thermal
- storage unit
- heat
- temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一套空心砖型蓄热式加热器设计与评估方法,包括:(一)加热器蓄热单元几何设计;(二)蓄热单元性能分析方法,包括蓄热单元压力损失分析方法和热应力分析方法;(三)加热器隔热层设计与分析方法;(四)加热器初步预热设计与分析方法。根据给定加热器的基准设计状态要求,采用本发明方法可以进行加热器方案几何设计,包括确定加热器的总高和最大直径、蓄热材料的高度和直径、隔热层的厚度等其它几何参数;采用本方法可以进一步评估蓄热式加热器的工作过程中的气动参数,即可以计算得到基准条件下的蓄热单元气流质量通量特性、基准运行状态下蓄热单元当地最大温差的沿程变化和预热过程蓄热阵温度径向分布等。
Description
技术领域
本发明涉及超燃冲压发动机地面实验技术设计领域,具体地说,涉及一种纯净空气空心砖蓄热式加热的超燃冲压发动机燃烧室实验系统设计评估方法。
背景技术
目前,国外蓄热式加热器的研制与应用主要集中在美国(如NASA、ASE公司和AEDC)、法国和日本(NAL)等。迄今为止,国外先后研制了30余套蓄热式加热系统,其中比较典型的有美国AEDC的9#风洞、NASA格林研究中心的高超声速风洞设备HTF、日本防卫厅研究开发局的RJTF、法国空间局的S4风洞等。目前,国内还没有用于高超声速风洞的空心砖蓄热式加热器。1998年大连北岛能源技术有限公司在国内首次实现了蓄热式技术燃烧高炉煤气在连续式轧钢加热炉上的应用,其蓄热体采用陶瓷小球。国内对蓄热式加热器的理论研究也有一些进展,黄河激等人对氧化铝多孔砖蓄热体的换热能力进行了简单估算,定义了蓄热体浮起危险系数(第三届高超声速科技学术会议,CSTAM-2010-0054,2010)。李龙飞等人提出蓄热阵的设计包括材料选择、蓄热阵几何设计、热性能和热应力分析。基于加热器使用温度选择适当的蓄热材料,以换热性能为目标,以蓄热量为标准,以热应力破坏为限制因素进行蓄热阵几何设计,难点和关键在于热应力分析(火箭推进38.2(2012):16-19)。罗飞腾等人针对空心砖型蓄热式加热器,开展了空心砖型蓄热单元初步设计研究,涉及换热性能、压降控制和热应力评估等方面(试验流体力学,2013,27(1):38-45)。截至目前,国内还没有一种完整、详细的空心砖型蓄热式加热器的总体设计评估方法。
发明内容
为了更好的开展空心砖蓄热式加热器的研制工作,本发明提出了出一套完整的空心砖型蓄热式加热器设计与评估方法,包括:(一)加热器蓄热单元几何设计;(二)蓄热单元性能分析方法,包括蓄热单元压力损失分析方法和热应力分析方法;(三)加热器隔热层设计与分析方法;(四)加热器初步预热设计与分析方法。
(一)蓄热单元几何设计
根据蓄热式加热器蓄热单元(图1)的特点(圆柱形空心砖蓄热单元、圆柱形通孔、通孔呈等边三角形分布等),设计步骤如下:
确定初始几何设计变量如下:蓄热床的外径Z与截面积A,气流通孔直径d,孔间距s,孔间距与孔径之比s/d,截面空隙率η,当量流体通径D,通孔个数N。其中d取值范围为5~10mm,s/d取值范围为1.4~2.4。
根据加热器的规模确定出蓄热单元截面流通面积A0,则当量流体通径蓄热体截面上气流通孔的个数
根据蓄热体截面上的等边三角形通孔分布形式,蓄热体截面的空隙率比表面积蓄热床的外径横截面积
通过上述步骤即可得到加热器全部几何参数。
(二)蓄热单元性能分析方法
1.蓄热单元压力损失分析方法
气流经过加热器蓄热单元的压力损失主要有局部损失和沿程损失两个方面。沿程摩擦压力损失的估算可按下式计算:
其中Cf为摩擦系数,与雷诺数相关;P为压力;为通孔内气流平均温度;R为气体常数;d为通孔内径;L为蓄热床通孔长度。局部压力损失可用下面实验经验公式计算:
2.蓄热单元热应力分析方法
蓄热阵内的圆柱型气流通道在截面上呈等边三角形均匀分布(如图3所示),这种几何构型由纯径向温度梯度的圆管模型来近似,适用于孔间距与孔径比值s/d>1.2的空心砖。加热器运行过程最大可允许的蓄热单元温度梯度和冷却率可由如下公式得到:
其中,a表示与蓄热材料物性参数相关的热扩散系数;Ia1和Ia3是与空心砖蓄热单元几何尺寸相关的参数,如下:
ΔTlim为Tm-T0,Tm为壁面厚度r0至ra范围内的平均温度,其定义如下:
T是指每一单孔截面上径向壁面温度T(r),在表面处To=T(ro)。
(三)隔热层设计
采用圆柱形隔热分析模型,图4为加热器圆柱段双层隔热的几何模型简图。采用圆筒壁的一维导热简化计算方法。蓄热体外表面到内层隔热层的传热:
内层隔热层外表面到外层隔热层的传热:
外壳与周围空气的大空间自然对流换热:
q3=hπD3(T3-Ta)
根据热平衡有
q1=q2=q3
其中qi为热流量;λ1和λ2为内层和外层隔热材料的热传导系数;T1为蓄热体外表面温度;T2为内层隔热层的外表面温度,也就是外层隔热层的内表面温度;T3为外层隔热层的外表面温度,假设等于外壳表面温度400K;Ta为周围环境温度并假定为288K;D1、D2、D3分别为蓄热体、内层隔热层和外层隔热层的外径;h为外壳与周围空气的自然对流换热系数。
由δi=(Di+1-Di)/2可得到每个隔热层厚度δ。求解双层隔热厚度时,采用固定内层厚度的方法进行设计,即首先设定内层厚度,然后利用单层厚度的优化设计方法对各物性参数进行计算,以内层外表面温度不大于外层保温材料允许使用温度作为边界条件,反复迭代,重新设计内层厚度和总厚度。
(四)加热器初步预热设计与分析
对于加热器预热分析,为了简化需要做出如下几点假设:(1)温度恒定的高温燃气从蓄热阵顶部进入,并均匀分散流过所有通孔;(2)当地流体温度与蓄热床温度相等,在蓄热阵横截面方向导热系数无穷大;(3)由于当地温度与周围环境温度的温差驱动,热量从蓄热床外围扩散损失;(4)忽略纵向热传导(即蓄热床长度方向的热传导)。
对于圆柱形壁面,当地隔热层总体热传导系数U可表示为:
其中Dbed为蓄热床外径;h为换热系数,包括自然对流和辐射;k为导热系数;rex为压力容器的外径;rin,rout分别为某隔热层的内表面和外表面半径;另外,下标1,2…表示同心轴的隔热层。下面为预热过程的简化传热方程:
T为当地温度;Tamb为初始温度,初次预热时为常温,二次预热时为上一试验周期结束后蓄热单元温度;为燃气质量流量;cp为燃气的定压比热容;z为与蓄热床顶端的距离。Cb为蓄热阵固体材料的比热容;M为蓄热阵的单位长度重量;t为时间变量。而且具有如下边界条件:1)当t=0时,T=Tamb+(Th-Tamb)e-bz,其中2)对于所有t>0,在z=0处均有T=Th,Th顶部进来的燃气温度。根据以上简化传热方程及边界条件,就可得到不同流量条件下预热之后的蓄热阵温度分布。
有益效果
本发明提出了一种空心砖型蓄热式加热器的总体设计评估方法,可用于空心砖型蓄热式加热器的总体设计和性能评估。
根据给定加热器的设计状态要求,采用本发明方法可以进行加热器方案设计,包括确定加热器的总高和最大直径、蓄热材料的高度和直径、隔热层的厚度等其它几何参数;采用本方法可以进一步评估蓄热式加热器的工作过程中的气动参数,即可以计算得到基准条件下的蓄热单元气流质量通量特性、基准运行状态下蓄热单元当地最大温差的沿程变化和预热过程蓄热阵温度径向分布等。本设计方法可为空心砖型蓄热式加热器方案设计和性能评估提供全套技术手段。
附图说明
下面结合附图和实施例对本发明一种空心砖蓄热式加热器的总体设计评估方法作进一步的详细说明。
图1为蓄热式加热器蓄热单元示意图。
图2为空心砖型蓄热单元单孔模型。
图3为空心砖型蓄热单元的单孔热应力分析模型。
图4为蓄热单元隔热计算的几何模型简图。
图5为基准条件下的蓄热单元气流质量通量特性(氧化铝)。
图6为基准条件下的蓄热单元气流质量通量特性(氧化锆)。
图7为基准运行状态下蓄热单元当地最大温差的沿程变化及最大允许值。
图8为加热器隔热层的径向截面视图。
图9为空心砖型蓄热式加热器总体结构示意图。
图10为第一调试阶段时高温燃气预热形成的蓄热阵温度分布评估
图11为第二调试阶段时高温燃气预热形成的蓄热阵温度分布评估。
图中:1.内环隔热层,2.外环隔热层,3.压力容器外壳,4.高温热衬,5.蓄热单元,6.预热燃烧器接口,7.辐射高温计,8.高温空气出口,9.温度测点,10.压力容器,11.隔热层,12.空心砖蓄热阵,13.底部支撑结构,14.空气进口。
具体实施方式
结合附图,对本发明实施例中的设计评估进行清楚、完整地描述,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。
实施例一
本实施例提供一种空心砖加热器蓄热单元设计及性能评估方法。
实施例中空心砖型蓄热式加热器设计的基准设计状态需求为:空气流量3.0kg/s、模拟总温2200K、模拟总压5.0MPa、有效运行时间至少50秒。表1给出了设计方案中的蓄热单元几何参数及采用本方法计算得到相关参数。蓄热材料采用高纯氧化铝砖和氧化锆砖。氧化铝材料配置于加热器下半部分较低温区域(1900K以下),氧化锆材料配置于加热器上半部分较高温区域(1900K-2400K),以满足2200K气流温度的加热需求。
基于空心砖型蓄热单元几何参数和氧化铝、氧化锆材料特性,在初始预热温度确定的前提下,对加热器性能进行计算与评估,在计算评估中取安全因子等于2。
表2为加热器设计运行状态及初始条件。根据以上初始条件,可以评估得到试验消耗热量大约占蓄热总量的14.3%。
表3列出了该基准状态下蓄热单元的几个重要参数,其中包括热应力限制参数与“漂浮限值”,带★上标分别为最大可允许的初始预热温度型斜率和当地冷却率两个热应力限制参数。
表1空心砖型蓄热式加热器蓄热单元几何参数
表2空心砖型蓄热加热器工作状态及初始条件
表3空心砖型蓄热单元的几个重要特性参数
图5和6分别显示了氧化铝、氧化锆蓄热阵中受限于热应力、“漂浮极限”的质量通量特性,星号*表示了对应于基准状态流量3.0kg/s时的实际气流质量通量,从图5和6可以看出,总体设计方案中蓄热阵实际气流质量通量远低于“漂浮限值”,可以确保加热器的安全运行。
图7显示了基准状态下0~50s运行时间内蓄热单元高度方向上当地最大温差的沿程变化以及最大允许温差。从评估结果来看,基准运行状态的当地最大温差均小于最大允许值,可以确保整体的热应力安全性;不同蓄热材质分界面处当地最大温差存在阶跃,这是因为材料特性差异引起的,但仍能满足安全性要求。
实施例二
本实施例为一种隔热层的校核方法。针对图8这一圆柱隔热结构布置进行校核计算。已知从里到外,各层分别为蓄热单元、高温热衬、内环隔热层、外环隔热层、不锈钢外壳,各层的直径分别如下:D1=450mm,D2=650mm,D3=850mm,D4=1050,D5=1150mm。采用多层圆筒隔热的一维稳态导热关系进行校核验证,可以确定出加热器中部等直圆柱段单位高度热损失功率为:
各层外表面温度为:
热面温度(即蓄热单元外表面温度)Tin按=2400K(高温段)和1900K(低温段)两种计算,经过隔热后不锈钢外壳温度为Tout=400K,周围环境温度按Tamb=288K算。代入各参数可计算得到如表5和6所示结果,各接触界面温度均满足相应隔热材料的使用温度要求。
表4最高热面温度2400K时加热器隔热结构布置及热工估算结果
表5最高热面温度1900K时加热器隔热结构布置及热工估算结果
由实施例1和实施例2可以大致确定出了空心砖型蓄热式加热器本体的基本规模,图9显示加热器本体的总体结构示意图,蓄热单元蓄热单元外径0.45m,高度5.0m,加热器本体高约7.0m,直径约1.2m。
本实施例为加热器初步预热评估方法。对于实施例一、二中的加热器设计方案,针对表6中的加热器两个运行状态,给出以下初始条件:a.从常温开始预热,即整个蓄热单元初始温度取288K;b.顶部进口燃气温度恒为1900K或2400K两种;c.预热过程分两个阶段进行,即小流量燃气预热阶段和大流量燃气预热阶段,小流量预热时为0.010kg/s,大流量预热时为0.100kg/s。图10和11分别显示了1900K和2400K时高温燃气预热形成的蓄热阵温度分布评估。其中小流量预热下达到准稳态温度分布条件大约需要3~4个小时,包括形成蓄热阵顶部恒温段在内的整个加热器预热时间大致在5个小时以内。
表6加热器蓄热材质垒积高度和目标状态要求
Claims (1)
1.一种空心砖型蓄热式加热器设计与评估方法,其特征包括:
(一)加热器蓄热单元几何设计,包括蓄热体截面的空隙率、比表面积、蓄热床的外径和横截面积等;
(二)蓄热单元性能分析方法,包括蓄热单元压力损失分析方法和热应力分析方法,可计算得到最大可允许的初始预热温度型斜率和当地冷却率两个热应力限制参数;
(三)加热器隔热层设计与分析方法,迭代计算隔热层各个隔热材料的厚度;
(四)加热器初步预热设计与分析方法。可得到不同流量条件下预热之后的蓄热阵温度分布;
采用本发明设计评估空心砖蓄热式加热器的步骤如下:
步骤一、根据给定加热器的基准设计状态要求,采用本发明方法可以进行加热器方案几何设计,包括确定加热器的总高和最大直径、蓄热材料的高度和直径、隔热层的厚度等其它几何参数;
步骤二、更进一步评估蓄热式加热器的工作过程中的气动参数,即可以计算得到加热器蓄热单元出口温度随时间的变化、蓄热单元内沿程气流温度和固体温度分布、预热过程蓄热阵温度径向分布等。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810189272.9A CN108595746A (zh) | 2018-03-08 | 2018-03-08 | 一种空心砖蓄热式加热器的总体设计评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810189272.9A CN108595746A (zh) | 2018-03-08 | 2018-03-08 | 一种空心砖蓄热式加热器的总体设计评估方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108595746A true CN108595746A (zh) | 2018-09-28 |
Family
ID=63625761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810189272.9A Pending CN108595746A (zh) | 2018-03-08 | 2018-03-08 | 一种空心砖蓄热式加热器的总体设计评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108595746A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112240719A (zh) * | 2020-10-20 | 2021-01-19 | 西安航天动力研究所 | 一种空心砖蓄热体结构以及蓄热式加热器 |
CN112378084A (zh) * | 2020-10-27 | 2021-02-19 | 西安航天动力研究所 | 一种空心砖蓄热阵的支撑结构 |
CN112902720A (zh) * | 2020-12-31 | 2021-06-04 | 中国航天空气动力技术研究院 | 一种陶瓷空心砖蓄热器的蓄热体设计方法 |
CN114199499A (zh) * | 2022-02-16 | 2022-03-18 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于大流量高温高压高速气体环境的蓄热元件 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104992042A (zh) * | 2015-08-10 | 2015-10-21 | 中国航天空气动力技术研究院 | 用于空心砖型蓄热式加热器的蓄热阵换热设计方法 |
CN105488292A (zh) * | 2015-12-14 | 2016-04-13 | 中广核工程有限公司 | 基于阀门仿真模型的高温阀门结构性能评估方法及系统 |
-
2018
- 2018-03-08 CN CN201810189272.9A patent/CN108595746A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104992042A (zh) * | 2015-08-10 | 2015-10-21 | 中国航天空气动力技术研究院 | 用于空心砖型蓄热式加热器的蓄热阵换热设计方法 |
CN105488292A (zh) * | 2015-12-14 | 2016-04-13 | 中广核工程有限公司 | 基于阀门仿真模型的高温阀门结构性能评估方法及系统 |
Non-Patent Citations (2)
Title |
---|
罗飞腾: "高温纯净空气风洞加热技术的应用与发展", 《世界科技研究与发展》 * |
罗飞腾: "高温风洞蓄热式加热器蓄热单元初步设计与分析", 《实验流体力学》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112240719A (zh) * | 2020-10-20 | 2021-01-19 | 西安航天动力研究所 | 一种空心砖蓄热体结构以及蓄热式加热器 |
CN112240719B (zh) * | 2020-10-20 | 2022-07-05 | 西安航天动力研究所 | 一种空心砖蓄热体结构以及蓄热式加热器 |
CN112378084A (zh) * | 2020-10-27 | 2021-02-19 | 西安航天动力研究所 | 一种空心砖蓄热阵的支撑结构 |
CN112378084B (zh) * | 2020-10-27 | 2022-01-07 | 西安航天动力研究所 | 一种空心砖蓄热阵的支撑结构 |
CN112902720A (zh) * | 2020-12-31 | 2021-06-04 | 中国航天空气动力技术研究院 | 一种陶瓷空心砖蓄热器的蓄热体设计方法 |
CN112902720B (zh) * | 2020-12-31 | 2022-07-05 | 中国航天空气动力技术研究院 | 一种陶瓷空心砖蓄热器的蓄热体设计方法 |
CN114199499A (zh) * | 2022-02-16 | 2022-03-18 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于大流量高温高压高速气体环境的蓄热元件 |
CN114199499B (zh) * | 2022-02-16 | 2022-04-22 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于大流量高温高压高速气体环境的蓄热元件 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108595746A (zh) | 一种空心砖蓄热式加热器的总体设计评估方法 | |
Tan et al. | Experimental study on liquid/solid phase change for cold energy storage of Liquefied Natural Gas (LNG) refrigerated vehicle | |
Mahdi et al. | Experimental study on the melting behavior of a phase change material in a conical coil latent heat thermal energy storage unit | |
He et al. | Employing spiral fins to improve the thermal performance of phase-change materials in shell-tube latent heat storage units | |
Zhang et al. | Design and heat transfer characteristics analysis of combined active and passive thermal protection system for hydrogen fueled scramjet | |
JP2012517789A5 (zh) | ||
CN106066235A (zh) | 超临界水窄通道自然循环实验装置及方法 | |
Sharma et al. | Solar tower based aluminum heat treatment system: Part I. Design and evaluation of an open volumetric air receiver | |
Aksoy et al. | Performance testing of a Fresnel/Stirling micro solar energy conversion system | |
CN111520760A (zh) | 一种冲击/气膜双层壁复合冷却方式的燃烧室火焰筒壁面结构 | |
CN108679843A (zh) | 一种用于高温风洞空心砖型蓄热式加热器隔热层设计 | |
Kumar et al. | Experimental study of multi tubular sensible heat storage system fitted with wire coil inserts | |
CN105097056A (zh) | 一种非能动余热排出热交换器模拟实验系统 | |
CN106248726A (zh) | 抗氧化涂层在500~2300℃区间热震/热疲劳性能和辐射特性测试装置 | |
Ji et al. | Experimental study on heat transfer performance of high temperature heat pipe with large length-diameter ratio for heat utilization of concentrated solar energy | |
Minzer et al. | Evaporation in parallel pipes––splitting characteristics | |
CN204720082U (zh) | 一种非能动余热排出热交换器模拟实验系统 | |
CN110359972B (zh) | 超临界co2为工质的熔盐电蓄热发电装置 | |
Vesely et al. | Design of experimental loop with supercritical carbon dioxide | |
CN102661972B (zh) | 燃气轮机隔热罩隔热性能的测定方法 | |
CN117072253A (zh) | 重型燃气轮机高温叶片热障涂层及设计、制造与评价方法 | |
Liu et al. | Experimental study on the heat transfer characteristics of fluoride salt in the new conceptual passive heat removal system of molten salt reactor | |
Yanhua et al. | Analysis of diffusion process and influence factors in the air ingress accident of the HTR-PM | |
CN112378084B (zh) | 一种空心砖蓄热阵的支撑结构 | |
CN205719474U (zh) | 超临界水窄通道自然循环实验装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180928 |