
[Sponsors] 
November 23, 2010, 06:11 
Implicit iterative GaussSeidel method

#1 
New Member
Yasin
Join Date: Nov 2010
Posts: 4
Rep Power: 7 
I have developed a numerical code. I want to solve the Navierstokes equations with implicit method. I used point by point GaussSeidel method. I applied the developed code for the flow around a cylinder. I use the finite volume method with body fitted mesh.the total mesh number is about 600 000. Also I use the collocated mesh and the dual time stepping method and LES for turbulence.
my questions: 1) is the gaussseidel method, unconditionally stable for my problem(even for large CFL~100 )? 2) what other methods do you suggest instead of GaussSeidel (direct methods or iterative methods)? 

November 24, 2010, 12:31 

#2  
New Member
Join Date: Jun 2009
Posts: 27
Rep Power: 9 
Quote:
(2) You can try conjugate gradient methods which are much faster and relatively easy to implement. 

November 24, 2010, 22:39 

#3 
New Member
Yasin
Join Date: Nov 2010
Posts: 4
Rep Power: 7 
Many thanks Daniel,
your help is useful. The used gaussseidel in my code converges very good( for 1020 iteration with underrelaxation factor about 0.65). I have used threelayer backward difference scheme for time discretization. It seems that this scheme should be unconditionally stable( I'm not sure!). So may be my time discretization is not suitable or not unconditionally stable. I will change it to Crank–Nicolson. If you have any helpful suggestion please let me know about it. Many thanks. 

November 28, 2010, 07:24 

#4 
New Member
Dennis Ernens
Join Date: Sep 2010
Posts: 9
Rep Power: 7 
The threelayer backward difference scheme is unconditionally stable, so it will not explode for high CFL numbers. Furthermore it is second order accurate (like CrankNicholson).
But, do not forget that you still want to compute physics hence you still need to set a time step appropriate for your problem. For the cylinder you want to resolve the oscillations of the wake for instance. For the linear system solve you can also look at Krylov subspace methods (like conjugate gradient suggested earlier), these are generally faster. If you are using a pressure correction method, you can also use a direct solver for the poisson term. 

November 28, 2010, 21:59 

#5 
New Member
Yasin
Join Date: Nov 2010
Posts: 4
Rep Power: 7 
Hi Dennis,
Many thanks for your reply. I know that I should use appropriate time steps. but I have used dual time stepping (Artificial compressibility method) in my solution and this method in each physical time need to be converged (with iterations). if we use the high CFL number in virtual (pseudo) time we can reduce the number of iterations. regarding to this issue, I need this unconditionally stable methods to increase the time steps in virtual time. I should mention that despite using the implicit method in my numerical code, my time steps is very small (about 0.00005 s). I don't know why for larger time steps my implicit solution is not stable. 

November 29, 2010, 10:18 

#6  
New Member
Join Date: Jun 2009
Posts: 27
Rep Power: 9 
Quote:
If you use PISO you still have to use small time steps. The momentum discretization is assumed to be frozen while doing multiple pressure corrections, which is true only at small time steps. More at page 152 of this document http://www.scribd.com/doc/27927305/N...outer_page_152 Edit: I see now you use Artificial copressibilty method, which I am not familiar with. Anyway implicit methods may still require small time steps for reasons like that. Last edited by dshawul; November 29, 2010 at 13:18. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Problem in fvschemes divSchemes cannot use Gauss linearUpwind  nico765  OpenFOAM Bugs  27  June 3, 2014 05:29 
About joint face in implicit method  quaintchewster  Main CFD Forum  0  November 11, 2010 08:04 
solution diverges when linear upwind interpolation scheme is used  subash  OpenFOAM  0  May 29, 2010 01:23 
Implicit method of FSI  cjtune  CDadapco  0  March 12, 2003 13:07 
implicit method on unstructured grid  Koo  Main CFD Forum  5  July 29, 2001 01:43 