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

Face and point ordering during interpolation

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 2 Post By Arnoldinho
  • 1 Post By Arnoldinho

Reply
 
LinkBack Thread Tools Search this Thread 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: 18
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
blue8803 and scleakey like this.
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: 18
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
blue8803 likes this.
Arnoldinho is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[Other] mesh airfoil NACA0012 anand_30 OpenFOAM Meshing & Mesh Conversion 13 March 7, 2022 18:22
[blockMesh] Internal walls of zero thickness anger OpenFOAM Meshing & Mesh Conversion 23 February 6, 2020 19:25
[Netgen] Import netgen mesh to OpenFOAM hsieh OpenFOAM Meshing & Mesh Conversion 32 September 13, 2011 06:50
[blockMesh] BlockMeshmergePatchPairs hjasak OpenFOAM Meshing & Mesh Conversion 11 August 15, 2008 08:36
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Meshing & Mesh Conversion 2 July 15, 2005 05:15


All times are GMT -4. The time now is 04:34.