发明名称 用于空心砖型蓄热式加热器的蓄热阵换热设计方法
摘要 本发明公开了用于空心砖型蓄热式加热器的蓄热阵换热设计方法,包括:建立蓄热阵的一维换热分析模型,将多个换热通道的蓄热阵简化为当量流通面积圆管模型;推导一维换热的控制方程与定解条件;构建控制方程的差分格式;求解蓄热阵和气体的温度分布;给定蓄热阵的热力学参数和工作条件,求解蓄热阵以及气体在不同的几何参数时的温度分布,以对比分析蓄热阵在不同的几何参数下的换热性能,进而确定出具备所须换热性能的蓄热阵的布局。本发明简便实用,仅需按照预定的步骤,代入蓄热阵的热力学参数、工作条件参数及几何参数,在几秒内即可获得一组参数对应的蓄热阵的换热性能。对比不同几何参数对应的蓄热阵换热性能,即可筛选出较佳的几何布局方案。
申请公布号 CN104992042A 申请公布日期 2015.10.21
申请号 CN201510486214.9 申请日期 2015.08.10
申请人 中国航天空气动力技术研究院 发明人 卢洪波;刘展;郭孝国
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京远大卓悦知识产权代理事务所(普通合伙) 11369 代理人 史霞
主权项 用于空心砖型蓄热式加热器的蓄热阵换热设计方法,其特征在于,包括下列步骤:步骤一、建立空心砖型蓄热式加热器的蓄热阵的一维换热分析模型,将具有多个换热通道的蓄热阵简化为当量流通面积圆管模型;步骤二、推导一维换热的控制方程与定解条件,控制方程为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mfrac><mrow><mo>&part;</mo><mi>T</mi></mrow><mrow><mo>&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>&theta;</mi><mo>-</mo><mi>T</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><mo>&part;</mo><mi>&theta;</mi></mrow><mrow><mo>&part;</mo><mi>&tau;</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mn>4</mn><mi>&eta;</mi></mrow><mrow><mi>d</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&eta;</mi><mo>)</mo></mrow></mrow></mfrac><mfrac><mi>h</mi><mrow><msub><mi>&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>&theta;</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000777758220000011.GIF" wi="654" he="340" /></maths>定解条件为<img file="FDA0000777758220000012.GIF" wi="487" he="174" />其中,θ(x,τ)、T(x,τ)分别为蓄热阵的温度和气体的温度,x、τ分别为轴向位置和运行时间,η、d分别为蓄热阵的孔隙率和单个换热通道的直径,ρ<sub>b</sub>、C<sub>b</sub>分别为蓄热阵材料的密度和比热,C<sub>p</sub>为气体的定压比热,h为对流换热系数,T<sub>in</sub>为入口处气体温度,<img file="FDA0000777758220000015.GIF" wi="113" he="71" />为起始时刻蓄热阵的轴向温度分布,当量流通面积圆管的流通面积A=M/(ρU),ρ、U和M分别为该当量流通面积圆管内气体的密度、流速和流量;步骤三、利用步骤二中的控制方程构建控制方程的差分格式,控制方程的差分格式由以下3个公式Eq.1、Eq.2和Eq.3组成,其中τ=mΔτ,x=pΔx,<img file="FDA0000777758220000013.GIF" wi="820" he="161" />θ<sub>m,p</sub>=θ(mΔτ,pΔx),T<sub>m,p</sub>=T(mΔτ,pΔx):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&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><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><mfrac><mrow><mn>1</mn><mo>-</mo><mi>E</mi></mrow><mrow><mn>1</mn><mo>+</mo><mi>E</mi></mrow></mfrac><msub><mi>&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><msub><mi>&theta;</mi><mrow><mi>m</mi><mo>,</mo><mi>p</mi></mrow></msub><mo>+</mo><msub><mi>&theta;</mi><mrow><mi>m</mi><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><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>&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>&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>&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>]]></math><img file="FDA0000777758220000014.GIF" wi="1748" he="465" /></maths>步骤四、利用步骤三的控制方程的差分格式求解蓄热阵和气体的温度分布,以蓄热阵的温度分布代表蓄热阵的换热性能,具体过程为:步骤(1)给定蓄热阵初始温度分布<img file="FDA0000777758220000025.GIF" wi="540" he="76" />步骤(2)给定气体入口温度T<sub>m,0</sub>=T<sub>in</sub>(p=0);步骤(3)由Eq.2知气体的初始温度分布<maths num="0003" id="cmaths0003"><math><![CDATA[<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>&theta;</mi><mrow><mn>0</mn><mo>,</mo><mi>p</mi></mrow></msub><mo>+</mo><msub><mi>&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><mrow><mo>(</mo><mi>m</mi><mo>=</mo><mn>0</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000777758220000021.GIF" wi="1022" he="141" /></maths>步骤(4)由Eq.1确定蓄热阵温度分布<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&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>&theta;</mi><mrow><mn>0</mn><mo>,</mo><mn>0</mn></mrow></msub><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>]]></math><img file="FDA0000777758220000022.GIF" wi="1054" he="138" /></maths>步骤(5)由Eq.3确定蓄热阵温度分布<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&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>&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>&theta;</mi><mrow><mn>1</mn><mo>,</mo><mn>0</mn></mrow></msub><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>]]></math><img file="FDA0000777758220000023.GIF" wi="1812" he="141" /></maths>步骤(6)由Eq.2确定气体温度<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>T</mi><mrow><mn>1</mn><mo>,</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mi>F</mi><mrow><mn>1</mn><mo>+</mo><mi>F</mi></mrow></mfrac><mrow><mo>(</mo><msub><mi>&theta;</mi><mrow><mn>1</mn><mo>,</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>&theta;</mi><mrow><mn>1</mn><mo>,</mo><mn>0</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>1</mn><mo>,</mo><mn>0</mn></mrow></msub><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>]]></math><img file="FDA0000777758220000024.GIF" wi="1013" he="142" /></maths>步骤(7)重复步骤(5)与步骤(6),求出τ=Δτ时刻蓄热阵和气体的温度分布θ<sub>1,p</sub>、T<sub>1,p</sub>(m=1,p≥2,且p为整数);步骤(8)重复步骤(4)~步骤(7)即可求出所有时刻蓄热阵和气体的温度分布;步骤五、给定蓄热阵的热力学参数和工作条件,重复执行步骤四,求解蓄热阵以及气体在不同的几何参数时的温度分布,以对比分析蓄热阵在不同的几何参数下的换热性能,进而确定出具备所须换热性能的蓄热阵的布局,其中,热力学参数包括ρ<sub>b</sub>、C<sub>b</sub>和h,工作条件包括T<sub>in</sub>、C<sub>p</sub>、<img file="FDA0000777758220000026.GIF" wi="108" he="71" />和A,几何参数包括η和d。
地址 100074 北京市丰台区云岗西路17号