《建立温度分布模型应用.doc》由会员分享，可在线阅读，更多相关《建立温度分布模型应用.doc（15页珍藏版）》请在三一文库上搜索。
1、2012年青岛理工大学数学建模竞赛承 诺 书我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C中选择一项填写）： A 所属学院和专业、年级（请填写完整的全名）： 环境学院 环境工程102 参赛队员：1. 姓
2、名 何康康 联系电话 13275424539 2. 姓名 郭皓 联系电话 18766215275 3. 姓名 马西涛 联系电话 13153242178 日期： 2012 年 6 月 10 日 题 目 气温分布问题研究 关 键 词 太阳辐射能 太阳高度角 热辐射定律 牛顿冷却公式 matlab软件 摘要 本文通过研究太阳辐射强度对温度的影响，从而对不同地点的温度进行探讨。 问题一：本文首先研究不同纬度下的太阳辐射强度的与纬度的关系，推导出地面吸收太阳辐射能的计算公式，并利用热平衡原理以及热辐射相关定律，总结出不同纬度地区春分日的最高温度的计算模型：，考虑海拔因素后，利用matlab软件得到该区域
3、温度分布图（伪彩色图、温度等高线图、温度曲面图）。 问题二：利用问题一得到的地面吸收太阳辐射能计算式，在冷空气影响下，考虑对流与辐射散失的热量，利用牛顿冷却公式和热平衡原理建立该情况下的温度模型，同样利用matlab软件拟合出该区域的温度分布图。 本文建立的温度模型计算简便，图像合理。一 模型假设1、 本文考虑的气温影响因子只考虑该地区所接收的太阳辐射能；2、 本文计算的温度是在晴朗天气得出的。二 符号说明太阳入射角：太阳光线与地面法线之间的夹角；太阳高度角：从地面某一观察点像太阳中心作一条射线，该射线在地面上有一投影线，这两线之间夹角；太阳赤纬角：午时太阳光线与地球赤道平面的夹角；太阳时角：
4、太阳午时，上午为负，下午为正，每小时相当于15度。（其余符号在文中指出）三 模型建立与求解问题一 空气温度也就是气温，是表示空气冷热程度的物理量。空气中的热量主要来源于太阳辐射，太阳辐射到达地面后，一部分被反射，一部分被地面吸收，使地面增热；地面再通过辐射、传导和对流把热传给空气，这是空气中热量的主要来源。所以说，影响气温的主要因素是太阳辐射强度，太阳光热并不能直接使气温升高。空气中的各种气体直接吸收阳光的热能只占据一小部分，地面吸收的热量占据较大部分。地面吸收了太阳辐射热量后，再通过辐射、对流、乱流等形式向空气中传导，这种向上输送热量的结果，是使气温增高的主要原因。所以说，地面温度的升高，主
5、要是吸收太阳辐射能量的结果。 本文从太阳辐射强度这个角度出发研究气温。太阳辐射强度取决于太阳高度角、日照长短和日地距离。地面所接收的能量主要是太阳辐射的能量，到达地面的能量主要取决于太阳高度角。 太阳高度角愈大，太阳辐射强度愈大。因为同一束光线照射时，直射时，照射面积最小，单位面积所获得的太阳辐射则多；反之，斜射时，照射面积大，单位面积上获得的太阳辐射则少。太阳高度角还与纬度有关，不同纬度地区太阳高度角不同，一日之中，太阳高度角正午大于早晚；夏季大于冬季；低纬地区大于高纬度地区。日地距离是指地球环绕太阳公转时，由于公转轨道呈椭圆形，日地之间的距离则不断改变。 地球上获得的太阳辐射强度与日地距离
6、的平方呈反比。地球位于近日点时，获得太阳辐射大于远日点。同时，太阳赤纬角也会影响太阳高度角的变化。根据资料，我们查找出太阳入射角与太阳高度角相关函数若在春分日，太阳光线直射赤道，即 公式变成 为了计算出日最高温，首先需要了解空气增温的过程。根据空气增温原理，太阳光照射到地面上，晒热了地面，地面吸收的热再放出去烘热空气，空气的温度主要是间接从地面得到的。虽然在中午的时候，太阳光照射地面最接近直射，地面和空气受热量强，但地面放出的热量，少于太阳所供应的热量，所以此时并非是热的时候。中午以后，地面温度仍能继续升高，一直等到地面放出的热量等于太阳所供应的热量时，地面温度才能升到最高；而近地面气温的升高
7、，必须从地面放出的热量得到，这一切都需要时间，因此夏天最热的时候不是中午，而是14时左右的时候。因此，当气温达到最高时，满足方程在这段时间，热量达到平衡状态，在此之后，由于支出的热量大于吸收的热量，气温开始下降。 接下来，我们研究太阳给地面所提供的热量，即。单位面积在单位时间内接受的太阳辐射能称之为太阳辐射强度，其单位为，由于大气层对太阳辐射的衰减作用，到达地面上的太阳辐射能有两种：一种是太阳直接辐射；另一种是由大气、灰尘、水滴等散射而来的，成为散射辐射。在天气晴朗时，打打地面的太阳辐射主要是直接辐射，在阴天时主要是散射辐射。在本文中考虑的是直接辐射。查找相关资料，大气层外的太阳辐射其中太阳常
8、数，为从元旦算起的天数。直射到地面的太阳辐射能其中为大气透过率。水平地面太阳直射辐射根据地面所吸收的太阳辐射能表示地面吸热系数根据式可得根据物体表面辐射热能力的大小与表面性质及表面绝对温度有关，根据斯蒂芬-玻尔兹曼定律积分整理后得：式成为斯蒂芬-玻尔兹曼定律，它揭示了黑体辐射能力与其表面温度的关系。但在自然界中，黑体是不存在的，它只是用来作为比较的标准。我们把地面看作是灰体，利用克西霍夫定律，得到公式：在同一温度下，灰体辐射能力与黑体辐射能力之比，定义为物体的黑度（或称为物体的发射率），用表示：即我们可以得到物体表面温度与辐射热量的关系式：利用和式，把地面散失即放出的热量近似认为等于其辐射的热
9、量，根据此关系，我们可以把公式化简为： 利用式，并查找相关数据，就可以得到地表温度其中太阳常数；地面热量吸收率；为从元旦至春分的天数，即；太阳时角的选取在各地的14时，这里所指的时间为太阳时，即；地面表面黑度，我们选用混凝土的表面黑度；为辐射系数，其值为5.67；大气透明度通过查找相关资料获得，在本文中，我们考虑的是晴朗空气，取。 根据以上数据及式，我们只要知道该地的纬度，就可以求出该地区的地表温度。已知表格数据：表一E90E92.5E95E97.5E100E102.5E105E107.5E110E112.5E115E117.5E120N5026802220178016301350113087
10、06859101150970740640N47.518401400150012009001090870670N4510001550150012301070980650N42.5115022202600142012001050530220N40130035002000137011806005N37.5256037003330217013408101200N3542004500330021009603300N32.549003800326031301300230195130N3048002880170013703000330021201250680130400480390根据上表提供的纬度，利用公
11、式计算得到不同纬度下的地表温度，如下表所示：表二纬度N30N32.5N35N37.5N40N42.5N45N47.5N50温度28.626.624.42219.516.713.610.47根据据算得出的结果，我们还绘制了纬度与温度关系的折线图，如下：很容易看出，随着纬度的增高，温度越来越低，符合实际情况。因为我们所处的位置在对流层下部，对流层下部的温度，随着海拔的增高，所吸收的地面辐射热量越少，空气越稀薄，温度越低。所以，我们考虑海拔对气温的影响其中，为垂直降温率，在正常气候条件下，高度每上升100米，气温下降，所以；为海拔。根据上式以及表二提供的数据，我们可以得出不同地区的最高气温，如下表三
12、所示（其中已把海拔换算成温度）表三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.2822N
13、35-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区域不同海拔下的温度分布。图一为伪彩色图，颜色越趋近于红色的温度越高，越蓝的区域温度越低。图二如下： 如图二所示，为温度的等高线图，比较清晰的分辨出温度的高低。本文主要借
14、助图三进行分析。图三 如图三（三维曲面）所示，此图的竖坐标为温度，横纵坐标为纬度和经度。根据图像分析，随着纬度的上升，曲面呈现下降趋势，即温度越低，符合我们的实际情况。在E90到E100，N30到N40区域中，该区域位于青藏高原附近，海拔高，空气稀薄，温度偏低，这从图像可反映出来。在E120，N30区域，海拔较低，处于平原或丘陵地带，温度较高。通过建立三维曲面温度分布图，我们可以容易的看到该区域的温度高低情况。问题二 由于冷空气入侵加速了空气流动，当地面吸收太阳所提供的能量等于地面放出能量时，温度达到最高，即由问题一式可得 地面放出的热量可以看作是地面散失的热量，即，而散失的热量在冷空气的影响
15、下，一方面与外界空气进行对流换热，另一方面还要与周围环境物体进行辐射换热，所以把散失的热量分为对流方式散失的热量和辐射方式散失的热量之和，即根据牛顿冷却公式可得 为对流传热系数，为辐射传热系数，单位为；为以对流方式散失的热量，为以辐射方式散失的热量，单位为；为地表温度，为空气温度。综合式，可得与问题一式相结合，可得查找相关文献结合式等式左边的未知参数与式相同，等式右边参数,根据已知的纬度，利用matlab软件求解上述方程。可求出冷空气入侵后的地表温度。纬度N30N32.5N35N37.5N40N42.5N45N47.5N50温度4.54.13.73.22.92.41.91.40.8考虑海拔对气
16、温的影响，为垂直降温率E90E92.5E95E97.5E100E102.5E105E107.5E110E112.5E115E117.5E120N50-15.82-12.52-9.88-2.98-10-5.98-4.42-3.31-4.66-6.1-5.02-3.64-3.04N47.5-9.64-7-7.6-5.8-4-5.14-3.82-2.62N45-4.1-7.4-7.1-5.48-4.52-3.98-2N42.5-4.5-10.92-13.2-6.12-4.8-3.9-0.781.08N40-4.9-18.1-9.1-5.32-4.18-0.72.87N37.5-12.16-19-16
17、.78-9.82-4.84-1.662.483.2N35-21.5-23.3-16.1-8.9-2.061.723.7N32.5-25.3-18.7-15.46-14.68-3.72.722.933.32N30-24.3-12.78-5.7-3.72-13.5-15.3-8.22-30.423.722.11.622.16 根据上表所得数据，利用matlab软件，编写程序模拟出冷空气入侵后的气温分布图，如下所示：伪彩色图，温度等高线图，三维曲面图：四 模型结果分析 通过比较问题一与问题二所得到的气温分布图，很容易看出纬度越高，气温越低，当冷空气入侵时，气温下降。根据问题一与问题二的温度分布曲面图
18、可知，原来的坡度大，当冷空气入侵时，坡度减小，这说明冷空气对高温地区影响系数大于低温地区。五 模型评价 优点：为得到春分日地区最高温度分布情况，本文经过推导得出的最高温度求解模型，方法较为合理，计算简便，用matlab软件模拟出的温度分布图清晰明了。可以运用此模型求解任何时候任何地点的气温分布。缺点：问题考虑的因素不够全面，比如气温只考虑太阳辐射强度的影响，忽略了气候，海陆等其他因素，模型的检验还需要验证。六 参考文献1张鹤飞，与太阳辐射量有关角度的定义和计算公式2高等学校教学用书环境工程原理第一章第五节3李业发 谢红 张有为，全玻真空管太阳能辐射量与其放置角度的关系4http:/ 黄鼎业 华
19、毅杰，混凝土表面对流换热系数测试研究6张建荣 刘照求，混凝土对流换热系数的风洞实验研究七 附录用matlab软件编写程序，模拟出温度分布图：问题一：A=30.0,90.0,-0.20;50,120,3.16;30.0,92.5,11.00;50,117.5,2.56;32.5,90.0,-2.80; 37.5,90,6.64;45,90,7.6;30.0,95.0,18.40;30.0,97.5,14.90;47.5,90,-0.64; 50,115,1.18;42.5,90,9.8;42.5,92.5,3.38;32.5,95.0,3.80;47.5,90,-0.64; 40,90,11.7
20、;40,95,-1.5;30.0,100.0,10.60;42.5,97.5,1.1;50,112.5,0.1; 40,100,7.5;40,105,11.28;35.0,90.0,-0.80;47.5,92.5,2;47.5,97.5,1.4; 50,110,1.54;30.0,102.5,8.80;37.5,92.5,-0.2;37.5,97.5,2.02;45,95,4.3; 47.5,102.5,3.2;35.0,95.0,-2.60;32.5,120.0,25.82;32.5,117.5,25.43; 50,105,1.78;50,107.5,2.89;45,100,4.6;35,1
21、00,4.6;35,105,11.8;35,110,18.64; 42.5,102.5,8.18;47.5,107.5,5;47.5,112.5,3.86;30.0,105.0,15.88;37.5,102.5,8.98; 40,110,12.42;40,115,15.9;40,120,19.47;32.5,102.5,7.82;37.5,107.5,13.96; 50,102.5,0.22;50,105,1.78;50,107.5,2.89;35,110,18.64;35,115,22.42; 30.0,107.5,21.10; 45,105,6.22;45,110,7.18;42.5,10
22、7.5,9.5;32.5,117.5,25.43; 32.5,120.0,25.82;50,92.5,-6.32;50,95,-3.68;47.5,120,6.38;50,90,-9.08; 30.0,110.0,24.52;30.0,112.5,27.82;42.5,112.5,10.4;42.5,117.5,13.52; 37.5,107.5,13.96;37.5,112.5,17.14;30.0,115.0,26.20;42.5,120,15.38; 32.5,120.0,25.82;30.0,120.0,26.26;37.5,120,24.4;x=A(:,1);y=A(:,2);z=A
23、(:,3);scatter(x,y,5,z)%散点图figureX,Y,Z=griddata(x,y,z,linspace(min(x),max(x),linspace(min(y),max(y),v4);pcolor(X,Y,Z);shading interp %伪彩色图figure,contourf(X,Y,Z) %等高线图figure,surf(X,Y,Z) %三维曲面问题二：求解方程的程序：x=*3.14/180%把相应纬度换算成弧度p=0.035,28.61,282.6-476.2857*cos(x)roots(p)温度分布图程序：A=30.0,90.0,-24.3;50,120,-
24、3.04;30.0,92.5,-12.78;50,117.5,-3.64;32.5,90.0,-25.3; 37.5,90,-12.16;45,90,-4.1;30.0,95.0,-5.7;30.0,97.5,-3.72;47.5,90,-9.64; 50,115,-5.02;42.5,90,-4.5;42.5,92.5,-10.92;47.5,90,-9.64; 40,90,-4.9;40,95,-18.1;30.0,100.0,-13.5;42.5,97.5,-13.2;50,112.5,-6.1; 40,100,-9.1;40,105,-5.32;35.0,90.0,-21.5;47.5
25、,92.5,-7;47.5,97.5,-7.6; 50,110,-4.66;30.0,102.5,-15.3;37.5,92.5,-19;37.5,97.5,-16.78;45,95,-7.4; 47.5,102.5,-5.8;35.0,95.0,-23.3;32.5,120.0,3.32;32.5,117.5,2.93; 50,105,-4.42;50,107.5,-3.31;45,100,-7.1;35,100,-16.1;35,105,-8.9;35,110,-2.06; 42.5,102.5,-6.12;47.5,107.5,-4;47.5,112.5,-5.14;30.0,105.0
26、,-8.22;37.5,102.5,-9.82; 40,110,-4.18;40,115,-0.7;40,120,2.87;32.5,102.5,-14.68;37.5,107.5,-4.84; 50,102.5,-5.98;50,105,-4.42;50,107.5,-3.31;35,110,-2.06;35,115,1.72; 30.0,107.5,-3;45,105,-5.48;45,110,-4.52;42.5,107.5,-4.8;32.5,117.5,2.93; 32.5,120.0,3.32;50,92.5,-12.52;50,95,-9.88;47.5,120,-2.62;50
27、,90,-15.82; 30.0,110.0,0.42;30.0,112.5,3.72;42.5,112.5,-3.9;42.5,117.5,-0.78; 37.5,107.5,-4.84;37.5,112.5,-1.66;30.0,115.0,2.1;42.5,120,1.08; 32.5,120.0,3.32;30.0,120.0,2.16;37.5,120,3.2;x=A(:,1);y=A(:,2);z=A(:,3);scatter(x,y,5,z)figureX,Y,Z=griddata(x,y,z,linspace(min(x),max(x),linspace(min(y),max(y),v4);pcolor(X,Y,Z);shading interpfigure,contourf(X,Y,Z) 参赛队员及联系电话：何康康 13275424539郭皓 18766215275马西涛 13153242178