Menu
▾
▴
#217 RenumberMesh cause fail to ggi pair
assigned
None
normal
major
always
none
foam-3.0
other
2014-05-07
2014-04-27
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
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.