CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

How to implement power law or hybrid schemes??

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

Like Tree5Likes
  • 1 Post By kamkari
  • 2 Post By deepsterblue
  • 1 Post By kamkari
  • 1 Post By arjun

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 7, 2011, 18:43
Default How to implement power law or hybrid schemes??
  #1
New Member
 
babak kamkari
Join Date: Dec 2010
Posts: 26
Rep Power: 15
kamkari is on a distinguished road
Hi every body
I am going to model a flow field in an enclosure by using power law scheme or hybrid scheme.
But in fvScheme I couldn’t find any of these schemes.
As far as I seen in openFoam solvers the most used scheme for divergence discritization are upwind scheme and linear scheme (also TVD, vanleer …). But these two schemes have some defects.
Upwind scheme: this scheme impose artificial dispersion (false diffusion) which Leeds to unrealistic solution.
Linear scheme (central differencing): causes oscillating answer for high Reynolds numbers and unrealistic solution.
Using power law scheme could overcome these infections.
Does anyone know how these schemes could be implemented?
Thanks for consideration
kamkari is offline   Reply With Quote

Old   May 8, 2011, 08:09
Default
  #2
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 620
Blog Entries: 6
Rep Power: 24
elvis will become famous soon enough
Hi,

http://openfoam-extend.svn.sourcefor.../viscoelastic/

http://openfoam-extend.svn.sourcefor...coelasticLaws/

http://openfoam-extend.svn.sourcefor...oelasticModel/

hope this goes in the direction you are looking for.
elvis is offline   Reply With Quote

Old   May 8, 2011, 11:07
Default
  #3
New Member
 
babak kamkari
Join Date: Dec 2010
Posts: 26
Rep Power: 15
kamkari is on a distinguished road
Hi Elvis
Thank you very much for your answer.
I think I couldn’t convey my problem clearly. Power law scheme is not related to viscoelastic modeling. It is on the most famous way for discretizing the convective terms in momentum equation (like second order upwind, central, quick, ...) . It is also one of the special options in Fluent software for discretizing the equations. I couldn’t find this scheme as one of the standard scheme in fvscheme file.
May be it exist in openFoam library with another name!!!!
mohsen cheraghi likes this.
kamkari is offline   Reply With Quote

Old   May 8, 2011, 12:40
Default
  #4
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 403
Rep Power: 25
deepsterblue will become famous soon enough
Could you point to a published manuscript that describes the scheme? Does it belong to the TVD/NVD family?

Edit: Never mind. Found the sources in Fluent's User Guide. This looks like a relatively straightforward implementation. Derive a class from upwind, compute the local Peclet number, and assign phi_f on a face-by-face basis. The scheme is still first-order, so I'm not sure if there's any real accuracy benefit to this, besides stability.
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon

Last edited by deepsterblue; May 8, 2011 at 12:50. Reason: Found source
deepsterblue is offline   Reply With Quote

Old   May 8, 2011, 13:13
Default
  #5
New Member
 
babak kamkari
Join Date: Dec 2010
Posts: 26
Rep Power: 15
kamkari is on a distinguished road
Following link is Fluent user guide which explain discretization schemes including power law



http://my.fit.edu/itresources/manual...ug/node992.htm
kamkari is offline   Reply With Quote

Old   May 8, 2011, 16:03
Default
  #6
New Member
 
babak kamkari
Join Date: Dec 2010
Posts: 26
Rep Power: 15
kamkari is on a distinguished road
May be it is better to ask which divergence scheme is closer to power law discrertization scheme.
Power law scheme is stable and doesn’t impose false diffusion to the solution. So it can lead to accurate results for convective terms in momentum equation.
amin_jalalian likes this.
kamkari is offline   Reply With Quote

Old   May 8, 2011, 18:34
Default
  #7
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by kamkari View Post
May be it is better to ask which divergence scheme is closer to power law discrertization scheme.
Power law scheme is stable and doesn’t impose false diffusion to the solution. So it can lead to accurate results for convective terms in momentum equation.

This is true only in 1D. Why? Because powerlaw scheme has basis on exact solution of transport equation in 1D. In 2D and 3D though, it loses its shine. Still I am big fan of hybrid schemes, they are much better than upwind type schemes.
mm.abdollahzadeh likes this.
arjun is offline   Reply With Quote

Old   May 9, 2011, 00:51
Default
  #8
New Member
 
babak kamkari
Join Date: Dec 2010
Posts: 26
Rep Power: 15
kamkari is on a distinguished road
Hi Arjun
Thank you for your reply,
What kind of div scheme do you recommend for applying hybrid scheme? Having both acuurecy and minimum diffusion
Thank you in advance
kamkari is offline   Reply With Quote

Old   May 9, 2011, 01:02
Default
  #9
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by kamkari View Post
Hi Arjun
Thank you for your reply,
What kind of div scheme do you recommend for applying hybrid scheme? Having both acuurecy and minimum diffusion
Thank you in advance

if you have access to Patankar's book, I think he explains how to impliment powerlaw and hybrid schemes.

He also wrote about how and why one should prefer hybrid schemes over other schemes. That is an interesting read.
arjun is offline   Reply With Quote

