Servicios Personalizados
Revista
Articulo
Indicadores
-
Citado por SciELO
-
Accesos
Links relacionados
-
Similares en SciELO
Compartir
Revista mexicana de física
versión impresa ISSN 0035-001X
Rev. mex. fis. vol.52 no.2 México abr. 2006
Investigación
A numerical study of stiffness effects on some high order splitting methods
J. SalcedoRuíz*, F.J. SánchezBernabe**
* Instituto Politécnico Nacional, Escuela Superior de Ingeniería Mecánica y Eléctrica, Unidad Culhuacán, Av. Santa Ana No. 1000, Col. San Francisco Culhuacán, 04430 México, D.F. México.
** Univ. Autónoma MetropolitanaIztapalapa, Depto. de Matemáticas, San Rafael Atlixco No. 186, Col. Vicentina, 09340 México, D.F México.
Recibido el 18 de marzo de 2005
Aceptado el 31 de enero de 2006
Abstract
In this paper we compare operator splitting methods of first, second, third and fourth orders that are applied to problems with stiff matrices. In order to efficiently solve the resultant subproblems is necessary to use implicit RungeKutta methods. It is known that, in this context, the precision order of operator splitting schemes is reduced. Furthermore, we propose a fifth order operator splitting method that is obtained by applying Richardson extrapolation to a fourth order method. All methods are tested with a model problem with matrices such that its condition number is taken up to 20,000. Our conclusion is that order reduction is more severe for low order operator splitting methods.
Keywords: Operator splitting; stiff matrix; Richardson extrapolation; implicit RungeKutta methods.
Resumen
En este trabajo se comparan métodos de descomposición de operadores de órdenes uno, dos, tres y cuatro, que se aplican a problemas cuyas matrices son de tipo rígido. A fin de poder resolver eficientemente los problemas intermedios que aparecen es necesario aplicar métodos de RungeKutta de tipo implícito. Se ha observado que en estas condiciones, el orden de precision de los esquemas de descomposición de operadores se reduce. Se propone además un método de descomposición de operadores de orden cinco que se obtiene al aplicar extrapolación de Richardson a un esquema de orden cuatro. Todos los metodos se aplican a un problema modelo con matrices cuyo número de condición se incrementa hasta 20,000. Se concluye que el fenomeno de reducción de orden es más severa para los métodos de orden bajo.
Descriptores: Descomposición de operadores; matriz rígida; extrapolación de Richardson; métodos de RungeKutta implícitos.
PACS: 02.70.Bf; 02.90.+p; 02.60.x
DESCARGAR ARTÍCULO EN FORMATO PDF
References
1. J.G. Verwer and B. Sportisse, A Note on Operator Splitting in a Stiff Linear Case (CWI Report: MASR9830, Amsterdam, 1998). [ Links ]
2. D. Goldman and T.J. Kaper, SIAM J. Num. Analysis 33 (1996) 349. [ Links ]
3. J. Stoer and R. Bulirsch, Introduction to Numerical Analysis, 3rd. ed. (SpringerVerlag, New York, 2002). [ Links ]
4. A. Sornborger and J. Stewart, Physical Review A 60 (1999) 1956. [ Links ]
5. E. Hairer and G. Wanner, Solving Ordinary Differential Equations II: Stiff and DifferentialAlgebraic Problems (SpringerVerlag, Berlin, 1996). [ Links ]
6. R. Glowinski Handbook of Numerical Analysis, Volume IX: Numerical Methods for Fluids (Part 3). P.G. Ciarlet and J.L. Lions ed. (Elsevier Science, Amsterdam, 2003). [ Links ]
7. G. Strang, SIAM J. Num. Analysis 5 (1968) 506. [ Links ]
8. J. Salcedo y F. Sánchez, Información Tecnológica, Revista Internacional: 14 (2003) 69. [ Links ]
9. S. Descombes, Mathematics of Computation 70 (2000) 1481. [ Links ]
10. K. Dekker and J.G. Verwer, Stability of RungeKutta methods for stiff nonlinear differential equations (NorthHolland, Amsterdam, 1984). [ Links ]
11. F.J. Sánchez, Computers Math. Applic. 36 (1998) 71. [ Links ]