|
[Sponsors] |
June 4, 2018, 13:46 |
GAMG crash
|
#1 |
Member
ff
Join Date: Feb 2010
Posts: 81
Rep Power: 16 |
Hi,
I got one case that GAMG solver is crashed. I only attached a piece of information here, because they are just repeating. [325] #1 Foam::sigFpe::sigHandler(int)[38] #1 addr2line failed [297] #2 ?[182] #1 addr2line failed [359] #1 Foam::sigFpe::sigHandler(int)[323] #3 Foam::PBiCGStab::solve(Foam::Field<float>&, Foam::Field<float> const&, unsigned char) constFoam::sigFpe::sigHandler(int) addr2line failed [87] #1 Foam::sigFpe::sigHandler(int)[139] #2 ?[109] #3 Foam::PBiCGStab::solve(Foam::Field<float>&, Foam::Field<float> const&, unsigned char) const addr2line failed [138] #1 Foam::sigFpe::sigHandler(int)Foam::sigFpe::sigHand ler(int)[261] #1 addr2line failed [68] #3 Foam::PBiCGStab::solve(Foam::Field<float>&, Foam::Field<float> const&, unsigned char) const addr2line failed addr2line failed addr2line failed [188] #1 Foam::sigFpe::sigHandler(int)[99] #2 ? addr2line failed [357] #1 Foam::sigFpe::sigHandler(int) addr2line failed addr2line failed Foam::sigFpe::sigHandler(int)[236] #1 Foam::sigFpe::sigHandler(int) addr2line failed addr2line failed [212] #1 Foam::sigFpe::sigHandler(int)[312] #1 Foam::sigFpe::sigHandler(int) addr2line failed [42] #1 Foam::sigFpe::sigHandler(int) addr2line failed addr2line failed [93] #1 Foam::sigFpe::sigHandler(int) addr2line failed [224] #2 ? addr2line failed [328] #1 Foam::sigFpe::sigHandler(int)[102] #2 addr2line failed addr2line failed [238] #1 Foam::sigFpe::sigHandler(int) addr2line failed [132] #2 ?[58] #1 Foam::sigFpe::sigHandler(int) addr2line failed addr2line failed [259] #2 ??[77] #1 Foam::sigFpe::sigHandler(int) addr2line failed [230] #2 ? addr2line failed [356] #1 addr2line failed [193] #2 ? addr2line failed [294] #2 ? addr2line failed [191] #1 Foam::sigFpe::sigHandler(int)[274] #1 Foam::sigFpe::sigHandler(int) addr2line failed [145] #1 Foam::sigFpe::sigHandler(int) addr2line failed [293] #1 Foam::sigFpe::sigHandler(int) addr2line failed [73] #2 ? addr2line failed [227] #2 Foam::sigFpe::sigHandler(int)[321] #3 Foam::PBiCGStab::solve(Foam::Field<float>&, Foam::Field<float> const&, unsigned char) const addr2line failed [220] #1 Foam::sigFpe::sigHandler(int) addr2line failed [167] #2 ? addr2line failed addr2line failed addr2line failed addr2line failed addr2line failed [290] #2 [190] #1 Foam::sigFpe::sigHandler(int) addr2line failed addr2line failed ? addr2line failed addr2line failed [381] #1 Foam::sigFpe::sigHandler(int) addr2line failed ? addr2line failed [140] #1 Foam::sigFpe::sigHandler(int)[329] #3 Foam::PBiCGStab::solve(Foam::Field<float>&, Foam::Field<float> const&, unsigned char) const[32] #1 Foam::sigFpe::sigHandler(int) addr2line failed addr2line failed [76] #1 addr2line failed [234] #1 Foam::s........................... I have tried to use directSolveCoarsest, by this setting, it is working and get the results. But problem is my simulation become 3 times longer, which is not affordable for me. I tried other solvers, but they all unstable. So does anyone have any suggestions if I can make GAMG working without using directSolveCoarsest option. |
|
June 4, 2018, 14:04 |
|
#2 |
Senior Member
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 16 |
You can only use PCG/CG instead of gauss seidel (directSolve) in GAMG. If you say is unstable you might want to change your settings
|
|
June 4, 2018, 14:42 |
|
#3 |
Senior Member
Join Date: Sep 2013
Posts: 353
Rep Power: 21 |
I'd like to raise the point that nearly always this boils down to improper boundary conditions, discretisation or a low quality mesh. Since you did not provide this information. What is the output does this give you:
Code:
checkMesh -allTopology -allGeometry |
|
June 5, 2018, 05:33 |
|
#4 |
Member
ff
Join Date: Feb 2010
Posts: 81
Rep Power: 16 |
Thanks very much for reply. I actually tried PCG/CG, but it is unstable and get divergence some points. The GAMG is only solver I found not divergence if it is not crash.
|
|
June 5, 2018, 05:37 |
|
#5 | |
Member
ff
Join Date: Feb 2010
Posts: 81
Rep Power: 16 |
Quote:
Thanks very much for reply. I used checkMesh and all checks are ok. My case is conjugate heat transfer, so in chtMultiRegionSimpleFoam. It is crashed at solid region. Actaully, if I switch off the radiation, GAMG working perfectly. But when radiation is on, GAMG could crash like this without any warning. Is it possible that qr radiation heat flux passed through coupled boundary and created some bad condition in solid region? |
||
June 5, 2018, 05:48 |
|
#6 |
Senior Member
Join Date: Sep 2013
Posts: 353
Rep Power: 21 |
Well in that case it is even more likely that something is wrong with your setup. You are only solving a laplace equation in the solid. This shouldn't crash. You might want to lower your relaxationFactor for h or e a bit. That should probably do the trick.
|
|
June 5, 2018, 06:09 |
|
#7 |
Member
ff
Join Date: Feb 2010
Posts: 81
Rep Power: 16 |
well, my case is tricky. My solid is just a round disc, I added single rotating frame in it to simulation rotating effect like Fluent moving solid. I tried reduce relaxation Factor, but if anything below than 0.99, it will take unaffordable number of iterations to reach converge. So I cannot reduce relaxation factor.
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Suppress twoPhaseEulerFoam energy | AlmostSurelyRob | OpenFOAM Running, Solving & CFD | 33 | September 25, 2018 18:45 |
chtMultiRegionSimpleFoam turbulent case | Aditya Patil | OpenFOAM Running, Solving & CFD | 6 | April 24, 2017 23:13 |
Wrong fluctuation of pressure in transient simulation | caitao | OpenFOAM Running, Solving & CFD | 2 | March 5, 2015 22:33 |
pimpleFoam: turbulence->correct(); is not executed when using residualControl | hfs | OpenFOAM Running, Solving & CFD | 3 | October 29, 2013 09:35 |
Differences between serial and parallel runs | carsten | OpenFOAM Bugs | 11 | September 12, 2008 12:16 |