|
[Sponsors] |
December 2, 2012, 16:35 |
GAMG issue
|
#1 |
Senior Member
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 334
Rep Power: 18 |
Dear OpenFOAMers,
I am solving a chtMultiRegionSimpleFoam case. On the air domain, for p_rgh, when I used GAMG solver, I got floating point exception right away, however, when I switched to PCG with DIC preconditioner, the case run fine - just converging slower. I am wondering what settings I can change to resolve the GAMG issue. Thanks! Pei-Ying ------------------------ GAMG setting used: p_rgh { solver GAMG; tolerance 1e-7; relTol 0.01; smoother GaussSeidel; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; maxIter 100; } -------------------- PCG that worked solver PCG; preconditioner DIC; tolerance 1e-7; relTol 0.01; ----------------- error messages: #0 Foam::errorprintStack(Foam::Ostream&) in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so" #2 in "/lib64/libc.so.6" #3 FoamdICPreconditioner::calcReciprocalD(Foam::Field <double>&, Foam::lduMatrix const&) in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so" #4 FoamdICPreconditionerdICPreconditioner(Foam::lduMa trix::solver const&, Foam::dictionary const&) in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so" #5 Foam::lduMatrixpreconditioner::addsymMatrixConstru ctorToTable<FoamdICPreconditioner>::New(Foam::lduM atrix::solver const&, Foam::dictionary const&) in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so" #6 Foam::lduMatrixpreconditioner::New(Foam::lduMatrix ::solver const&, Foam::dictionary const&) in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so" #7 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so" #8 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so" #9 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so" #10 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so" #11 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libfiniteVolume.so" #12 Foam::fvMatrix<double>::solve() in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/bin/chtMultiRegionSimpleFoam" #13 in "/home/hsieh/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/bin/chtMultiRegionSimpleFoam" #14 __libc_start_main in "/lib64/libc.so.6" #15 at /home/abuild/rpmbuild/BUILD/glibc-2.15/csu/../sysdeps/x86_64/elf/start.S:116 Floating point exception |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
pimpleFoam: turbulence->correct(); is not executed when using residualControl | hfs | OpenFOAM Running, Solving & CFD | 3 | October 29, 2013 09:35 |
GAMG solver causes troubles | colinB | OpenFOAM Running, Solving & CFD | 5 | June 13, 2013 08:48 |
Interfoam blows on parallel run | danvica | OpenFOAM Running, Solving & CFD | 16 | December 22, 2012 03:09 |
rhoPorousMRFSimpleFoam/compressor simulation blows up | mpat | OpenFOAM Running, Solving & CFD | 2 | August 10, 2012 10:52 |
Differences between serial and parallel runs | carsten | OpenFOAM Bugs | 11 | September 12, 2008 12:16 |