24小时热门版块排行榜    

查看: 253  |  回复: 0

droiyton

新虫 (初入文坛)

[交流] 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;
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 droiyton 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见