方腔流--格子玻尔兹曼方法，基于MATLAB，带画图程序。

最新推荐文章于 2023-01-30 23:27:45 发布

QJ Meng 01

最新推荐文章于 2023-01-30 23:27:45 发布

阅读量6.2k

点赞数 6

分类专栏： MATLAB 数值模拟 LBM 文章标签： LBM MATLAB 格子玻尔兹曼方腔流

本文链接：https://blog.csdn.net/weixin_37783345/article/details/102727293

版权

还是老规矩先宣传一下QQ群：格子玻尔兹曼救星：293267908。

% To prove the credibility of fluid part of solver.
%方腔流体--速度计算。
%来自书--莫哈默德中文 A.A.Mohamad P141

%Qiujie Meng 2019.10.24.

%     ^y  
%   7 3 6   
%    \|/       
%   4-1-2   --->   x  
%    /|\   
%   8 5 9  
 
% Authors: Qiujie Meng, Copyright 2019-2021
clear;close all;clc;
 

% GENERAL FLOW CONSTANTS
Lx= 100;%Lx=  0~x_max
Ly= 100;
m=100; % m=  1~x_max
n=100;
nu= 0.01;   % kinematic viscosity
dt=1;
dx=1;
dy=1;
x=1:dx:100;
y=1:dy:100;
rhoo=5;
u0     = 0.1  ;      % maximum velocity of Poiseuille inflow
Re     = u0*Lx/nu ;      % Reynolds number
omega  = 1. / (3*nu+1./2.);      % relaxation parameter
maxT   = 40000 ;   % total number of iterations
tPlot  = 5;        % cycles
cs=1/(3)^0.5;

% D2Q9 LATTICE