发明名称 一种预测平流层气球的白天和夜晚温度的新方法
摘要 一种预测平流层气球的白天和夜晚温度的新方法,它有五大步骤:步骤一、提出预测平流层气球的白天和夜晚温度的新方法的假设条件;步骤二、根据平流层气球的使用环境和工程热力学的基本知识,确定平流层气球的能量交换机理和各能量交换项;步骤三、针对步骤二提出的各能量交换项分别建立其热分析模型,获得具体各能量交换项具体的参数化表达形式;步骤四、根据能量平衡原理,分别建立平流层气球在白天和夜晚的能量平衡方程;步骤五、将步骤三得到的各能量交换项代入步骤四的能量平衡方程,再利用数值求解方法中的简单迭代法即求得平流层气球在白天和夜晚的温度。本发明仅需要囊布材料参数、悬浮高度和时间,就能得到平流层气球的白天和夜晚温度。
申请公布号 CN103471723A 申请公布日期 2013.12.25
申请号 CN201310406617.9 申请日期 2013.09.09
申请人 北京航空航天大学 发明人 熊峻江;白江波;陈勇;满孜郁
分类号 G01J5/58(2006.01)I 主分类号 G01J5/58(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 1.一种预测平流层气球的白天和夜晚温度的新方法,其特征在于:该方法具体步骤如下:步骤一、提出预测平流层气球的白天和夜晚温度的新方法的假设条件;假设条件包括:(1)当平流层气球受载时,实际的形状类似水滴形,为了简化模型,将平流层气球简化为半径为R的圆形,平流层气球的特征尺寸L用直径来表示;(2)由于平流层气球的囊布非常薄,在微米量级,因此,设囊布厚度方向的温度相同,不考虑囊布内层与外层之间的热传递;(3)不考虑平流层气球内部氦气的损失,设平流层气球的重量是恒定的;(4)由于地球半径R<sub>e</sub>为6371.23km,而太阳距地球的距离为1.5×10<sup>8</sup>km,因此,设太阳光到达地球上空时是平行光;(5)设地球为热平衡体;(6)不考虑太阳辐射强度随太阳辐射角和大气消光系数的变化;(7)不考虑二氧化碳和臭氧红外辐射的影响;步骤二、根据平流层气球的使用环境和工程热力学的基本知识,确定平流层气球的能量交换机理和各能量交换项;平流层气球通过辐射和对流来进行能量交换,包括:平流层气球上半球吸收的太阳直接辐射项Q<sub>1</sub>(t)、平流层气球吸收的太阳散射辐射项Q<sub>2</sub>(t)、平流层气球下半球吸收的地球反射辐射项Q<sub>3</sub>(t)、平流层气球下半球吸收的地球红外辐射项Q<sub>4</sub>(t)、平流层气球对外的红外辐射项Q<sub>5</sub>(t)、平流层气球吸收的内部反射红外辐射项Q<sub>6</sub>(t)、平流层气球与外部大气之间的对流换热项Q<sub>7</sub>(t)和平流层气球与内部气体之间的对流换热项Q<sub>8</sub>(t);步骤三、针对步骤二提出的各能量交换项分别建立其热分析模型,获得具体各能量交换项具体的参数化表达形式;平流层气球的微元的面积表示为<img file="FDA00003792174800011.GIF" wi="1198" he="61" />式中,dA为平流层气球上的微元面积,θ为球坐标系的纬度或太阳辐射角或太阳光与微元法向之间的夹角,<img file="FDA00003792174800012.GIF" wi="42" he="44" />为球坐标系的经度;平流层气球上半球上的微元所吸收的太阳直接辐射为dQ<sub>1</sub>(t)=α<sub>1</sub>τI<sub>0</sub>cosθdA(2)式中,α<sub>1</sub>为平流层气球的囊布对太阳直接辐射的吸收率,τ为大气透明系数,I<sub>0</sub>为太阳常数;平流层气球上的微元所吸收的太阳散射辐射为dQ<sub>2</sub>(t)=α<sub>1</sub>I<sub>1</sub>dA(3)式中,I<sub>1</sub>为太阳散射辐射强度;平流层气球下半球上的微元所吸收的地球反射辐射为<maths num="0001"><![CDATA[<math><mrow><mi>d</mi><msub><mi>Q</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><msub><mi>&rho;</mi><mi>g</mi></msub><mi>&tau;</mi><msub><mi>I</mi><mn>0</mn></msub><msup><mrow><mo>(</mo><mfrac><msub><mi>R</mi><mi>e</mi></msub><mrow><msub><mi>R</mi><mi>e</mi></msub><mo>+</mo><mi>H</mi></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mi>cos</mi><mi></mi><mi>&theta;</mi><mi>cos</mi><mi>&phi;dA</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,ρ<sub>g</sub>为地面反射率平均值,R<sub>e</sub>为地球半径,H平流层气球的悬浮高度,φ为太阳光方向与地球中心和平流层气球中心连线方向之间的夹角;平流层气球下半球上的微元所吸收的地球红外辐射为<maths num="0002"><![CDATA[<math><mrow><mi>d</mi><msub><mi>Q</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&alpha;</mi><mn>2</mn></msub><mi>&tau;</mi><msub><mi>I</mi><mn>0</mn></msub><msubsup><mi>R</mi><mi>e</mi><mn>2</mn></msubsup><mfrac><mrow><mn>1</mn><mo>-</mo><msub><mi>&rho;</mi><mi>g</mi></msub></mrow><mrow><mn>4</mn><msup><mrow><mo>(</mo><msub><mi>R</mi><mi>e</mi></msub><mo>+</mo><mi>H</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mi>cos</mi><mi>&theta;dA</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,α<sub>2</sub>为平流层气球的囊布对红外辐射的吸收率;平流层气球上的微元对外的红外辐射为<maths num="0003"><![CDATA[<math><mrow><mi>d</mi><msub><mi>Q</mi><mn>5</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mn>2</mn><mi>&epsiv;&sigma;</mi><msubsup><mi>T</mi><mn>1</mn><mn>4</mn></msubsup><mi>dA</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,ε为平流层气球的囊布对红外辐射的发射率,σ为斯蒂芬-波尔兹曼常数,T<sub>1</sub>为平流层气球的囊布温度;平流层气球上的微元所吸收的内部反射红外辐射为<maths num="0004"><![CDATA[<math><mrow><mi>d</mi><msub><mi>Q</mi><mn>6</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>&alpha;</mi><mn>2</mn></msub><mi>&epsiv;&sigma;</mi><msubsup><mi>T</mi><mn>1</mn><mn>4</mn></msubsup></mrow><mrow><mn>1</mn><mo>-</mo><mi>r</mi></mrow></mfrac><mi>dA</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,r为平流层气球的反射率;平流层气球上的微元与外部大气之间的对流换热为dQ<sub>7</sub>(t)=h<sub>1</sub>(T<sub>1</sub>-T<sub>2</sub>)dA(8)式中,h<sub>1</sub>为平流层气球与外部大气之间的对流换热系数,T<sub>2</sub>为外部大气的温度;平流层气球上的微元与内部氦气之间的对流换热为dQ<sub>8</sub>(t)=h<sub>2</sub>(T<sub>1</sub>-T<sub>3</sub>)dA(9)式中,h<sub>2</sub>为平流层气球与内部氦气之间的对流换热系数,T<sub>3</sub>为内部氦气的温度;式(3)中的太阳散射辐射强度I<sub>1</sub>表示为<maths num="0005"><![CDATA[<math><mrow><msub><mi>I</mi><mn>1</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>I</mi><mn>0</mn></msub><mi>cos</mi><mi>&theta;</mi><msup><mi>cos</mi><mn>2</mn></msup><mfrac><mi>&theta;</mi><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(4)中的太阳光方向与地球中心和平流层气球中心连线方向之间的夹角φ表示为<maths num="0006"><![CDATA[<math><mrow><mi>&phi;</mi><mo>=</mo><mfrac><mrow><mi>&pi;</mi><mo>|</mo><mn>12</mn><mo>-</mo><mi>t</mi><mo>|</mo></mrow><mn>12</mn></mfrac><mo>,</mo><mrow><mo>(</mo><mn>6</mn><mi>a</mi><mo>.</mo><mi>m</mi><mo>.</mo><mo>&le;</mo><mi>t</mi><mo>&le;</mo><mn>18</mn><mi>p</mi><mo>.</mo><mi>m</mi><mo>.</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,t表示时间,从早上6点到下午18点为一个白天的周期;式(8)和式(9)中的h<sub>1</sub>和h<sub>2</sub>分别表示为<maths num="0007"><![CDATA[<math><mrow><msub><mi>h</mi><mn>1</mn></msub><mo>=</mo><mi>N</mi><msub><mi>u</mi><mn>1</mn></msub><mfrac><msub><mi>&lambda;</mi><mn>1</mn></msub><mi>L</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msub><mi>h</mi><mn>2</mn></msub><mo>=</mo><mi>N</mi><msub><mi>u</mi><mn>1</mn></msub><mfrac><msub><mi>&lambda;</mi><mn>1</mn></msub><mi>L</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,Nu<sub>1</sub>为强迫对流的努塞尔数,λ<sub>1</sub>为强迫对流的热传导系数,Nu<sub>2</sub>为自由对流的努塞尔数,λ<sub>2</sub>为自由对流的热传导系数;在平流层气球的热分析中,太阳常数I<sub>0</sub>为1367W/m<sup>2</sup>,大气透明系数τ为0.9,地面反射率平均值ρ<sub>g</sub>为0.3,斯蒂芬-波尔兹曼常数σ为5.6696×10<sup>-8</sup>W(m<sup>2</sup>·K<sup>4</sup>),平流层气球的囊布对红外辐射的发射率与吸收率相等,即ε=α<sub>2</sub>;分别对式(2)至式(9)进行积分,并代入相关参数,即得到各能量交换项的总量:<maths num="0009"><![CDATA[<math><mrow><msub><mi>Q</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1230.3</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>,</mo><mrow><mo>(</mo><mn>6</mn><mi>a</mi><mo>.</mo><mi>m</mi><mo>.</mo><mo>&le;</mo><mi>t</mi><mo>&le;</mo><mn>18</mn><mi>p</mi><mo>.</mo><mi>m</mi><mo>.</mo><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0010"><![CDATA[<math><mrow><msub><mi>Q</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1139.2</mn><mi>&pi;</mi><msub><mi>&alpha;</mi><mn>1</mn></msub><msup><mi>R</mi><mn>2</mn></msup><mo>,</mo><mrow><mo>(</mo><mn>6</mn><mi>a</mi><mo>.</mo><mi>m</mi><mo>.</mo><mo>&le;</mo><mi>t</mi><mo>&le;</mo><mn>18</mn><mi>p</mi><mo>.</mo><mi>m</mi><mo>.</mo><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0011"><![CDATA[<math><mrow><msub><mi>Q</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>366.8</mn><mi>&pi;</mi><msub><mrow><msup><mi>R</mi><mn>2</mn></msup><mi>&alpha;</mi></mrow><mn>1</mn></msub><msup><mrow><mo>(</mo><mfrac><msub><mi>R</mi><mi>e</mi></msub><mrow><msub><mi>R</mi><mi>e</mi></msub><mo>+</mo><mi>H</mi></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mi>cos</mi><mfrac><mrow><mi>&pi;</mi><mo>|</mo><mn>12</mn><mo>-</mo><mi>t</mi><mo>|</mo></mrow><mn>12</mn></mfrac><mo>,</mo><mrow><mo>(</mo><mn>6</mn><mi>a</mi><mo>.</mo><mi>m</mi><mo>.</mo><mo>&le;</mo><mi>t</mi><mo>&le;</mo><mn>18</mn><mi>p</mi><mo>.</mo><mi>m</mi><mo>.</mo><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow></math>]]></maths><img file="FDA00003792174800036.GIF" wi="1839" he="175" /><maths num="0012"><![CDATA[<math><mrow><msub><mi>Q</mi><mn>5</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Integral;</mo><mi>A</mi></munder><mn>2</mn><mi>&epsiv;&sigma;</mi><msubsup><mi>T</mi><mn>1</mn><mn>4</mn></msubsup><mi>dA</mi><mo>=</mo><mn>2</mn><mi>A&epsiv;A&sigma;</mi><msubsup><mi>T</mi><mn>1</mn><mn>4</mn></msubsup><mo>=</mo><mn>45.4</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>&alpha;</mi><mn>2</mn></msub><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>8</mn></mrow></msup><msubsup><mi>T</mi><mn>1</mn><mn>4</mn></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0013"><![CDATA[<math><mrow><msub><mi>Q</mi><mn>6</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Integral;</mo><mi>A</mi></munder><mfrac><mrow><msub><mi>&alpha;</mi><mn>2</mn></msub><mi>&epsiv;&sigma;</mi><msubsup><mi>T</mi><mn>1</mn><mn>4</mn></msubsup></mrow><mrow><mn>1</mn><mo>-</mo><mi>r</mi></mrow></mfrac><mi>dA</mi><mo>=</mo><mfrac><mrow><mi>A</mi><msub><mi>&alpha;</mi><mn>2</mn></msub><mi>&epsiv;&sigma;</mi><msubsup><mi>T</mi><mn>1</mn><mn>4</mn></msubsup></mrow><mrow><mn>1</mn><mo>-</mo><mi>r</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mn>22.7</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>&alpha;</mi><mn>2</mn></msub><mi>&epsiv;</mi><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>8</mn></mrow></msup><msubsup><mi>T</mi><mn>1</mn><mn>4</mn></msubsup></mrow><mrow><mn>1</mn><mo>-</mo><mi>r</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>19</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0014"><![CDATA[<math><mrow><msub><mi>Q</mi><mn>7</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Integral;</mo><mi>A</mi></munder><msub><mi>h</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>-</mo><msub><mi>T</mi><mn>2</mn></msub><mo>)</mo></mrow><mi>dA</mi><mo>=</mo><mi>A</mi><msub><mi>h</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>-</mo><msub><mi>T</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><mn>4</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>h</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>-</mo><msub><mi>T</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>20</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0015"><![CDATA[<math><mrow><msub><mi>Q</mi><mn>8</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Integral;</mo><mi>A</mi></munder><msub><mi>h</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>-</mo><msub><mi>T</mi><mn>3</mn></msub><mo>)</mo></mrow><mi>dA</mi><mo>=</mo><mi>A</mi><msub><mi>h</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>-</mo><msub><mi>T</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>=</mo><mn>4</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>h</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>-</mo><msub><mi>T</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>21</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,步骤三中所述的微元是指平流层气球上截取的微小单元,面积无限小;步骤四、根据能量平衡原理,分别建立平流层气球在白天和夜晚的能量平衡方程;对于白天,平流层气球的囊布的总能量为△Q<sub>film</sub>(t)=Q<sub>1</sub>(t)+Q<sub>2</sub>(t)+Q<sub>3</sub>(t)+Q<sub>4</sub>(t)-Q<sub>5</sub>(t)+Q<sub>6</sub>(t)-Q<sub>7</sub>(t)-Q<sub>8</sub>(t)(22)平流层气球的内部氦气的总能量为△Q<sub>gas</sub>(t)=Q<sub>8</sub>(t)(23)当平流层气球的能量交换达到平衡状态时,△Q<sub>film</sub>(t)=0,△Q<sub>gas</sub>(t)=0,代入式(22)和式(23),得到关于平流层气球的囊布和内部氦气在白天的能量平衡方程:<maths num="0016"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>Q</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>Q</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>Q</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>Q</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Q</mi><mn>5</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>Q</mi><mn>6</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Q</mi><mn>7</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Q</mi><mn>8</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mn>8</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>24</mn><mo>)</mo></mrow></mrow></math>]]></maths>对于夜晚,没有太阳直接辐射、太阳散射辐射和地球反射辐射,即Q<sub>1</sub>(t)=0,Q<sub>2</sub>(t)=0,Q<sub>3</sub>(t)=0,因此,平流层气球的囊布和内部氦气在夜晚的能量平衡方程为<maths num="0017"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>Q</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Q</mi><mn>5</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>Q</mi><mn>6</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Q</mi><mn>7</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Q</mi><mn>8</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mn>8</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>25</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤五、将步骤三得到的各能量交换项代入步骤四的能量平衡方程,再利用数值求解方法中的简单迭代法即求得平流层气球在白天和夜晚的温度;将式(13)至式(21)分别代入到式(24)和式(25)中,得<maths num="0018"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mrow><mo>(</mo><mfrac><mrow><mn>22.7</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>&alpha;</mi><mn>2</mn></msub><mi>&epsiv;</mi><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>8</mn></mrow></msup></mrow><mrow><mn>1</mn><mo>-</mo><mi>r</mi></mrow></mfrac><mo>-</mo><mn>45.4</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>&alpha;</mi><mn>2</mn></msub><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>8</mn></mrow></msup><mo>)</mo></mrow><msubsup><mi>T</mi><mn>1</mn><mn>4</mn></msubsup><mo>-</mo><mn>4</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>h</mi><mn>1</mn></msub><msub><mi>T</mi><mn>1</mn></msub><mo>+</mo><mo>[</mo><mn>366.8</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>&alpha;</mi><mn>1</mn></msub><msup><mrow><mo>(</mo><mfrac><msub><mi>R</mi><mi>e</mi></msub><mrow><msub><mi>R</mi><mi>e</mi></msub><mo>+</mo><mi>H</mi></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mi>cos</mi><mfrac><mrow><mi>&pi;</mi><mo>|</mo><mn>12</mn><mo>-</mo><mi>t</mi><mo>|</mo></mrow><mn>12</mn></mfrac></mtd></mtr><mtr><mtd><mo>+</mo><mn>2369.5</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>+</mo><mn>214.0</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>&alpha;</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><mfrac><msub><mi>R</mi><mi>e</mi></msub><mrow><msub><mi>R</mi><mi>e</mi></msub><mo>+</mo><mi>H</mi></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>h</mi><mn>1</mn></msub><msub><mi>T</mi><mn>2</mn></msub><mo>]</mo><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>T</mi><mn>1</mn></msub><mo>-</mo><msub><mi>T</mi><mn>3</mn></msub><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced></math>]]></maths>(26)<maths num="0019"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mrow><mo>(</mo><mfrac><mrow><mn>22.7</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>&alpha;</mi><mn>2</mn></msub><mi>&epsiv;</mi><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>8</mn></mrow></msup></mrow><mrow><mn>1</mn><mo>-</mo><mi>r</mi></mrow></mfrac><mo>-</mo><mn>45.4</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>&alpha;</mi><mn>2</mn></msub><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>8</mn></mrow></msup><mo>)</mo></mrow><msubsup><mi>T</mi><mn>1</mn><mn>4</mn></msubsup><mo>-</mo><mn>4</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>h</mi><mn>1</mn></msub><msub><mi>T</mi><mn>1</mn></msub><mo>+</mo><mo>[</mo><mn>214.0</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>&alpha;</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><mfrac><msub><mi>R</mi><mi>e</mi></msub><mrow><msub><mi>R</mi><mi>e</mi></msub><mo>+</mo><mi>H</mi></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup><msub><mi>h</mi><mn>1</mn></msub><msub><mi>T</mi><mn>2</mn></msub><mo>]</mo><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>T</mi><mn>1</mn></msub><mo>-</mo><msub><mi>T</mi><mn>3</mn></msub><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced></math>]]></maths>(27)由(26)和式(27)可知,当平流层气球的能量交换达到平衡状态时,白天和夜晚的平流层气球的囊布温度与内部氦气温度相等;通过数值求解方法中的简单迭代法对方程(26)和(27)分别进行求解,即得到平流层气球在白天和夜晚的囊布温度与内部氦气温度。
地址 100191 北京市海淀区学院路37号
您可能感兴趣的专利