El metodo de diferencias nitas aplicado al estudio en la distribucion de temperaturas en una barra unidimensional sometida a una fuente de calor
Jose Angel Neria Perez Bajo la instruccion de la doctora Lucia Medina G. Gilberto Gomez C. y Mauricio V. Noviembre de 2011
Resumen.
mas real.
Se considera una barra unidimensional de longitud L la cual esta sumergida en un bath ter- 3 Estudio analitico. mico a 0o , posteriormente se le saca de ese bath La ecuacion que debe cumplir la distribucion y se coloca de tal modo que solo sus extremos de temperaturas en cualquier tiempo y en estan en bath termico a 0o y se le coloca una cualquier punto de la barra esta dada como fuente de calor de la forma q(x, t) = xet . Se comparan las soluciones analitica y numerica, la segunda provista por el uso del metodo de diferencias nitas u 2u = k 2 2 + q(x, t). (1) t x
Los extremos de la barra se mantienen o Usualmente al estudiar sistemas sicos surgen a 0 , entonces la solucion debe cumplir las escuaciones diferenciales, ya sea ordinarias o conciones (a la frontera) parciales, las cuales no siempre son faciles de resolver, pero al hacer ciertas asumpciones sou(0, t) = u(L, t) = 0. bre el sistema dichas ecuaciones pueden ser reducidas a otras cuya solucion analitica es cono- Como al principio del estudio la barra pero cida, pero hacer estas simplicaciones implica manece a 0 se tiene la condicion inicial estudiar una caricatura del sistema pero, bajo dichas simplicaciones es de gran interes tanto teorico como practico saber si dado un metodo numerico se halla la misma solucion al probu(x, 0) = 0. (2) lema simplicado, si lo hace, entonces con algunas modicaciones el metodo puede ser us- De aqui se deduce que la distribucion de ado para resolver la ecuacion sin las simpli- temperaturas es de la forma caciones hechas, es decir, estudiar un sistema 1
Introduccion.
Como esta ecuacion se cumple para toda x en el intervalo [0.L] entonces, el coeciente de n u(x, t) = An (t) sin x . (3) sin n x es el que se debe de anular, por n=1 L L tanto se debe cumplir la ecuacion Por otro lado, la fuente q(x,t) puede ser escrita de la forma q(x, t) = bn (t) sin n=1 n x . L
2L (1)n et . (7) n Entonces, aplicando la ortogonalidad de la para cada numero natural n. n x sobre el intervalo [0,L] se La ecuacion (7) es de primer orden de la forma funcion sin L obtiene An (t) + 2 An (t) = bn = n A + hA = wet . bn (t) = 2 e t L
L 0 x sin
n x dx. L
Cuya solucion al aplicar la condicion (2) A= w (et ceht ) h1
bn (t) = Por lo cual
2L(1)n t e . n
(4) entonces,
A n (t) = 2L3 (1)n q(x, t) = 2L t (1)n n e n=1 sin x n L (5) n(L2 k 2 n2 2 )
k 2 n2 2 t t L2 . e e
(8)
Substituyendo (3) y (5) en (1) se obtiene An (t) + 2 An (t) bn n n=1 n sin x = 0, L con
Por tanto, u(x,t) esta dada como
u(x, t) = n et en t sin n=1 con n = 2L3 (1)n , n(L2 (nk)2 ) n =
n x . L (9)
n =
kn . L
(6)
nk . L Observese que n presenta indeterminaciones cuando L = nk.
3.1
Comportamiento de la solucion. pseudocodigo para hallar la solucion a esta
ecuacion, a saber,
Debido a como se escogio a la fuente q(x, t) se espera que la temperatura de la barra se perturbe de su estado original, pero despues regrese a ese mismo, a saber a 0o . Tomando e 2.7, se observa la fuente molesta a la barra por un tiempo menor a 20s, por ejemplo, a una posicion ja y t=10s, q(x, t) < 104 .
ui,j+1 = (1 2s)ui,j + s(ui+1,j + ui1,j ) +q(x(i), t(j)), (11)
De aqui se deduce que un tiempo de estudio con razonable es tomar t 10. Por otro lado, u(x, t) comienza siendo nula en toda la barra, comienza a crecer, alcanza un maximo y despues decae hacia 0o . Para observar esto k2 s= 2. gracamente se tomaran los valores particuh lares Sujeto a las condiciones de frontera k = 1, L = , 2
(12)
con lo cual u(x, t) toma la forma
u(0, j) = f (ih) 0 i m, u(0, j) = u(n, j) = 0 0 < j.
u(x, t) = n n=1 y
Donde h es el numero de nodos en el que se divide el intervalo [0,L], f es la condicion 2t inicial a la cual esta sujeta la barra, en este e4n et sin(2nx) . caso f = 0 y es el size de la particion que se (10) le dara al intervalo temporal. Se tomaran 5 nodos y = 0.01,con lo cual el perl de temperaturas se ve de la forma (1)n . n(4n2 1)
n =
Con ayuda de Matlab, y truncando la suma hasta 10 terminos, el perl de temperaturas se Como se puede notar, los resultados numericos ve de la siguiente forma. y los obtenidos por medio de la solucion analitica son muy similares, pero se observo un comportamiento extrano al intentar usar mas nodos y mas small, por ejemplo, con 10 nodos 4 Estudio numerico. y = 0.03 se vio el siguiente comportamiento Se usara el metodo de diferencias nitas, de la solucion. Esto debido a que el error crece el cual, en [1] esta esquematizado por el al aumentar el size de las particiones. 3
Discusion de resultados.
Bibliograa.
[1] N. H. Asmar,Partial Dierential Equations With Fourier Series And Boundary Value Problems, Second Edition, Prentice Hall, New Jersey, 2005. [2] W. Cheney and D. Kincaid, Numerical Mathematics and Computing, Fifth Edition, University of Texas, Austin, 2004. [3] S. Kottwitz, LaTeX: Beginners Guide, Packt Publishing, 2011.
