Math-Model(三)高斯羽烟模型计算气体扩散浓度

原创
2019/03/28 23:01
阅读数 4.2K

高斯羽烟模型

已经是数学建模退役老油子了,偶尔发点模型计算方法出来vanvan
对高斯羽烟模型方程:

式中c为污染物浓度(单位:kg/m3)
Q为源强(单位:kg/s)
u为泄漏高度的平均风速(单位:m/s)
y、z分别用浓度标准偏差表示的y轴及z轴上的扩散参数
H为泄漏有效高度(单位:m)

扩散系数系数

稳定度 a b c d
A 0.527 0.865 0.28 0.90
B 0.371 0.866 0.23 0.85
C 0.209 0.897 0.22 0.80
D 0.123 0.905 0.20 0.76
E 0.098 0.902 0.15 0.73
F 0.065 0.902 0.12 0.67

Matlab 仿真代码:

clc;
clear all;

[x,y]=meshgrid(0:20:5000,-2000:20:2000);

Q=1000;%源强
%sigma_y=ax^b 水平扩散系数
a=0.123;
b=0.905;
sigy=a*x.^b;
%sigma_z=cx^d 垂直扩散系数
c=0.2;
d=0.76;
sigz=c*x.^d;
u=1;%风速
H=200;%泄漏有效高度
z=0;%高度

Qpi=2*Q./(2*pi*u.*sigy.*sigz+eps);
ex1=exp(-0.5*(y./(sigy+eps)).^2);
ex2=exp(-0.5*((z-H)./(sigz+eps)).^2)+exp(-0.5*((z+H)
C=Qpi.*ex1.*ex2;

contour(x,y,C,'ShowText','on')


展开阅读全文
打赏
0
0 收藏
分享
加载中
更多评论
打赏
0 评论
0 收藏
0
分享
返回顶部
顶部