主权项 |
用于空心砖型蓄热式加热器的蓄热阵换热设计方法,其特征在于,包括下列步骤:步骤一、建立空心砖型蓄热式加热器的蓄热阵的一维换热分析模型,将具有多个换热通道的蓄热阵简化为当量流通面积圆管模型;步骤二、推导一维换热的控制方程与定解条件,控制方程为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mfrac><mrow><mo>∂</mo><mi>T</mi></mrow><mrow><mo>∂</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>θ</mi><mo>-</mo><mi>T</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><mo>∂</mo><mi>θ</mi></mrow><mrow><mo>∂</mo><mi>τ</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mn>4</mn><mi>η</mi></mrow><mrow><mi>d</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>η</mi><mo>)</mo></mrow></mrow></mfrac><mfrac><mi>h</mi><mrow><msub><mi>ρ</mi><mi>b</mi></msub><msub><mi>C</mi><mi>b</mi></msub></mrow></mfrac><mrow><mo>(</mo><mi>T</mi><mo>-</mo><mi>θ</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>θ</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>θ</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>θ</mi><mrow><mi>m</mi><mo>,</mo><mi>p</mi></mrow></msub><mo>+</mo><msub><mi>θ</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>θ</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>θ</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>θ</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>θ</mi><mrow><mn>0</mn><mo>,</mo><mi>p</mi></mrow></msub><mo>+</mo><msub><mi>θ</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>θ</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>θ</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>θ</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>θ</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>θ</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>θ</mi><mrow><mn>1</mn><mo>,</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>θ</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。 |