Accessing X and Y coordinates of Face centre

R_21 · May 8, 2015, 16:00

Hi Foamers,
I have been trying to access the x- and y- co-ordinate of a mesh cell face closest to a particular set of points. I have written the following set of code referring to various posts on this Forum.
"
label actMeshID = mesh.findCell(point(4.0,0.75,1.0));
const vectorField& ActPoint = mesh.Cf().boundaryField()[actMeshID];
const scalarField xa = ActPoint.component(vector::X);
const scalarField ya = ActPoint.component(vector::Y);

"

I am getting this sort of error while running:
"--> FOAM FATAL ERROR:
hanging pointer at index 669 (size 6), cannot dereference

From function PtrList:

perator[] const
in file /home/rafath/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/PtrListI.H at line 134.
"
Can Someone please help Please!!!
Thanks in advance

fumiya · May 9, 2015, 02:26

This code is an example:

Code:

vector probePoint(4.0,0.75,1.0);
label probeCell = mesh.findCell(probePoint);
const labelList& f = mesh.cells()[probeCell]; // face label list of probeCell

forAll(f, facei)
{
    label nFace = f[facei];
    Info<< mesh.Cf()[nFace] << endl;
    // mesh.Cf()[nFace][0]: x
    // mesh.Cf()[nFace][1]: y
}

I hope this will be of help
Fumiya

R_21 · May 9, 2015, 05:15

Thanks a lot Fumiya,
That really helped my case.

May 8, 2015, 16:00	Accessing X and Y coordinates of Face centre	#1
R_21 New Member rafath Join Date: Jun 2014 Location: mumbai Posts: 24 Rep Power: 11	Hi Foamers, I have been trying to access the x- and y- co-ordinate of a mesh cell face closest to a particular set of points. I have written the following set of code referring to various posts on this Forum. " label actMeshID = mesh.findCell(point(4.0,0.75,1.0)); const vectorField& ActPoint = mesh.Cf().boundaryField()[actMeshID]; const scalarField xa = ActPoint.component(vector::X); const scalarField ya = ActPoint.component(vector::Y); " I am getting this sort of error while running: "--> FOAM FATAL ERROR: hanging pointer at index 669 (size 6), cannot dereference From function PtrList:perator[] const in file /home/rafath/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/PtrListI.H at line 134. " Can Someone please help Please!!! Thanks in advance

May 9, 2015, 02:26		#2
fumiya Senior Member Fumiya Nozaki Join Date: Jun 2010 Location: Yokohama, Japan Posts: 266 Blog Entries: 1 Rep Power: 18	This code is an example: Code: vector probePoint(4.0,0.75,1.0); label probeCell = mesh.findCell(probePoint); const labelList& f = mesh.cells()[probeCell]; // face label list of probeCell forAll(f, facei) { label nFace = f[facei]; Info<< mesh.Cf()[nFace] << endl; // mesh.Cf()[nFace][0]: x // mesh.Cf()[nFace][1]: y } I hope this will be of help Fumiya Luttappy and micheall like this. __________________ [Personal] Blog: https://caefn.com/ Twitter: https://twitter.com/FumiyaNozaki SlideShare: http://www.slideshare.net/fumiyanozaki96 Last edited by fumiya; May 9, 2015 at 06:26.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[Netgen] Import netgen mesh to OpenFOAM	hsieh	OpenFOAM Meshing & Mesh Conversion	32	September 13, 2011 05:50
[blockMesh] BlockMeshmergePatchPairs	hjasak	OpenFOAM Meshing & Mesh Conversion	11	August 15, 2008 07:36
fluent add additional zones for the mesh file	SSL	FLUENT	2	January 26, 2008 11:55
[blockMesh] Axisymmetrical mesh	Rasmus Gjesing (Gjesing)	OpenFOAM Meshing & Mesh Conversion	10	April 2, 2007 14:00
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues	michele	OpenFOAM Meshing & Mesh Conversion	2	July 15, 2005 04:15

May 9, 2015, 05:15		#3
R_21 New Member rafath Join Date: Jun 2014 Location: mumbai Posts: 24 Rep Power: 11	Thanks a lot Fumiya, That really helped my case.