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

CFD and corners

Register Blogs Community New Posts Updated Threads Search

Like Tree4Likes
  • 1 Post By cfdnewbie
  • 2 Post By sbaffini
  • 1 Post By sbaffini

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 11, 2013, 09:23
Arrow CFD and corners
  #1
Senior Member
 
James
Join Date: May 2013
Posts: 116
Rep Power: 13
Tensian is on a distinguished road
Hi everyone,

I have a very stupid question so I guess somebody could answer it easily. Why CFD gives bad results when itīs applied over corners? My teacher always told me that if the geometry has very marked corners then probably CFD will not provide accurate results, but why?
And why CFd provides bad results when cells have high skewness or aspect ratios?

Thanks in advance, I am a beginner...
Tensian is offline   Reply With Quote

Old   September 11, 2013, 13:40
Default
  #2
Senior Member
 
ghost82's Avatar
 
Rick
Join Date: Oct 2010
Posts: 1,016
Rep Power: 27
ghost82 will become famous soon enough
It's not very clear to me what your teacher says: I could imagine that main problems could be related to:
- mesh, as you wrote: skewness
- flow separation in turbulence regime: but there are some specific models to predict separation.

Daniele
ghost82 is offline   Reply With Quote

Old   September 11, 2013, 13:48
Arrow
  #3
Senior Member
 
James
Join Date: May 2013
Posts: 116
Rep Power: 13
Tensian is on a distinguished road
Hi Daniele,

Thanks for the reply, but you have not answered my question. Why a CFD code(for example FLUENT or OpenFOAM) provides bad results (or even do not achieve convergence) when the geometry have very marked corners? Is it due to numerical schemes (finite volume method) or is it a problem related to conservation equations? Or perharps is it due to the meshing produced to discretize this corner (skewed elements)?

Thanks
Tensian is offline   Reply With Quote

Old   September 11, 2013, 14:08
Default
  #4
Senior Member
 
cfdnewbie
Join Date: Mar 2010
Posts: 557
Rep Power: 20
cfdnewbie is on a distinguished road
The statement that codes give bad results at corners is an overly general one, and as ghost82 pointed out, there are a number of possiblities why such problems may occur. That would depend on the scheme, the way the BCs are enforced, the type of equation solved, the flow situation, the meshing etc...

It is a very broad statement and not true in general, so it is difficult to answer what your prof had in mind.
psakievich likes this.
cfdnewbie is offline   Reply With Quote

Old   September 16, 2013, 06:34
Default
  #5
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,191
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
More often than not, domain corners imply some form of singularity which, for example, has to be carefully treated in higher order methods in order to mantain the expected order of convergence (e.g., spectral methods).

Mesh skewness is a totally separate issue, which might be present or not near a domain corner. Very roughly speaking, in FV the issue for the skewness and mesh non orthogonality in general is related to the alignment of face normals and face centroids with respect to the cell centers.
FMDenaro and aero_head like this.
sbaffini is offline   Reply With Quote

Old   October 8, 2013, 00:29
Default
  #6
Senior Member
 
Ashwani
Join Date: Sep 2013
Location: Hyderabad
Posts: 154
Rep Power: 13
AshwaniAssam is on a distinguished road
One more reason can be fluent,etc, solver is based on cell centred calculation, which then is extrapolated on the corner. This would be depending on the Mesh size at the corner points.
AshwaniAssam is offline   Reply With Quote

Old   October 8, 2013, 06:40
Default
  #7
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,865
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
It is a problem due to singularity appearing in the representation of the continuous geometry over a discrete mesh representation.
FMDenaro is offline   Reply With Quote

Old   October 8, 2013, 06:45
Default
  #8
Senior Member
 
Ashwani
Join Date: Sep 2013
Location: Hyderabad
Posts: 154
Rep Power: 13
AshwaniAssam is on a distinguished road
Thank you. I have not come across singularity problem in CFD much. Is singularity related to zero volume cells which gets created while generating mesh.
AshwaniAssam is offline   Reply With Quote

Old   October 9, 2013, 03:03
Default
  #9
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,191
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
The singularity is related to the fact that you are fixing your solution over a non regular boundary. The consequence is that some derivatives of the solution might not exist (which might become an issue for high order methods).

The issue is described, for example, in Boyd: Chebyshev and Fourier Spectral Methods, Chapter 19. However, nearly all the books on spectral methods have at least a small description of this issue.
AshwaniAssam likes this.
sbaffini is offline   Reply With Quote

Old   October 9, 2013, 22:25
Default
  #10
Senior Member
 
Julien de Charentenay
Join Date: Jun 2009
Location: Australia
Posts: 231
Rep Power: 18
julien.decharentenay is on a distinguished road
Send a message via Skype™ to julien.decharentenay
Hi Tensian,

First, good on you for posting your very excellent question and trying to understand the ins and outs of what your professor has been teaching/telling you.

As pointed out in some of the responses, your question is stated in very general terms and would benefit from being more specific. As a counter-example, sharp corners in bluff body are well defined separation points which are easier to predict in CFD.

Below may be an explanation of what your teacher may be referring to:

Fluent and openFoam are general purpose CFD software using the finite volume method. As you may be aware, the finite volume method solves the variable at the cell center, and uses projection of the variable to the faces to calculate the flux used in the conservation equation.

Typically boundary layer meshes are applied near walls, resulting in cells with large aspect ratio (thin, but long and large). At sharp corners, the mesh would look as follows (depending on what software you may be using):



In such case, non-orthogonality issues associated with the face (shown in red) and interpolation of the cell centered variable to the face can be problematic.
__________________
---
Julien de Charentenay
julien.decharentenay is offline   Reply With Quote

Old   August 12, 2020, 13:31
Default No pictures shown
  #11
Member
 
Saeed Pashazanousi
Join Date: Mar 2020
Location: Iran
Posts: 66
Rep Power: 6
pashazanousi is on a distinguished road
hi...thanks for your good reply....But The pictures not available...
if it's possible for you upload it again....Thank you
__________________
Best regards

Saeed Pashazanousi
Urmia University
Email: st_s.pashazanousi@urmia.ac.ir
pashazanousi is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
gambit 64 bit?? artemiss1984 ANSYS Meshing & Geometry 23 November 21, 2016 02:07
[ANSYS Meshing] ANSYS 14.0 Workbench ICEM CFD Interactive meshing Bill Tuer ANSYS Meshing & Geometry 4 May 30, 2013 11:51
Corners Michael1408 Main CFD Forum 1 January 22, 2012 17:19
[ICEM] beginner questions for ICEM CFD johnp ANSYS Meshing & Geometry 14 August 10, 2010 18:45
ANSYS ICEM CFD smoothing mesh corner chamfered Hans CFX 4 January 26, 2006 17:00


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