[kaifu]

Hi Foamers

I have a source term Gamma on RHS of energy equation, which of course is a type of volScalarField. However this source term only appears in near wall cells and it is zero elsewhere. In physics, this source term should only come from the wall, which should be a surfaceScalarField. The flux of Gamma on the wall surface should be:

Gamma*volume of cell

However in energy equation, every term should be written in form of volScalarField (may interpolate to surface later). In such situation, when using central differencing scheme, the surface flux on boundary (wall) may be smaller than it should be. And what is worse, the interpolated value Gamma on the other side of surface will not be zero. (The near wall cells are all cubics) This is definitely not acceptable in physics. I have considered using upwind differencing scheme for Gamma. It is easy to make the flux of Gamma to be zero on the other side of cells (not wall surface) but the problem is how could I interpolate Gamma to the wall boundary? Thx.

// Kai

[kaifu]

Some updates: My mistake. All the explicit term in energy equation will be treated like

Gamma*volume of cell

which means there is no interpolation of Gamma onto cell faces. Only the implicit one will be interpolated.

// Kai