Old   May 9, 2011, 10:09
Default
  #10
Senior Member
 
Ben K
Join Date: Feb 2010
Location: Ottawa, Canada
Posts: 140
Rep Power: 19
benk is on a distinguished road
I'm also interested in knowing how to implement power law/hybrid schemes, since these are the schemes that I've studied before and are well documented in CFD books. A good summary of these schemes is available in the book Versteeg and Malalasekera (An intro to computational fluid dynamics), chapter 5. QUICK is documented in Versteeg (and available in openFOAM) but they say that it may not be stable in certain circumstances.

When I first started to use openFOAM I was surprised not to see these schemes as options, but I also saw schemes that I didn't recognize so I figured that they might be implemented in other ways via blending in some way...still not sure.
benk is offline   Reply With Quote

Old   May 9, 2011, 11:29
Default
  #11
Senior Member
 
Mirko Vukovic
Join Date: Mar 2009
Posts: 159
Rep Power: 17
mirko is on a distinguished road
Quote:
Originally Posted by kamkari View Post
Hi every body
I am going to model a flow field in an enclosure by using power law scheme or hybrid scheme.
But in fvScheme I couldn’t find any of these schemes.
As far as I seen in openFoam solvers the most used scheme for divergence discritization are upwind scheme and linear scheme (also TVD, vanleer …). But these two schemes have some defects.
Upwind scheme: this scheme impose artificial dispersion (false diffusion) which Leeds to unrealistic solution.
Linear scheme (central differencing): causes oscillating answer for high Reynolds numbers and unrealistic solution.
Using power law scheme could overcome these infections.
Does anyone know how these schemes could be implemented?
Thanks for consideration
When I was starting out in this field, Patankar's book was also my first reference, and my first code used the hybrid/power-law scheme.

In my second reference, much to my surprise, Ferziger and Peric strongly discourage its use (Computational Methods for Fluid Dynamics p. 64). I learned a lot from reading articles by Prof. Brian Leonard. For example, see (sorry about fonts -- cutting and pasting issues)
Why you should not use ‘Hybrid’, ‘Power-Law’ or related exponential schemes for convective modelling—there are much better alternatives

International Journal for Numerical Methods in Fluids

Special Issue: Numerical Methods for Thermal Problems
Volume 20, Issue 6, pages 421–442, 30 March 1995


As others have mentioned, the power-law/hybrid scheme is good for sourceless 1D convection-diffusion problems.

Once you add sources and move to 2D, you are better off using other methods. Here I defer to others to suggest them.

One of Fluent's basic tutorials (if you have access to it) demonstrates it quite nicely: fluid flow in an bend with a hot fluid added at a second point in the bend. The hybrid solution diffuses the added fluid much more than a second order scheme does.

Mirko
mirko is offline   Reply With Quote

Old   September 30, 2015, 15:09
Default
  #12
Member
 
amin jalalian
Join Date: Mar 2014
Posts: 34
Rep Power: 12
amin_jalalian is on a distinguished road
Quote:
Originally Posted by kamkari View Post
Hi Arjun
Thank you for your reply,
What kind of div scheme do you recommend for applying hybrid scheme? Having both acuurecy and minimum diffusion
Thank you in advance
Hi
see this
http://cfd.direct/openfoam/user-guid...s/#x20-1100276
amin_jalalian is offline   Reply With Quote

Old   September 30, 2015, 15:23
Default
  #13
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by mirko View Post
In my second reference, much to my surprise, Ferziger and Peric strongly discourage its use (Computational Methods for Fluid Dynamics p. 64)
The reason is that these schemes are derived from 1d solution and thus does not go well in 2d andf three d.
arjun is offline   Reply With Quote

Old   January 13, 2018, 07:35
Default
  #14
Senior Member
 
A. Min
Join Date: Mar 2015
Posts: 305
Rep Power: 12
alimea is on a distinguished road
Quote:
Originally Posted by arjun View Post
This is true only in 1D. Why? Because powerlaw scheme has basis on exact solution of transport equation in 1D. In 2D and 3D though, it loses its shine. Still I am big fan of hybrid schemes, they are much better than upwind type schemes.
Hi arjun
I agree with you that "hybrid" is the best choice for div terms. but why isn't there hybrid scheme in openfoam?
if it is available what is the name of that?
my foam is 2.3.1. is it in new versions? what should I do for download and install it?

I heared that limitedLinear is hybrid! is it correct? if it is correct, there is a problem here: it has limiter and I don't want to use limiter! because it dissipates a lot.

thanks
alimea is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Non-Newtonian Power Law for Viscosity John FLUENT 16 September 12, 2015 07:18
Non-Newtonian Flow: Power law (2D BFS) hellouch FLUENT 3 April 21, 2009 04:03
about one seventh power law in CFX Eric CFX 3 January 10, 2007 08:16
power law flow through porous media Francesca Chiusa FLUENT 1 November 3, 2006 15:30
Non newtonian power law model in pipe annulus Poly FLUENT 0 March 24, 2005 16:12


All times are GMT -4. The time now is 16:05.