[LESlie]

Hi Foamers,

I am working on needle transient motion of injectors and I have implemented my own class with layer removal based on movingTopoFvMesh. The injector mesh is converted from a fluent mesh. What I did is to find the patch, say NEEDLE with

label patchID = boundaryMesh().findPatchID("NEEDLE");

Info << "patchID = " << patchID << endl;

const polyPatch& cPatch = boundaryMesh()[patchID];

then put the faces of this polyPatch as extrusion faces for the mesh modifier, and also mark the points on these faces for the motion, so that only these points will be moved.

The code is working in serial mode (attached pictures), however failed in parallel.

Serial log file
Create mesh for time = 0.0002

Selecting dynamicFvMesh needleTopoFvMesh
Initial time:0.0002 Initial curMotionVel_0 0 1)
Time = 0.0002

Adding zones and modifiers to the mesh

patchID = 1
cPatch =
type wall;
nFaces 20732;
startFace 3529376;

number of face in zone1 20732

Adding mesh zones.
Adding 1 mesh modifiers
patchID = 1
number of points marked 20986

.... reading data....

time:0.0002000017749 curMotionVel_0 0 1)
No topology change
Executing mesh motion
Execution time for mesh.update() = 1.67 s

..... flow calculation .....

Total Volume change = -3.598473946e-15
ExecutionTime = 31.52 s


Parallel log file
Selecting dynamicFvMesh needleTopoFvMesh
[0] Initial time:0.0002 Initial curMotionVel_0 0 1)
Time = 0.0002
Adding zones and modifiers to the mesh
[1] Initial time:0.0002 Initial curMotionVel_0 0 1)
[2] Initial time:0.0002 Initial curMotionVel_0 0 1)
[3] Initial time:0.0002 Initial curMotionVel_0 0 1)

patchID = 1
cPatch =
type wall;
nFaces 0;
startFace 879157;

number of face in zone1 0

Adding mesh zones.
Adding 1 mesh modifiers
patchID = 1
number of points marked 0

.... reading data....

time:0.0002000017749 curMotionVel_0 0 1)
No topology change
Executing mesh motion
Execution time for mesh.update() = 1.67 s

..... flow calculation .....

Total Volume change = 0
ExecutionTime = 14.72 s


Apparently when running in parallel (4 procs), polyPatch& cPatch = boundaryMesh()[patchID] failed to create the polyPatch. I am running on 2.2.x. I would really appreciate some help on this!

Have a nice day!