[samiam1000]

Dear All,

what I am trying to do has been discussed, yet.

Anyway, I am finding it a bit difficult to find a solution, since there are many threads about this question and I can not cope with them all.

So, I would like to ask you, for help.

What I wanna do is to edit the buoyantSimpleFoam solver, in order to fix the value of the temperature in a sub-volume of my domain.

I have a mesh whit a volume patch which contains the cell I want to fix the T (or h) value.

The problem I have to face first, are the following:
a. in the solver we do not have T, but h. This is not a problem, since for a PIG it stands h = C_p*T, with C_p = const.

Now, the question is: what are the next step?

Following this tutorial I created my own solver buoyantSimpleFoam_Mod.

I edited the file `files' in the Make directory and then I started editing the hEqn.H and the createFields files.

Following this thread, I figured out that my hEqn.H file should looks like this:

{
    fvScalarMatrix hEqn
    (
        fvm::div(phi, h)
      - fvm::Sp(fvc::div(phi), h)
      - fvm::laplacian(turbulence->alphaEff(), h)
     ==
        fvc::div(phi/fvc::interpolate(rho)*fvc::interpolate(p))
      - p*fvc::div(phi/fvc::interpolate(rho))
    );

    hEqn.relax();
    hEqn.setValues(hCells, hVals);
    hEqn.solve();

    thermo.correct();
}

where hCells is a vector of cells and hVals is a vector containing the values that I want to impose at each cell of the hCells vector.

One question that comes to my mind is: where should I define these 2 vectors? And how?
Also, is it possible to insert the label of the volume patch, instead of the label of each cell?

Fuerthermore, I guess that I have to change the file fvSolution in the case/system directory and I think that it is enough to add the following line:

 
   
hCells	    0;
sVals      1000;

Is that correct? And how can I say to OF that hCells is the vector that contains all the cells in the patch I am interested in?

Could anyone help, please?

Thanks a lot,
Samuele