Resumen
En este artículo se presenta un nuevo método mimético de diferencias finitas para resolver la ecuación no estática de difusión. Éste usa el esquema de Crank-Nicholson para obtener aproximaciones en tiempo y discretizaciones miméticas de segundo orden, para los operadores gradiente y divergencia, en el espacio. La convergenica de este nuevo método es analizada usando el teorema de equivalencia de Lax-Friedrichs. Este análisis es desarrollado para el caso unidimensional. Además del estudio teórico, se dan pruebas prácticas que evidencian que el esquema mimético tipo Crank-Nicholson es mejor que el esquema tradicional de diferencias finitas ya que arroja tasas de convergencia cuadráticas, errores de truncamiento de segundo orden y mejores aproximaciones a la solución exacta.
Citas
Castillo, J.E.; Grone, R.D. (2003) “A matrix analysis approach to higher-order approximations for divergence and gradients satisfying a global conservation law”, SIAM Journal on Matrix Analysis and Applications 25(1): 128–142.
Freites-Villegas, M.; Guevara-Jordan, J.M.; Rojas, O.; Castillo, J.E.; Rojas S. (2004) “A mimetic finite difference scheme for solving the steady state diffusion equation with singular sources”, VII International Congress of Numerical Methods in Engineering and Science, San Cristobal, Venezuela.
Castillo, J.E.; Yasuda, M. (2005) “Linear system arising for second order mimetic divergence and gradient operators”, Journal of Mathematical Modeling and Algorithm 4(1): 67–82.
Shashkov, M. (1996) “Conservative Finite–Differences Methods on General Grids”, Symbolic and Numerical Computation Series, CRC Press, Boca Raton FL.
Guevara-Jordan, J.M.; Rojas, S.; Freites-Villegas, M.; Castillo J.E. (2007) “Convergence of a mimetic finite difference method for static diffusion equation”, Advances in Difference Equations, Volume 2007, Article ID 12303, 12 pages.
Mannarino, I. (2007) “Un método mimético de diferencias finitas para la ecuación no estática de difusión”, Master thesis, Universidad Central de Venezuela, Caracas.
Mannarino, I.; Quintana, Y.; Guevara-Jordan, J.M. (2007) “A numerical study of mimetic scheme for the unsteady heat equation”, submitted to FACYT Review (Faraute).