solve(fvm::ddt(rhoH)+ divrhoHFlux+ divPFlux- fvc::laplacian(turbulence->alphaEff(), hs) - fvc::div(sigmaDotU)==fvc::ddt(p)+ combustion->Sh()
combustion->Sh()