| 查看: 311 | 回复: 0 | |||
[交流]
Boltzmann方法求解决、、、
|
|
用Boltzmann方法加入我的多孔介质模型后，生成的速度矢量图总是和我的多孔介质不匹配。 请问有没有大神知道什么原因。。。 global Cd global Vr kn=50; [A,Cd]=model(kn); omega=1.0; tau=1.0/omega; m=50; n=50; deltaU=1e-7; %----------------------------------------- t=0; % 初始化密度 rho=ones(m,n); % 初始化速度(u,v) u=zeros(m,n); v=zeros(m,n); % 由已知条件给f赋初值 f=repmat(1.0/9.0,[m,n,9]); feq=f; % 设置边界 bound=A; % bound=zeros(m,n); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % bound=rand(m,n)>0.8; %设置孔隙率 Vr=0.6; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % A(:,[1 n])=1; balk=find(A); imshow(1-A);figure; |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有258人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有114人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复