数学建模社区-数学中国
标题:
谁能用MATLAB模拟高斯气体分布
[打印本页]
作者:
核科学
时间:
2010-5-5 18:11
标题:
谁能用MATLAB模拟高斯气体分布
高斯气体分布
9 A; m7 X& E. C4 H( r; F$ |
数值模拟
! H" X3 x$ e5 Y0 V) w
作出动态分布图
作者:
923050019
时间:
2010-8-6 22:56
function y=shuzhimoni
- k: z. }/ [+ h2 a
%syms x y z t;
( n# U+ f' k5 G5 g: p8 z; o2 y
%Q=10;%初始浓度定为10mg/m^3;
6 n( F$ P" @: O% c6 p, M
%1.32*10^-5=1.32*10^-5;%扩散系数
& {( [ ]9 o( A7 ]: S' f5 C& ~& d
%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
6 [1 V* H; G0 ~
%统一并申明一下单位
- }$ J" ]! v; F3 g
%C (mg/m^3)
& `9 |, n3 P# ~% R" d2 [
%Q (mg/m^3)
6 o& N" h+ z* d$ G7 m+ C+ G
%1.32*10^-5 (m^2/s)
* l, P, u% Z6 O2 ` ?
%t (s)
2 e8 N! b: w! Y+ A4 l" Y
%x,y,z (m)
1 l6 f m: S+ e% {0 }6 N
%申明一下扩散的条件如下：
$ R; D {- f+ o
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
( H8 i9 S" G" a- d# T
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
3 l- {+ }0 n& p- A6 S. k
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
4 {4 e0 [) R8 K
for i=1:4
6 l+ f/ T* q; ~7 o) E2 f7 V* ^
[x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
: U( B: c+ G6 I8 q
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));
' k' f2 }1 S; f7 }! V
xslice=[-0.25,1.75];
- H& K) _; ^, h& ^
yslice=[-0.25,1.25];
: m3 K$ y% o: Z5 l/ S. a* P0 V6 Z' J" x
zslice=0.25;
9 O, J3 b; x: r7 W& a, @' G% K
subplot(2,2,i);
% H* j* _7 n2 {0 v7 h( y
slice(x,y,z,C,xslice,yslice,zslice);
* D4 E P) i, K3 }5 O5 l( G% U
xlabel('x坐标');
8 _& y' @3 {4 j k
ylabel('y坐标');
' `6 Z2 @% }5 S2 c
zlabel('z坐标');
0 g0 \! V, I1 |4 _* D
title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
8 e+ S& o" a# y3 h
colormap hsv
# N8 A$ C& N& \; S
colorbar
0 S5 k8 z0 K) \: i) J
end
" Q8 ^; S8 {0 {! f- N( `# M9 z0 u0 G& V
作者:
13ansys
时间:
2014-8-11 14:24
看不懂什么意思 啊
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5