[Parisa_Khiabani]

Hello everyone,
The TEqn.H in my solver is as below:

surfaceScalarField kappaf = twoPhaseProperties->kappaf();

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

TEqn.solve();

in system/fvSchemes, I defined:
div (rhophiCpf,T) Gauss upwind;
However, after one iteration, I get a fatal error as: keyword div(rho*phi*Cpf,T) is undefined in dictionary.
I also change the initial form of div (rhophiCpf,T) to div(rho*phi*Cpf,T), but the same error happened. That's weird because I did define the divergence.
I really appreciate if you guys can help me.

Best,
Parisa