[Oden]

Hi All, I am trying to add a temperature dependant cooling source fvModel in a coldRoom simulation case. Below is my code:

coolingModel
{
type coded;
active yes;

codedCoeffs
{
name coolingSource;
selectionMode all; // Apply to entire mesh or specify a cellZone
// cellZone coolingZone; // Uncomment if applying to a cell zone
fields (e); // Energy field to which the source is applied

code
#{
#include "fvCFD.H"

// Access the mesh and time objects
const fvMesh& mesh = this->mesh();
const Time& runTime = mesh.time();

// Source term addition function for scalar equations
void addSup(fvMatrix<scalar>& eqn)
{
const volScalarField& T = mesh.lookupObject<volScalarField>("T");
const scalarField& V = mesh.V();
scalarField& source = eqn.source();

scalar Tmax = 280.15; // Upper temperature limit (5°C)
scalar Q = -100; // Cooling rate in W/m³

forAll(source, celli)
{
if (T[celli] > Tmax)
{
source[celli] += Q * V[celli];
}
}

Info << "Applied cooling source term at time " << runTime.timeName() << endl;
}

// Source term addition function for vector equations (if needed)
void addSup(fvMatrix<vector>& eqn)
{
// Do nothing for vector equations
}

// Correct function (called after the solution)
void correct()
{
// No corrections needed
}
#};
}
}

As much as the fvModel is executed during simulation, I don't see its impact on the temperature field as the temperature in my cold room keeps rising beyond the set threshold (Tmax) until it settles at ambient temperature.

Any ideas on how I can better this approach?

[dasith0001]

Hi,

Just double check if you can implement a negative value for the heat source?

you can certainly decrease of power say 100 to 0 W/m^3 over specified period of time.

hope this helps