cyclicGgi and AMG solvers

cliffoi · July 25, 2014, 09:41

I am trying to use cyclicGgi boundaries in combination with AMG solvers in foam-extend 3.1 and receive the following error:

Code:

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

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

FOAM aborting

I've tried this with GAMG and the AMG solvers in lduSolvers and the result is the same. Is there an incompatibility between this combination of BCs and solver.

Thanks advance
Ivor

hjasak · July 27, 2014, 03:56

No incompatibility: this must work and we have tutorials + test loop items using it.

Did you set up the zones correctly? If so, can you share the case - I'd like to have a look.

Hrv

cliffoi · July 30, 2014, 09:24

Hi Hrv,
I would have thought it should work as well. The mesh works fine using AMI on OF-2.2.1, I'm also fairly certain the zone definitions are correct since there are no issues using the standard conjugate gradient solver with incomplete Cholesky preconditioning. The only thing I think it could be is the mesh quality. I'm trying to use cfMesh for the meshing and am still working through some issues capturing the feature edges and boundary layers correctly.
Unfortunately the geometry is proprietary and it's not mine to give away. I will try to reproduce the error using modified geometry and then I can perhaps send it to you.

Best Regards
Ivor

July 25, 2014, 09:41	cyclicGgi and AMG solvers	#1
cliffoi Member Ivor Clifford Join Date: Mar 2009 Location: Switzerland Posts: 94 Rep Power: 17	I am trying to use cyclicGgi boundaries in combination with AMG solvers in foam-extend 3.1 and receive the following error: Code: --> FOAM FATAL ERROR: Error in interface update: incorrect size of zone fields Field size = 656 Zone size = 451 From function ggiGAMGInterfaceField::updateInterfaceMatrix in file matrices/lduMatrix/solvers/GAMG/interfaceFields/ggiGAMGInterfaceField/ggiGAMGInterfaceField.C at line 114. FOAM aborting I've tried this with GAMG and the AMG solvers in lduSolvers and the result is the same. Is there an incompatibility between this combination of BCs and solver. Thanks advance Ivor

July 27, 2014, 03:56		#2
hjasak Senior Member Hrvoje Jasak Join Date: Mar 2009 Location: London, England Posts: 1,907 Rep Power: 33	No incompatibility: this must work and we have tutorials + test loop items using it. Did you set up the zones correctly? If so, can you share the case - I'd like to have a look. Hrv __________________ Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk

July 30, 2014, 09:24		#3
cliffoi Member Ivor Clifford Join Date: Mar 2009 Location: Switzerland Posts: 94 Rep Power: 17	Hi Hrv, I would have thought it should work as well. The mesh works fine using AMI on OF-2.2.1, I'm also fairly certain the zone definitions are correct since there are no issues using the standard conjugate gradient solver with incomplete Cholesky preconditioning. The only thing I think it could be is the mesh quality. I'm trying to use cfMesh for the meshing and am still working through some issues capturing the feature edges and boundary layers correctly. Unfortunately the geometry is proprietary and it's not mine to give away. I will try to reproduce the error using modified geometry and then I can perhaps send it to you. Best Regards Ivor