[MartinB]

Hi Jeff,

add this snippet to an appropriate place in your code, for example in createFields.H:

    label TRefCell = 0;
    scalar TRefValue = 0.0;
    setRefCell(T, mesh.solutionDict().subDict("PISO"), TRefCell, TRefValue);

And in your fvSolution:

PISO
{
    . . .
    //TRefCell	0; // <--- definition via cell label
    TRefPoint       (0 0 0); // <--- definition via coordinates
    TRefValue      1000;
    . . .
}

Martin

[ozzythewise]

Update to this old thread:

Upon playing with this it didn't do it quite as i wanted. I'm not sure what it was doing but it wasn't keeping my temperature at 1 value, it was changing in time. So what I did is made my TEqn.H file read as follows:

/**************************
solve
(
    fvm::ddt(T)
    + fvm::div(phi, T)
    - fvm::laplacian(DT, T)
);
**************************/

{
    fvScalarMatrix TEqn
    (
        fvm::ddt(T)
        + fvm::div(phi, T)
        - fvm::laplacian(DT,T)
    );

    point TRefPoint(2, 0, 0);   // metres
    label TRefCell = mesh.findCell(TRefPoint);

    Info<< "TRefCell = " << TRefCell << nl
        << "TRefValue = " << TRefValue << nl
        << endl;

    // create a labelList containing only the label TRefCell
    labelList fixedCells(1, TRefCell);

    // create a field having as many values as the labelList,
    // specifically the values we want to set
    scalarField fixedValues(fixedCells.size(), TRefValue);

    // set the values.  This eliminates the specified cells from the matrix
    TEqn.setValues(fixedCells, fixedValues);

    TEqn.solve();
}

Hopefully this works for everyone else!
Jeff

[samiam1000]

Hi there,

I am writing since I am facing the problem you had and discussed about in this thread.

I am running a simulation and - in a certain area of my domain - I want to impose a fixed temperature. Let's say that I want to impose a boundary condition in a volume and not on a surface.

The solver I am using is buoyantSimpleFoam and I do not know which is the file that I need to edit.

Thanks a lot,

Samuele