|
[Sponsors] |
November 14, 2011, 04:33 |
Bananas!!! error in GGI and GAMG
|
#1 |
Member
Join Date: Nov 2009
Posts: 65
Rep Power: 16 |
Dear all,
I simulated 3-D vertical wind turbine rotors using pimpleDyMFoam with GGI on OF1.6-ext (latest Git ver.). Some cases worked fine but another case failed with FATAL error message like below at p iteration stage which uses GAMG solver. In this error message, what "bananas!!!" means ? Thanks in advance, waku ---- error ----- [waku@ensis11 pimpleDyMFoam_VAWT2_NACA0018]$ pimpleDyMFoam [snip: FOAM headers] Create time Create dynamic mesh for time = 0 Selecting dynamicFvMesh mixerGgiFvMesh Initializing the GGI interpolator between master/shadow patches: statorInterface/rotorInterface void mixerGgiFvMesh::addZonesAndModifiers() : Zones and modifiers already present. Skipping. Mixer mesh: origin: (0 0 1) axis : (0 0 1) rpm : -191 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type RASModel Selecting RAS turbulence model kOmegaSST kOmegaSSTCoeffs { alphaK1 0.85034; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.85616; gamma1 0.5532; gamma2 0.4403; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; c1 10; } Reading field rAU if present Starting time loop Courant Number mean: 0 max: 0.0365714 velocity magnitude: 5 Time = 5e-05 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 7.55501e-06, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 9.75065e-08, No Iterations 3 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.70441e-06, No Iterations 1 --> FOAM FATAL ERROR: Bananas!!! From function ggiGAMGInterfaceField::updateInterfaceMatrix in file matrices/lduMatrix/solvers/GAMG/interfaceFields/ggiGAMGInterfaceField/ggiGAMGInterfaceField.C at line 113. FOAM aborting --------------------------------- |
|
November 14, 2011, 05:40 |
|
#2 |
Senior Member
Pawel Sosnowski
Join Date: Mar 2009
Location: Munich, Germany
Posts: 105
Rep Power: 18 |
||
November 14, 2011, 06:46 |
|
#3 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 29 |
Pawel, I had the same though reading the thread title, but this error is unrelated in my opinion.
Edit: The offending piece of code is: Code:
110 // New treatment. HJ, 26/Jun/2011 111 if (zonePnf.size() != faceCells.size()) 112 { 113 FatalErrorIn("ggiGAMGInterfaceField::updateInterfaceMatrix") 114 << "Bananas!!!" 115 << abort(FatalError); 116 } Last edited by akidess; November 14, 2011 at 06:49. Reason: extended answer |
|
November 14, 2011, 09:21 |
|
#4 |
Member
Join Date: Nov 2009
Posts: 65
Rep Power: 16 |
Thanks Anton and Pawel
I know the banana method mentioned by Pawel. Anton's comment gives me useful suggestions. I'll check my mesh. Best regards waku Last edited by waku2005; November 15, 2011 at 01:41. |
|
November 15, 2011, 03:50 |
|
#5 | |
Member
Join Date: Nov 2009
Posts: 65
Rep Power: 16 |
Dear all,
I checked my mesh but no error was found in my data with checkMesh command, however, I found zonePf.size() had non zero value even faceCell.size() was zero during p calculation stage of the below if statement. Does the way of GGI come to include something new after jun/2011 update ? Any suggestions are highly welcomed. waku Quote:
|
||
November 17, 2011, 15:01 |
|
#6 |
New Member
Sinisa Majer
Join Date: Oct 2009
Location: Zagreb, CROATIA
Posts: 8
Rep Power: 16 |
Hi waku2005,
I've had a same problem several weeks ago, and I solved it by replacing GAMG solver for p with PCG. After few iterations I replaced it again with GAMG and after that everything worked just fine. I cannot explain the reason for this behavior, but it worked for me. Regards, S |
|
November 17, 2011, 17:41 |
|
#7 |
Member
Join Date: Nov 2009
Posts: 65
Rep Power: 16 |
Dear Sinisa
Thanks a lot ! I've just replace GAMG to PCG solver for a workaround. I'll do the same with your suggestion. Regards waku |
|
November 21, 2011, 22:22 |
|
#8 |
Member
Join Date: Nov 2009
Posts: 65
Rep Power: 16 |
Dear Sinisa
Sorry to my late response. I've just test PCG solver for several time steps from startTime, and then changed to the pressure solver to GAMG but failed to continue my run with the same "Bananas!" error. Regards waku |
|
November 22, 2011, 05:47 |
|
#9 |
New Member
Sinisa Majer
Join Date: Oct 2009
Location: Zagreb, CROATIA
Posts: 8
Rep Power: 16 |
Sorry,but I have no idea where is the problem. This worked for me and I didn't try to find out why.
Perhaps you should iterate more time steps with PCG (but this looks like a 'blind man's attempt'). S |
|
November 22, 2011, 06:26 |
|
#10 |
Member
Join Date: Nov 2009
Posts: 65
Rep Power: 16 |
Thanks a lot!
I'll continue to find where is the problem. waku at 24/11/2011 I've just found that the parallel run with GAMG solver works ! Last edited by waku2005; November 23, 2011 at 23:00. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
GGI in OpenFOAM-1.5-dev | philippose | OpenFOAM Running, Solving & CFD | 14 | November 13, 2011 14:55 |
Performance of GGI case in parallel | hannes | OpenFOAM Running, Solving & CFD | 26 | August 3, 2011 03:07 |
Increasing Iter. number of GAMG in pimpleDyMFoam with GGI | waku2005 | OpenFOAM | 0 | July 3, 2011 20:21 |
GAMG does not work using the regionCouple-BC | herbert | OpenFOAM Running, Solving & CFD | 1 | December 17, 2009 10:01 |