[satkinson]

Hi All,

I have been trying to understand the EEqn.h using http://cfd.direct/openfoam/energy-equation/ and the file itself. From the file;

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)
==
rho*(U&g)
+ radiation->Sh(thermo)
+ fvOptions(rho, he)
);

I fully understand the LHS and the process that is happening there from the website. I am a little confused why I, the Identity tensor gets swapped in it becomes U in equation 12.

Following this I want to know what rho*(U&g) is in physical terms. I understand that rho*g gives joules, which makes sense, but what is rho*U kg/m2.s.

Thanks,

satkinson