[yzf1215]

Hello, world. I have costumed the enthalpy in hEqn.H which is included rhoSimpleFoam. The problem is, I inserted a scalarField Q in it. But it is a steady-state solver, but I don't know the relation of scalarField Q and the heat(unit Joule) I need to add per second per volume. And, what should I set the Q's unit?

The hEqn.H I costumed is:
{
fvScalarMatrix hEqn
(
fvm::div(phi, h)
- fvm::Sp(fvc::div(phi), h)
- fvm::laplacian(turbulence->alphaEff(), h)
- Q
==
fvc::div(phi/fvc::interpolate(rho), rho/psi, "div(U,p)")
- (rho/psi)*fvc::div(phi/fvc::interpolate(rho))
);

hEqn.relax();

hEqn.solve();

thermo.correct();
}

[yzf1215]

OK, I have solved this, the unit of Q is J/(M^3*s)

[zhuangli]

Originally Posted by yzf1215

OK, I have solved this, the unit of Q is J/(M^3*s)

Hi！boy，
How did you solve it, could you give me some advice?
anything would be appreciate!
lizhuang
bestWish！