Aproximación numérica de la distribución de temperatura en una barra metálica utilizando un esquema de diferencias finitas

Resumen

En el presente trabajo de investigación se desarrolla el problema de conducción bidimensional del calor usando el método numérico de diferencias finitas, utilizando los esquemas de diferencias finitas progresivas, retrógradas y Crank-Nicolson. En la parte matemática, se desarrollan las series de Taylor para una, dos y tres variables independientes, se obtienen los esquemas de diferencias finitas progresivas, retrógradas y Crank-Nicolson. Se deduce la ecuación de conducción bidimensional del calor que depende de la geometría del dominio y se le impone condiciones iniciales y cie frontera, de acuerdo a estas condiciones el problema se torna con un cierto grado de dificultad para su solución analítica, entonces se puede utilizar métodos de aproximación numérica como el método de diferencias finitas. Para hallar la solución aproximada se han utilizado los esquemas de diferencias progresivas y Crank-Nicolson

Para el contraste de los resultados, se realiza la comparación de las soluciones aproximadas obtenidas con los esquemas de diferencias y la solución analítica, obtenida por el método de separación de variables utilizando los mismos parámetros de entrada, llegando a determinar que el error de aproximación es muy pequeño, concluyendo así que los esquemas de diferencias utilizados resultan eficientes en su aplicación. Se ha implementado computacionalmente los esquemas de diferencias finitas progresivas y de Crank-Nicolson para el problema de conducción bidimensional del calor en una placa metálica rectangular con condiciones de frontera de Dirichlet y condiciones iniciales dadas, logrando determinar el comportamiento de la temperatura para diferentes tiempos