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

TVD Schemes in OF

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 14, 2010, 10:06
Default TVD Schemes in OF
  #1
Member
 
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 8
joern is on a distinguished road
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.

joern is offline   Reply With Quote

Old   January 29, 2010, 09:50
Default
  #2
Member
 
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 8
joern is on a distinguished road
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 -.-
joern is offline   Reply With Quote

Old   June 19, 2010, 23:20
Default
  #3
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 14
santiagomarquezd will become famous soon enough
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.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   June 20, 2010, 16:54
Default
  #4
Member
 
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 8
joern is on a distinguished road
Quote:
Originally Posted by santiagomarquezd View Post
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 View Post
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 View Post
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?
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.

Quote:
Originally Posted by santiagomarquezd View Post
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.
joern is offline   Reply With Quote

Old   June 28, 2010, 21:28
Default
  #5
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 14
santiagomarquezd will become famous soon enough
Hello Joern, thanks for your reply,

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.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   July 7, 2010, 15:34
Default
  #6
Member
 
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 8
joern is on a distinguished road
Quote:
Originally Posted by santiagomarquezd View Post
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.
joern is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


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


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