基于高斯模型的放射性物質擴散模型-資料下載頁
2025-06-20 12:45本頁面
【正文】 濃度的預測模型。此問題可以利用模型二進行求解，將模型二中的替換成,其中為風速的向量表示，為泄漏點放射性氣體自身擴散速度的向量表示。根據(jù)假設，風速沿軸正方向，恒為正。對上風向公里處放射性氣體濃度計算時，方向沿軸負方向，為負；對下風向公里處放射性氣體濃度計算時，方向沿軸正方向，為正。因此此問題轉化為求和兩處的放射性氣體濃度，即風速為時，上風和下風公里處放射性物質濃度的預測模型為： （30）我們根據(jù)模型分別模擬下風向和上風向的放射性物質擴散仿真圖，如圖五、圖六：（程序分別見附錄二、三）圖五：核泄漏下風向擴散圖圖六：核泄漏上風向擴散圖七 參考文獻[1] 卓金武，魏永生，秦建，李必文．Matlab在數(shù)學建模中的應用 [M] ．北京：北京航空航天大學出版社，2011．[2] 程勇，于林，姚安林．采用高斯模型分析輸氣管道泄漏后氣體的擴散 [J] ．內蒙古石油化工，2010 ，14：4951．[3] 張斌才，趙軍．大氣污染擴散的高斯煙羽模型及其GIS集成研究[J] ．環(huán)境監(jiān)測管理與技術，2008，20（5）：1719．八 附錄附錄一：%繪制當沒有風影響時放射性物質的擴散圖形function []=expand(k,Q,C,t,r,ws,wa,wd)R=sqrt((4)*k*t*log((C/((1r)*Q)).*(4*pi*k*t)^()))for i=0:24:t/3600if (wd*(i/24))=R R=Rwd。 R0=R。t=0::*pi。p=0::2*pi。wsx=ws*cos(wa)*i。wsy=ws*sin(wa)*i。[theta,phi]=meshgrid(t,p)。x=R*sin(theta).*cos(phi)+wsx。y=R*sin(theta).*sin(phi)+wsy。z=R*cos(theta)。hold onsurfc(x,y,z) endendR0=R0%繪制其它相關信息view(20,70) %觀察角度title(39。核泄漏無風影響擴散39。)。%圖形標題xlabel(39。x39。)。 %x軸標記ylabel(39。y39。)。 %y軸標記zlabel(39。z39。)。 %z軸標記axis([10,15,10,20,0,12]) %xyz軸顯示范圍余量x0=ws*cos(wa)*1000000/3600y0=ws*sin(wa)*1000000/3600R1=sqrt((x0)^2+(y0)^2)End附錄二：Q=2。 %輸入源強k=2。 %輸入風速s=。u=k+s。d=1。 %步長Z=。 x=10:d:300。 %下風向距離y=70:d:70。 %橫風向距離[x,y]=meshgrid(x,y)。by0=*x.*(1+*x).^(1/2)。bz0=*x.*(1+*x).^(1/2)。tempy1=y.*y./by0./by0./2。tempy2=.^(tempy1)。c=Q/pi/u*((by0.*bz0).^(1)).*tempy2Cs=2000。 %輸入求解的條數(shù)contour(x,y,c,Cs)。shading interp。colorbar。grid。title(39。核泄漏模擬圖一39。)附錄三：Q=2。 %輸入源強k=2。 %輸入風速s=。u=ks。d=1。 %步長Z=。 x=300:d:10。 %下風向距離y=70:d:70。 %橫風向距離[x,y]=meshgrid(x,y)。by0=*x.*(1+*x).^(1/2)。bz0=*x.*(1+*x).^(1/2)。tempy1=y.*y./by0./by0./2。tempy2=.^(tempy1)。c=Q/pi/u*((by0.*bz0).^(1)).*tempy2Cs=2000。 %輸入求解的條數(shù)contour(x,y,c,Cs)。shading interp。colorbar。grid。title(39。核泄漏模擬圖二39。)
教學課件相關推薦
鄂ICP備17016276號-1