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

advice needed in solving linear system

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 3, 2001, 07:31
Default advice needed in solving linear system
  #1
yf yap
Guest
 
Posts: n/a
dear sirs:

in solving a 2 dimensional poisson equation, i obtained a set of linear equations in the following form:

C1*phi[i,j] + C2*phi[i+1,j] + C3*phi[i-1,j] + C4*phi[i,j+1] + C5*phi[i,j-1] = B[i,j] -[1]

C1 is always the largest. for all the equations,

C1 = C2 + C3 + C4 + C5 ----------- [2]

when this set of equations is solved using SOR, it is found that the solutions diverge slowly.

Questions: 1) how to avoid divergent in solving all those equations?

2) does [2] still imply strongly diagonally dominance ?

any advice would be greatly appreciated, thanks.

regards, yfyap
  Reply With Quote

Old   February 3, 2001, 21:16
Default Re: advice needed in solving linear system
  #2
John C. Chien
Guest
 
Posts: n/a
(1). First of all, the coefficients C1, C2, C3, C4 and C5 will all have effect on the results. (2). I would just create a five-point template problem, and check out the SOR result. You could do the hand calculations to look at the result of each iteration. Naturally, with the appropriate boundary conditions applied. (3). If the phi values are known at the boundary, then you should be able to get the solution quickly. If gradient type of b.C.'s is used, it will take longer to converge. But I guess, you can try different under-relaxation parameter values. But you still need to know the behavior of the coefficients.
  Reply With Quote

Old   February 4, 2001, 20:24
Default Re: advice needed in solving linear system
  #3
John C. Chien (Alter Ego)
Guest
 
Posts: n/a
Check out the book:

Templates for the Solution of Linear Systems in the resource section. I found it useful. You might also try some different solvers like CGSTAB, maybe even start with a brute force algorithm which always converges to see if you have set up your problem correctly.
  Reply With Quote

Old   February 5, 2001, 00:32
Default Re: advice needed in solving linear system
  #4
Abhijit Tilak
Guest
 
Posts: n/a
hi,

u will run in to trouble if C1,C2,... turn out to be negative(there are caes where they may go negative). just check that. Consider using Line Sor(LSOR) or ADI that will acclerate convergence. BTW from ur past postings looks like ur posting ur classroom course assignment questions on this forum eh? u may also consider under-relaxing as jhon said (SUR) hope that helps abhijit
  Reply With Quote

Old   February 5, 2001, 01:21
Default Re: advice needed in solving linear system
  #5
yf yap
Guest
 
Posts: n/a
thanks, John. regards, yfyap
  Reply With Quote

Old   February 5, 2001, 01:21
Default Re: advice needed in solving linear system
  #6
yf yap
Guest
 
Posts: n/a
dear abhijit, your valuable advice is greatly appreciated. please note that i'm not posting my assignments in this forum. for the previous postings, i feel that it is necessary for me to decribe my situations properly so that experts like your guys can have a better view of the situation. thanks. regards, yfyap
  Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Forces in OF15 richard OpenFOAM Running, Solving & CFD 180 July 9, 2018 11:54
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM Running, Solving & CFD 8 March 28, 2016 12:08
Velocity blows up suddenly after 30,000+ iterations lordvon OpenFOAM Running, Solving & CFD 15 October 19, 2015 14:52
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 06:24
Error log vw.cfd OpenFOAM 6 August 7, 2009 06:44


All times are GMT -4. The time now is 16:41.