[Asghari_M]

Dear Foamers;

In my problem, I have a few zones and cell zones.
I need to select a cell located on center of a zone after problem solution.
How can I access that cell ?
Then, for example I want to access some values (quantities) such as temperature for that cell. How can I do these operations;
Thx in advance for any answer and help.

M.Asghari

[Asghari_M]

Dear OF users,
For example, in Chtmultiregionfoam, I'd like to write T(nx/2, ny/2) for heater Zone in a .dat text file (nx, ny are x and y mesh numbers for heater zone,respectively ).
Can anybody help me ?

[fumiya]

You can get the label list of cells in the particular cellZone.
I hope my slide(p.7) will be of help:
http://www.slideshare.net/fumiyanoza...ogramming-tips

[Asghari_M]

Hi Fumiya.
I don't have any username ,pass for slideshare.net.
Can you send your slide my email address, please?
My address is: badzohreh.h57@gmail.com
Best Regards;

[fumiya]

Hi,

This is an example:

label cellZoneID = mesh.cellZones().findZoneID("your_cellZone_name");
const labelList& cellLabels = mesh.cellZones()[cellZoneID];

// Loop over the cells that are in the specified cellZone and output temperature
forAll(cellLabels, i)
{
    label cellI = cellLabels[i];       // label of a cell that is in your cellZone
    Info<< mesh.C()[cellI] << " " << T[cellI] << endl; // Write cell center coordinate and temperature
}

I hope this will be of some help.
Fumiya

[Asghari_M]

Dear Fumiya;
Thanks for your suggested code.
However, I don't know how to use ur suggested code.
1-Where should I use it in?
2-Really, I have geometrical coordinate of a point in a specified zone and I want to access its cell number and then I want to obtain temperature of that cell Id (lable),
What's the relation between your suggested sub-code and my purpose?
3-What 's the meaning of Info command and its supplements in following your suggested sub-code?

    Info<< mesh.C()[cellI] << " " << T[cellI] << endl; // T: your variable name
}

Thanks in advance for your help.

[fumiya]

Hi,
My code in #5 outputs temperature in each cell that is in your cellZone.

1. What is your solver?
2. If you know the geometrical coordinate of a point in a specified zone,
you can use findCell(const point & location) command to find the cell containing the point.

vector probePoint(x,y,z);
label probeCell = mesh.findCell(probePoint);

You can access the temperature at the cell center by T[probeCell].

3. Info command output to the terminal. If you want to output to your text file,
you can use OFstream. You might want to consult the following thread:
http://www.cfd-online.com/Forums/ope...scii-file.html

Hope this helps,
Fumiya

[Asghari_M]

Dear Fumiya;
Thanks for ur kindly helps.
My solver is OpenFoam 2.3.1.
I am not familiar with these commands & symbols and where they are written.

vector probePoint(x,y,z); label probeCell = mesh.findCell(probePoint);

How can I download OFstream? Is it a part of OF2.3.1?
I'm reading ur suggested thread. It seems to be interesting & helpful.

Best Regards.

M.Asghari