[shackman287]

I'm performing mesh morphing with the tetFEM solver in foam-extend-4.0. As expected it truly is far more robust than FVM (cell-centered) laplacian solvers. Now I want to try some more complex motions and define explicitly the pointMotion (now motionU in tetFem solver), but I can't seem to access the mesh.boundaryMesh[patchI].points.blah(), like I normally do. In fact even within the deformedTetFemGeom utility when I try to report the points, it just lists the FVM patch points, not the added FEM (mid) points.

I see some functions within tetPolyMesh libs, but really have no idea how I can access these functions (calcLocalPoints?) at the solver level to get a list of the FEM points on a given patch, so I can use to prescribe the motionU for a given patch.

HJasak Help! Your my only hope.... unless others also know

::FoamOn

[shackman287]

So I figured it out tinkering around in the tetFiniteElement src.

Here are my solution steps for anyone trying to access the FEM points for a given patch.

1. Hijacked the deformedTetFemGeom.C file to basically create a "tetPolyMesh tetMesh(mesh); " object and call member functions from that object " Info<< tetMesh.getPoints() <<endl; "

2. Created "pointField tetPolyMesh::getPoints()" member function in tetPolyMesh.C/.H to call a "boundaryTriPoints" member function.

3. Created "boundaryTriPoints" function in tetPolyBoundaryMesh.C/.H
to provide a pointField using the 'operator[](patchid).localPoints'

Hope this helps as an easy example of bringing "lower src code" functionality up to the executable level (deformedTetFemGeom.C).


::FoamOn