[Bufacchi]

Hi Foamers,

I'm using OpenFOAM 2.1.x.

I noticed that in the file reactingOneDim.C, in the energy equation:

fvScalarMatrix TEqn
(
fvm::ddt(rhoCp, T_)
- fvm::laplacian(K_, T_)
==
chemistrySh_
+ fvc::div(phiQr)
+ fvc::div(phiGas)
);

the components ddt(rhoCp, T_), laplacian(K_, T_) and chemistrySh_ all have the same dimension = [1 -1 -3 0 0 0 0] or J/s/m3, but the components div(phiQr) and div(phiGas) have a different dimension = [1 1 -3 0 0 0 0] or J/s/m.

Does anyone know the reason for that?
Am I missing something?

Regards,

Paulo