CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Face and point ordering during interpolation

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   November 28, 2011, 14:32
Default Face and point ordering during interpolation
  #1
Senior Member
 
Arne Stahlmann
Join Date: Nov 2009
Location: Hanover, Germany
Posts: 209
Rep Power: 8
Arnoldinho is on a distinguished road
Hi all,

I have two questions concerning the interpolation from faces to points and getting the corresponding point coordinates.
To move the boundary of a movingMesh case, I have calculated the displacement/motion at the face centers on the moving boundary patch (stored in volScalarField hbVol.boundaryField()[patchi]).

As I need the point displacement/motion at the points and not the faces, I used primitivePatchInterpolation faceToPointInterpolate on the face values. The point values are now used to modify the pointdisplacement/motionU file at the boundary patch within the solver, by doing e.g.

Quote:
pointVectorField& pointDisplacement
(
const_cast<pointVectorField&>(mesh.objectRegistry: :lookupObject<pointVectorField>("pointDisplacement "))
);

vectorField &pDisp = refCast<vectorField>(pointDisplacement.boundaryFie ld()[patchi]);

pDisp += newPointHeight; // newPointHeight from faceToPointInterpolate
My questions are:
1. How do I know or can make sure that faceToPointInterpolate is giving the same order of the points that is used in pointDisplacement file?

2. How can I get the x,y,z coordinates of these boundary patch points, and again make sure that the point order is right?
Something like mesh.boundaryMesh()[patchi].points()[pointI].x() does not give me the right points.

Arne
Arnoldinho is offline   Reply With Quote

Old   November 29, 2011, 05:54
Default
  #2
Senior Member
 
Arne Stahlmann
Join Date: Nov 2009
Location: Hanover, Germany
Posts: 209
Rep Power: 8
Arnoldinho is on a distinguished road
Ok, to answer my questions:

2. Point coordinates on the patch (e.g. for x) are found using mesh.boundaryMesh()[patchi].localPoints()[pointI].x().

1. I have checked the point order (sorry, not ordering!) by looping 1. over mesh.boundaryMesh()[patchi].localPoints() and 2. over the point list received from face to point interpolation. Result: Order is equal!

Greetings,
Arne
Arnoldinho is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Internal walls of zero thickness anger OpenFOAM Native Meshers: blockMesh 21 March 19, 2015 10:21
mesh airfoil NACA0012 anand_30 OpenFOAM Meshing & Mesh Conversion 12 December 12, 2011 05:16
Import netgen mesh to OpenFOAM hsieh Open Source Meshers: Gmsh, Netgen, CGNS, ... 32 September 13, 2011 05:50
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36
Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15


All times are GMT -4. The time now is 15:19.