CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Fortran code-lattice boltzmann-boundary condition

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   October 26, 2012, 04:20
Default Fortran code-lattice boltzmann-boundary condition
  #1
New Member
 
Amir sadeghi pour marvi
Join Date: Dec 2011
Posts: 25
Rep Power: 5
Prince of Persia is on a distinguished road
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:
Prince of Persia is offline   Reply With Quote

Old   October 26, 2012, 10:23
Default
  #2
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,098
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
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.
flotus1 is offline   Reply With Quote

Old   October 26, 2012, 11:37
Default
  #3
New Member
 
Amir sadeghi pour marvi
Join Date: Dec 2011
Posts: 25
Rep Power: 5
Prince of Persia is on a distinguished road
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?
Prince of Persia is offline   Reply With Quote

Old   October 27, 2012, 11:42
Default
  #4
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,098
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
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.
flotus1 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 19:00.