[Kummi]

Hello Foamers, Hope everyone doing good. Here in the attachment, the image contains nodal points at boundary and internal cells for 1D problem.
POINT 0 - Boundary node
POINTS 1,2,3 - Internal nodes

In order to extract variables at the boundary node (@ POINT 0), the below code was implemented, compiled good and end with reasonable results.

forAll (regionMesh().boundaryMesh(), patchI)
{
forAll(T_.boundaryField()[patchI], faceI)
{
scalar tK_ = K_.boundaryField()[patchI][faceI]; //VARIABLE K_ extracted from boundary POINT 0
scalar tT_ = T_.boundaryField()[patchI][faceI]; //VARIABLE T_ extracted from boundary POINT 0
QB_[cellI] = tK_* (tT_ - Tliq) / DX.value(); //Heat evaluation based on K_ & T_
}
}

However, I would like to have some clue as how to extract the variables at the internal nodes (POINTS 1,2,3). Correct me if I'm wrong anywhere please.
Kindly someone share some ideas.
Thank you

[Jagan G]

I have similar issues with 3D model.
Cell-centred calculation for internal cells is optable in my case.

[random_ran]

Originally Posted by Kummi

Hello Foamers, Hope everyone doing good. Here in the attachment, the image contains nodal points at boundary and internal cells for 1D problem.
POINT 0 - Boundary node
POINTS 1,2,3 - Internal nodes

In order to extract variables at the boundary node (@ POINT 0), the below code was implemented, compiled good and end with reasonable results.However, I would like to have some clue as how to extract the variables at the internal nodes (POINTS 1,2,3). Correct me if I'm wrong anywhere please.
Kindly someone share some ideas.
Thank you

Not sure this can help you or not.

You can access the `k` value by referencing the cellId`

`scalar k = turbulence->k()()[cellId];`

Then interpolated it to the patch face. After that, interpolated it to
the node of your interest. I never did this before, but I speculate
there's mature functions to let you do this within native OpenFOAM.

[Kummi]

Thank you for your response, Jaganathan and Ran.

Accessing the cellId may be required for specific purpose (for example, to imply the source term). Near the boundary, the cell ids are accessed through the boundary condition. So, I have created a patch based on boundaryMesh and obtained the values at the boundary.

For obtaining values at internal cells, I tried accessing neighbor cells, as in code below.

const Foam::List<int> neigbhorcellI = regionMesh().cellCells(cellI);//gives the neighbours of cell cellI
forAll(neigbhorcellI, I) //Loop inside internal cells
{
QB_[cellI] = K_[I]* ((T_[I] - Tliq) /DX.value());
break ;
}

But the problem here is, the code doesn't loop over all the cells and as a result, the results are not good as expected.
Kindly share your thoughts.

Thank you