CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

GAMG issue

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   December 2, 2012, 16:35
Default GAMG issue
  #1
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 272
Rep Power: 9
phsieh2005 is on a distinguished road
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
phsieh2005 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


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 07: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 09:52
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 11:16


All times are GMT -4. The time now is 13:34.