
[Sponsors] 
March 13, 2006, 11:27 
Hello,
where can I find refer

#1 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 27 
Hello,
where can I find references for the limitedLinear scheme to cite in a paper? Thanks in advance, Alberto
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

March 13, 2006, 14:04 
There is no printed reference

#2 
Super Moderator
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 17 
There is no printed reference but it is Henry's application of the Sweby limiter to central differencing.


March 13, 2006, 18:57 
Thanks Mattijs.
Alberto

#3 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 27 
Thanks Mattijs.
Alberto
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

September 7, 2006, 09:42 
After looking into the code I

#4 
New Member
Steffen Jahnke
Join Date: Mar 2009
Posts: 14
Rep Power: 9 
After looking into the code I found in limitedLinear.H:
... return max(min(twor, 1), 0); ... with "twor" beeing equal to 2*r for default setup, e.g. "div(phi,B) Gauss limitedLinear 1;". For me this limiter looks very similar to the minmod limiter instead of the factor 2: minmod: return max(min(r,1),0); The Sweby limiter reads max(0,max(min(beta*r,1),min(beta,r))). I wondering why this modified limiter is used because due to the factor 2 it lies on the opposite limit (lower left) compared to the minmod in the Sweby diagram. 

February 27, 2008, 13:22 
Hello everyone.
As a newcom

#5 
New Member
nicasch
Join Date: Mar 2009
Posts: 23
Rep Power: 9 
Hello everyone.
As a newcomer I have two concrete questions concerning two divschemes:  Gauss limitedLinearV 1 and  Gauss interfaceCompression. I suppose that the interfaceCompression scheme does some blending between UD and CD, based on the blending factor (limiter) calculated according to either quartic or quadratic formula. Is this true? Or where is the limiter for this scheme used? In limitedLinearV, I suppose that it is a kind of limited CD, but how is the scheme limited? Where is the limiter for this scheme used? Is there any written text about these two schemes, or can anyone shortly explain me how they work? Many thanks in advance, best regards. 

March 1, 2013, 07:55 

#6 
New Member
Simon Tornros
Join Date: Aug 2012
Posts: 2
Rep Power: 0 
(sorry for bumping this old thread)
Why does the limiter in limitedLinear look like max(min(twoByk_*r, 1), 0); with twoByk_ = 2/k_ where 0 < k <1 when Sweby limiter is max(0,min(beta*r,1), min(r,beta)) with 1 < beta <2 ? From what I can see this is not the same. Last edited by tsimon; March 1, 2013 at 08:00. Reason: miscalculation.. 

March 12, 2013, 11:18 

#7  
New Member
Robin Debroux
Join Date: Oct 2012
Posts: 22
Rep Power: 6 
Quote:
I'm working with limitedLinear and I'm also trying to understand the limiter. Have you found any informations about it? Thank you. 

May 13, 2013, 10:18 

#8 
Member
Join Date: Mar 2013
Posts: 94
Rep Power: 5 
Hi to all,
I'm working with this limiter too and I would know if anybody found some information about it. Thank to all 

May 21, 2014, 22:02 

#9 
Member
Eric Robertson
Join Date: Jul 2012
Posts: 95
Rep Power: 6 
I think the twoByk_ term is just a calibrated term (i.e. a construct) that best suits what the limitedLinear scheme in OF is attempting to accomplish. I think it's to help place even more weight on higherorder differencing (in this case, central differencing vs upwind).
If you were to go through Jasak's thesis (pp. 9899) and go back to theory, one could see that this method is actually pretty effective in enforcing either central differencing or upwinding based on the userinput (k). One could actually show that, for a given ratio of successive gradients ("r"), inputting a low value of k (i.e. 0.1) would give more weight to central differencing, whereas a larger value (i.e. 0.8) places more emphasis on upwinding. I've also tried with beta values (1 < beta < 2) and the behavior is actually quite different. As such, I will have to conduct further studies. Also, just hypothetically speaking, I don't think there would be any difference between: Code:
max(min(twoByk_*r, 1), 0) Code:
max(0,min(twoByk_*r,1), min(r,twoByk_)) 

May 22, 2014, 00:56 

