#217 RenumberMesh cause fail to ggi pair

assigned
None
normal
major
always
none
foam-3.0
other
2014-05-07
2014-04-27
Pekka Kanninen
No

When the computation case includes ggi interface which face number is over than ~50 000 and renumberMesh is used. Error occur when the pressure solver is GAMG. Same model works fine if it build without running renumberMesh.

Other issue related to renumberMesh when the mixingPlane interface is used. The renumbering generate two coordinates systems types. From boundary file:

ŒŒŒŒŒŒŒ coordinateSystemŒŒŒŒŒŒŒŒ
ŒŒŒŒŒŒŒ {
ŒŒŒŒŒŒŒŒŒŒŒ typeŒŒŒŒŒŒŒŒŒŒŒ cylindrical;
ŒŒŒŒŒŒŒŒŒŒŒ nameŒŒŒŒŒŒŒŒŒŒŒ mixingCS;
ŒŒŒŒŒŒŒŒŒŒŒ typeŒŒŒŒŒŒŒŒŒŒŒ cylindrical;
ŒŒŒŒŒŒŒŒŒŒŒ originŒŒŒŒŒŒŒŒŒ (0 0 0);
ŒŒŒŒŒŒŒŒŒŒŒ e1ŒŒŒŒŒŒŒŒŒŒŒŒŒ (1 0 0);
ŒŒŒŒŒŒŒŒŒŒŒ e3ŒŒŒŒŒŒŒŒŒŒŒŒŒ (0 0 1);
ŒŒŒŒŒŒŒ }

it does not affected to foam solver but it cause paraview error when run "pvbatch --use-offscreen-rendering"
---- additional_information ----
Attached zip file includes modified tutorial. Error can be repeat by this model using AllRun.py script. Example case out put:

Starting time loop

Time = 0.001

smoothSolver:ΠSolving for Ux, Initial residual = 0.9999999944, Final residual = 0.09473563285, No Iterations 3
smoothSolver:ΠSolving for Uy, Initial residual = 0.9999999992, Final residual = 0.07470442025, No Iterations 3
smoothSolver:ΠSolving for Uz, Initial residual = 1, Final residual = 0.08002172049, No Iterations 3

--> FOAM FATAL ERROR:
Error in interface update: incorrect size of zone fields
Field size = 125000 Zone size = 71762

From function ggiGAMGInterfaceField::updateInterfaceMatrix
in file matrices/lduMatrix/solvers/GAMG/interfaceFields/ggiGAMGInterfaceField/ggiGAMGInterfaceField.C at line 114.

FOAM aborting

Discussion

  • Pekka Kanninen
    Pekka Kanninen
    2014-05-07

    Actually same error occur without renumbering in the larger model. Model, where is over 200 000 face on the ggi patch cause the same GAMG error. Same model works, for example, PCG pressure solver. In other words, renumberMesh is not the key factor on this issue.

    --> FOAM FATAL ERROR:
    Error in interface update: incorrect size of zone fields
    Field size = 2 Zone size = 1

    ŒŒŒ From function ggiGAMGInterfaceField::updateInterfaceMatrix
    ŒŒŒ in file matrices/lduMatrix/solvers/GAMG/interfaceFields/ggiGAMGInterfaceField/ggiGAMGInterfaceField.C at line 114.