Temperature Dependent Thermal Properties in LaplacianFoam

'm solving the heat equation using the finite volume method in OpenFOAM and the way it is set up right now is to average the thermal properties of the temperature field. What I want to achieve is to implement the thermal properties that correspond to the temperature of each cell/mesh. Is there a way to do that through the solver itself or in the transportProperties file.

Attached is the code I'm using. I just replaced the source term with the volumetric power that is calculated by another code.

Code:

Info<< "\nCalculating temperature distribution\n" << endl;

while (simple.loop(runTime))
{
Info<< "Time = " << runTime.timeName() << nl << endl;

while (simple.correctNonOrthogonal())
{
solve
(
fvm::ddt(T) - fvm::laplacian(DT, T) == volPower/(RHO * CP)
);

}

#include "write.H"

Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< "  ClockTime = " << runTime.elapsedClockTime() << " s"
<< nl << endl;
}

Info<< "End\n" << endl;

return 0;

This is the script that reads the thermal properties from the transportProperties file.

Code:

Info<< "Reading field T\n" << endl;

volScalarField T
(
IOobject
(
"T",
runTime.timeName(),
mesh,
IOobject::AUTO_WRITE
),
mesh
);

volScalarField volPower
(
IOobject
(
"volPower",
runTime.timeName(),
mesh,
IOobject::AUTO_WRITE
),
mesh
);

IOdictionary transportProperties
(
IOobject
(
"transportProperties",
runTime.constant(),
mesh,
IOobject::NO_WRITE
)
);

Info<< "Reading diffusivity DT\n" << endl;

dimensionedScalar DT
(
transportProperties.lookup("DT")
);

Info<< "Reading density RHO\n" << endl;

dimensionedScalar RHO
(
transportProperties.lookup("RHO")
);

Info<< "Reading thermal conductivity CP\n" << endl;

dimensionedScalar CP
(
transportProperties.lookup("CP")

This is how the transportProperties file is structured

Code:

DT              DT  [0 2 -1 0 0 0 0] 1.7620248e-5; // Diffusion Coefficent

RHO             RHO [1 -3 0 0 0 0 0] 1.52082743e3;  // Density

CP              CP  [0 2 -2 -1 0 0 0] 707.36193;  // Heat Capacity

I want DT and CP to be dependant on the temperature of each cell

OpenFOAM 05-01-22, 10:52 a.m. Akram_
To average the thermal properties of the temperature field and to achieve is to implement the thermal properties that correspond to the temperature of each cell/mesh driving directions.
14-01-22, 10:18 a.m. marfavor

I should express that I like this article so a ton and moreover using it in regular. I believe you would continue sharing better than average articles here. run 3
21-02-22, 1:42 p.m. run3donlineaz

21-03-22, 7:31 p.m. geebranz