用于空心砖型蓄热式加热器的蓄热阵换热设计方法
【技术领域】
[0001] 本发明涉及一种蓄热阵的换热设计方法，尤其涉及一种用于空心砖型蓄热式加热 器的蓄热阵换热设计方法。
【背景技术】
[0002] 高超声速风洞是探宄高超声速气动力学与气动热力学问题和确认各类高超声速 飞行器设计方法正确性所不可或缺的地面设备，需复现飞行器的真实工作条件。在高超声 速条件下，风洞试验气体必须进行加热，以保证风洞能够正常工作、甚至复现飞行器实际工 作的总温。进入喷管前的空气温度若不够高，经喷管的膨胀加速降温后，空气将凝结，形成 气-液两相流，相应地风洞不能正常工作。此外，某些特种风洞试验如超燃冲压发动机性能 试验要求进入喷管前的空气温度达到飞行器真实飞行的总温，以确保试验数据尽可能地真 实和可信。因此加热器是高超声速风洞的关键部件之一，直接关系到风洞能否正常工作、乃 至复现飞行器实际工作的总温。
[0003] 国内外广泛采用的加热方式主要有三大类：电弧加热、燃烧加热与蓄热式加热，但 只有蓄热式加热方式能获得与真实空气完全一致的气流条件，应用于常规高超声速风洞， 所获得的试验数据（特别是超燃冲压发动机地面性能试验数据）最为可靠。目前，国内尚 无高温（T> 1000K)、蓄热式加热的纯净空气风洞，亟需研制与建设这类地面试验设施，以 拓宽地面风洞模拟范围并提高超燃冲压发动机等特种试验的数据可靠性。高温蓄热式加热 器主要有两种结构：鹅卵石状加热器和空心砖加热器。与鹅卵石状结构相比，多孔空心砖结 构可缓解蓄热结构的应力集中问题，提高加热器工作温度，并有效减少粉尘量，同时圆柱形 的气体通道也有效地降低了流阻，进而大幅度降低了浮起事故的发生概率。然而如何快速 设计出换热性能较优的空心砖型蓄热式加热器蓄热阵布局方案，目前尚无实用的换热设计 方法，本发明就是在这样的背景下提出的。
【发明内容】
[0004] 针对上述技术问题，本发明提供了一种简便实用的用于空心砖型蓄热式加热器的 蓄热阵换热设计方法。
[0005] 本发明提供的技术方案为：
[0006] 用于空心砖型蓄热式加热器的蓄热阵换热设计方法，包括下列步骤：
[0007] 步骤一、建立空心砖型蓄热式加热器的蓄热阵的一维换热分析模型，将具有多个 换热通道的蓄热阵简化为当量流通面积圆管模型；
[0008] 步骤二、推导一维换热的控制方程与定解条件，
[0009] 控制方程：
[0010] 其中，0 (X，T)、T(X，T)分别为蓄热阵的温度和气体的温度，X、T分别为轴向位 置和运行时间，n、d分别为蓄热阵的孔隙率和单个换热通道的直径，pb、cb分别为蓄热阵 材料的密度和比热，cp为气体的定压比热，h为对流换热系数，tin为入口处气体温度，妒(工） 为起始时刻蓄热阵的轴向温度分布，单个当量流通面积圆管的流通面积A=MAPU)，p、U和M分别为该当量流通面积圆管内气体的密度、流速和流量；
[0011] 步骤三、利用步骤二中的控制方程构建控制方程的差分格式，控制方程的 差分格式由以下3个公式Eq.l、Eq. 2和Eq. 3组成，其中T=mAT，x=pAx，
0 m，p= 0(mAt，pAx)，Tm，p=T(mAt，pAx):
[0012]
[0013] 步骤四、利用步骤三的控制方程的差分格式求解蓄热阵和气体的温度分布，以蓄 热阵的温度分布代表蓄热阵的换热性能，具体过程为：
[0014] 步骤（1)给定蓄热阵初始温度分珣
[0015] 步骤⑵给定气体入口温度1^。=Tin(p= 0);
[0016] 步骤（3)由Eq. 2知气体的初始温度分布
[0020]步骤（5)由Eq. 3确定蓄热阵温度分布[0021]
[0017]
[0018]
[0019]
[0022] 步骤（6)由Eq. 2确定气体温度】
[0023] 步骤（7)重复步骤（5)与步骤（6)，求出t=At时刻蓄热阵和气体的温度分布 0u'Drn= 1，p彡2,且p为整数）；
[0024] 步骤（8)重复步骤（4)~步骤（7)即可求出所有时刻蓄热阵和气体的温度分布；
[0025] 步骤五、给定蓄热阵的热力学参数和工作条件，重复执行步骤四，求解蓄热阵以及 气体在不同的几何参数时的温度分布，以对比分析蓄热阵在不同的几何参数下的换热性 能，进而确定出具备所须换热性能的蓄热阵的布局，其中，热力学参数包括Pb、Cb和h，工作 条件包括Tin、Cp、p〇)和A，几何参数包括n和d。
[0026] 优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤五 中，重复执行步骤四的过程中，先求解蓄热阵在预热阶段的温度分布，再求解蓄热阵在吹风 阶段的温度分布，并假设蓄热阵在预热阶段的结束时刻的温度分布与吹风阶段的起始时刻 的温度分布相同。
[0027] 优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤五 中，在预热阶段，T=0时识〇0为环境温度，x= 0为蓄热阵的顶端，Tin为预热混合燃气温 度；在吹风阶段：T= 〇时识00为预热结束时刻的蓄热阵温度，X= 〇为蓄热阵的底端，Tin 为待加热的试验气体温度，与环境温度相同。
[0028] 优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤二 中，控制方港
推导过程为：
[0029] 在当量流通面积圆管模型中，气流方向上的微元控制体dx内，在dt时 间内，气体从控制边界吸收的热量为h( 9 -T)(NJrd?dx) ?dt，气体增加的内能为
，气体从控制边界带走的净能量为
蓄热阵含有换热通道的数目为N，Cv为气体的定容比热，则有
[0030]
,
[0033] 由理想气体状态方程知
蓄热器工作过程中，
0忽略不计，又PU=M/A，进一步整理可得：
[0034]
[0035] 优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤二 中，控制方彳
> 的推导过程为：
[0036] 在dT时间内，dx长度蓄热阵减小的内能，
，蓄热阵通 过控制边界向气体传递的能量为h( 0 -T)(N31d?dx) ?dT，则有
[0037]
[0038]
[0039]
[0040] 优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤三 中，利用步骤二中的控制方程构建控制方程的差分格式的具体过程为：
[0041] 用差分代替微分，并令t=mAt、x=pAx，则
由步骤（2)所述的控制方程可得：
[0042]
[0043] 其中，At时间内的平均{
Ax长 度上的平均彳1
代入并整理可得：
[0046]又
中整理可得：
[0047]
[0048] 联立式Eq. 1与Eq. 2可得式Eq. 3 :
[0049]
[0050] 优选的是，所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法中，步骤一 中，假设蓄热阵在横截面上的温度近似一致，每个换热通道为圆柱形通道，且每个换热通道 内气体在横截面上的温度均匀，气体与蓄热阵之间的热量只通过热对流方式传递，忽略蓄 热阵向外的热损失，从而将具有多个换热通道的蓄热阵简化为当量流通面积圆管模型。
[0051] 本发明所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法与现有技术相 比的优势如下：
[0052] (1)本发明思路明确、简便实用，仅需按照预定的步骤，代入蓄热阵的热力学参数、 工作条件参数及几何参数，在几秒内即可获得一组参数对应的蓄热阵的换热性能。
[0053] (2)对比不同几何参数对应的蓄热阵换热性能，即可筛选出较佳的几何布局方案。
[0054] (3)换热设计中结合蓄热阵的工作特点，采用预热与吹风两工作阶段联合分析，吹 风阶段的初始条件摆脱了人为的假设，使计算结果更为精确。
【附图说明】
[0055] 图1为本发明所述的用于空心砖型蓄热式加热器的蓄热阵换热设计方法的流程 图；
[0056] 图2(a)为空心砖型蓄热式加热器的蓄热阵的结构示意图；图2(b)为本发明建立 的当量流通面积圆管模型；
[0057] 图3(a)为设计工况下预热阶段蓄热阵的温度特性曲线，图3(b)为设计工况下吹 风阶段蓄热阵的温度特性曲线，其中，设计工况条件如下：预热混合燃气温度和流量分别为 1200K和lkg/s，环境温度为300K，预热结束时刻蓄热阵底端温度为500K，试验空气流量和 压力分别为l〇kg/s和 12MPa，d= 7mm，q= 0? 3 ;
[0058] 图4(a)为设计工况下吹风阶段出口处气体温度特性曲线，d= 7mm，n取值范 围为0. 2~0. 4,图4(b)为设计工况下吹风阶段出口处气体温度特性曲线，d取值范围为 5~10mm，n取值0.3,其中，设计工况条件如下：预热混合燃气温度和流量分别为1200K和 lkg/s，环境温度为300K，预热结束时刻蓄热阵底端温度为500K，试验空气流量和压力分别 为 10kg/s和 12MPa。
【具体实施方式】
[0059] 下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照说明书文 字能够据以实施。
[0060] 如图1所示，本发明提供一种用于空心砖型蓄热式加热器的蓄热阵换热设计方 法，包括下列步骤：
[0061] 步骤一、建立空心砖型蓄热式加热器的蓄热阵的一维换热分析模型，将单个换热 通道简化为当量流通面积圆管模型。
[0062] 假设空心砖型蓄热式加热器的蓄热阵在横截面上的温度可近似一致，每个换热通 道为圆柱形通道，且每个换热通道的气体在横截面上的温度分布均匀，气体与蓄热阵之间 的热量只通过热对流方式传递，忽略蓄热阵向外的热损失，于是具有多个换热通道的蓄热 阵可简化为当量流通面积圆管模型的一维换热分析模型。当量流通面积圆管内气体的密 度P、流速U、流量M与蓄热阵内气体的对应量相同，当量流通面积圆管的流通面积A=M/ (PU)。若蓄热阵含有换热通道的数目为N，则气体等效流通面积也可表示为A=NJrd2/4。 [0063] 步骤二、推导一维换热的控制方程与定解条件。
[0064] 对于气体流动而言，取图2(b)虚线所示的微元控制体为研宄对象，则dT时 间内，微元控制体内气体内能的增量为
气体从控制边界进 出带走的净能量戈
气体从控制边界吸收的热量为h(0_T) (NJrd?dx) ?dt，由能量守恒知，即气体从控制边界吸收的热量=微元控制体内气体增加 的内能+气体从控制边界带走的净能量，即r