No self-introduction
Loading
Profile
Follow
codes (1)
Two dimensional simulation of cavity flow with LBM
no vote
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
quocyn 0 0
1
No more~
Follows0
Fans0