CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

using petsc for large coupled linear systems

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

LinkBack Thread Tools Display Modes
Old   January 31, 2011, 13:22
Default using petsc for large coupled linear systems
New Member
Join Date: Jan 2011
Location: canada
Posts: 1
Rep Power: 0
sd_vakip is on a distinguished road
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?

Many Thanks,
sd_vakip is offline   Reply With Quote


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
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

All times are GMT -4. The time now is 07:21.