errors about div(phi,h) when move 2.0.1solver to openfoam2.2.1
 

Hi ,Foamers
I moves fireFoam of 2.0.1 to 2.2.1,then rewrite hEqn according to the EEqn.H in fireFoam of 2.2.1 and other necessary codes.But ,when I wmake my rewrited fireFoam,it gives me errors:undefine div(phi,h) in fvscheme.
The following is codes of EEqn.H:
{
volScalarField& he = thermo.he();
fvScalarMatrix EEqn
(
fvm::ddt(rho, he) + fvm::div(phi, he)
+ fvc::ddt(rho, K) + fvc::div(phi, K)
+ (
he.name() == "e"
? fvc::div
(
fvc::absolute(phi/fvc::interpolate(rho), U),
p,
"div(phiv,p)"
)
: -DpDt
)
- fvm::laplacian(turbulence->alphaEff(), he)
==
fvOptions(rho, he)
);
EEqn.relax();
fvOptions.constrain(EEqn);
EEqn.solve();
fvOptions.correct(he);
thermo.correct();
Info<< "T min/max = " << min(T).value() << ", " << max(T).value() << endl;
}
There is div(phi,he) insteading of div(phi,h),So It should be no need to define (phi,h) in fvshceme,but define div(phi,he).However,I see the fvsheme in smallpoolfire2D of 2.2.1,there indeed is div(phi,h).In fact,there is the term of "div(phi,he)" in EEqn.H,why define div(phi,h) in fvsheme?I think I have some thing not understand.Can some kinds foamers tell me?I am very anxious!I have report to my teacher soon!

Could someones kinds help me about the problem?