# Flux Limiter and Interpolation Scheme

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

 March 23, 2015, 15:54 Flux Limiter and Interpolation Scheme #1 New Member   Paolo Capobianchi Join Date: Sep 2014 Posts: 11 Rep Power: 2 Dear All, I have been working for a while on the implementation of a new FLux Limiter in OpneFOAM, and after some attempts it seems that I found they way to do that. I have compiled it and it seems that it works properly (at least it allowed me to get reasonable results). However, I have a doubt, tthus I would ask if someone could be so kind to help me. I know that a flux limiter can be seen as a function that allows a "blended" scheme to behave like an UD scheme or like a HR scheme in different region of the flow field depending on the particular conditions in the flow field itself. In principle (I guess), the HR scheme could be whatever you want, but in my case I need to use the QUICK scheme of Leonard with a particular limiter. Is it possible to specify this scheme in the field "interpolationScheme" into the "fvScheme" of the case? (I tried but it seems that it does not work. Probably I used wrong settings) Does anyone know what is the default HR scheme? Maybe the linear? Thank you very much in avantage! Paolo

 March 24, 2015, 00:10 #2 Senior Member   Tushar Chourushi Join Date: Jul 2009 Location: IIT-Indore, India Posts: 318 Blog Entries: 1 Rep Power: 8 I think GAMMA scheme is mostly used for OpenFOAM work - Best Luck! __________________ Mistakes Are Inevitable, You Can Only Grow Stronger If You Lean From Them

 March 24, 2015, 08:28 #3 New Member   Paolo Capobianchi Join Date: Sep 2014 Posts: 11 Rep Power: 2 Hi Tushar, Thanks a lot for the reply! In practice, my problem is the following. The flux at the cell face can be calculated as: (phi)f = (phi)_UD + psi(r)*[(phi)_HR + (phi)_UD] where (phi)_HR is the flux at the face "f" evaluated whit a Higher Order scheme, which in my case should be the QUICK scheme, and psi(r) is the flux limiter that I have added on my version of OpenFOAM. by considering what you said, I am not allowed to use a QUICK scheme for the interpolation at the cell face? However I am a bit confised. In fact, I had have a read to to an article in OpenFOAM Wiki, and I found the this: The surfaceInterpolation classes implemented in the finite volume library perform interpolation from volume fields to face fields, a critical calculation in the discretization process for the finite volume method, particularly for co-located meshes, as is used by OpenFOAM. OpenFOAM implements dozens of different schemes for this operation. This kind of interpolation is denoted: Notation used by Marupio in this wiki 1 Schemes The schemes are selected in the fvSchemes file, and loaded using runTime selection. There are more than fifty schemes available in OpenFOAM-extend, including: QUICK - quadratic upwind interpolation; linear - central differencing scheme; upwind - upwind differencing scheme; and skewCorrected - skew upwind scheme. and seems that would be possible to select the QUICK scheme, but I do not understand how to do that. Does anyone have any suggestion? Thanks a lot! Paolo

 March 25, 2015, 00:31 #4 Senior Member   Tushar Chourushi Join Date: Jul 2009 Location: IIT-Indore, India Posts: 318 Blog Entries: 1 Rep Power: 8 Refer the following http://www.openfoam.org/docs/user/fvSchemes.php check out section 4.4.5, you can apply the above said scheme as: Code: `div(phi,U) Gauss QUICK;` I think this will work For more information refer the cfd-online post (I guess you know this): Numerical schemes: under the hood - Best Luck! __________________ Mistakes Are Inevitable, You Can Only Grow Stronger If You Lean From Them Last edited by Tushar@cfd; March 25, 2015 at 01:58. Reason: Added line

March 25, 2015, 07:06
#5
New Member

Paolo Capobianchi
Join Date: Sep 2014
Posts: 11
Rep Power: 2
Hi Tushar,

Thanks again for helping me! I am afraid that when you specify "Gauss QUICK" on the fvSchemes dictionary, you are actually specifying the QUICK limiter. I have been looking for a long time the QUICK scheme inside OpenFOAM and the only files I found about it are actually inside the following foleder:

src/fiiteVolume/interpolation/surfaceInterpolation/limitedSchemes/QUICK

on which there are two files: QUICK.H and QUICK.C

I am pretty sure that those files define the QUICK flux limiter and not the quadratic interpolation scheme. In fact as you can see here:
Quote:

http://en.wikipedia.org/wiki/Flux_limiter

My goal is to specify a limiter in this way:

PHP Code:
``` div(phi,U)        Gauss My_Limiter;  ```
and use the quadratic interpolation for the face flux as a Higher Order scheme.

Best Wishes,
Paolo

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post yeharav OpenFOAM Programming & Development 2 May 15, 2014 09:21 vkrastev OpenFOAM Running, Solving & CFD 1 May 29, 2012 09:04 glaspina OpenFOAM Programming & Development 0 January 26, 2012 05:42 mchurchf OpenFOAM Running, Solving & CFD 3 November 17, 2009 00:48 alberto OpenFOAM Running, Solving & CFD 0 March 13, 2008 22:18

All times are GMT -4. The time now is 23:39.