Fortran code-lattice boltzmann-boundary condition

Prince of Persia · October 26, 2012, 05:20

Hi everyone
anybody can help me to write boundary condition fortran code.
I have boundary between cylinder and neighboring lattice point.

rw is the point at the cylinder surface , rp is the neighboring point inside the cylinder , rl is the neighboring site in the fluid area.
Yu-Mei-Lue-shyy model:

this method uses following equation:

flotus1 · October 26, 2012, 11:23

Isn't this just a bounce-back boundary condition with interpolation?
Where exactly do you see the problem for the code implementation?

Ah now I think I see: The main problem will be to provide the link information, since you need the populations from neighboring collision sites.
One way to do this is to store the link information (the number of the neighboring site for each link) in a separate field. This way, r_l' can easily be found.

Prince of Persia · October 26, 2012, 12:37

thanks for your reply.
it's easy to store the number of the neighboring sites.

my problem is the delta.

How can I determine it?

flotus1 · October 27, 2012, 12:42

This might be done with a little geometry.
Basically, you need to determine the intersection between a straight line (defined by the collision site and the link vector) and a circle (defined by the cylinder).
From the intersection point and the position of the collision site it should be straightforward to determine the cut length Delta.

October 26, 2012, 05:20	Fortran code-lattice boltzmann-boundary condition	#1
Prince of Persia New Member Amir sadeghi pour marvi Join Date: Dec 2011 Posts: 25 Rep Power: 14	Hi everyone anybody can help me to write boundary condition fortran code. I have boundary between cylinder and neighboring lattice point. rw is the point at the cylinder surface , rp is the neighboring point inside the cylinder , rl is the neighboring site in the fluid area. Yu-Mei-Lue-shyy model: this method uses following equation:

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
An error has occurred in cfx5solve:	volo87	CFX	5	June 14, 2013 18:44
How to set the Heat Flux boundary condition at Outlet	creddy_trddc	CFX	3	September 21, 2011 08:44
External Radiation Boundary Condition (Two sided wall), Grid Interface	CFD XUE	FLUENT	0	July 8, 2010 07:49
vorticity boundary condition	bearcharge	Main CFD Forum	0	May 14, 2010 12:32
Airfoil boundary condition	Frank	Main CFD Forum	1	April 21, 2008 19:36

October 26, 2012, 11:23		#2
flotus1 Super Moderator Alex Join Date: Jun 2012 Location: Germany Posts: 3,427 Rep Power: 49	Isn't this just a bounce-back boundary condition with interpolation? Where exactly do you see the problem for the code implementation? Ah now I think I see: The main problem will be to provide the link information, since you need the populations from neighboring collision sites. One way to do this is to store the link information (the number of the neighboring site for each link) in a separate field. This way, r_l' can easily be found.

October 26, 2012, 12:37		#3
Prince of Persia New Member Amir sadeghi pour marvi Join Date: Dec 2011 Posts: 25 Rep Power: 14	thanks for your reply. it's easy to store the number of the neighboring sites. my problem is the delta. How can I determine it?

October 27, 2012, 12:42		#4
flotus1 Super Moderator Alex Join Date: Jun 2012 Location: Germany Posts: 3,427 Rep Power: 49	This might be done with a little geometry. Basically, you need to determine the intersection between a straight line (defined by the collision site and the link vector) and a circle (defined by the cylinder). From the intersection point and the position of the collision site it should be straightforward to determine the cut length Delta.