【物理应用】基于Matlab实现LBM-D2Q9模型粗糙界面流动

2024-06-17 06:18:43

1 内容介绍

采用D2Q9模型，反弹边界，粗糙界面采用规则矩形

2 仿真代码

%%lbm方法模拟粗糙界面流动程序，

%%采用D2Q9模型，反弹边界，粗糙界面采用规则矩形

clear all;

clc;

%++++++++++++++++++++++++++++++++++++++++++++++++++++++++??????????????????????+++++++++++++++++++++++++++++++++++++++++

nx=40; ny=250;

bound=zeros(nx,ny);bound([1 nx],:)=1;A=1:38;

left_corner=zeros(1,5);right_corner=zeros(1,5);

for i=1:5

bound_point=40*i+(0:10);

bound(1:6,bound_point)=1;

bound(nx-(1:5),bound_point)=1;

left_corner(i)=40*i;

right_corner(i)=40*i+10;

A=[A ((40*i+1):(40*i+9)) (40*i+12):(40*(i+1)-2)];

end

A=[A A(end)+1:250];

B=bwboundaries(bound,8);

B1=B{1};B2=B{2};

B1(1:249,:)=[];B2(end-248:end,:)=[];

[~,tempi,~]=intersect(B1(:,2),A);b1=B1(tempi,:);

[~,tempi,~]=intersect(B2(:,2),A);b2=B2(tempi,:);

boundary=find(bound==1);

imshow(1-bound);

figure;

%+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++????????????+++++++++++++++++++++++++++++++++++++++++

Pin=1.001;Pout=1;

Rb=0.7;

tau= 1;

dx=1;dy=1;dt=1;c=dx/dt;t1=4/9; t2=1/9; t3=1/36; c_squ=c^2/3;

rouin=Pin/c_squ; rouout=Pout/c_squ;

DENSITY=rouout*ones(nx,ny);DENSITY(:,1)=rouin*ones(nx,1);

DENSITY(boundary)=rouout;

UX = zeros(nx,ny);UY= zeros(nx,ny);UX0=zeros(nx,ny);UY0=zeros([nx ny]);DENSITY0=zeros(nx,ny);

3 运行结果

4 参考文献

[1]刘立军, 胡文彬, 梅红岩. 基于变精度粗糙集模型的Matlab实现[J]. 计算机工程与应用, 2007.

[2]徐辉, 陶文铨. 熵格子Boltzmann方法模拟高Reynolds数流动[J]. 工程热物理学报, 2009(1):3.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。

【物理应用】基于Matlab实现LBM-D2Q9模型粗糙界面流动相关推荐

基于matlab的股票投资,基于matlab的股票估价模型系统.doc
基于matlab的股票估价模型系统.doc 东海科学技术学院毕业论文(设计)题目基于MATLAB的股票估计模型系统系机电工程学生姓名专业班级指导教师起止日期3浙江海洋学院本科生毕业论文1基于MATLA ...
matlab 摄动波浪理论,基于MATLAB的三维海浪模型数值仿真_齐宁.pdf
ISSN1009-3044 E-mail:eduf@ 第9卷第25期 (2013年09月) ComputerKnowledgeandTechnology电脑知识与技术 ComputerKnowledg ...
如何用Matlab建立信用风险模型,基于Matlab计算的KMV模型在商业银行信用风险管理中的实践应用.pdf...
科技论坛 2014.6 基于Matlab计算的KMV模型在商业银行信用风险管理中的应用李园 (天津大学管理与经济学部,天津,300072) 摘要:社会经济的不断发展,金融行业也处于不断发展之中, ...
基于Matlab人脸识别签到系统(GUI界面)
文件大小:5.3M 代码行数:298行(主程序) 开发环境:Matlab2016.2017.2018.2020.2021 点击下载:点击下载简要概述:基于Matlab人脸识别签到系统(GUI界面) ...
基于Matlab的虹膜识别系统(GUI界面)
点击查看:基于Matlab的虹膜识别系统(GUI界面) 文件大小:56M 操作系统:Windows10旗舰版开发工具:Matlab2016.2018.2019.2020 开发语言:.m 代码注释:
基于MATLAB的数字图像处理系统GUI界面设计
基于MATLAB的数字图像处理系统GUI界面设计图像读入从图形文件中读入图像 imread Syntax: A = imread(filename, fmt) filename:指定的灰度或彩色图 ...
D2Q9模型--基于MATLAB的LBM代码
哇咔咔又扒代码啦,感谢我的难兄老田同志贡献了自己的账号.这个代码上面没署名,不知道是哪个朋友写的,侵权立刻删哈. 还是老规矩先宣传一下QQ群群: 格子玻尔兹曼救星:293267908. % -- ...
如何使用matlab仿真物理实验,基于Matlab的大学物理实验模拟.pdf
维普资讯浙江科技学院学报,第 2O卷第 2期,2008年 6月 JournalofZhejiangUniversityofScienceandTechnology Vo1．20No．2,June20 ...
基于matlab的硅晶体模型,基于Matlab的图像处理技术识别硅太阳电池的缺陷
第 44 卷第 7 期 2010 年 7 月上海交通大学学报 JOURNAL OF SHANGHAI J IAOTON G UNIVERSITY Vol. 44 No. 7 Jul. ...

最新文章

热门文章