CN104992042B - 用于空心砖型蓄热式加热器的蓄热阵换热设计方法 - Google Patents

用于空心砖型蓄热式加热器的蓄热阵换热设计方法 Download PDF

Info

Publication number
CN104992042B
CN104992042B CN201510486214.9A CN201510486214A CN104992042B CN 104992042 B CN104992042 B CN 104992042B CN 201510486214 A CN201510486214 A CN 201510486214A CN 104992042 B CN104992042 B CN 104992042B
Authority
CN
China
Prior art keywords
mrow
msub
mfrac
accumulation
heat
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.)
Active
Application number
CN201510486214.9A
Other languages
English (en)
Other versions
CN104992042A (zh
Inventor
卢洪波
刘展
郭孝国
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Academy of Aerospace Aerodynamics CAAA
Original Assignee
China Academy of Aerospace Aerodynamics CAAA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Academy of Aerospace Aerodynamics CAAA filed Critical China Academy of Aerospace Aerodynamics CAAA
Priority to CN201510486214.9A priority Critical patent/CN104992042B/zh
Publication of CN104992042A publication Critical patent/CN104992042A/zh
Application granted granted Critical
Publication of CN104992042B publication Critical patent/CN104992042B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了用于空心砖型蓄热式加热器的蓄热阵换热设计方法，包括：建立蓄热阵的一维换热分析模型，将多个换热通道的蓄热阵简化为当量流通面积圆管模型；推导一维换热的控制方程与定解条件；构建控制方程的差分格式；求解蓄热阵和气体的温度分布；给定蓄热阵的热力学参数和工作条件，求解蓄热阵以及气体在不同的几何参数时的温度分布，以对比分析蓄热阵在不同的几何参数下的换热性能，进而确定出具备所须换热性能的蓄热阵的布局。本发明简便实用，仅需按照预定的步骤，代入蓄热阵的热力学参数、工作条件参数及几何参数，在几秒内即可获得一组参数对应的蓄热阵的换热性能。对比不同几何参数对应的蓄热阵换热性能，即可筛选出较佳的几何布局方案。

Description

用于空心砖型蓄热式加热器的蓄热阵换热设计方法
技术领域
本发明涉及一种蓄热阵的换热设计方法，尤其涉及一种用于空心砖型蓄热式加热器的蓄热阵换热设计方法。
背景技术
高超声速风洞是探究高超声速气动力学与气动热力学问题和确认各类高超声速飞行器设计方法正确性所不可或缺的地面设备，需复现飞行器的真实工作条件。在高超声速条件下，风洞试验气体必须进行加热，以保证风洞能够正常工作、甚至复现飞行器实际工作的总温。进入喷管前的空气温度若不够高，经喷管的膨胀加速降温后，空气将凝结，形成气-液两相流，相应地风洞不能正常工作。此外，某些特种风洞试验如超燃冲压发动机性能试验要求进入喷管前的空气温度达到飞行器真实飞行的总温，以确保试验数据尽可能地真实和可信。因此加热器是高超声速风洞的关键部件之一，直接关系到风洞能否正常工作、乃至复现飞行器实际工作的总温。
国内外广泛采用的加热方式主要有三大类：电弧加热、燃烧加热与蓄热式加热，但只有蓄热式加热方式能获得与真实空气完全一致的气流条件，应用于常规高超声速风洞，所获得的试验数据(特别是超燃冲压发动机地面性能试验数据)最为可靠。目前，国内尚无高温(T＞1000K)、蓄热式加热的纯净空气风洞，亟需研制与建设这类地面试验设施，以拓宽地面风洞模拟范围并提高超燃冲压发动机等特种试验的数据可靠性。高温蓄热式加热器主要有两种结构：鹅卵石状加热器和空心砖加热器。与鹅卵石状结构相比，多孔空心砖结构可缓解蓄热结构的应力集中问题，提高加热器工作温度，并有效减少粉尘量，同时圆柱形的气体通道也有效地降低了流阻，进而大幅度降低了浮起事故的发生概率。然而如何快速设计出换热性能较优的空心砖型蓄热式加热器蓄热阵布局方案，目前尚无实用的换热设计方法，本发明就是在这样的背景下提出的。
发明内容
针对上述技术问题，本发明提供了一种简便实用的用于空心砖型蓄热式加热器的蓄热阵换热设计方法。
本发明提供的技术方案为：
用于空心砖型蓄热式加热器的蓄热阵换热设计方法，包括下列步骤：
步骤一、建立空心砖型蓄热式加热器的蓄热阵的一维换热分析模型，将具有多个换热通道的蓄热阵简化为当量流通面积圆管模型；
步骤二、推导一维换热的控制方程与定解条件，
控制方程为定解条件为
其中，θ(x，τ)、T(x，τ)分别为蓄热阵的温度和气体的温度，x、τ分别为轴向位置和运行时间，η、d分别为蓄热阵的孔隙率和单个换热通道的直径，ρb、Cb分别为蓄热阵材料的密度和比热，Cp为气体的定压比热，h为对流换热系数，Tin为入口处气体温度，为起始时刻蓄热阵的轴向温度分布，单个当量流通面积圆管的流通面积A＝M/(ρU)，ρ、U和M分别为该当量流通面积圆管内气体的密度、流速和流量；
步骤三、利用步骤二中的控制方程构建控制方程的差分格式，控制方程的差分格式由以下3个公式Eq.1、Eq.2和Eq.3组成，其中τ＝mΔτ，x＝pΔx，θm，p＝θ(mΔτ，pΔx)，Tm，p＝T(mΔτ，pΔx)：
步骤四、利用步骤三的控制方程的差分格式求解蓄热阵和气体的温度分布，以蓄热阵的温度分布代表蓄热阵的换热性能，具体过程为：
步骤(1)给定蓄热阵初始温度分布
步骤(2)给定气体入口温度Tm，0＝Tin(p＝0)；
步骤(3)由Eq.2知气体的初始温度分布
步骤(4)由Eq.1确定蓄热阵温度分布
步骤(5)由Eq.3确定蓄热阵温度分布
步骤(6)由Eq.2确定气体温度
步骤(7)重复步骤(5)与步骤(6)，求出τ＝Δτ时刻蓄热阵和气体的温度分布θ1，p、T1，p(m＝1，p≥2，且p为整数)；
步骤(8)重复步骤(4)～步骤(7)即可求出所有时刻蓄热阵和气体的温度分布；
步骤五、给定蓄热阵的热力学参数和工作条件，重复执行步骤四，求解蓄热阵以及气体在不同的几何参数时的温度分布，以对比分析蓄热阵在不同的几何参数下的换热性能，进而确定出具备所须换热性能的蓄热阵的布局，其中，热力学参数包括ρb、Cb和h，工作条件包括Tin、Cp、和A，几何参数包括η和d。
优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤五中，重复执行步骤四的过程中，先求解蓄热阵在预热阶段的温度分布，再求解蓄热阵在吹风阶段的温度分布，并假设蓄热阵在预热阶段的结束时刻的温度分布与吹风阶段的起始时刻的温度分布相同。
优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤五中，在预热阶段，τ＝0时为环境温度，x＝0为蓄热阵的顶端，Tin为预热混合燃气温度；在吹风阶段：τ＝0时为预热结束时刻的蓄热阵温度，x＝0为蓄热阵的底端，Tin为待加热的试验气体温度，与环境温度相同。
优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤二中，控制方程的推导过程为：
在当量流通面积圆管模型中，气流方向上的微元控制体dx内，在dτ时间内，气体从控制边界吸收的热量为h(θ-T)(Nπd·dx)·dτ，气体增加的内能为气体从控制边界带走的净能量为蓄热阵含有换热通道的数目为N，Cv为气体的定容比热，则有
整理可得：
由理想气体状态方程知：蓄热器工作过程中，可忽略不计，又ρU＝M/A，进一步整理可得：
优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤二中，控制方程的推导过程为：
在dτ时间内，dx长度蓄热阵减小的内能为蓄热阵通过控制边界向气体传递的能量为h(θ-T)(Nπd·dx)·dτ，则有
整理可得：
优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤三中，利用步骤二中的控制方程构建控制方程的差分格式的具体过程为：
用差分代替微分，并令τ＝mΔτ、x＝pΔx，则 由步骤(2)所述的控制方程可得：
其中，Δτ时间内的平均值Δx长度上的平均值代入并整理可得：
又令并整理可得：
联立式Eq.1与Eq.2可得式Eq.3：
优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤一中，假设蓄热阵在横截面上的温度近似一致，每个换热通道为圆柱形通道，且每个换热通道内气体在横截面上的温度均匀，气体与蓄热阵之间的热量只通过热对流方式传递，忽略蓄热阵向外的热损失，从而将具有多个换热通道的蓄热阵简化为当量流通面积圆管模型。
本发明所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法与现有技术相比的优势如下：
(1)本发明思路明确、简便实用，仅需按照预定的步骤，代入蓄热阵的热力学参数、工作条件参数及几何参数，在几秒内即可获得一组参数对应的蓄热阵的换热性能。
(2)对比不同几何参数对应的蓄热阵换热性能，即可筛选出较佳的几何布局方案。
(3)换热设计中结合蓄热阵的工作特点，采用预热与吹风两工作阶段联合分析，吹风阶段的初始条件摆脱了人为的假设，使计算结果更为精确。
附图说明
图1为本发明所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法的流程图；
图2(a)为空心砖型蓄热式加热器的蓄热阵的结构示意图；图2(b)为本发明建立的当量流通面积圆管模型；
图3(a)为设计工况下预热阶段蓄热阵的温度特性曲线，图3(b)为设计工况下吹风阶段蓄热阵的温度特性曲线，其中，设计工况条件如下：预热混合燃气温度和流量分别为1200K和1kg/s，环境温度为300K，预热结束时刻蓄热阵底端温度为500K，试验空气流量和压力分别为10kg/s和12MPa，d＝7mm，η＝0.3；
图4(a)为设计工况下吹风阶段出口处气体温度特性曲线，d＝7mm，η取值范围为0.2～0.4，图4(b)为设计工况下吹风阶段出口处气体温度特性曲线，d取值范围为5～10mm，η取值0.3，其中，设计工况条件如下：预热混合燃气温度和流量分别为1200K和1kg/s，环境温度为300K，预热结束时刻蓄热阵底端温度为500K，试验空气流量和压力分别为10kg/s和12MPa。
具体实施方式
下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。
如图1所示，本发明提供一种用于空心砖型蓄热式加热器的蓄热阵换热设计方法，包括下列步骤：
步骤一、建立空心砖型蓄热式加热器的蓄热阵的一维换热分析模型，将单个换热通道简化为当量流通面积圆管模型。
假设空心砖型蓄热式加热器的蓄热阵在横截面上的温度可近似一致，每个换热通道为圆柱形通道，且每个换热通道的气体在横截面上的温度分布均匀，气体与蓄热阵之间的热量只通过热对流方式传递，忽略蓄热阵向外的热损失，于是具有多个换热通道的蓄热阵可简化为当量流通面积圆管模型的一维换热分析模型。当量流通面积圆管内气体的密度ρ、流速U、流量M与蓄热阵内气体的对应量相同，当量流通面积圆管的流通面积A＝M/(ρU)。若蓄热阵含有换热通道的数目为N，则气体等效流通面积也可表示为A＝Nπd2/4。
步骤二、推导一维换热的控制方程与定解条件。
对于气体流动而言，取图2(b)虚线所示的微元控制体为研究对象，则dτ时间内，微元控制体内气体内能的增量为气体从控制边界进出带走的净能量为气体从控制边界吸收的热量为h(θ-T)(Nπd·dx)·dτ，由能量守恒知，即气体从控制边界吸收的热量＝微元控制体内气体增加的内能+气体从控制边界带走的净能量，即
整理可得：
由理想气体状态方程知：蓄热器工作过程中，气体的压力随时间变化非常缓慢，相对而言，该项是个小量，可忽略不计，又ρU＝M/A，进一步整理可得：
对于蓄热阵而言，在dτ时间内，dx长度蓄热阵减少的内能为蓄热阵通过控制边界向气体传递的能量为：h(θ-T)(Nπd·dx)·dτ，由能量守恒知，蓄热阵减少的内能全部通过控制边界散去，即dτ→0，整理可得：
对应的定解条件为预热阶段：τ＝0时为环境温度，x＝0为蓄热阵的顶端，Tin为预热混合燃气温度；吹风阶段：τ＝0时为预热结束时刻的蓄热阵温度，x＝0为蓄热阵的底端，Tin为待加热的试验气体温度、与环境温度相同。
步骤三、利用步骤二中的控制方程构建控制方程的差分格式，具体过程为：
用差分代替微分，并令τ＝mΔτ、x＝pΔx，则 由步骤二所述的控制方程可得：
其中，Δτ时间内的平均值Δx长度上的平均值代入并整理可得：
又令并整理可得：
联立式Eq.1与Eq.2可得式Eq.3：
步骤四、利用步骤三的控制方程的差分格式求解蓄热阵和气体的温度分布，以蓄热阵的温度分布代表蓄热阵的换热性能，具体过程为：
步骤(1)给定蓄热阵初始温度分布
步骤(2)给定气体入口温度Tm，0＝Tin(p＝0)；
步骤(3)由Eq.2知气体初始温度分布
步骤(4)由Eq.1确定蓄热阵温度
步骤(5)由Eq.3确定蓄热阵温度
步骤(6)由Eq.2确定气体温度
步骤(7)重复步骤(5)与步骤(6)即可求出τ＝Δτ时刻蓄热阵和气体的温度分布θ1，p、T1，p(m＝1，p＝2，3，4…)；(即m＝1，p≥2，且p为整数)；
步骤(8)重复步骤(4)～(7)即可求出所有时刻蓄热阵和气体的温度分布。
步骤五、给定蓄热阵的热力学参数和工作条件，重复执行步骤四，求解蓄热阵以及气体在不同的几何参数时的温度分布，以对比分析蓄热阵在不同的几何参数下的换热性能，进而确定出具备所须换热性能(即换热性能优秀的)的蓄热阵的布局，其中，热力学参数包括ρb、Cb和h，工作条件包括Tin、Cp、和A，几何参数包括η和d。
本发明提供的方法思路清晰、简便实用，可用于高超声速风洞蓄热式加热器和工业加热炉的蓄热阵换热性能的快速分析，并筛选换出热性能较优的蓄热阵布局方案。
优选地，步骤五中，重复执行步骤四的过程中，先求解蓄热阵在预热阶段的温度分布，再求解蓄热阵在吹风阶段的温度分布，并假设蓄热阵在预热阶段的结束时刻的温度分布与吹风阶段的起始时刻的温度分布相同，从而摆脱蓄热阵在吹风阶段的初始时刻温度分布的人工假设，使计算结果更准确。
实施例
图3(a)为设计工况下预热阶段蓄热阵的温度特性曲线，图3(b)为设计工况下吹风阶段蓄热阵的温度特性曲线，其中，设计工况条件如下：预热混合燃气温度和流量分别为1200K和1kg/s，环境温度为300K，预热结束时刻蓄热阵底端温度为500K，试验空气流量和压力分别为10kg/s和12MPa，d＝7mm，η＝0.3。图3(a)和图3(b)中，吹风阶段的初始时刻蓄热阵的温度分布与预热阶段的终止时刻蓄热阵的温度分布相同。
图4(a)为设计工况下吹风阶段出口处气体温度特性曲线，d＝7mm，η取值范围为0.2～0.4，图4(b)为设计工况下吹风阶段出口处气体温度特性曲线，d取值范围为5～10mm，η取值0.3，其中，设计工况条件如下：预热混合燃气温度和流量分别为1200K和1kg/s，环境温度为300K，预热结束时刻蓄热阵底端温度为500K，试验空气流量和压力分别为10kg/s和12MPa。图4(a)中，给定了d值为7mm，随着η值增加，在相同时刻出口处气体(即待加热的试验气体)的温度更高，即提高蓄热阵的孔隙率有助于提高换热性能。但当η值取值为0.4时，在吹风持续到60s时曲线斜率增加，出口处气体的温度下降速度加快，当出口处气体的温度下降至一定温度以下，蓄热阵就需要重新预热，因此，η值取值过大会导致蓄热阵的工作时间缩短。在图4(b)中，η取值0.3，d取值范围为5～10mm，随着d值增加，在相同时刻出口处气体的温度更低，即减小d值有助于提高换热性能。当d值过大(比如取值10mm)时，曲线的斜率较大，出口处气体的温度下降过快，即蓄热阵的工作时间缩短。
通过给定其他的热力学参数和工作条件，分别改变几何参数d和η，可以得出不同几何参数下蓄热阵的换热性能，进而设计具有合适的d和η的蓄热阵。
尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。

Claims (7)

1.用于空心砖型蓄热式加热器的蓄热阵换热设计方法，其特征在于，包括下列步骤：
步骤一、建立空心砖型蓄热式加热器的蓄热阵的一维换热分析模型，将具有多个换热通道的蓄热阵简化为当量流通面积圆管模型；
步骤二、推导一维换热的控制方程与定解条件，
控制方程为定解条件为
其中，θ(x,τ)、T(x,τ)分别为蓄热阵的温度和气体的温度，x、τ分别为轴向位置和运行时间，η、d分别为蓄热阵的孔隙率和单个换热通道的直径，ρb、Cb分别为蓄热阵材料的密度和比热，Cp为气体的定压比热，h为对流换热系数，Tin为入口处气体温度，为起始时刻蓄热阵的轴向温度分布，当量流通面积圆管的流通面积A＝M/(ρU)，ρ、U和M分别为该当量流通面积圆管内气体的密度、流速和流量；
步骤三、利用步骤二中的控制方程构建控制方程的差分格式，控制方程的差分格式由以下3个公式Eq.1、Eq.2和Eq.3组成，其中τ＝mΔτ，x＝pΔx，θm,p＝θ(mΔτ,pΔx)，Tm,p＝T(mΔτ,pΔx)：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mi>E</mi> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>E</mi> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> </mrow> </mfrac> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>E</mi> <mi>q</mi> <mn>.1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mi>F</mi> <mrow> <mn>1</mn> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>F</mi> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>E</mi> <mi>q</mi> <mn>.2</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mi>E</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>F</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mi>E</mi> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>F</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mi>E</mi> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mi>F</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mi>E</mi> <mi>F</mi> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>E</mi> <mi>q</mi> <mn>.3</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
步骤四、利用步骤三的控制方程的差分格式求解蓄热阵和气体的温度分布，以蓄热阵的温度分布代表蓄热阵的换热性能，具体过程为：
步骤(1)给定蓄热阵初始温度分布
步骤(2)给定气体入口温度Tm,0＝Tin(p＝0)；
步骤(3)由Eq.2知气体的初始温度分布
<mrow> <msub> <mi>T</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mi>F</mi> <mrow> <mn>1</mn> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>F</mi> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>T</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>=</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
步骤(4)由Eq.1确定蓄热阵温度分布
<mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mn>1</mn> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mi>E</mi> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mrow> <mn>1</mn> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>+</mo> <msub> <mi>T</mi> <mrow> <mn>0</mn> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>E</mi> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> </mrow> </mfrac> <msub> <mi>&amp;theta;</mi> <mrow> <mn>0</mn> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>,</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> <mo>=</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
步骤(5)由Eq.3确定蓄热阵温度分布
<mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mn>1</mn> <mo>,</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>E</mi> <mo>)</mo> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>F</mi> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>&amp;theta;</mi> <mrow> <mn>0</mn> <mo>,</mo> <mn>1</mn> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mi>E</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>F</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>T</mi> <mrow> <mn>0</mn> <mo>,</mo> <mn>1</mn> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mi>E</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>F</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>T</mi> <mrow> <mn>1</mn> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mi>E</mi> <mi>F</mi> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>&amp;theta;</mi> <mrow> <mn>1</mn> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>,</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> <mo>=</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
步骤(6)由Eq.2确定气体温度
步骤(7)重复步骤(5)与步骤(6)，求出τ＝Δτ时刻蓄热阵和气体的温度分布θ1,p、T1,p(m＝1，p≥2，且p为整数)；
步骤(8)重复步骤(4)～步骤(7)即可求出所有时刻蓄热阵和气体的温度分布；
步骤五、给定蓄热阵的热力学参数和工作条件，重复执行步骤四，求解蓄热阵以及气体在不同的几何参数时的温度分布，以对比分析蓄热阵在不同的几何参数下的换热性能，进而确定出具备所须换热性能的蓄热阵的布局，其中，热力学参数包括ρb、Cb和h，工作条件包括Tin、Cp、和A，几何参数包括η和d。
2.如权利要求1所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法，其特征在于，步骤五中，重复执行步骤四的过程中，先求解蓄热阵在预热阶段的温度分布，再求解蓄热阵在吹风阶段的温度分布，并假设蓄热阵在预热阶段的结束时刻的温度分布与吹风阶段的起始时刻的温度分布相同。
3.如权利要求2所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法，其特征在于，步骤五中，在预热阶段，τ＝0时为环境温度，x＝0为蓄热阵的顶端，Tin为预热混合燃气温度；在吹风阶段：τ＝0时为预热结束时刻的蓄热阵温度，x＝0为蓄热阵的底端，Tin为待加热的试验气体温度，与环境温度相同。
4.如权利要求1至3中任一项所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法，其特征在于，步骤二中，控制方程的推导过程为：
在当量流通面积圆管模型中，气流方向上的微元控制体dx内，在dτ时间内，气体从控制边界吸收的热量为h(θ-T)(Nπd·dx)·dτ，气体增加的内能为气体从控制边界带走的净能量为蓄热阵含有换热通道的数目为N，Cv为气体的定容比热，U为当量流通面积圆管模型内气体的流速，则有
<mrow> <mi>h</mi> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>-</mo> <mi>T</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>N</mi> <mi>&amp;pi;</mi> <mi>d</mi> <mo>&amp;CenterDot;</mo> <mi>d</mi> <mi>x</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>d</mi> <mi>&amp;tau;</mi> <mo>=</mo> <mfrac> <mo>&amp;part;</mo> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>{</mo> <mfrac> <mrow> <msup> <mi>N&amp;pi;d</mi> <mn>2</mn> </msup> </mrow> <mn>4</mn> </mfrac> <msub> <mi>&amp;rho;UC</mi> <mi>p</mi> </msub> <mi>T</mi> <mo>}</mo> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>&amp;tau;</mi> <mo>+</mo> <mfrac> <mo>&amp;part;</mo> <mrow> <mo>&amp;part;</mo> <mi>&amp;tau;</mi> </mrow> </mfrac> <mo>{</mo> <mfrac> <mrow> <msup> <mi>N&amp;pi;d</mi> <mn>2</mn> </msup> </mrow> <mn>4</mn> </mfrac> <mi>&amp;rho;</mi> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mi>v</mi> </msub> <mi>T</mi> <mo>)</mo> </mrow> <mo>}</mo> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>&amp;tau;</mi> <mo>,</mo> </mrow>
整理可得：
<mrow> <mfrac> <mi>d</mi> <mn>4</mn> </mfrac> <mo>&amp;lsqb;</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mrow> <mo>(</mo> <msub> <mi>&amp;rho;C</mi> <mi>v</mi> </msub> <mi>T</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <mi>&amp;tau;</mi> </mrow> </mfrac> <mo>+</mo> <msub> <mi>&amp;rho;UC</mi> <mi>p</mi> </msub> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>T</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>&amp;rsqb;</mo> <mo>=</mo> <mi>h</mi> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>-</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
由理想气体状态方程知：蓄热器工作过程中，可忽略不计，又ρU＝M/A，进一步整理可得：
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>T</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mn>4</mn> <mi>A</mi> <mi>h</mi> </mrow> <mrow> <msub> <mi>dMC</mi> <mi>p</mi> </msub> </mrow> </mfrac> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>-</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
5.如权利要求1至3中任一项所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法，其特征在于，步骤二中，控制方程的推导过程为：
在dτ时间内，dx长度蓄热阵减小的内能为蓄热阵通过控制边界向气体传递的能量为h(θ-T)(Nπd·dx)·dτ，蓄热阵含有换热通道的数目为N，则有
<mrow> <mfrac> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>4</mn> <mi>&amp;eta;</mi> </mrow> </mfrac> <msub> <mi>&amp;rho;</mi> <mi>b</mi> </msub> <msub> <mi>C</mi> <mi>b</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>&amp;delta;</mi> <mi>&amp;theta;</mi> <mo>=</mo> <mi>h</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>-</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>d</mi> <mi>&amp;tau;</mi> <mo>,</mo> </mrow>
整理可得：
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;theta;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>&amp;tau;</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mn>4</mn> <mi>&amp;eta;</mi> </mrow> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mfrac> <mi>h</mi> <mrow> <msub> <mi>&amp;rho;</mi> <mi>b</mi> </msub> <msub> <mi>C</mi> <mi>b</mi> </msub> </mrow> </mfrac> <mrow> <mo>(</mo> <mi>T</mi> <mo>-</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
6.如权利要求1至3中任一项所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法，其特征在于，步骤三中，利用步骤二中的控制方程构建控制方程的差分格式的具体过程为：
用差分代替微分，并令τ＝mΔτ、x＝pΔx，则 由步骤(2)所述的控制方程可得：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mfrac> <mrow> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>x</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mn>4</mn> <mi>A</mi> <mi>h</mi> </mrow> <mrow> <msub> <mi>dMC</mi> <mi>p</mi> </msub> </mrow> </mfrac> <mover> <msub> <mrow> <mo>(</mo> <mrow> <mi>&amp;theta;</mi> <mo>-</mo> <mi>T</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>x</mi> </mrow> </msub> <mo>&amp;OverBar;</mo> </mover> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>&amp;tau;</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mn>4</mn> <mi>&amp;eta;</mi> </mrow> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mi>&amp;eta;</mi> </mrow> <mo>)</mo> </mrow> </mrow> </mfrac> <mfrac> <mi>h</mi> <mrow> <msub> <mi>&amp;rho;</mi> <mi>b</mi> </msub> <msub> <mi>C</mi> <mi>b</mi> </msub> </mrow> </mfrac> <mover> <msub> <mrow> <mo>(</mo> <mrow> <mi>T</mi> <mo>-</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>&amp;tau;</mi> </mrow> </msub> <mo>&amp;OverBar;</mo> </mover> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>
其中，Δτ时间内的平均值Δx长度上的平均值代入并整理可得：
<mrow> <mfrac> <mrow> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>x</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mi>A</mi> <mi>h</mi> </mrow> <mrow> <msub> <mi>dMC</mi> <mi>p</mi> </msub> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow>
<mrow> <mfrac> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>&amp;tau;</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mi>&amp;eta;</mi> <mi>h</mi> </mrow> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;eta;</mi> <mo>)</mo> <msub> <mi>d&amp;rho;</mi> <mi>b</mi> </msub> <msub> <mi>C</mi> <mi>b</mi> </msub> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>,</mo> </mrow>
又令并整理可得：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mi>E</mi> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>E</mi> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> </mrow> </mfrac> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>E</mi> <mi>q</mi> <mn>.1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mi>F</mi> <mrow> <mn>1</mn> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>F</mi> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>E</mi> <mi>q</mi> <mn>.2</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>
联立式Eq.1与Eq.2可得式Eq.3：
<mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>E</mi> <mo>)</mo> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>F</mi> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mi>E</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>F</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>p</mi> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mi>E</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>F</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mi>E</mi> <mi>F</mi> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mi>E</mi> <mo>+</mo> <mi>F</mi> </mrow> </mfrac> <msub> <mi>&amp;theta;</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>p</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mi>E</mi> <mi>q</mi> <mn>.3.</mn> </mrow>
7.如权利要求1至3中任一项所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法，其特征在于，步骤一中，假设蓄热阵在横截面上的温度近似一致，每个换热通道为圆柱形通道，且每个换热通道内气体在横截面上的温度均匀，气体与蓄热阵之间的热量只通过热对流方式传递，忽略蓄热阵向外的热损失，从而将具有多个换热通道的蓄热阵简化为当量流通面积圆管模型。
CN201510486214.9A 2015-08-10 2015-08-10 用于空心砖型蓄热式加热器的蓄热阵换热设计方法 Active CN104992042B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510486214.9A CN104992042B (zh) 2015-08-10 2015-08-10 用于空心砖型蓄热式加热器的蓄热阵换热设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510486214.9A CN104992042B (zh) 2015-08-10 2015-08-10 用于空心砖型蓄热式加热器的蓄热阵换热设计方法

Publications (2)

Publication Number Publication Date
CN104992042A CN104992042A (zh) 2015-10-21
CN104992042B true CN104992042B (zh) 2018-03-30

Family

ID=54303856

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510486214.9A Active CN104992042B (zh) 2015-08-10 2015-08-10 用于空心砖型蓄热式加热器的蓄热阵换热设计方法

Country Status (1)

Country Link
CN (1) CN104992042B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108595746A (zh) * 2018-03-08 2018-09-28 西北工业大学 一种空心砖蓄热式加热器的总体设计评估方法
CN108827051A (zh) * 2018-07-09 2018-11-16 西北工业大学 一种用于超高温蓄热式加热器的氧化铝蓄热体
CN112378084B (zh) * 2020-10-27 2022-01-07 西安航天动力研究所 一种空心砖蓄热阵的支撑结构
CN112902720B (zh) * 2020-12-31 2022-07-05 中国航天空气动力技术研究院 一种陶瓷空心砖蓄热器的蓄热体设计方法
CN112762444B (zh) * 2020-12-31 2022-08-12 中国航天空气动力技术研究院 一种通气过程初始时刻蓄热体温度分布曲线的优化方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100980782B1 (ko) * 2007-12-27 2010-09-10 한국생산기술연구원 다항식 회귀 방법을 통한 잠열 방출량 연산 장치 및 방법
CN102819643B (zh) * 2012-08-14 2016-01-20 中冶南方工程技术有限公司 热风炉传热蓄热仿真模型的建模方法及其应用

Also Published As

Publication number Publication date
CN104992042A (zh) 2015-10-21

Similar Documents

Publication Publication Date Title
CN104992042B (zh) 用于空心砖型蓄热式加热器的蓄热阵换热设计方法
CN107218133A (zh) 一种预冷吸气式发动机用高效紧凑预冷换热器
Zhang Heat and mass transfer in a cross-flow membrane-based enthalpy exchanger under naturally formed boundary conditions
Zhang et al. Heat and moisture transfer in application scale parallel-plates enthalpy exchangers with novel membrane materials
CN105973052A (zh) 错流导管换热器中的管
CN103592142A (zh) 适用于大型发电设备换热器传热与阻力特性的测试系统
CN109506876B (zh) 一种大气边界层环境风洞中温度层结模拟装置及方法
CN115290288B (zh) 一种模拟Ma8环境的超高温纯净空气发生系统及方法
CN115618758A (zh) 一种电弧加热高超风洞前室及其气动设计方法
CN114154438A (zh) 一种碱金属热管冷态启动三阶段计算方法
Li et al. Experimental insights into thermal performance of a microtube precooler with drastic coolant properties variation and precooling impacts on turbojet engine operation
Yang et al. Mass transfer study on high-pressure membrane dehumidification applied to aircraft environmental control system
Wang et al. Research on the cooling performance of the discontinuous transpiration surface structure for the leading edge of a hypersonic vehicle
CN104480288B (zh) 预防局部热处理时换热管产生凹痕的装置及方法
CN110705075A (zh) 一种静压箱孔板的优化设计方法
CN101606036A (zh) 热量回收方法及设备
CN110940526B (zh) 一种发动机试车台用单股高温热气与冷空气的掺混装置
CN205279811U (zh) 一种带壳体冷却的管翅式气－液换热器
CN109030789A (zh) 盐浴式煤油加热与裂解模拟试验方法
CN106871496B (zh) 室内换热器和空调器
CN101303205B (zh) 充气热管在防止露点腐蚀中的应用
CN109050984A (zh) 一种皱褶式主动冷却热防护承力一体化结构
CN205279810U (zh) 一种管翅式气－液换热器
CN206584652U (zh) 一种换热器的流量分配演示机构和制冷循环机构的流量分配演示机构
CN208420910U (zh) 盐浴式煤油加热与裂解模拟试验系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant