A numerical study of stiffness effects on some high order splitting methods



Título del documento: A numerical study of stiffness effects on some high order splitting methods
Revista: Revista mexicana de física
Base de datos: PERIÓDICA
Número de sistema: 000461155
ISSN: 0035-001X
Autores: 1
2
Instituciones: 1Instituto Politécnico Nacional, Escuela Superior de Ingeniería Mecánica y Eléctrica, México, Distrito Federal. México
2Universidad Autónoma Metropolitana, Departamento de Matemáticas, Iztapalapa, Distrito Federal. México
Año:
Periodo: Abr
Volumen: 52
Número: 2
Paginación: 129-134
País: México
Idioma: Inglés
Tipo de documento: Artículo
Enfoque: Analítico, teórico
Resumen en español 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 Runge-Kutta 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
Resumen en inglés 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 Runge-Kutta 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
Disciplinas: Física y astronomía
Palabras clave: Física,
Descomposición de operadores,
Matriz rígida,
Extrapolación de Richardson,
Métodos de Runge-Kutta implícitos
Keyword: Physics,
Operator splitting,
Stiff matrix,
Richardson extrapolation,
Implicit Runge-Kutta methods
Texto completo: Texto completo (Ver HTML) Texto completo (Ver PDF)