[cute]

Dear members,

I am running modified version of cavity test case with dimension of 1 mm^3 and grid size of 100 x 100 x 100. I want to access all values at z = 95. I can write C++ code but don't know how to do it in openFoam. Could you please help? I searched the forum but did not find related information.

    for (i=0; i<200; i++)
   {
       for (j=0; j<200; j++)
           cout<< T(i,j,95)<<endl;
   }

Thanks.

[Tobi]

Hi,

I think you want to have access to the cell values of all cells within (x, y, 95). There are different options. The best one would to make a cell zone of all these cells first (you will get a cell list that contains the ID's of the cells, this will save time).

Otherwise you can search for these cells for each loop. What you have to do is, get the cells (mesh.C()), search for the components and done. But you should more or less define somehow a range. Lets say:

//- Check the doxygen for the correct Types
const cellField& pF = mesh.C();

//- Search the cells with 94.5 < z < 95.5
List<label> cellIDs;

forAll(pF, ID)
{
       //- z coordinate
       const scalar zCoor = pF[ID].z();

       if ( zCoor < 95.5 && zCoor > 94.5)
       {
             cellIDs.append(ID);
       }
}


//- Print all values from Temperature field
const scalarField& T_ = T.internalField();

forAll(cellIDs, ID)
{
     Info<<  T_[ID].value() << endl;
}

Something like that should do the stuff. It is easier to have a pre-defined cellZone. Then it is much shorter.

Hope it will help you.

[cute]

Thanks for the reply. I will implement as soon as I will be back. But one more thing, If I also need to include boundary value at z=95 what I need to change?

[Tobi]

Then you have to do the same stuff with the boundary conditions. Check for the boundary and search for the face-center and get the ID.

[cute]

Thanks, I will try this afternoon.

[cute]

I am getting the following error:

error: ‘cellField’ does not name a type
     const cellField& pF = mesh().C();

whats wrong here, could not understand?


(2) For creating a list at z=95, why do we need

zCoor < 95.5 && zCoor > 94.5

can we use ZCoor == 95?

[Tobi]

As I told before, you have to check the types and functions yourself. Therefore you can use doxygen (I will not present the solution).

You can take 95 but if the cell centers are not exactly at that position this will not work. You should more like check if the cell goes through this fixed z-line. So it is up to you what you finally do.

For example: http://cpp.openfoam.org/v4/a00886.html

You will see that the first line has to be:

const volVectorField foobar = mesh.C()

As you can see further, the boundaries are also included (should be because it is a volVectorField).

[cute]

Thanks, it is volVectorField.

Now I am getting different error:

error: request for member ‘value’ in ‘(& X_)->Foam::Field<double>::<anonymous>.Foam::List<double>::<anonymous>.Foam::UList<T>::operator[]<double>(ID)’, which is of non-class type ‘const double’
            if (X_[ID].value() > 0.0)

what is the reason for this?

[cute]

Got it, after deleting
" .value() "

[cute]

Thaks Tobi for your help. It is running now.

[cute]

This simulation gets stuck for parallel processing .

If I remove the above code, simulation runs and creates data files. Do i need to add something in above code for parallel processing? I added the above code before while-loop (solution-loop).

[cute]

After doing some work on code, I realized that forAll LOOP does not work in parallel.

What should I add or modify in above code for parallel processing?
Need help plz.

[cute]

I am not able to find the answer yet. Some of them in the forum also reported this issue but no solution.

forAll loop is not running in parallel and my simulation gets stuck.

Any hint would be appreciated.