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

核科学

高斯气体分布
" o+ Y+ X! F: a数值模拟
$ c/ A% q; b2 k6 K- e/ n# q作出动态分布图

13ansys

看不懂什么意思  啊

923050019

function y=shuzhimoni
%syms x y z t;
: q3 u( c5 k, x3 X9 U) D%Q=10;%初始浓度定为10mg/m^3;
%1.32*10^-5=1.32*10^-5;%扩散系数
%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
* I7 {& W! G& L" k%统一并申明一下单位
4 C9 c# ?' ?9 K) L0 W6 s1 X%C (mg/m^3)
%Q (mg/m^3)
3 ^& h' F/ M' R4 q%1.32*10^-5 (m^2/s)
( q, m/ I/ t' A%t (s)
%x,y,z (m)
! s4 y+ t/ n$ c%申明一下扩散的条件如下：
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
/ ]/ h" `! [# J+ V( e0 [4 @randnum=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);
& H" u- C. n' o3 h    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];
/ C" h- S2 H, {3 E  K# P    yslice=[-0.25,1.25];
/ P% \) [- k5 O4 Z' M    zslice=0.25;
( ?7 l1 n- Z6 T) G4 \" L    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
3 f* {- _0 |" a    xlabel('x坐标');
    ylabel('y坐标');
    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
* `- B0 k& [! R/ e* B    colorbar
end
. {4 g& {/ D2 n: w: W4 s4 X