CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

Wrong treatment of time derivatives in matrices (Jasak Presentation)

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

LinkBack Thread Tools Display Modes
Old   February 1, 2017, 16:44
Default Wrong treatment of time derivatives in matrices (Jasak Presentation)
Senior Member
Tobi's Avatar
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,662
Blog Entries: 6
Rep Power: 31
Tobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Dear community,

today I got an interesting presentation that was held last year Hrvoje Jasak with the title Pressure-Velocity Coupling in FOAM - Consistent Derivation for Steady and Transient Flow Solvers. Here it is described that the included time term in the matrices:

fvVectorMatrix UEqn
lead to wrong or lets say time-step depended results. Therefore, I came across this problem again (already forgot it). At the moment I did not read everything but I will do that. It seems that the included ddt term makes troubles in both, transient and steady-state solvers and therefore the ddtCorr() function takes care of removing the troubles somehow.

Can someone provide some information about that and does anyone know if this will be re-implemented in a correct and consistent way in future like Jasak suggested? However, there is one question left because based on the fact that the developers introduced the ddtCorr() function, they know that there is a problem (more or less), so why do they patch it instead of fix it. Maybe there is some special meaning of doing it in that way or there is no other possible way (but again, Jasak demonstrated it).

In addition I want ask kindly if this ddtCorr is related to the Rhie Chow Interpolation because something is mentioned here: DdtPhiCorr However, here it is said that this stuff is consistent while Jasak mentioned in his presentation, that it is not consistent. In addition I will check out the paper that is mentioned in the above mentioned thread.

Thanks in advance.
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote



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
decomposePar problem: Cell 0contains face labels out of range vaina74 OpenFOAM Pre-Processing 27 November 2, 2015 18:04
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
Unstabil Simulation with chtMultiRegionFoam mbay101 OpenFOAM Running, Solving & CFD 13 December 28, 2013 14:12
pisoFoam with k-epsilon turb blows up - Some questions Heroic OpenFOAM Running, Solving & CFD 26 December 17, 2012 04:34
AMI speed performance danny123 OpenFOAM 19 October 24, 2012 07:44

All times are GMT -4. The time now is 13:48.