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

核科学

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

923050019

function y=shuzhimoni
, L9 Z7 e; A9 Y! X2 Q0 e# E%syms x y z t;
%Q=10;%初始浓度定为10mg/m^3;
( j6 C# a( ^4 h2 B%1.32*10^-5=1.32*10^-5;%扩散系数
3 }8 C. B& t& }5 s%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
% L+ _3 x: u" c5 G7 L%统一并申明一下单位
7 f1 d/ ^6 Q$ N/ z% o%C (mg/m^3)
* Z$ W; _! ~3 g" m; Y7 L( T  h/ ~%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
%t (s)
5 `4 J4 _- e" a%x,y,z (m)
" v' |7 v" g+ Y) [: f7 _1 s%申明一下扩散的条件如下：
0 Z9 ^0 e; j( ^( y$ k/ ]%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
) z- u  y6 O4 L# K3 g7 c# }# q%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
' U# Y4 g! W, R0 nrandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
for i=1:4
    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
7 m, ~- W# ]0 W$ W2 d    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];
7 y  f- l, E+ X! p    yslice=[-0.25,1.25];
( ?8 p) M+ N5 M( @- q/ n7 o/ ]9 n    zslice=0.25;
    subplot(2,2,i);
% h+ r, x  L! L: Q    slice(x,y,z,C,xslice,yslice,zslice);
; h' r  ^/ [, F    xlabel('x坐标');
' w( F: v5 d2 G! m' k7 h' ~    ylabel('y坐标');
    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
7 H; }& P+ `2 S( g+ a! D, q    colormap hsv
4 N- e1 A/ K) W5 I* g    colorbar
" j8 o1 h( \1 s0 hend

13ansys

看不懂什么意思  啊