放射性气体扩散的预估模型数学建模优秀论文

1A题：放射性气体扩散的预估模型一：题目设有一座核电站遇自然灾害发生泄漏，浓度为p0的放射性气体以匀速排出，速度为mkg/s，在无风的情况下，匀速在大气中向四周扩散,速度为sm/s.1）请你建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。2）当风速为km/s时，给出核电站周边放射性物质浓度的变化情况。3）当风速为km/s时，分别给出上风和下风L公里处，放射性物质浓度的预测模型。4）将你建立的模型应用于福岛核电站的泄漏，计算出福岛核电站的泄漏对我国东海岸，及美国西海岸的影响。计算所用数据可以在网上搜索或根据具体情况自己模拟。二：摘要本论文关于核泄漏核电站周边不同距离地区、不同时段放射性物质浓度的预测模型，根据“泄漏放射性物质质量守恒定律”和“热传导定律”单位时间通过单位法向面积的流量与浓度梯度成正比。要探究风速对放射性物质浓度分布的影响:须考虑到P-G-T方法,Pasquill把风速和辐射状况作为划分稳定度等级的指标。利用常规气象资料把大气的扩散能力划分为六个稳定度等级，从A到F（极不稳定—稳定）。还考虑到城市:污染源多种多样,下垫面粗糙热岛效应。使得微气象特征和大气扩散规律与平原有着显著不同，城市中的高层建筑物、体形大的建筑物和构筑物，都能造成气流在小范围内产生涡流，阻碍气流运动，减小平均风速，降低了近地层风速梯度，并使风向摆动很大，近地层风场变得很不规则。关于问题四，在结合模型一、模型二的条件下，在参阅整理大量的气象、地理新闻质料，日本核事故期间核泄露产生放射性物质在低层先向南再向东扩散，据中央气象台预报，核泄漏期间日本近地面以西北风为主，核辐射物质的辐射量非常微量远低于当年切尔诺贝利灾害带给亚洲的辐射量不会影响到公众健康，不必担心福岛核电站事故产生的辐射物造成的危害。关键字:放射性气体扩散泄露放射物质质量守恒热传导定律P-G-T模型热岛效应2三：符号说明与名词解释t—气体扩散时间，气体由泄露源泄漏时刻t=0x,y,z—以泄漏源为坐标原点，空间任意一点的坐标C—空间中任一点的气体浓度k—气体扩散系数Q—气体由扩散源扩散时施放的气体总量—平均风速y—用浓度标准偏差表示y轴上的扩散参数z—用浓度标准偏差表示的z轴上的扩散参数H—气体扩散的有效高度x—下风方向到泄漏点源的距离y—侧风方向离泄漏源点的距离z—垂直向上方向离泄漏源点的距离k-风速s-放射性气体排除后向四周扩散的速度四：问题重述设有一座核电站遇自然灾害发生泄漏，浓度为p0的放射性气体以匀速排出，速度为mkg/s，在无风的情况下，匀速在大气中向四周扩散,速度为sm/s.5）请你建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。6）当风速为km/s时，给出核电站周边放射性物质浓度的变化情况。7）当风速为km/s时，分别给出上风和下风L公里处，放射性物质浓度的预测模型。8）将你建立的模型应用于福岛核电站的泄漏，计算出福岛核电站的泄漏对我国东海岸，及美国西海岸的影响。3五:问题分析(1):对于问题一的分析先研究在理想情况下，发射性气体在大气中的扩散问题。描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。放射性物质向四周扩散，形成圆形区域。无穷空间由瞬时点源导致的扩散过程，用二阶偏微分方程描述放射性物质浓度的变化。模型假设:放射性物质在无穷空间扩散，不受大地和风的影响；扩散服从热传导定律。模型建设：放射性物质浓度C（x,y,z,t）的变化规律热传导定律：单位时间通过单位法向面积的流量与浓度梯度成正比Ω内放射性物质的该变量曲面积分的奥氏公式：sVqnddivqdVqkgradC初始条件C(x,y,z,0)=Q(x,y,z)Q:代表放射元素释放的总量代表单位强度点源函数解得：对任意t,C的等值面是球面x2+y2+z2=R2;RCgradCkq流量通过],[ttttttsdtdnqQ12[(,,,)(,,,)]VQCxyztCxyzttdV12QQ222222[()]CCCCkdivgradCktxyz222222,,,,0CCCCkxyzttxyz222432(,,,)(4)xyzktQCxyztekt4仅当t,对任意点(x,y,z),C0nvΩ浓度扩散示意图上面的模型是最理想化浓度预测模型，因为它没有考虑风向地形微气象特征城市的扩散以及平原的扩散。但是这个模型远远不能解决实际问题，在实际中影响放射性物质扩散的因数很多。（2）对问题二的研究:当风速为km/s时，给出核电站周边放射性物质浓度的变化情况，这个就要解决实际问题考虑到风速的影响本文还考虑到城市以及山谷特殊气象条件下的扩散模式。a：城市扩散模型城市扩散模型:下垫面粗糙,产生热岛效应。风向与线源垂直时连续排放的无限长线源在横风向产生的浓度处处相等。则：2222,,0,expexp22zyyzqHyxyHdyu对上式积分得到:222,0,0,exp22zzqHxHub：山区的扩散模型：山风和谷风的方向是相反的，但比较稳定。在山风与谷风的转换期，风方向是不稳定的，山风和谷风均有机会出现，时而山风，时而谷风。这时若有大量污染物排入山谷中，由于风向的摆动，污染物不易扩散，在山谷中停留时间很长，有可能造成严重的大气污染。山区由于受到复杂地形的热力和动力因子影响，流场均匀和定常的假定以及正态分布假设均难以成立,因此应取向不稳定方向提及后的扩散参在垂直方向仍为正态分布。q1QS5由：得：对地面浓度：z=0；得高斯模式：污染物扩散所波及的垂直范围都处于同一温度层结中，即具有同一稳定度。实际中难以实现。c:封闭式扩散：上部有逆温，使污染物的垂直扩散受到限制，地面的反射作用，利用“像源法”进行处理只能在地面和逆温层地之间进行扩散。相当于两镜面之间无穷次全反射。实源在两个镜子里分别形成n个像计算公式:222/220/2,exp2exp2zWWzzxzAxzQuAxdydz222,exp22zzQzxzuW2,2zQxzuW6简化的计算公式，分三种情况处理：当：Dxx时，放射物质的扩散尚未受到上部逆温层的影响有：2.15zDH根据值查P—G曲线可得Dx值地面轴线的浓度：当2Dxx时，放射性物质经过两界面多次反射，达到某一距离X后，在Z方向的浓度分布将趋于均匀，但Y方向浓度分布仍为正态分布：方向的分布函数：22,()exp2yycxyAx则：220expDyyQuAxdydz即：22,exp2yyQycxyuD设泄漏源有效高度为H，取其在地面投影为坐标原点，x轴指向风向。考虑地面反射作用，可得到烟羽模型的浓度分布算式为：222222()()222(,,,)()2yzzyzHzHyzQCxyzteee模型二参数的确定有风时的扩散参数，(,)yz的确定采用Briggs给出一套扩散参数，如表1和表2所示。表1Briggs扩散参数（开阔平原田野）大气稳定度yz22(,0,0,)exp2zyzQHcxHu0111DdZDD7A21)0001.01(22.0xxx2.0B21)0001.01(16.0xxx12.0C21)0001.01(11.0xx21)0002.01(08.0xxD21)0001.01(08.0xx21)0015.01(06.0xxE21)0001.01(06.0xx)0003.01(03.0xxF21)0001.01(04.0xx)0003.01(016.0xx表2Briggs扩散参数（城市）大气稳定度xyA-B21)0004.01(32.0xx21)0001.01(24.0xxC21)0004.01(22.0xxx20.0D21)0004.01(16.0xx21)0003.01(14.0xxE-F21)004.01(11.0xx21)0015.01(08.0xx根据国家标准(GB/T13201--1991))制定地方大气污染物排放标准的技术方法的规定，划分大气稳定度的级别，共分为6级A~F，A为极不稳定；F为极稳定。大气稳定度的级别参考表地面风速1()ms白天太阳辐射阴天的白天或夜间有云的夜晚强中弱薄云遮天或低云≧0.5云量≦0.43AA-BBD--2-3BBCDEF3-5BB-CCDDE5-6CC-DDDDD6CDDDDD模型二求解：根据以上数据我们利用Ｃ语言程序编程计算程序如下：#includestdio.h#includemath.h#definee2.718281828459#definePI3.14159268voidqiuxy(charwending,doublex,doubledy,doubledz){switch(wending){case'A':dy=0.22*x/(sqrt(1+0.0001*x)),dz=0.2*x;break;case'B':dy=0.16*x/(sqrt(1+0.0001*x)),dz=0.12*x;break;case'C':dy=0.11*x/(sqrt(1+0.0001*x)),dz=0.08*x/(sqrt(1+0.0002*x));break;case'D':dy=0.08*x/(sqrt(1+0.0001*x)),dz=0.06*x/(sqrt(1+0.0015*x));break;case'E':dy=0.06*x/(sqrt(1+0.0001*x)),dz=0.03*x/(1+0.0003*x);break;case'F':dy=0.04*x/(sqrt(1+0.0001*x)),dz=0.016*x/(1+0.0003*x);break;}}voidmain(){doublec1,c2,x,z,t,p0,H,dx,dy,dz,a,s,k,I,Ts;chardaqi;printf(请按以下格式输入：\n);printf(下风口L处坐标（米）：\n);scanf(%f,x);printf(地点的纵坐标（米）：\n);printf(0\n);printf(地点的铅直坐标（米）：\n);scanf(%f,z);printf(扩散时间（秒）：\n);scanf(%f,t);printf(地面反射系数：\n);scanf(%f,a);printf(放射性气体扩散速度（m/s）：\n);scanf(%f,s);printf(风速（m/s）：\n);scanf(%f,k);printf(逆温的影响（Q\C）：\n);scanf(%f,I);printf(初始浓度（千克/立方米）：\n);scanf(%f,p0);printf(有效高度（米）：\n);scanf(%f,H);printf(大气稳定度等级：\n);scanf(%f,daqi);qiuxy(daqi,x,dy,dz);dx=dz;9c1=p0/(2*PI*(s-k)*dy*dz)*pow(e,-(0.693*x)/(3600*(Ts*24)*(k+s))-(8*e-5*pow(I,0.6)*x)/(k+s))*(pow(e,-((z-H)+(0.01*x/(s-k)))*((z-H)+(0.01*x/(s-k)))/2*dz*dz)+a*pow(e,-(z+H+(0.01*x/(s-k)))*(z+H+(0.01*x/(s-k)))/(2*dz*dz)));printf(下风口L处放射性物质浓度=%f\n,c1);c2=p0/(2*PI*(s-k)*dy