Two dimensional simulation of cavity flow with LBM
2016-08-23
0 0 0
no vote
Other
Earn points
Application background
This code is about simple LBM cavity driven flow. It is writen in F language.
! computer code for lid-driven cavity
parameter (n=1000,m=40)
real f(0:8,0:n,0:m)
real feq(0:8,0:n,0:m),rho(0:n,0:m)
real w(0:8), cx(0:8),cy(0:8)
real u(0:n,0:m), v(0:n,0:m)
integer i
open(2,file='uvfield')
open(3,file='uvely')
open(4,file='vvelx')
open(8,file='timeu')
!
uo=0.20
sumvelo=0.0
rhoo=5.00
dx=1.0
dy=dx
dt=1.0
alpha=0.02
Re=uo*m/alpha
print *, "Re=", Re
omega=1.0/(3.*alpha+0.5)
mstep=4000
w(0)=4./9.
!
do i=1,4
w(i)=1./9.
end do
!
do i=5,8
w(i)=1./36.
end do
!
cx(0)=0
cx(1)=1
cx(2)=0
cx(3)=-1
cx(4)=0
cx(5)=1
cx(6)=-1
cx(7)=-1
cx(8)=1
cy(0)=0
cy
fortran
lbm
模拟
二维
空腔
Related Source Codes
Deform secondary development
0
0
no vote
unifiber-VUMAT
0
0
no vote
Dir command simulating DOS
0
0
no vote
3D Hashin failure criterion subroutine, Umat, for
0
0
no vote
Matlab two-dimensional photonic crystal transfer m
0
0
no vote
No comment