June 12, 2012, 09:51
Join Date: May 2012
Location: Calumet, MI
Rep Power: 5
Originally Posted by grtabor
Its to do with the discretisation of the
U & fvc::grad(T)
term. In the FVM we discretise this using Gauss theorem, converting it into fluxes into and out of the domain. If we discretise the term written like this then it is not guaranteed to preserve conservation of the quantity being advected (here essentially the energy); i.e. your derivative may be loosing or gaining energy. If you do some manipulation of the term to get it into the form
where phi is essentially the same as U, then when we apply Gauss we get something which preserves continuity at the _numerical_ level, rather than the mathematical level. This would be a strong conservative implementation.
Forgive my lack of mathematical background, but is there someone that would be able to show
the manipulation steps such that
I've taken a few whacks at it and my lack of FVM, I think, is making me miss some obvious connection.
Many thanks to whomever is able to assist me!
Edit: Nevermind. Figured it out. Odd way, though it works.
Last edited by iamed18; June 15, 2012 at 09:47.
Reason: added definition of [MATH]\phi[/MATH], figured it out.