Sobre la solución numérica de la ecuación de onda

Abstract

La solución numérica de ecuaciones diferenciales parciales que evolucionan en el tiempo es un área de trabajo en constante desarrollo. En este trabajo se aborda la solución computacional de la ecuación de onda bajo dos métodos para problemas que involucran el tiempo: el de Newmark y el de diferencias finitas (DF). El método de Newmark tiene una alta precisión y una excelente tasa de convergencia, comparado con el de DF. Por su parte, el método de DF es de fácil implementación. Con el ánimo de comparar estos dos métodos se han puesto en funcionamiento dos problemas típicos-test en FORTRAN: el de una membrana cuadrada totalmente fija en sus bordes con una velocidad inicial en el centro y una viga simplemente empotrada con una velocidad inicial en uno de sus extremos. Cada uno de estos problemas son implementados, espacialmente, con el método de los elementos finitos y, temporalmente, con Newmark y DF. Los resultados muestran que Newmark permite utilizar pasos de tiempo mayores que DF, pero presentan mayor oscilación numérica. Con estos resultados se espera obtener datos iniciales para comparar con otros métodos que serán implementados posteriormente.