Incompatible Field of operation
 

I want to include pressure in my energy eqn in the form below. It compiles fine but whenever I run it the case, it gives an error
"incompatible fields of operation
[T] + [p_rgh]"

I guess it means I can only solve for one field per scalar matrix. How can I go around this to include pressure in the energy eqn? Thanks.

fvScalarMatrix TEqn
(
fvm::ddt(rhoCp,T)
+ fvm::div(rhoPhiCpf,T)
+ fvm::div(phi,p_rgh)
- fvm::laplacian(kappa,T)
);

Including it explicitly might work, but I don't think that is what you are after.

Thanks. it works when I solve it explicitly but according to your statement, what should I look out for?