|January 31, 2011, 13:22||
using petsc for large coupled linear systems
Join Date: Jan 2011
Posts: 1Rep Power: 0
Hello to all,
I am in trouble using petsc.
I have solved Free-Convection flow in annulus on a 90x100 grid using PGMRES-ILUT of Saad's SPARSKIT. I am using a coupled algorithm with PWIM for face velocities to solve continuty equation with momentum and energy equations simultaneously. to It worked well up to 90x100 cells (90x100x4 equations) with non-dimensionalized time step t=5. To study the grid independency, I used 180x200 grid with the same solver settings and B.C.'snd t=5. However, it does not converge on the fine grid. I reduced the time step size up to t=0.001, but it did not work. Up to now, I have no idea about this issue: it is caused by the solver or the problem is too taugh to solve.
Therefore, I used PETSc in order to utilize more options on both preconditioners as well as solver and also, to use its parallel processing benefit.
I tests almost all preconditioners and ksp solver of petsc and they can not find a solution (do not converege) even though, for 90x100 cells!
1- How could I use the same preconditioner and solver (i.e. ILUT and GMRES with fill-in=100 and droptol=0.01) by petsc to get the same result for 90x100 cells?
2- Do you have any idea about what makes the solver to diverge when I use 180x200?
|Thread||Thread Starter||Forum||Replies||Last Post|
|Inversion of large sparse linear systems||pankos||Main CFD Forum||4||June 24, 2010 06:24|
|solution diverges when linear upwind interpolation scheme is used||subash||OpenFOAM||0||May 29, 2010 01:23|
|IdeasUnvToFoam Bug amp Fix||benru||OpenFOAM Bugs||42||November 13, 2009 08:59|
|Water pump OpenFOAM 15 ANSYS CFX 110 comparation||waynezw0618||OpenFOAM Running, Solving & CFD||39||March 5, 2009 13:57|
|Non linear coupled boundaies||tamer||Main CFD Forum||0||June 12, 2003 02:46|