高斯羽烟模型
已经是数学建模退役老油子了，偶尔发点模型计算方法出来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')
