完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
问题描述:
为了实现《液氯运输泄漏事故扩散风险分析》
液氯运输泄漏事故扩散风险分析.pdf
(1.05 MB, 下载次数: 160
)
论文中的一个改进的高斯烟团模型,移动点源的高斯烟团模型。 主要公式: 主要公式中涉及到的一个式子: 绘图最后结果: 而我的结果却是不管车速v设的有多大,浓度值最大点的横坐标总是21附近,只受风速的影响不受车速的影响,为了便于运算我的积分时间仅为10s:。 我的程序如下:%% clc clear %% %根据论文中设置参数 t=10; % 10秒 H=1; % 泄露点距离地面的高度为1m u=2; % 平均风速u=2m/s v=12.5; % 车速v=45km/h=12.5m/s z=1.5; %% for x=1:130 %deta_x、deta_y、deta_z由大气稳定度决定 E级 deta_x=(0.06* x *((1+0.0001*x).^(-0.5))); deta_y=deta_x; deta_z=(0.03 * x * ((1+0.0003*x).^(-0.5))); for y=-100:100 s(x,y+101)=integral(@(tba) (((6500000.*(1-0.007.*(tba.^0.5)))./(((2.*pi).^1.5).*deta_x.*deta_y.*deta_z)).* (exp(-((x-u*t+u*tba-v*tba)).^2) ./(2.*(deta_x.^2))).*(exp(-((y).^2)./(2.*(deta_y.^2)))).*(exp(-((z-H).^2)./(2.*(deta_z.^2)))+exp(-((z+H).^2)./(2.*(deta_z.^2))))),0,10); end end x=1:130; y=-100:100; s=s';%转置 s=flipud(s);%上下翻转 save test_movecontinue1_5.mat load test_movecontinue1_5.mat %移动点源持续泄露 figure surf(x,y,s);%画三维图 [cs, h]=contour(x,y,s,'LevelList',[6 90 300 2500]); %1.5m时画出值为等值线 grid minor,xlabel( ' x-axis ' ),ylabel( ' y-axis '),zlabel( ' z-axis ' ); title( '1.5m高度的扩散情况' ) ; |
|
相关推荐
|
|
你正在撰写答案
如果你是对答案或其他答案精选点评或询问,请使用“评论”功能。
1652 浏览 1 评论
200圆!求助大佬给一份VSG并网和离网模式的simulink仿真
1928 浏览 0 评论
MATLAB(3)--矩阵的引用(sub2ind、ind2sub、reshape函数使用)
2751 浏览 0 评论
3084 浏览 0 评论
4096 浏览 1 评论
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2024-11-26 09:14 , Processed in 0.635390 second(s), Total 74, Slave 55 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191 工商网监 湘ICP备2023018690号