Number of iteration difference between PCG and GAMG

hr_nt · November 22, 2022, 06:52

I ran two interFOAM multiphase simulation in which one had the p_rgh solver set as PCG and one in which it was GAMG. The preconditioner for both was set to DIC. I observed that the number of iterations taken to converge to 1e-8 is higher for PCG (about 1000) and lower for GAMG (about 40). Why does this occur ?

dlahaye · November 22, 2022, 07:36

GAMG requires less iterations because it iterates on various levels of accuracy. See literature on Multigrid Methods, e.g., book of Moukalled-Darwish-Mangani or https://en.wikipedia.org/wiki/Multigrid_method.

I suggest using GAMG as preconditioner for CG to combine best of two worlds (zat least in case that number of processors is one or small).

November 22, 2022, 06:52	Number of iteration difference between PCG and GAMG	#1
hr_nt New Member HR Join Date: Nov 2022 Posts: 17 Rep Power: 4	I ran two interFOAM multiphase simulation in which one had the p_rgh solver set as PCG and one in which it was GAMG. The preconditioner for both was set to DIC. I observed that the number of iterations taken to converge to 1e-8 is higher for PCG (about 1000) and lower for GAMG (about 40). Why does this occur ?

November 22, 2022, 07:36		#2
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 802 Blog Entries: 1 Rep Power: 18	GAMG requires less iterations because it iterates on various levels of accuracy. See literature on Multigrid Methods, e.g., book of Moukalled-Darwish-Mangani or https://en.wikipedia.org/wiki/Multigrid_method. I suggest using GAMG as preconditioner for CG to combine best of two worlds (zat least in case that number of processors is one or small). joshwilliams and hr_nt like this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
LES, Courant Number, Crash, Sudden	Alhasan	OpenFOAM Running, Solving & CFD	5	November 22, 2019 03:05
[solids4Foam] How to calculate drag coeff when using solids4Foam	amuzeshi	OpenFOAM CC Toolkits for Fluid-Structure Interaction	15	November 7, 2019 13:50
Large number of iterations using PCG in pressure equation	mgg	OpenFOAM Running, Solving & CFD	3	January 8, 2019 03:23
p_rgh initial residual no change with different settings	manuc	OpenFOAM Running, Solving & CFD	3	June 26, 2018 16:53
should Courant number always be kept below 1?	wc34071209	OpenFOAM Running, Solving & CFD	16	March 9, 2014 20:31