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

Boundary condition of cell-centered FVM

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 11, 2010, 12:04
Default Boundary condition of cell-centered FVM
  #1
New Member
 
Edward
Join Date: Jul 2009
Posts: 12
Rep Power: 0
jackson is on a distinguished road
Is it possible to make the velocity be zero exactly at solid boundary(viscous no slip BC) by the cell-centerd finite volume method, i knew finite difference method and vertex based finite volume method can do that ,but not clear of whether cell centered FVM could or not
jackson is offline   Reply With Quote

Old   May 11, 2010, 16:10
Default
  #2
New Member
 
Jiannan
Join Date: Mar 2009
Posts: 10
Rep Power: 8
jntan is on a distinguished road
Quote:
Originally Posted by jackson View Post
Is it possible to make the velocity be zero exactly at solid boundary(viscous no slip BC) by the cell-centerd finite volume method, i knew finite difference method and vertex based finite volume method can do that ,but not clear of whether cell centered FVM could or not
hi jackson,my answer is a big YES coz I've tried that way and it work just fine.
Good luck!
jntan is offline   Reply With Quote

Old   May 11, 2010, 18:26
Default
  #3
New Member
 
Rich
Join Date: Aug 2009
Location: Montrose, Colorado, USA
Posts: 11
Rep Power: 7
Rich is on a distinguished road
Jackson:

I agree with Jiannan, it does work. I've also used it. For a viscous flow, you can simply set the velocity in the next-to-boundary cells to zero.

For an inviscid flow, where the velocity in those cells is finite and always parallel to the surface, you have to be more formal. I've found that you can derive governing equations that apply specifically to the boundary cells. In an inviscid flow, the "V dot dS" term in the integral form of the governing equations is always zero for the cell wall that forms the boundary. Rederive your equations with that term set to zero for that wall and apply them to the boundary cells, and those cells operate just as they should.


Rich
Rich is offline   Reply With Quote

Old   May 11, 2010, 21:02
Default
  #4
New Member
 
Edward
Join Date: Jul 2009
Posts: 12
Rep Power: 0
jackson is on a distinguished road
Set the velocity in the next-to-boundary cells to zero is not strictly right way? Is there a more accurate method? For my code, i give some special treatments for no slip B.C:
1.gradient computation in the cell next to wall boundary(Gauss Green Method,the face velocity on wall is set zero )
2.i think the flux on wall is the same as inviscid flux,so no more special treatment
3.viscous flux calculation contains face velocity and it is set zero exactly

Then the code ran, but i found the velocity on boundary face is not zero by Tecplot.....................

IS THERE SOME PROBLEMS OF MY TREATMENT?
THANK YOU FOR YOUR HELP

Quote:
Originally Posted by Rich View Post
Jackson:

I agree with Jiannan, it does work. I've also used it. For a viscous flow, you can simply set the velocity in the next-to-boundary cells to zero.

For an inviscid flow, where the velocity in those cells is finite and always parallel to the surface, you have to be more formal. I've found that you can derive governing equations that apply specifically to the boundary cells. In an inviscid flow, the "V dot dS" term in the integral form of the governing equations is always zero for the cell wall that forms the boundary. Rederive your equations with that term set to zero for that wall and apply them to the boundary cells, and those cells operate just as they should.


Rich
jackson is offline   Reply With Quote

Old   May 12, 2010, 00:09
Default
  #5
Super Moderator
 
praveen's Avatar
 
Praveen. C
Join Date: Mar 2009
Location: Bangalore
Posts: 248
Blog Entries: 6
Rep Power: 9
praveen is on a distinguished road
Quote:
Originally Posted by jackson View Post
Set the velocity in the next-to-boundary cells to zero is not strictly right way? Is there a more accurate method? For my code, i give some special treatments for no slip B.C:
1.gradient computation in the cell next to wall boundary(Gauss Green Method,the face velocity on wall is set zero )
2.i think the flux on wall is the same as inviscid flux,so no more special treatment
3.viscous flux calculation contains face velocity and it is set zero exactly

Then the code ran, but i found the velocity on boundary face is not zero by Tecplot.....................

IS THERE SOME PROBLEMS OF MY TREATMENT?
THANK YOU FOR YOUR HELP
You are doing it correctly.

For cell-centered scheme, boundary conditions are imposed through appropriate definition of boundary fluxes.

Dont worry if tecplot shows velocity on boundary not satisfying the boundary condition. Tecplot does not know your boundary condition, so it is doing some interpolation to get the boundary values.
praveen is offline   Reply With Quote

Old   May 13, 2010, 15:55
Default
  #6
New Member
 
Rich
Join Date: Aug 2009
Location: Montrose, Colorado, USA
Posts: 11
Rep Power: 7
Rich is on a distinguished road
Quote:
Originally Posted by jackson View Post
Set the velocity in the next-to-boundary cells to zero is not strictly right way? Is there a more accurate method? For my code, i give some special treatments for no slip B.C:
1.gradient computation in the cell next to wall boundary(Gauss Green Method,the face velocity on wall is set zero )
2.i think the flux on wall is the same as inviscid flux,so no more special treatment
3.viscous flux calculation contains face velocity and it is set zero exactly

Then the code ran, but i found the velocity on boundary face is not zero by Tecplot.....................

IS THERE SOME PROBLEMS OF MY TREATMENT?
THANK YOU FOR YOUR HELP
No, I think your treatment is correct. I was just pointing out what would be necessary IF you were working with an INviscid system. Sorry I confused you.


Rich
Rich is offline   Reply With Quote

Old   May 13, 2010, 16:19
Default
  #7
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 369
Rep Power: 10
arjun is on a distinguished road
Quote:
Originally Posted by jackson View Post
Is it possible to make the velocity be zero exactly at solid boundary(viscous no slip BC) by the cell-centerd finite volume method, i knew finite difference method and vertex based finite volume method can do that ,but not clear of whether cell centered FVM could or not

that probably is interpolation issue.
arjun is offline   Reply With Quote

Old   May 14, 2010, 08:23
Default
  #8
New Member
 
Edward
Join Date: Jul 2009
Posts: 12
Rep Power: 0
jackson is on a distinguished road
i output the cell-centered value for tecplot cell centered mode,tecplot auto interpolate them to vertexes, i changed it to cell centered display mode and found it is not closed to zero. so i think it is not interpolation's problem. can someone offer me some paper about
the implement of BC?
Quote:
Originally Posted by praveen View Post
You are doing it correctly.

For cell-centered scheme, boundary conditions are imposed through appropriate definition of boundary fluxes.

Dont worry if tecplot shows velocity on boundary not satisfying the boundary condition. Tecplot does not know your boundary condition, so it is doing some interpolation to get the boundary values.
jackson 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
SymmetryPlane Boundary Condition raytracer OpenFOAM Running, Solving & CFD 0 July 11, 2008 17:04
Airfoil boundary condition Frank Main CFD Forum 1 April 21, 2008 18:36
Boundary Condition When using FVM? leaf Main CFD Forum 2 December 5, 2007 04:53
Rotating interpolated boundary condition hani OpenFOAM Running, Solving & CFD 0 July 4, 2006 07:09
How to visualise cell centered FVM mesh? assabari Main CFD Forum 1 March 1, 2001 08:39


All times are GMT -4. The time now is 14:13.