#11 
Member
Eric Robertson
Join Date: Jul 2012
Posts: 95
Rep Power: 6 

May 22, 2014, 04:23 

#12 
Senior Member
Tushar Chourushi
Join Date: Jul 2009
Location: IITIndore, India
Posts: 319
Blog Entries: 1
Rep Power: 10 
Hello All,
I have some doubts with the representation of the following schemes in OpenFOAM. I asked this question here because it seems many of you are working on these. (1) Code:
default limitedLinear 1.0 phi; (2) Code:
default limitedVanLeer 2.0 3.0; 

July 17, 2014, 04:52 

#13  
New Member
Marvin
Join Date: Jan 2014
Posts: 11
Rep Power: 4 
Quote:
Thanks for any help 

July 23, 2014, 00:51 

#14  
Senior Member
Tushar Chourushi
Join Date: Jul 2009
Location: IITIndore, India
Posts: 319
Blog Entries: 1
Rep Power: 10 
Quote:
It is TVD for any of the value ranging from 0 to 1. But, the order of accuracy decrease as the k value decrease (i.e., less than 1). For k =1, it is most accurate. I hope it cleared your doubt.  Best Luck! 

July 24, 2014, 05:20 

#15 
New Member
Marvin
Join Date: Jan 2014
Posts: 11
Rep Power: 4 
Tushar,
Looking at Versteeg p. 168: Sweby (1984) has given necessary and sufficient conditions for a scheme to be TVD in terms of the r − ψ relationship: • If 0<r<1 the upper limit is ψ(r)=2r, so for TVD schemes ψ®≤2r • If r≥1 the upper limit is ψ(r)=2, so for TVD schemes ψ(r)≤2 but twor is greater 2 which wouldn't satisfy the first criteria. Can you give some reference for you're statement? Thanks Marvin 

July 24, 2014, 08:26 

#16 
Senior Member
Tushar Chourushi
Join Date: Jul 2009
Location: IITIndore, India
Posts: 319
Blog Entries: 1
Rep Power: 10 
Marvin,
Did you read it thoroughly before commenting? Anyways, If you reread the same page no. 168169 of the mentioned ebook it will clear your doubt. I have attached a file of the same book with marked lines of page no. 168, which clears explains the TVD region. Please find attached note with this comment.  Best Regards! 

September 8, 2014, 04:38 

#17  
New Member
Marvin
Join Date: Jan 2014
Posts: 11
Rep Power: 4 
Quote:
(compare http://www.openfoam.org/docs/user/fvSchemes.php) > For k = 1 it is most stable. For k = 0 it is most accurate. Sorry I didn't put that right earlier. Further do I still not see why limitedLinear is TVD. sourcecode: return max(min(twoByk_*r,1),0); (1) If I draw that function (1) for k = 0.5 that is definitely not in the shaded region marked in verteeg for k <1. It demands that twoByk_ < 2. 

September 30, 2014, 20:23 

#18 
New Member
Gecamp
Join Date: Oct 2010
Posts: 13
Rep Power: 8 
Hi Marvin,
I might be wrong here, but I think that the way the limiter has been implemented allows the scheme to behave as a TVD. In fact, with k=1, your limiter will be ψ(r) = 2*r. (as long as 2*r < 1). Conversely, with k = 0 (which means == SMALL), your ψ(r) becomes somethingBIG*r. This is most likely >> 1 and 'the switch' will choose ψ(r) = 1 as limiter. When such thing occurs, you are using pure CD to interpolate phi values on the internal faces. Since I believe that the value r = (Phi_C  Phi_U)/(Phi_D  Phi_C) changes during the run, I expect the limiter ψ(r) to do the same. By choosing k = 1, one just increases the possibility to include in the interpolation process some UD diffusivity, I think. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Help!explicit multistage scheme with upwind scheme  D .T.  Main CFD Forum  1  June 4, 2003 03:42 
References  Lego  Main CFD Forum  0  February 26, 2002 10:57 
References on DNS/LES  J. Xia  Main CFD Forum  5  April 6, 2001 11:28 
references  sundar  Main CFD Forum  1  June 21, 2000 12:47 
References for MacCormack scheme  L Srinivasa Mohan  Main CFD Forum  1  March 10, 2000 14:45 