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

核科学

高斯气体分布
数值模拟
: p- S& ?3 T& l% u3 H作出动态分布图

923050019

function y=shuzhimoni
: G: ~  l) S9 K7 x' T6 K8 ^) ?%syms x y z t;
* a( M/ C7 d5 ^/ `5 a* t2 B, c. H%Q=10;%初始浓度定为10mg/m^3;
; l1 @# j+ f0 C( f3 J( p1 [%1.32*10^-5=1.32*10^-5;%扩散系数
9 Z+ C* Q0 u1 s( B( N. Y%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)
& k" H9 J# {( E% f$ E. I%Q (mg/m^3)
# I3 ^. G, p& n5 x. [' m%1.32*10^-5 (m^2/s)
%t (s)
+ B# l0 ^4 a3 [9 Q! Q+ ]1 q( _%x,y,z (m)
4 S3 U. C- F/ K2 T%申明一下扩散的条件如下：
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
" v7 Z+ `1 s5 L- ^; `, {%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
5 F6 d7 `! I. A+ G& ~) X: ofor i=1:4
& S1 N7 E% L- X    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
# q) ^! s* v3 t6 ?  D. f    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));
. ]+ a' w; s0 G2 t% E    xslice=[-0.25,1.75];
: m6 Q: ?/ D9 |2 s, F    yslice=[-0.25,1.25];
$ \' p( ]! |7 a* d    zslice=0.25;
7 {: f, t6 P% O8 h4 \9 ~    subplot(2,2,i);
* \( R  w6 R: ]4 n    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
/ |; _) s* i8 B' d( y0 R    ylabel('y坐标');
# p/ R# l% T. w* }' F. D6 \: c    zlabel('z坐标');
- r* o# b0 a8 A" L% U, z% L    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
; ?& [/ G! X2 [/ l0 j2 u    colormap hsv
    colorbar
; l3 d/ ]# p8 F+ c- `  gend
" ?- L  m  R; [) {3 b9 c# v

13ansys

看不懂什么意思  啊