CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Wiki > Stone's method

Stone's method

From CFD-Wiki

(Difference between revisions)
Jump to: navigation, search
(References: wikilink)
Line 29: Line 29:
#{{reference-paper|author=Acosta, J.M.|year=2001|title=Numerical Algorithms for Three Dimensional Computational Fluid Dynamic Problems|rest=PhD Thesis, Polytechnic University of Catalunia}}
#{{reference-paper|author=Acosta, J.M.|year=2001|title=Numerical Algorithms for Three Dimensional Computational Fluid Dynamic Problems|rest=PhD Thesis, Polytechnic University of Catalunia}}
-
 
+
== External link ==
 +
*[http://en.wikipedia.org/wiki/Stone_method Wikipedia article on Stome method]
----
----
<i> Return to [[Numerical methods | Numerical Methods]] </i>
<i> Return to [[Numerical methods | Numerical Methods]] </i>

Revision as of 07:33, 15 February 2006

Contents

Introduction

We have seen that LU decomposition is an excellent general purpose linear equation solver. The biggest disadvantage of these type of solver is they fail to take advantage of coefficient matrix. In preconditioned iterative methods, if the preconditioner matrix M is a good approximation of coefficient matrix A then the convergence is faster. This brings us to idea of using approximate factorization LU of A as the iteration matrix M. A version of incomplete lower-upper decomposition method was proposed by Stone (1968). This method is also called the strongly implicit procedure or SIP. This method is designed for equation system arising from discretisation of partial differential equations. This method does not apply to generic system of equations.

Algorithm

For the linear system Ax = b
calculate Incomplete LU factorization of matrix A
Ax = (M-N)x = (LU-N)x = b
Mx(k+1) = Nx(k)+b , with ||M|| >> ||N||
Mx(k+1) = LUx(k+1) = c(k)
LUx(k) = L(Ux(k+1)) = Ly(k) = c(k)
set a guess
k = 0, x(k)
r(k)=b - Ax(k)
while ( ||r(k)||2  \ge \varepsilon ) do
evaluate new right hand side
c(k) = Nx(k) + b
solve Ly(k) = c(k) by forward substitution
y(k) = L-1c(k)
solve Ux(k+1) = y(k) by back substitution
x(k+1) = U-1y(k)
end while

References

  1. Ferziger, J.H. and Peric, M. (2001), Computational Methods for Fluid Dynamics, ISBN 3540420746, 3rd Rev. Ed., Springer-Verlag, Berlin..
  2. Acosta, J.M. (2001), "Numerical Algorithms for Three Dimensional Computational Fluid Dynamic Problems", PhD Thesis, Polytechnic University of Catalunia.

External link


Return to Numerical Methods

My wiki