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