[iamed18]

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

fvc::div(phi, T)

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.

Gavin

Forgive my lack of mathematical background, but is there someone that would be able to show the manipulation steps such that

where ?

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.

[Mohammad Jam]

Hi

How to implement .oldTime(); after writing volScalarField in chtMultiregionFoam?

When i do that it gives this error:
class Foam::PtrList<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >’ has no member named ‘oldTime’