[robingilbert]

Hi i want to add a heat source inside a small region in a fluid domain. i added a heat source term in TEqn.H and it works fine for the whole region. But i want to impose a condition that the heat source gets activated inside a small control volume. i have been reading up a lot but cant figure out how to go about doing it. if atleast someone can tell me where to start, I will figure it out. how do i give an if-else condition to check if it is solving the domain i need to add the source term.

i tried the following in the createFields.H

volScalarField x = mesh.C().component(vector::X);
volScalarField y = mesh.C().component(vector::Y);
volScalarField z = mesh.C().component(vector::Z);

forAll(x,celli)
{
    if(x[celli] >= 1 && y[celli] >= 1 && z[celli] >= 2.2 && x[celli] <= 1.5 && y[celli] <= 1.5 && z[celli] <= 2.5)
    {
       dimensionedScalar Q
    (
    "Q",
    dimensionSet(1,-1,-3,0,0,0,0),
    scalar(5000)
    );
    }else
    {
       dimensionedScalar Q
    (
    "Q",
    dimensionSet(1,-1,-3,0,0,0,0),
    scalar(1000)
    );
    }
}

unfortunately the Q value is not set by this method.

Thanks.