建立温度分布模型应用

2012年青岛理工大学数学建模竞赛承诺书我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C中选择一项填写）：A所属学院和专业、年级（请填写完整的全名）：环境学院环境工程102参赛队员：1.姓名何康康联系电话132754245392.姓名郭皓联系电话187662152753.姓名马西涛联系电话[phone]日期：2012年6月10日题目气温分布问题研究关键词太阳辐射能太阳高度角热辐射定律牛顿冷却公式matlab软件摘要本文通过研究太阳辐射强度对温度的影响，从而对不同地点的温度进行探讨。问题一：本文首先研究不同纬度下的太阳辐射强度的与纬度的关系，推导出地面吸收太阳辐射能的计算公式，并利用热平衡原理以及热辐射相关定律，总结出不同纬度地区春分日的最高温度的计算模型：40scξCcosωcos365n/360cos0.0331τρG100T，考虑海拔因素后，利用matlab软件得到该区域温度分布图（伪彩色图、温度等高线图、温度曲面图）。问题二：利用问题一得到的地面吸收太阳辐射能计算式，在冷空气影响下，考虑对流与辐射散失的热量，利用牛顿冷却公式和热平衡原理建立该情况下的温度模型t-t9.55v3.06t-t0.035cosωcos365n/360cos0.0331τρGwwsc，同样利用matlab软件拟合出该区域的温度分布图。本文建立的温度模型计算简便，图像合理。一模型假设1、本文考虑的气温影响因子只考虑该地区所接收的太阳辐射能；2、本文计算的温度是在晴朗天气得出的。二符号说明太阳入射角i：太阳光线与地面法线之间的夹角；太阳高度角θ：从地面某一观察点像太阳中心作一条射线，该射线在地面上有一投影线，这两线之间夹角；太阳赤纬角σ：午时太阳光线与地球赤道平面的夹角；太阳时角ω：太阳午时0ω，上午为负，下午为正，每小时相当于15度。（其余符号在文中指出）三模型建立与求解问题一空气温度也就是气温，是表示空气冷热程度的物理量。空气中的热量主要来源于太阳辐射，太阳辐射到达地面后，一部分被反射，一部分被地面吸收，使地面增热；地面再通过辐射、传导和对流把热传给空气，这是空气中热量的主要来源。所以说，影响气温的主要因素是太阳辐射强度，太阳光热并不能直接使气温升高。空气中的各种气体直接吸收阳光的热能只占据一小部分，地面吸收的热量占据较大部分。地面吸收了太阳辐射热量后，再通过辐射、对流、乱流等形式向空气中传导，这种向上输送热量的结果，是使气温增高的主要原因。所以说，地面温度的升高，主要是吸收太阳辐射能量的结果。本文从太阳辐射强度这个角度出发研究气温。太阳辐射强度取决于太阳高度角、日照长短和日地距离。地面所接收的能量主要是太阳辐射的能量，到达地面的能量主要取决于太阳高度角θ。太阳高度角愈大，太阳辐射强度愈大。因为同一束光线照射时，直射时，照射面积最小，单位面积所获得的太阳辐射则多；反之，斜射时，照射面积大，单位面积上获得的太阳辐射则少。太阳高度角还与纬度有关，不同纬度地区太阳高度角不同，一日之中，太阳高度角正午大于早晚；夏季大于冬季；低纬地区大于高纬度地区。日地距离是指地球环绕太阳公转时，由于公转轨道呈椭圆形，日地之间的距离则不断改变。地球上获得的太阳辐射强度与日地距离的平方呈反比。地球位于近日点时，获得太阳辐射大于远日点。同时，太阳赤纬角σ也会影响太阳高度角θ的变化。根据资料1，我们查找出太阳入射角i与太阳高度角θ相关函数1icoscosωcosδosδsinδinδsinθ若在春分日，太阳光线直射赤道，即0δ公式变成2cosicosωcossinθ为了计算出日最高温，首先需要了解空气增温的过程。根据空气增温原理4，太阳光照射到地面上，晒热了地面，地面吸收的热再放出去烘热空气，空气的温度主要是间接从地面得到的。虽然在中午的时候，太阳光照射地面最接近直射，地面和空气受热量强，但地面放出的热量，少于太阳所供应的热量，所以此时并非是热的时候。中午以后，地面温度仍能继续升高，一直等到地面放出的热量等于太阳所供应的热量时，地面温度才能升到最高；而近地面气温的升高，必须从地面放出的热量得到，这一切都需要时间，因此夏天最热的时候不是中午，而是14时左右的时候。因此，当气温达到最高时，满足方程3QQ放出吸收在这段时间，热量达到平衡状态，在此之后，由于支出的热量大于吸收的热量，气温开始下降。接下来，我们研究太阳给地面所提供的热量，即吸收Q。单位面积在单位时间内接受的太阳辐射能称之为太阳辐射强度，其单位为2W/m，由于大气层对太阳辐射的衰减作用，到达地面上的太阳辐射能有两种：一种是太阳直接辐射；另一种是由大气、灰尘、水滴等散射而来的，成为散射辐射。在天气晴朗时，打打地面的太阳辐射主要是直接辐射，在阴天时主要是散射辐射。在本文中考虑的是直接辐射。查找相关资料3，大气层外的太阳辐射4360n/3650.033cos1GQscon其中太阳常数2scW/m1367G，n为从元旦算起的天数。直射到地面的太阳辐射能5τQQon1其中τ为大气透过率。水平地面太阳直射辐射6sinθQQ12根据地面所吸收的太阳辐射能7ρQQ2吸收ρ表示地面吸热系数根据71式可得9cosωcos360n/3650.033cos1τρGQ8cosωτρcosQQQscon放出吸收放出根据物体表面辐射热能力的大小与表面性质及表面绝对温度有关，根据斯蒂芬-玻尔兹曼定律210dλ1-eλCdλIE0/λλc-5100λ02积分整理后得：11100TCTδE40400式11成为斯蒂芬-玻尔兹曼定律，它揭示了黑体辐射能力与其表面温度的关系。但在自然界中，黑体是不存在的，它只是用来作为比较的标准。我们把地面看作是灰体，利用克西霍夫定律，得到公式：12100TCAEAE40b00在同一温度下，灰体辐射能力与黑体辐射能力之比，定义为物体的黑度（或称为物体的发射率），用ξ表示：13EEξ0即我们可以得到物体表面温度与辐射热量的关系式：14100TξCE40利用9和14式，把地面散失即放出的热量近似认为等于其辐射的热量，根据此关系，我们可以把公式化简为：15ξCcosωcos365n/360cos0.0331τρG100T40sc利用15式，并查找相关数据，就可以得到地表温度TK其中太阳常数2scW/m1367G；地面热量吸收率5.0ρ；n为从元旦至春分的天数，即80n；太阳时角ω的选取在各地的14时，这里所指的时间为太阳时，即030ω；地面表面黑度ξ，我们选用混凝土的表面黑度88.0ξ；0C为辐射系数，其值为5.6742KmW/；大气透明度τ通过查找相关资料获得，在本文中，我们考虑的是晴朗空气，取0.8τ。根据以上数据及15式，我们只要知道该地的纬度，就可以求出该地区的地表温度C273-Tt01。已知表格数据：表一根据上表提供的纬度，利用公式15计算得到不同纬度下的地表温度，如下表所示：表二纬度N30N32.5N35N37.5N40N42.5N45N47.5N50温度28.626.624.42219.516.713.610.47根据据算得出的结果，我们还绘制了纬度与温度关系的折线图，如下：E90E92.5E95E97.5E100E102.5E105E107.5E110E112.5E115E117.5E120N502680222017801630135011308706859101150970740640N47.518401400150012009001090870670N4510001550150012301070980650N42.5115022202600142012001050530220N40130035002000137011806005N37.5256037003330217013408101200N3542004500330021009603300N32.549003800326031301300230195130N3048002880170013703000330021201250680130400480390很容易看出，随着纬度的增高，温度越来越低，符合实际情况。因为我们所处的位置在对流层下部，对流层下部的温度，随着海拔的增高，所吸收的地面辐射热量越少，空气越稀薄，温度越低。所以，我们考虑海拔对气温的影响)100Ar(-tt1其中，r为垂直降温率，在正常气候条件下，高度每上升100米，气温下降C0.60，所以6.0r；A为海拔。根据上式以及表二提供的数据，我们可以得出不同地区的最高气温，如下表三所示（其中已把海拔换算成温度）表三E90E92.5E95E97.5E100E102.5E105E107.5E110E112.5E115E117.5E120N50-9.08-6.32-3.68-2.78-1.10.221.782.891.540.11.182.563.16N47.5-0.6421.43.253.865.186.38N457.64.34.66.227.187.729.7N42.59.83.381.18.189.510.413.5215.38N4011.7-1.57.511.2812.4215.919.47N37.56.64-0.22.028.9813.9617.1421.2822N35-0.8-2.64.611.818.6422.4224.4N32.5-2.83.87.047.8218.825.2225.4325.82N30-0.211.3218.414.910.68.815.8821.124.5227.8226.225.7226.26下面，我们通过运用matlaB软件，对表三的数据进行编程处理（程序见附录），得到了该区域春分日最高温度分布图，如下所示：图一图一所示的是东经90到120度，北纬30到50区域不同海拔下的温度分布。图一为伪彩色图，颜色越趋近于红色的温度越高，越蓝的区域温度越低。图二如下：如图二所示，为温度的等高线图，比较清晰的分辨出温度的高低。本文主要借助图三进行分析。图三如图三（三维曲面）所示，此图的竖坐标为温度，横纵坐标为纬度和经度。根据图像分析，随着纬度的上升，曲面呈现下降趋势，即温度越低，符合我们的实际情况。在E90到E100，N30到N40区域中，该区域位于青藏高原附近，海拔高，空气稀薄，温度偏低，这从图像可反映出来。在E120，N30区域，海拔较低，处于平原或丘陵地带，温度较高。通过建立三维曲面温度分布图，我们可以容易的看到该区域的温度高低情况。问题二由于冷空气入侵加速了空气流动，当地面吸收太阳所提供的能量等于地面放出能量时，温度达到最高，即由问题一3式可得放出吸收QQ地面放出的热量可以看作是地面散失的热量，即散失Q，而散失的热量在冷空气的影响下，一方面与外界空气进行对流换热，另一方面还要与周围环境物体进行辐射换热，所以把散失的热量分为对流方式散失的热量和辐射方式散失的热量之和，即16QQQ辐射对流散失根据牛顿冷却公式5t-tAaQwfc可得t-taQw117t-ta