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

Sharp corners

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 30, 2012, 07:46
Default Sharp corners
  #1
Member
 
Ford Prefect's Avatar
 
Join Date: Mar 2009
Posts: 72
Rep Power: 8
Ford Prefect is on a distinguished road
Hello,

What is the most common way of dealing with sharp corners in a FDM/FVM solver if we have a node-centered, colocated, approach (i.e. the corner coincides with a computational node).

I have implemented the boundary condition in such a way that node A "feel" a vertical wall at C. Node B "feel" a horizontal wall at C. (This is used with the normal momentum equation when calculating the pressure boundary condition)

Is this correct? Are there any drawbacks in this approach?

Cheers!
Attached Images
File Type: jpg corner.jpg (3.9 KB, 12 views)
__________________
"Trying is the first step to failure." - Homer Simpson
Ford Prefect is offline   Reply With Quote

Old   January 30, 2012, 10:40
Default
  #2
Senior Member
 
Andrew
Join Date: Mar 2009
Location: Washington, DC
Posts: 194
Rep Power: 8
mettler is on a distinguished road
Roache's book, "Computational Fluid Dynamics", has a good section on defining corners.
mettler is offline   Reply With Quote

Old   January 30, 2012, 18:38
Default
  #3
Member
 
Ford Prefect's Avatar
 
Join Date: Mar 2009
Posts: 72
Rep Power: 8
Ford Prefect is on a distinguished road
Quote:
Originally Posted by mettler View Post
Roache's book, "Computational Fluid Dynamics", has a good section on defining corners.
Thank you. Could you please give some indication what is written in that section with regard to my original question?

Cheers!
__________________
"Trying is the first step to failure." - Homer Simpson
Ford Prefect is offline   Reply With Quote

Old   January 31, 2012, 09:50
Default
  #4
Senior Member
 
Andrew
Join Date: Mar 2009
Location: Washington, DC
Posts: 194
Rep Power: 8
mettler is on a distinguished road
Roache's section is based on the streamfunction-vorticity set of equations, and one suggestion is to use the value of the wall being used in the computation. That means that, in your diagram, computing from left to right, the value at C would be based on node A first, since the vertical wall would be the first BC, and the when the horizontal wall is the BC it would be based on B.
So, when 0<x<c them point C is based on A, and when c<x<end C is based on B. Make sense?

I used that method doing a backwards facing step with good results.
mettler is offline   Reply With Quote

Old   January 31, 2012, 13:11
Default
  #5
Member
 
Ford Prefect's Avatar
 
Join Date: Mar 2009
Posts: 72
Rep Power: 8
Ford Prefect is on a distinguished road
Quote:
Originally Posted by mettler View Post
Roache's section is based on the streamfunction-vorticity set of equations, and one suggestion is to use the value of the wall being used in the computation. That means that, in your diagram, computing from left to right, the value at C would be based on node A first, since the vertical wall would be the first BC, and the when the horizontal wall is the BC it would be based on B.
So, when 0<x<c them point C is based on A, and when c<x<end C is based on B. Make sense?

I used that method doing a backwards facing step with good results.
Hey Mettler,

Ok, so this is basically the same approach as I use, except I do not change the value of C, instead I keep it at zero and add it as a source term in points A and B. Alternatively I modify my sparse coefficient matrix to incorporate the boundary conditions.

Now my problem is that I can solve cavity flow with my code to very nice accuracy. If I solve flow without internal obstacles I get good results as well. However, if I solve using an internal obstacle I get some (very minor) asymmetries in the solution even though I should not (I think). Flow over a cylinder in 2d gives ok values I guess but I want the asymmetries to disappear. I thought this was related to the treatment of corners. Well back to debugging I guess - Well unless anyone has any additional comments.

Cheers!
__________________
"Trying is the first step to failure." - Homer Simpson
Ford Prefect is offline   Reply With Quote

Old   January 31, 2012, 13:30
Default
  #6
Senior Member
 
Andrew
Join Date: Mar 2009
Location: Washington, DC
Posts: 194
Rep Power: 8
mettler is on a distinguished road
when I did the backwards facing step the value of C , vorticity in my example, definitely changed. It was updated each time, and also after each iteration the values were averaged to further take into account nearby node(s) - C in this example. Then, thru the iteration again, which would depend on the nearby/updated node values.

This works are low Re. I tried to incorporate this into a scheme with high Re that would let me see separation at a point, but I couldn't get that to work
mettler 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
[ICEM] Hole near sharp trailing edge of airplane/wing geometry jlichtwa ANSYS Meshing & Geometry 2 September 21, 2010 15:19
SnappyHexMesh and sharp corners madad2005 OpenFOAM 2 September 7, 2009 09:26
Issue with SimpleFoam with Sharp corners kdneroorkar OpenFOAM Running, Solving & CFD 13 September 15, 2008 18:12
Flow near sharp corners Harish Main CFD Forum 4 February 21, 2007 22:55
How to implement ghost cell wall for sharp angle? Qu Kun Main CFD Forum 3 August 2, 2004 07:29


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