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

核科学

高斯气体分布
: N+ x9 E0 _7 ^; |( w数值模拟
. @0 s7 |0 R5 i- A作出动态分布图

13ansys

看不懂什么意思  啊

923050019

function y=shuzhimoni
%syms x y z t;
  H. C3 q: L6 k. Q0 ]# d; S%Q=10;%初始浓度定为10mg/m^3;
4 s; F1 r) }1 N" d2 E% j) }%1.32*10^-5=1.32*10^-5;%扩散系数
( g/ I2 G: v) I# 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);
8 b" v5 m, F2 a2 S* s+ |2 {9 Y%统一并申明一下单位
%C (mg/m^3)
0 l; X2 ?  n5 H; x) J# h! T% w%Q (mg/m^3)
, B: z/ A) \9 Z7 {3 a- b%1.32*10^-5 (m^2/s)
%t (s)
9 u- ?# p9 F. [%x,y,z (m)
%申明一下扩散的条件如下：
! v; T  h7 e. ~% i+ F" w%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
- x5 T  F  ~# A8 ]$ Q4 X%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
6 R9 k, h' P9 vrandnum=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);
    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];
4 k, G6 l7 {; I    yslice=[-0.25,1.25];
    zslice=0.25;
  r" i/ x" Z0 r7 B; z( ~    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
$ A7 S. y* Q8 y, e! D/ h/ @8 @    xlabel('x坐标');
    ylabel('y坐标');
    zlabel('z坐标');
2 f* z' L8 e; [! d. w. U6 Q- k    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
    colorbar
end
4 I. }# t  J$ D2 j/ |