CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - DDt (https://www.cfd-online.com/Forums/openfoam-solving/58580-ddt.html)

Quote:

Originally Posted by grtabor (Post 196218)

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

$\vec{U}\cdot\nabla\vec{U}=\nabla\cdot(\phi\vec{U})$ where $\phi=\vec{U}\cdot\vec{S}_f$ ?

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.

.oldTime(); in chtMultiregionFoam

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’

Re: How to write material derivative of ln(T) in the consitutive equation

Hi Dear Former's
I'm going to Add temperature Equation to viscoelasticFluidFoam solver.
there is material derivative term of log(T)(screenshot attached) in the constitutive equation(I.e in this case Olroyd-B)

kindly any could guide me how to erite the that term in OldoydB.C file

I write it like fvm::ddt(tau_,logT) + fvc::div(phi()*tau_,logT)

Note:change to fvc and fvm respectively but getting same error's.(scrrenshot Attached)

regards
Idrees khan

How to write material derivative of ln(T) in the consitutive equation