基于高斯模型的放射性物質(zhì)擴(kuò)散模型-資料下載頁(yè)
2025-06-20 12:45本頁(yè)面
【正文】 濃度的預(yù)測(cè)模型。此問(wèn)題可以利用模型二進(jìn)行求解，將模型二中的替換成,其中為風(fēng)速的向量表示，為泄漏點(diǎn)放射性氣體自身擴(kuò)散速度的向量表示。根據(jù)假設(shè)，風(fēng)速沿軸正方向，恒為正。對(duì)上風(fēng)向公里處放射性氣體濃度計(jì)算時(shí)，方向沿軸負(fù)方向，為負(fù)；對(duì)下風(fēng)向公里處放射性氣體濃度計(jì)算時(shí)，方向沿軸正方向，為正。因此此問(wèn)題轉(zhuǎn)化為求和兩處的放射性氣體濃度，即風(fēng)速為時(shí)，上風(fēng)和下風(fēng)公里處放射性物質(zhì)濃度的預(yù)測(cè)模型為： （30）我們根據(jù)模型分別模擬下風(fēng)向和上風(fēng)向的放射性物質(zhì)擴(kuò)散仿真圖，如圖五、圖六：（程序分別見(jiàn)附錄二、三）圖五：核泄漏下風(fēng)向擴(kuò)散圖圖六：核泄漏上風(fēng)向擴(kuò)散圖七 參考文獻(xiàn)[1] 卓金武，魏永生，秦建，李必文．Matlab在數(shù)學(xué)建模中的應(yīng)用 [M] ．北京：北京航空航天大學(xué)出版社，2011．[2] 程勇，于林，姚安林．采用高斯模型分析輸氣管道泄漏后氣體的擴(kuò)散 [J] ．內(nèi)蒙古石油化工，2010 ，14：4951．[3] 張斌才，趙軍．大氣污染擴(kuò)散的高斯煙羽模型及其GIS集成研究[J] ．環(huán)境監(jiān)測(cè)管理與技術(shù)，2008，20（5）：1719．八 附錄附錄一：%繪制當(dāng)沒(méi)有風(fēng)影響時(shí)放射性物質(zhì)的擴(kuò)散圖形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%繪制其它相關(guān)信息view(20,70) %觀察角度title(39。核泄漏無(wú)風(fēng)影響擴(kuò)散39。)。%圖形標(biāo)題xlabel(39。x39。)。 %x軸標(biāo)記ylabel(39。y39。)。 %y軸標(biāo)記zlabel(39。z39。)。 %z軸標(biāo)記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。 %輸入源強(qiáng)k=2。 %輸入風(fēng)速s=。u=k+s。d=1。 %步長(zhǎng)Z=。 x=10:d:300。 %下風(fēng)向距離y=70:d:70。 %橫風(fēng)向距離[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。 %輸入源強(qiáng)k=2。 %輸入風(fēng)速s=。u=ks。d=1。 %步長(zhǎng)Z=。 x=300:d:10。 %下風(fēng)向距離y=70:d:70。 %橫風(fēng)向距離[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。)
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1