# TVD Schemes in OF

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

 January 14, 2010, 10:06 TVD Schemes in OF #1 Member   Joern Bader Join Date: Mar 2009 Posts: 33 Rep Power: 9 Hi, I read the OpenFOAM Code parts for the TVD-Schemes and i have some questions about it. In literature the TVD-schemes 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: 9 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 20, 2010, 16:54
#4
Member

Join Date: Mar 2009
Posts: 33
Rep Power: 9
Quote:
 Originally Posted by santiagomarquezd 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).
not i do so, OF does.

Quote:
 Originally Posted by santiagomarquezd 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.
Thats right, r has to be calculated every timestep but most literature, i found, either used a complete explicit TVD sheme or ignored the time index in their text.
When i worte that i just meant the syntax.

Quote:
 Originally Posted by santiagomarquezd 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?
TVD does a mix of high and low schemes by using factors z and 1-z. So you can get values from full-high- to full-low-scheme.

Quote:
 Originally Posted by santiagomarquezd 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.
my literature was:
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:995-1011,
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: 430
Rep Power: 16

Quote:
 fade? thats what TVD does. TVD does a mix of high and low schemes by using factors z and 1-z. So you can get values from full-high- to full-low-scheme.
That's that I supposed, I think in FOAM framework it is called blending (HRV PhD Ts. chap. 3.4.3).

Quote:
 H. K. Versteeg, W. Malalasekera. An Introduction to Computational Fluid Dynamics: the fnite volume method. Pearson Edudacion Limited, 2. Aufage, 2007
Except for this book I have the literature you suggested. I've read previous edition, but Ola Widlund told me that it is improved now. Do you have an electronic copy of it (my mail is in the contact information page)?

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Research Scientist
Research Center for Computational Methods (CIMEC) - CONICET/UNL
Tel: 54-342-4511594 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

Join Date: Mar 2009
Posts: 33
Rep Power: 9
Quote:
 Originally Posted by santiagomarquezd Except for this book I have the literature you suggested. I've read previous edition, but Ola Widlund told me that it is improved now. Do you have an electronic copy of it (my mail is in the contact information page)? Regards.
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 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 Shuo Main CFD Forum 5 September 19, 2008 21:47 Shuo Main CFD Forum 0 August 19, 2008 21:51 Ritesh Main CFD Forum 4 June 16, 2005 10:12 Sudeep FLUENT 0 October 14, 2000 21:31 Matt Umbel Main CFD Forum 1 January 25, 2000 05:21

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