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

核科学

高斯气体分布
/ d9 Z% q( ^$ x' h( |1 I数值模拟
4 t6 [* p! Z1 }& q; |4 L& ]作出动态分布图

13ansys

看不懂什么意思  啊

923050019

function y=shuzhimoni
, r9 e8 F5 k! |& T0 a%syms x y z t;
  I3 x8 b3 ^* E& `- E5 J' l( i%Q=10;%初始浓度定为10mg/m^3;
%1.32*10^-5=1.32*10^-5;%扩散系数
, r7 m. Y* u, m+ T5 e. K8 }; ]3 `%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
%统一并申明一下单位
%C (mg/m^3)
. S" Q& y) \( T0 g! R$ A%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
%t (s)
%x,y,z (m)
3 s; l- ]4 v7 |8 S%申明一下扩散的条件如下：
) v4 ~& `' P$ F/ }% j%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
" x. }: B# h: o; Z% o' q% b0 c8 V. W%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
5 I* `: j$ h6 y! g/ d9 e/ v/ J, ]3 @randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
! {8 d; E; V5 ~+ d1 ifor i=1:4
    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
9 ]" Y* d! ~1 }9 L* f0 J( c    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));
    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
    zslice=0.25;
" r3 N  a( Y: \; L3 K    subplot(2,2,i);
  D9 t5 \2 K* M    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
    ylabel('y坐标');
    zlabel('z坐标');
6 Y1 j7 N/ U, Q( u: S    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
6 O. M, o$ {' b  r4 c, U    colorbar
end