主权项 |
一种用于大口径太阳望远镜主镜镜面温度测量的非接触式测量装置,其特征在于:由主镜面板(1),空气喷管(2),环境温度传感器(3),空气射流温度传感器(4),总辐射传感器(5),风速仪(6)和数据处理模块(7)七个部分组成;其中,环境温度传感器(3)安装在主镜附近,用于监测环境温度;空气射流温度传感器(4)安装在空气喷管(2)出口,用于监测空气射流温度;总辐射传感器(5)安装在望远镜机架上,并正对太阳,用于检测太阳辐照度;风速仪(6)安装在主镜附近,用于检测望远镜周边风速;数据处理模块(7)与上述所有传感器相连,用于对上述所有传感器实时数据的记录、处理与显示;该装置的数据处理过程如下所述:首先,定义参数;环境参数定义:环境温度传感器(3)在第i时刻测得的环境温度为<img file="FDA0000744007990000011.GIF" wi="89" he="81" />(单位:℃);空气射流温度传感器在第i时刻测得的空气射流温度为<img file="FDA0000744007990000012.GIF" wi="64" he="81" />(单位:℃);总辐射传感器(5)在第i时刻测得的太阳辐照度为q<sup>i</sup>(单位:W/m<sup>2</sup>),风速仪(6)在第i时刻测得的环境空气流速为<img file="FDA0000744007990000013.GIF" wi="69" he="81" />(单位:m/s);主镜面板参数定义:面板材料热传导系数、比热容和密度分别为:k(单位:W/m·℃),C<sub>p</sub>(单位:J/Kg·℃),ρ(单位:Kg/m<sup>3</sup>),主镜厚度L<sub>thickness</sub>(单位:m),空间域步长<img file="FDA0000744007990000014.GIF" wi="278" he="133" />(N为空间离散化点数),时间域步长△t(单位:s);环境空气参数定义:主镜直径L(单位:m),环境空气动力学黏度系数ν<sub>1</sub>(查表,单位:m<sup>2</sup>·s),环境空气热传导系数λ<sub>1</sub>(查表,单位:W/m·℃),环境空气普朗特数Pr<sub>1</sub>(查表,无量纲数);射流空气参数定义:空气喷管(2)内径D(单位:m),喷管出口与镜面面板间距H(单位:m),射流区半径r(单位:m),空气射流动力学黏度系数ν<sub>2</sub>(查表,位:m<sup>2</sup>·s),空气射流热传导系数λ<sub>2</sub>(查表,单位:W/m·℃),空气射流普朗特数Pr<sub>2</sub>(查表,无量纲数);然后,计算热传导参数;主镜镜面对流传热系数h<sub>n</sub>计算:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>Re</mi><mi>L</mi><mi>i</mi></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>u</mi><mi>L</mi><mi>i</mi></msubsup><mi>L</mi></mrow><msub><mi>ν</mi><mn>1</mn></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000744007990000021.GIF" wi="1169" he="147" /></maths>当Re<sub>L</sub>≤5×10<sup>5</sup>时,使用方程(2);当Re<sub>L</sub>>5×10<sup>5</sup>时,使用方程(3);<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>h</mi><mi>n</mi></msub><mo>=</mo><mn>0.664</mn><msup><mrow><mo>(</mo><msubsup><mi>Re</mi><mi>L</mi><mi>i</mi></msubsup><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><msup><msub><mi>Pr</mi><mn>1</mn></msub><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mfrac><msub><mi>λ</mi><mn>1</mn></msub><mi>L</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000744007990000022.GIF" wi="1307" he="133" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>h</mi><mi>n</mi></msub><mo>=</mo><mn>0.037</mn><mo>[</mo><msup><mrow><mo>(</mo><msubsup><mi>Re</mi><mi>L</mi><mi>i</mi></msubsup><mo>)</mo></mrow><mrow><mn>4</mn><mo>/</mo><mn>5</mn></mrow></msup><mo>-</mo><mn>871</mn><mo>]</mo><msup><msub><mi>Pr</mi><mn>1</mn></msub><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mfrac><msub><mi>λ</mi><mn>1</mn></msub><mi>L</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000744007990000023.GIF" wi="1396" he="128" /></maths>主镜面板背侧对流传热系数h<sub>C</sub>计算:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>h</mi><mi>c</mi></msub><mo>=</mo><mn>2</mn><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>u</mi><mi>D</mi></msub><mi>D</mi></mrow><msub><mi>ν</mi><mn>2</mn></msub></mfrac><mo>)</mo></mrow><mn>0.5</mn></msup><msup><msub><mi>Pr</mi><mn>2</mn></msub><mn>0.42</mn></msup><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mn>0.005</mn><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>u</mi><mi>D</mi></msub><mi>D</mi></mrow><msub><mi>ν</mi><mn>2</mn></msub></mfrac><mo>)</mo></mrow><mn>0.5</mn></msup></mrow><mo>)</mo></mrow><mn>0.5</mn></msup><mfrac><msub><mi>λ</mi><mn>2</mn></msub><mi>D</mi></mfrac><mfrac><mrow><mn>1</mn><mo>-</mo><mn>1.1</mn><mi>D</mi><mo>/</mo><mi>r</mi></mrow><mrow><mn>1</mn><mo>+</mo><mn>0.1</mn><mrow><mo>(</mo><mi>H</mi><mo>/</mo><mi>D</mi><mo>-</mo><mn>6</mn><mo>)</mo></mrow><mi>D</mi><mo>/</mo><mi>r</mi></mrow></mfrac><mfrac><mi>D</mi><mi>r</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000744007990000024.GIF" wi="1644" he="208" /></maths>最后,将上述参数代入递推方程:<img file="FDA0000744007990000025.GIF" wi="1923" he="846" />其中,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mi>τ</mi><mo>=</mo><mfrac><mrow><mi>k</mi><mi>Δ</mi><mi>t</mi></mrow><mrow><mi>ρ</mi><msub><mi>C</mi><mi>p</mi></msub><mi></mi><msup><mi>Δx</mi><mn>2</mn></msup></mrow></mfrac><mo>,</mo><mi>B</mi><mo>=</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mn>2</mn><mi>τ</mi><mo>)</mo></mrow><mo>,</mo><mi>C</mi><mo>=</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mn>2</mn><mi>τ</mi><mo>+</mo><mn>2</mn><mi>τ</mi><mfrac><mrow><mi>Δ</mi><mi>x</mi><mo>·</mo><msub><mi>h</mi><mi>C</mi></msub></mrow><mi>k</mi></mfrac><mo>)</mo></mrow><mo>,</mo><mi>D</mi><mo>=</mo><mn>2</mn><mi>τ</mi><mfrac><mrow><mi>Δ</mi><mi>x</mi><mo>·</mo><msub><mi>h</mi><mi>C</mi></msub></mrow><mi>k</mi></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000744007990000026.GIF" wi="1695" he="154" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>E</mi><mo>=</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mn>2</mn><mi>τ</mi><mo>+</mo><mn>2</mn><mi>τ</mi><mfrac><mrow><mi>Δ</mi><mi>x</mi><mo>·</mo><msub><mi>h</mi><mi>n</mi></msub></mrow><mi>k</mi></mfrac><mo>)</mo></mrow><mo>,</mo><mi>F</mi><mo>=</mo><mn>2</mn><mi>τ</mi><mfrac><mrow><mi>Δ</mi><mi>x</mi><mo>·</mo><msub><mi>h</mi><mi>n</mi></msub></mrow><mi>k</mi></mfrac><mo>,</mo><mi>G</mi><mo>=</mo><mn>2</mn><mi>τ</mi><mfrac><mrow><mi>Δ</mi><mi>x</mi></mrow><mi>k</mi></mfrac><mo>,</mo><mi>τ</mi><mo>=</mo><mfrac><mi>k</mi><mrow><msub><mi>ρC</mi><mi>p</mi></msub></mrow></mfrac><mo>×</mo><mi>Δ</mi><mi>t</mi><mo>/</mo><msup><mi>Δx</mi><mn>2</mn></msup><mo>,</mo><mi>Δ</mi><mi>x</mi><mo>=</mo><mfrac><msub><mi>L</mi><mrow><mi>t</mi><mi>h</mi><mi>i</mi><mi>c</mi><mi>k</mi><mi>n</mi><mi>e</mi><mi>s</mi><mi>s</mi></mrow></msub><mi>N</mi></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000744007990000027.GIF" wi="1922" he="157" /></maths>经过一夜恒温,认为望远镜镜体各处温度与环境温度相同,可作为初始条件,通过递推公式(5),即可求得每一时刻主镜镜面温度<img file="FDA0000744007990000028.GIF" wi="118" he="82" /> |