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

High order convective schemes

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

Like Tree1Likes
  • 1 Post By pbohorquez

Reply
 
LinkBack Thread Tools Display Modes
Old   November 16, 2011, 09:01
Default High order convective schemes
  #1
Senior Member
 
Amir's Avatar
 
Amir
Join Date: May 2009
Location: Shiraz, Iran
Posts: 739
Blog Entries: 1
Rep Power: 14
Amir is on a distinguished road
Hi,

Can anyone suggest a reliable high order scheme for convective term of a tensor equation in which convective term is much more dominated? I examine upwind, limitedLinear, Gamma, SFCD, QUICK and Minmod but the results change with increasing grid resolution and it becomes unreasonable for my simple 2D geometry.

Thanks in advance,
__________________
Amir
Amir is offline   Reply With Quote

Old   November 16, 2011, 11:23
Default
  #2
Member
 
Patricio Bohorquez
Join Date: Mar 2009
Location: Jaén, Spain
Posts: 94
Rep Power: 8
pbohorquez is on a distinguished road
Hilary and Henry Weller implemented new polynomial-fit higher-order schemes (quadraticLinearFit, quadraticLinearPureUpwindFit, among others) in OF versions >= 1.6. For further details on the discretisation schemes see their papers:

http://www.met.reading.ac.uk/~sws02hs/publications/
http://www.met.reading.ac.uk/~sws02hs/AtmosFOAM/

However, spatial discretisation in OF is formally based on Taylor series including only the first two terms of the expansion. So I think that we are limited to 2nd order accuray in space. Recall that higher order schemes are difficult to be implemented in arbitrary polyhedral meshes.

Regards,
patricio.
pbohorquez is offline   Reply With Quote

Old   November 16, 2011, 13:14
Default
  #3
Senior Member
 
Amir's Avatar
 
Amir
Join Date: May 2009
Location: Shiraz, Iran
Posts: 739
Blog Entries: 1
Rep Power: 14
Amir is on a distinguished road
Quote:
Originally Posted by pbohorquez View Post
Hilary and Henry Weller implemented new polynomial-fit higher-order schemes (quadraticLinearFit, quadraticLinearPureUpwindFit, among others) in OF versions >= 1.6. For further details on the discretisation schemes see their papers:

http://www.met.reading.ac.uk/~sws02hs/publications/
http://www.met.reading.ac.uk/~sws02hs/AtmosFOAM/

However, spatial discretisation in OF is formally based on Taylor series including only the first two terms of the expansion. So I think that we are limited to 2nd order accuray in space. Recall that higher order schemes are difficult to be implemented in arbitrary polyhedral meshes.

Regards,
patricio.
Dear Patricio,

Thanks for your reply; certainly I'll check these schemes but before I take a look over these paper, do you know about their order or stability?
BTW, we have also SFCD as a second order scheme with good stability but it is said that it's very diffusive; do you have any experience with? because my convective term is very large so I prefer more diffusive schemes.

Bests,
__________________
Amir

Last edited by Amir; November 16, 2011 at 16:09.
Amir is offline   Reply With Quote

Old   November 16, 2011, 15:01
Default
  #4
Member
 
Patricio Bohorquez
Join Date: Mar 2009
Location: Jaén, Spain
Posts: 94
Rep Power: 8
pbohorquez is on a distinguished road
Hi Amir, the schemes are high order with exponent very close to 2.00. Curious, for convective dominated problem on simple geometries I usually prefer low diffusive WENO schemes (order >=5). Best wishes, patricio
pbohorquez is offline   Reply With Quote

Old   November 16, 2011, 16:29
Default
  #5
Senior Member
 
Amir's Avatar
 
Amir
Join Date: May 2009
Location: Shiraz, Iran
Posts: 739
Blog Entries: 1
Rep Power: 14
Amir is on a distinguished road
Quote:
Originally Posted by pbohorquez View Post
Hi Amir, the schemes are high order with exponent very close to 2.00. Curious, for convective dominated problem on simple geometries I usually prefer low diffusive WENO schemes (order >=5). Best wishes, patricio
Really!? I thought such high order schemes may lead to instability in convection dominated cases! Is there any WENO schemes in OpenFOAM?
__________________
Amir
Amir is offline   Reply With Quote

Old   November 17, 2011, 07:02
Default
  #6
Member
 
Patricio Bohorquez
Join Date: Mar 2009
Location: Jaén, Spain
Posts: 94
Rep Power: 8
pbohorquez is on a distinguished road
Quote:
Originally Posted by Amir View Post
Really!? I thought such high order schemes may lead to instability in convection dominated cases! Is there any WENO schemes in OpenFOAM?
Hi Amir, WENO works with success without introducing stability difficulties (Shu, 2009). Unfortunately, OF does not implement this family of numerical schemes that will improve the result of many solvers (e.g. interFoam)

C.-W. Shu, High order weighted essentially non-oscillatory schemes for convection dominated problems , SIAM Review, v51 (2009), pp.82-126.
songwukong likes this.
pbohorquez is offline   Reply With Quote

Old   November 17, 2011, 09:11
Default
  #7
Senior Member
 
Amir's Avatar
 
Amir
Join Date: May 2009
Location: Shiraz, Iran
Posts: 739
Blog Entries: 1
Rep Power: 14
Amir is on a distinguished road
Quote:
Originally Posted by pbohorquez View Post
Hi Amir, WENO works with success without introducing stability difficulties (Shu, 2009). Unfortunately, OF does not implement this family of numerical schemes that will improve the result of many solvers (e.g. interFoam)

C.-W. Shu, High order weighted essentially non-oscillatory schemes for convection dominated problems , SIAM Review, v51 (2009), pp.82-126.
Thanks for your time,
__________________
Amir
Amir is offline   Reply With Quote

Old   February 6, 2013, 18:01
Default
  #8
Member
 
dw
Join Date: Jul 2012
Posts: 32
Rep Power: 5
1/153 is on a distinguished road
Hi, what's your experience, among these higher order schemes, quadraticLinearFit, quadraticLinearPureUpwindFit, MUSCL, etc. which works better if I want to use higher order schemes for better accuracy? Say my flow is convection dominated.

Thanks a lot!
1/153 is offline   Reply With Quote

Old   May 16, 2013, 03:15
Default
  #9
Member
 
Join Date: Aug 2011
Posts: 74
Rep Power: 5
idefix is on a distinguished road
Hello,

I hope it´s correct to ask this question here.
In Jasak´s PhD thesis I found the following sentences:

Quote:
For good accuracy, it is necessary for the order of the discretisation to
be equal to or higher than the order of the equation that is being discretised
(page 77)

Can anyone explain me why it is so?

Thanks
Idefix
idefix 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
second order schemes marine OpenFOAM 66 October 26, 2014 12:10
when to change discretization schemes from first order to second order gyanesh FLUENT 0 April 23, 2011 05:23
Higher order schemes RameshK Main CFD Forum 0 July 15, 2010 22:35
CFL condition for higher order schemes Shyam Main CFD Forum 2 February 14, 2008 15:24
High order FVM scheme .... cfd_newbie FLUENT 0 February 4, 2008 14:43


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