谁能用MATLAB模拟高斯气体分布

核科学

高斯气体分布
数值模拟
作出动态分布图

923050019

function y=shuzhimoni
%syms x y z t;
( p0 @2 R; Z' z* B. Y%Q=10;%初始浓度定为10mg/m^3;
%1.32*10^-5=1.32*10^-5;%扩散系数
1 p3 R) g( x3 u) x' `%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
7 _& R0 I5 b7 v/ U$ K%统一并申明一下单位
%C (mg/m^3)
7 k1 c& V( Z5 j* z( o# R6 V% \%Q (mg/m^3)
' a$ }3 L+ X* O+ ?( `%1.32*10^-5 (m^2/s)
; ?7 A. k9 _. I8 C1 x5 s%t (s)
5 `: L: L3 o% U% S6 f3 H%x,y,z (m)
9 H& g7 y1 x; J& l' W; O( o2 e%申明一下扩散的条件如下：
1 Q( u1 X1 N- p, c# }%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
9 E, I$ ^' K% w* a%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
. S. d" Q/ H7 A8 i/ orandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
for i=1:4
6 ~2 l- N6 l# t5 I" o! m    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
7 s5 I8 ?+ |- U: p( n4 P' z7 _- a# V% i5 j    C=10/(4*pi.*1.32*10^-5.*randnum(i))^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*randnum(i));
* ^8 X4 j& m5 _; B+ i! y    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
    zslice=0.25;
    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
    ylabel('y坐标');
8 F0 X! |' n2 p# t  L% d    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
5 X  P( H, x$ T6 h9 p5 {- R! W    colorbar
end
1 c' e9 r0 o5 c/ A' y

13ansys

看不懂什么意思  啊