diff --git a/applications/solvers/multiphase/compressibleInterFoam/TEqn.H b/applications/solvers/multiphase/compressibleInterFoam/TEqn.H index 135fb5ee8..bf7a4bfc2 100644 --- a/applications/solvers/multiphase/compressibleInterFoam/TEqn.H +++ b/applications/solvers/multiphase/compressibleInterFoam/TEqn.H @@ -4,8 +4,8 @@ fvm::ddt(rho, T) + fvm::div(rhoPhi, T) - fvm::Sp(contErr, T) - fvm::laplacian(turbulence.alphaEff(), T) + ( - fvc::div(fvc::absolute(phi, U), p)() - contErr/rho*p - + (fvc::ddt(rho, K) + fvc::div(rhoPhi, K))() - contErr*K + fvc::div(fvc::absolute(phi, U), p)()() - contErr/rho*p + + (fvc::ddt(rho, K) + fvc::div(rhoPhi, K))()() - contErr*K ) *( alpha1()/mixture.thermo1().Cv()()