[KangX1]

Dear Foamers,

I was using pimpleDyMFoam to simulate some dynamic mesh case.
My case can run in serial but not in parallel.

Time = 0.0001

Inverting RBF motion matrix
volume continuity errors : volume = 63.99998804, max error = 0, sum local = 0, global = 0
smoothSolver:  Solving for Ux, Initial residual = 0.9999999697, Final residual = 1.940639766e-06, No Iterations 1
DICPCG:  Solving for p, Initial residual = 1, Final residual = 0.009113619622, No Iterations 10
DICPCG:  Solving for p, Initial residual = 0.1448696186, Final residual = 0.001360322043, No Iterations 12
DICPCG:  Solving for p, Initial residual = 0.03288304611, Final residual = 0.0003125233115, No Iterations 29
time step continuity errors : sum local = 1.062426731e-09, global = -8.417019943e-11, cumulative = -8.417019943e-11
Moving mesh time step continuity errors : sum local = 1.062426731e-09, global = -8.417019943e-11, cumulative = -1.683403989e-10
DICPCG:  Solving for p, Initial residual = 0.01141814687, Final residual = 0.0001126015515, No Iterations 40
DICPCG:  Solving for p, Initial residual = 0.004746542347, Final residual = 4.588051877e-05, No Iterations 92
DICPCG:  Solving for p, Initial residual = 0.004053370021, Final residual = 0.0003471630065, No Iterations 2
time step continuity errors : sum local = 1.223790396e-09, global = 1.079705363e-12, cumulative = -1.672606935e-10
Moving mesh time step continuity errors : sum local = 1.223790396e-09, global = 1.079705363e-12, cumulative = -1.661809881e-10
smoothSolver:  Solving for Ux, Initial residual = 6.683720577e-05, Final residual = 4.632345119e-10, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.001679268542, Final residual = 1.645379057e-05, No Iterations 51
DICPCG:  Solving for p, Initial residual = 0.001234871242, Final residual = 1.192388694e-05, No Iterations 11
DICPCG:  Solving for p, Initial residual = 0.0003613903132, Final residual = 3.483720583e-06, No Iterations 59
time step continuity errors : sum local = 1.236435531e-11, global = -2.547561024e-12, cumulative = -1.687285492e-10
Moving mesh time step continuity errors : sum local = 1.236435531e-11, global = -2.547561024e-12, cumulative = -1.712761102e-10
DICPCG:  Solving for p, Initial residual = 0.000356150859, Final residual = 3.559849006e-06, No Iterations 62
DICPCG:  Solving for p, Initial residual = 0.0004029486792, Final residual = 3.729478687e-06, No Iterations 12
DICPCG:  Solving for p, Initial residual = 0.0001090817749, Final residual = 7.626294024e-06, No Iterations 2
time step continuity errors : sum local = 2.706125338e-11, global = 2.824445213e-12, cumulative = -1.68451665e-10
Moving mesh time step continuity errors : sum local = 2.706125338e-11, global = 2.824445213e-12, cumulative = -1.656272198e-10
ExecutionTime = 16.72 s  ClockTime = 17 s

Courant Number mean: 5.566444905e-07 max: 0.0001664760881 velocity magnitude: 0.01569961157
forces output:
    forces(pressure, viscous)((0 0 0) (0 0 0))
    moment(pressure, viscous)((0 0 0) (0 0 0))

forceCoeffs output:
    Cd = 0
    Cl = 0
    Cm = 0

Time = 0.0002

[0] 
[0] 
[0] --> FOAM FATAL ERROR: 
[0] face 3013 area does not match neighbour by 0.12569% -- possible face ordering problem.
patch: procBoundary0to1 my area: 1.24799e-06 neighbour area: 1.24642e-06 matching tolerance: 0.0001
Mesh face: 551216 vertices: 3((0.00971941 -0.0168322 -0.0854581) (0.00981709 -0.0153605 -0.0854373) (0.00850303 -0.0162224 -0.0843145))
Rerun with processor debug flag set for more information.
[0] 
[0]     From function processorPolyPatch::calcGeometry()
[0]     in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 217.
[0] 
FOAM parallel run exiting
[0] 
[1] --> FOAM FATAL ERROR: 
[1] face 3013 area does not match neighbour by 0.12569% -- possible face ordering problem.
patch: procBoundary1to0 my area: 1.24642e-06 neighbour area: 1.24799e-06 matching tolerance: 0.0001
Mesh face: 711111 vertices: 3((0.00971941 -0.0168322 -0.0854581) (0.00850303 -0.0162224 -0.0843145) (0.00981098 -0.0153605 -0.0854373))
Rerun with processor debug flag set for more information.
[1] 
[1]     From function processorPolyPatch::calcGeometry()
[1]     in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 217.
[1] 
FOAM parallel run exiting
[1] 
[3] 
[3] 
[3] --> FOAM FATAL ERROR: 
[3] face 128 area does not match neighbour by 0.0344233% -- possible face ordering problem.
patch: procBoundary3to2 my area: 2.55386e-07 neighbour area: 2.55474e-07 matching tolerance: 0.0001
Mesh face: 548575 vertices: 3((0.00972533 -0.00624834 -0.0735605) (0.00928942 -0.00670668 -0.0732511) (0.00881222 -0.00615683 -0.0732511))
Rerun with processor debug flag set for more information.
[3] 
[3]     From function processorPolyPatch::calcGeometry()
[3]     in file [2] 
[2] 
[2] --> FOAM FATAL ERROR: 
[2] face 128 area does not match neighbour by 0.0344233% -- possible face ordering problem.
patch: procBoundary2to3 my area: 2.55474e-07 neighbour area: 2.55386e-07 matching tolerance: 0.0001
Mesh face: 708031 vertices: 3((0.00972533 -0.00624834 -0.0735605) (0.00881186 -0.00615683 -0.0732511) (0.00928942 -0.00670668 -0.0732511))
Rerun with processor debug flag set for more information.
[2] meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 217.
[3] 
FOAM parallel run exiting
[3] 

[2]     From function processorPolyPatch::calcGeometry()
[2]     in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 217.
[2] 
FOAM parallel run exiting

I have checked many threads in the forums. They both mentioned the cyclic boundary condition problem. But in my case, I don't have this boundary condition. Here's my boundary dict.

FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

4
(
    wing
    {
        type            wall;
        nFaces          26996;
        startFace       2442294;
    }
    front
    {
        type            wall;
        nFaces          10638;
        startFace       2469290;
    }
    back
    {
        type            wall;
        nFaces          10613;
        startFace       2479928;
    }
    side
    {
        type            empty;
        nFaces          42549;
        startFace       2490541;
    }
)

// ************************************************************************* //

If the problem is caused by the mesh, what can I do to check it? I have used the checkMesh. It seems all fine...

BTW, I am using Foam-extend 3.1

Thanks,
Lida

[KangX1]

can anyone help??