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

New AMG Solvers in OF141dev SVN

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

Like Tree4Likes
  • 1 Post By philippose
  • 3 Post By hjasak

Reply
 
LinkBack Thread Tools Display Modes
Old   October 27, 2007, 14:05
Default Hello and a Good day :-)! I
  #1
Senior Member
 
Philippose Rajan
Join Date: Mar 2009
Location: Germany
Posts: 530
Rep Power: 16
philippose will become famous soon enough
Hello and a Good day :-)!

I just finished compiling the latest SVN version of OpenFOAM-1.4.1-dev, and was generally looking around the sources, and found that there are new AMG solvers available....

Where could I get more information regarding what the advantages / disadvantages of, and the situations in which the following solvers can be used:

a. fpeAmgSolver
b. mpeAmgSolver
c. rreAmgSolver
d. amgSymSolver

And also... the following solvers:

d. bicgStabSolver
e. cgSolver
f. deflationSolver
g. gmresSolver

In addition, what do the "AMG Policies" do? There too.. there seem to be a couple of options. I remember finding a paper by Hrv some time ago with a performance comparison of some of these solvers, but I cant seem to find that paper now!

Have a great weekend!

Philippose
mm.abdollahzadeh likes this.
philippose is offline   Reply With Quote

Old   October 27, 2007, 17:34
Default A paper, of course: http:/
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
A paper, of course:

http://powerlab.fsb.hr/ped/kturbo/OpenFOAM/papers/CFD2007CanadaPrecon.pdf

and

http://powerlab.fsb.hr/ped/kturbo/OpenFOAM/papers/CFD2007CanadaExtrapolation.pdf

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   October 27, 2007, 19:04
Default Hello Hrv :-)! Thanks a lot
  #3
Senior Member
 
Philippose Rajan
Join Date: Mar 2009
Location: Germany
Posts: 530
Rep Power: 16
philippose will become famous soon enough
Hello Hrv :-)!

Thanks a lot for those two papers.... the first of those two was the one I had seen earlier, and the one I was referring to in the post :-)!

And... nice to hear from you after a long time too !

Shall read through those two tomorrow... and see if its applicable to the simulations I work on.

Regards,
Philippose
philippose is offline   Reply With Quote

Old   October 28, 2007, 11:24
Default I just compiled the OpenFOAM-1
  #4
New Member
 
Qiang Li
Join Date: Mar 2009
Location: Shenyang, LiaoNing, P. R. China
Posts: 14
Rep Power: 8
frank is on a distinguished road
I just compiled the OpenFOAM-1.4.1-dev- 26oct version, got a error. I guess GAMG solver setting changes.

Hrv, please give me a hint!

"

PBiCG: Solving for alpha, Initial residual = 0.00018123, Final residual = 1.09717e-20, No Iterations 1
Dispersed phase volume fraction = 0.275 Min(alpha) = 0 Max(alpha) = 0.55
PBiCG: Solving for alpha, Initial residual = 5.20764e-06, Final residual = 4.23583e-22, No Iterations 1
Dispersed phase volume fraction = 0.275 Min(alpha) = -2.85747e-101 Max(alpha) = 0.55


--> FOAM FATAL IO ERROR : Cannot read control minIter

file: -1.4.1-dev/run/twoPhaseEulerFoam/bed/system/fvSolution::p from line 29 to line 40.

From function void lduMatrix::solver::readControl
(
const dictionary& controlDict,
T& control,
const word& controlName
)
in file matrices/lduMatrix/lduMatrix/lduMatrixTemplates.C at line 57.

FOAM aborting

Aborted (core dumped)
"

Li Q
frank is offline   Reply With Quote

Old   October 28, 2007, 11:59
Default This problem have been solved
  #5
New Member
 
Qiang Li
Join Date: Mar 2009
Location: Shenyang, LiaoNing, P. R. China
Posts: 14
Rep Power: 8
frank is on a distinguished road
This problem have been solved through system/controlDict
solvers
{
p GAMG
{
minIter 0; //
maxIter 500; //
tolerance 1e-8;
relTol 0;

smoother DIC;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;

scaleCorrection true;
directSolveCoarsest false;

cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
};

Li Q
frank is offline   Reply With Quote

Old   October 28, 2007, 16:35
Default Simply, I don't like the idea
  #6
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Simply, I don't like the idea of hidden optional parameters, so I've made them all required. minIter tells you how many solver sweeps will be forced without checking the convergence tolerance: usually zero.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   November 21, 2007, 11:24
Default Hi, ... it seems that the
  #7
Senior Member
 
Fabian Braennstroem
Join Date: Mar 2009
Posts: 407
Rep Power: 10
braennstroem is on a distinguished road
Hi,

... it seems that the last post got lost....

I wondered how the the GAMG solvers are implemented in 1.4.1-dev? As I read the old AMG solver were transfering there data at every level in a parallel setup. Is it the same procedure for the new solver? Isn't it possible to do the interchange just at the finest level and pass the value to the coarse one; a lot of transfer time could be saved?

Fabian
braennstroem is offline   Reply With Quote

Old   November 26, 2007, 05:29
Default Hi,Prof. Hrvoje Jasak Your AM
  #8
New Member
 
Harry
Join Date: Mar 2009
Posts: 27
Rep Power: 8
harry is on a distinguished road
Hi,Prof. Hrvoje Jasak
Your AMG solver sounds amazing. I plan to use it in my specific multiphase CFD code, but do not want to completely give up my own in-house code. This is possible or not? I mean have you planed to release it as a independent solver of linear equations,in particular for CFD, as done by hype and PETSc?
Regards,
Harry
harry is offline   Reply With Quote

Old   July 27, 2009, 14:52
Default
  #9
Senior Member
 
John Deas
Join Date: Mar 2009
Posts: 160
Rep Power: 8
johndeas is on a distinguished road
Would it be possible to have an example case of the fastest one from the paper, in order to test the performance easily, and quickly use it on own cases ? For example, the rreAmg ?
johndeas 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
OF141dev undefined reference to METIS symbols on SUSE Enterprise matthias OpenFOAM Installation 5 August 29, 2008 10:55
OF141dev installation Woes chegdan OpenFOAM Installation 13 July 18, 2008 17:16
Compiling OF141dev with GCC430 matthias OpenFOAM Installation 5 June 23, 2008 13:39
Installation of OF141dev on redhatEL 64bits problems with Allmake young OpenFOAM Installation 3 May 1, 2008 00:42
Solvers DB CD-adapco 3 December 6, 2005 09:26


All times are GMT -4. The time now is 08:55.