数学建模社区-数学中国
标题:
谁能用MATLAB模拟高斯气体分布
[打印本页]
作者:
核科学
时间:
2010-5-5 18:11
标题:
谁能用MATLAB模拟高斯气体分布
高斯气体分布
3 P& n9 N2 {0 E8 E! z3 E. Z
数值模拟
% U& V5 b* q' }
作出动态分布图
作者:
923050019
时间:
2010-8-6 22:56
function y=shuzhimoni
9 x; A8 ]; S1 d/ p4 L5 q5 ]! ?
%syms x y z t;
! Y5 Z S, {8 r% O4 {
%Q=10;%初始浓度定为10mg/m^3;
7 K: H6 q, I0 B* m4 O- m
%1.32*10^-5=1.32*10^-5;%扩散系数
' h5 _9 X4 y9 I# T& U( ]9 g
%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
}1 V! N, ]' H2 C: b
%统一并申明一下单位
/ K* s( T$ ] I, ]$ x/ I
%C (mg/m^3)
8 n) |" j# Q( Q* G* b
%Q (mg/m^3)
b+ o7 H- C, c+ z r6 l
%1.32*10^-5 (m^2/s)
! m5 _' |7 g) P
%t (s)
# A/ H; Z! ?( P2 A
%x,y,z (m)
9 n% j1 k6 h) {/ c0 k _
%申明一下扩散的条件如下：
7 `. [2 Y7 t1 d6 b- j8 H
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
$ e3 g% W T. B6 U/ o) C/ R5 Q( m
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
( q' K2 l& S# p3 K) b; x
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
) J9 K0 {: i# F/ U$ c, O% ^9 V
for i=1:4
* ]4 Y7 c: D3 w" b' h
[x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
9 o+ s7 B6 k# B7 B( c, S- @$ j+ [
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));
8 h$ R" L1 E, ]) H _ P
xslice=[-0.25,1.75];
2 Q) W3 V/ ?7 g
yslice=[-0.25,1.25];
/ }& y! v$ r! m0 |& r
zslice=0.25;
* {8 |, ^2 T4 |# A2 i
subplot(2,2,i);
; p# I! N$ a/ p) H- Z' x
slice(x,y,z,C,xslice,yslice,zslice);
6 q. x& O- `: P+ [- s2 U) w# G
xlabel('x坐标');
1 s- Y- ~7 ?" L
ylabel('y坐标');
# W9 G& S5 ~' Y1 }0 r. q
zlabel('z坐标');
; ^/ ?6 l' w% _( _3 Z0 i
title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
7 M+ k* l- z2 B5 B: a, X" ~, B- x
colormap hsv
$ K: c& L$ I* J+ A$ I: {
colorbar
3 {' x# A% k6 D6 n
end
; u. a4 \3 M8 m# Y i; O X; U
作者:
13ansys
时间:
2014-8-11 14:24
看不懂什么意思 啊
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5