One question about the discretization of governing equation in OpenFOAM
I found there is an extra term in some solvers. I am not sure about the reason. For example in the TEqn.H of buoyantBoussinesqSimpleFoam. The scalar matrix is constructed via
- fvm::Sp(fvc::div(phi), T)
- fvm::laplacian(kappaEff, T)
The second term seems not necessary here. And it also make the discretization become nonconservative. I wonder why this is necessary. In contrast in the UEqn.H in the same solver
Such term does not appear. Can anyone explain why we need this?
Thanks a lot.
The line is added to have values of T bounded in the algorithms like SIMPLE where we have partial convergence of pressure. Check this: http://www.cfd-online.com/Forums/ope...silon-eqn.html
|All times are GMT -4. The time now is 18:11.|