
[Sponsors] 
January 14, 2010, 10:06 
TVD Schemes in OF

#1 
Member
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 10 
Hi,
I read the OpenFOAM Code parts for the TVDSchemes and i have some questions about it. In literature the TVDschemes are derived in a time independend manner. For a value U it is: U_f = (U_f)_low + limiter(r)*((U_f)_high  (U_f)_low)
with (upwind direction) r = (U_P  U_W)/(U_E  U_P)In OpenFOAM r is calculated like Darwish and Moukalled recommend. In my literature these terms are time independend. In OF the calculation is time dependend (clear), but the factor r is calculated explicit but it is used for an implicit fade. In other words: r is calculated for a time t^n but it is used in a fade of (U_f)_high and (U_f)_low at the time t^(n+1). Is this still consistent to literature? Why? Why does this scheme still satisfy the TVD condition? I suppose its still a TVD Scheme, because the limiter is time independed and still satisfy Swebys diagram. Is this true? I hope someone can help, because i need to explain the reason, why this way of calculation is consistent. Thx. 

January 29, 2010, 09:50 

#2 
Member
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 10 
Can anyone help me?
i need help. OF uses the TVD schemes in an implicit way, but the factor r for the limiter is calculated explicit. why is this kind of scheme still TVD? In literature i only find full explicit or full implicit schemes . 

June 19, 2010, 23:20 

#3 
Senior Member
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 438
Rep Power: 17 
Hello Joern, I started to study NVD/TVD limiters in FOAM recently, I have some comments and questions about this topic.
1. As you said r is calculated explicitly from values of U_i^n. Next you assemble time terms and solve the system for U_i^(n+1). 2. Why you propose, "In my literature these terms are time independent"? TVD criterion have to be accomplished every timestep, so in each timestep it is neccessary to recalculate r an then psi(r) to limite the flux. r is different across the time because the changes in U. 3. About this paragraph: "In OF the calculation is time dependend (clear), but the factor r is calculated explicit but it is used for an implicit fade. In other words: r is calculated for a time t^n but it is used in a fade of (U_f)_high and (U_f)_low at the time t^(n+1)." What is fading? 4. What literature do have? I used the book from Hirsch, Jasak's Thesis and original papers. What literature have been helpful for you?, specially in implicit TVD. Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D. Research Scientist Research Center for Computational Methods (CIMEC)  CONICET/UNL Tel: 543424511594 Int. 7032 Colectora Ruta Nac. 168 / Paraje El Pozo (3000) Santa Fe  Argentina. http://www.cimec.org.ar 

June 20, 2010, 16:54 

#4  
Member
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 10 
Quote:
Quote:
When i worte that i just meant the syntax. Quote:
TVD does a mix of high and low schemes by using factors z and 1z. So you can get values from fullhigh to fulllowscheme. Quote:
M. S. Darwish, F. Moukalled. TVD schemes for unstructured grids. International Journal of Heat and Mass Transfer 46, S. 599611, 2003 J. H. Ferziger, M. Peric. Computational Methods for Fluid Dynamics. Springer Verlag, 2. Au age, 1997 P. K. Sweby. High resolution schemes using flux limiters for hyperbolic conservation laws. SIAM J. Numer. Analysis, 21:9951011, 1984 H. K. Versteeg, W. Malalasekera. An Introduction to Computational Fluid Dynamics: the fnite volume method. Pearson Edudacion Limited, 2. Aufage, 2007 and some other papers i think. 

June 28, 2010, 21:28 

#5  
Senior Member
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 438
Rep Power: 17 
Hello Joern, thanks for your reply,
Quote:
Quote:
Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D. Research Scientist Research Center for Computational Methods (CIMEC)  CONICET/UNL Tel: 543424511594 Int. 7032 Colectora Ruta Nac. 168 / Paraje El Pozo (3000) Santa Fe  Argentina. http://www.cimec.org.ar 

July 7, 2010, 15:34 

#6 
Member
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 10 
Sry, i have no electronic copy. i found the book in our bib but not in the part for math. i found it within books for mechanical engeneering.


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
TVD schemes for acoustics  Shuo  Main CFD Forum  5  September 19, 2008 21:47 
TVD schemes: 2 questions  Shuo  Main CFD Forum  0  August 19, 2008 21:51 
TVD Schemes  Ritesh  Main CFD Forum  4  June 16, 2005 10:12 
TVD Schemes  Sudeep  FLUENT  0  October 14, 2000 21:31 
TVD Schemes  Matt Umbel  Main CFD Forum  1  January 25, 2000 05:21 