[rapierrz]

Hi friends,

I want to add a new face-flux field to openfoam solver.

I impliment below code but it has error.

surfaceScalarField phiEnergy
(
IOobject
(
"phiEnergy",
runTime.timeName(),
mesh
),
rho*Cp*(mesh.Sf() & fvc::interpolate(U))
);

fvScalarMatrix TEqn
(
fvm::ddt(rho,T)
+ fvm::div(phiEnergy,T)
- fvm::laplacian(K,T)
);

please Help me

[ghas]

Hi,

It didn't work because you have tried to multiply a volScalarField (rho*Cp) by a surfaceScalarField (mesh.Sf() & fvc::interpolate(U)).

To resolve the problem, you can define the face-flux as:

surfaceScalarField phiEnergy
(
IOobject
(
"phiEnergy",
runTime.timeName(),
mesh
),
(mesh.Sf() & fvc::interpolate(rho*Cp*U))
);

Best regards,

Ghassan

[rapierrz]

Thank you Ghassan,

it's work correctly.