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

为什么80%的码农都做不了架构师？>>>   

高斯羽烟模型

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

转载于:https://my.oschina.net/VenusV/blog/3029463