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

Implicit iterative Gauss-Seidel method

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 23, 2010, 06:11
Default Implicit iterative Gauss-Seidel method
  #1
New Member
 
Yasin
Join Date: Nov 2010
Posts: 4
Rep Power: 6
dearboy is on a distinguished road
I have developed a numerical code. I want to solve the Navier-stokes equations with implicit method. I used point by point Gauss-Seidel 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 gauss-seidel method, unconditionally stable for my problem(even for large CFL~100 )?

2) what other methods do you suggest instead of Gauss-Seidel (direct methods or iterative methods)?
dearboy is offline   Reply With Quote

Old   November 24, 2010, 12:31
Default
  #2
New Member
 
Join Date: Jun 2009
Posts: 27
Rep Power: 8
dshawul is on a distinguished road
Quote:
Originally Posted by dearboy View Post
I have developed a numerical code. I want to solve the Navier-stokes equations with implicit method. I used point by point Gauss-Seidel 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 gauss-seidel method, unconditionally stable for my problem(even for large CFL~100 )?

2) what other methods do you suggest instead of Gauss-Seidel (direct methods or iterative methods)?
(1)I am not an expert but I think that the method used for iterative solution is not affected by the CFL number. The CFL number is used to control the stability of the time discretization scheme (i.e for explicit or crank-nicholson). Gauss seidel will always converge for a diagonally dominant matrix.
(2) You can try conjugate gradient methods which are much faster and relatively easy to implement.
dshawul is offline   Reply With Quote

Old   November 24, 2010, 22:39
Default
  #3
New Member
 
Yasin
Join Date: Nov 2010
Posts: 4
Rep Power: 6
dearboy is on a distinguished road
Many thanks Daniel,
your help is useful.
The used gauss-seidel in my code converges very good( for 10-20 iteration with under-relaxation factor about 0.65). I have used three-layer 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.
dearboy is offline   Reply With Quote

Old   November 28, 2010, 07:24
Default
  #4
New Member
 
Dennis Ernens
Join Date: Sep 2010
Posts: 9
Rep Power: 6
Dennis is on a distinguished road
The three-layer backward difference scheme is unconditionally stable, so it will not explode for high CFL numbers. Furthermore it is second order accurate (like Crank-Nicholson).

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.
Dennis is offline   Reply With Quote

Old   November 28, 2010, 21:59
Default
  #5
New Member
 
Yasin
Join Date: Nov 2010
Posts: 4
Rep Power: 6
dearboy is on a distinguished road
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.
dearboy is offline   Reply With Quote

Old   November 29, 2010, 10:18
Default
  #6
New Member
 
Join Date: Jun 2009
Posts: 27
Rep Power: 8
dshawul is on a distinguished road
Quote:
Originally Posted by dearboy View Post
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.
Hi,
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.
dshawul is offline   Reply With Quote

Reply

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
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 CD-adapco 0 March 12, 2003 13:07
implicit method on unstructured grid Koo Main CFD Forum 5 July 29, 2001 01:43


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