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

核科学

高斯气体分布
数值模拟
2 ?0 `1 r! K, e  I  F8 u作出动态分布图

923050019

function y=shuzhimoni
5 N7 ]8 H  \  p- h' Y: }0 Z  ^%syms x y z t;
3 _) A0 d: f. X' L%Q=10;%初始浓度定为10mg/m^3;
%1.32*10^-5=1.32*10^-5;%扩散系数
6 |, t% m% z# u  W( r! c, r% q%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
& w( P- }4 g* p( A%统一并申明一下单位
1 b" E  L' l! p; U%C (mg/m^3)
%Q (mg/m^3)
4 M, e7 j  Q; S: u5 S! `4 s%1.32*10^-5 (m^2/s)
, L- {* Z1 P% Z! [' M8 P6 h%t (s)
%x,y,z (m)
( @' U& B( c3 [7 g/ _! o# V6 Z%申明一下扩散的条件如下：
; \% I3 \2 g6 H# ~; c) j8 e%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
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);
    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));
$ O1 m7 Z' k1 x3 ^- W    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
! o  r) o1 d) T% j# u( T" Q    zslice=0.25;
    subplot(2,2,i);
, u$ v6 y  A) j7 s# @    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
" Y5 I/ r; r. T7 r4 V( a    ylabel('y坐标');
    zlabel('z坐标');
8 p4 A9 ?7 V$ S, ]    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
0 }4 {  _$ |6 r" n7 \$ K! {    colormap hsv
    colorbar
end

13ansys

看不懂什么意思  啊