[S.Samridh]

Hello,

I have been trying to add div(sigmaDotU), energy term, which defines mechanical energy dissipation in fluid. Checking various threads, i haven't found any energy solver being used for icoFoam solver.
The rhoCentralFoam has this equation being solved, but has many functions in tow, to complete energy equation there.

I need help to add this(code below) specifically to icoFoam, and I am not sure about the header files that should be included or the proper framing of this code, any reference example would be a great help (I know how to add Temperature solver in icoFoam, but as energy dissipation term wouldn't need initial conditions, i am confused if it's added in createfields.H)

From rhoCentralFoam

// --- Solve energy
202 surfaceScalarField sigmaDotU
203 (
204 "sigmaDotU",
205 (
206 fvc::interpolate(muEff)*mesh.magSf()*fvc::snGrad(U )
207 + fvc::dotInterpolate(mesh.Sf(), tauMC)
208 )
209 & (a_pos*U_pos + a_neg*U_neg)
210 );
211
212 solve
213 (
214 fvm::ddt(rhoE)
215 + fvc::div(phiEp)
216 - fvc::div(sigmaDotU)
217 );

Out of above written code, I need to solve div(sigmaDotU) in icoFoam.

Thank you,
Samridh