Práctica de Métodos de la Física Matemática.
Área de Física Teórica. Dpto. Física. UEx

[Graphics:Images/DifusMetExpli2D.nb_gr_1.gif]: Santos Bravo Yuste. email:  santos@unex.es

Resolución numérica mediante el método explícito de la ecuación difusiva bidimensional
                             [Graphics:Images/DifusMetExpli2D.nb_gr_2.gif]= [Graphics:Images/DifusMetExpli2D.nb_gr_3.gif]+[Graphics:Images/DifusMetExpli2D.nb_gr_4.gif] .

En este cuaderno de Mathematica vamos a resolver la ecuación difusiva bidimensional anterior en un cuadrado de lado L donde los puntos situados sobre una cruz centrada en (L/2,L/2) se fijan a temperatura 1.  Esta cruz tiene brazos de tamaño  L/2  que se disponen paralelos a los lados del cuadrado.    
[Graphics:Images/DifusMetExpli2D.nb_gr_5.gif]
Introducimos las condiciones iniciales  y las  condiciones de contorno:  u(x=0,t)=u(x=1,t)=0

Empezamos asignando temperatura cero a todos los puntos (x_j,y_l)  en el instante inicial (t=0): u(x,y,0)=0.
  Recordad que x_m=m Δx,   e    y_l=l Δy.  Aquí tomamos jmax=lmax=16 como el valor máximo de j y l.  Podemos usar otro valor para jmax y lmax, pero hay que ser prudentes y no usar un número demasiado grande si no queremos hacer trabajar demasiado  al ordenador y no queremos aburrirnos esperando la salida de los resultados.

[Graphics:Images/DifusMetExpli2D.nb_gr_6.gif]

Asignamos temperatura u=1 para todo instante  t  a los puntos sobre una cruz centrada de tamaño L/2 por L/2

[Graphics:Images/DifusMetExpli2D.nb_gr_7.gif]

Asignamos temperatura u=-1 en el instante  t=0 a los puntos situados en el lado superior (j=0)

[Graphics:Images/DifusMetExpli2D.nb_gr_8.gif]

Asignamos, para cualquier instante posterior a t=0,  la temperatura u=-1 a los puntos situados en el borde superior (j=0) y temperatura 0 a  los  situados en los otros tres lados:

[Graphics:Images/DifusMetExpli2D.nb_gr_9.gif]

Comprobamos que las condiciones iniciales y de contorno son las deseadas.

[Graphics:Images/DifusMetExpli2D.nb_gr_10.gif]
[Graphics:Images/DifusMetExpli2D.nb_gr_11.gif]

Comprobamos que las condiciones de contorno son las deseadas. Para ello mostramos la solución en un instante tm  posterior a t=0.

[Graphics:Images/DifusMetExpli2D.nb_gr_12.gif]
[Graphics:Images/DifusMetExpli2D.nb_gr_13.gif]

Representamos las condiciones iniciales  en una gráfica tridimensional.

[Graphics:Images/DifusMetExpli2D.nb_gr_14.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_15.gif]

Fórmula de integración del método explícito.  La instrucción   u[j_,l_,m_]:=u[j,l,m]=...., hace que una vez evaluado u[j.l,m], este valor quede en memoria y no necesite ser evaluado de nuevo.  Esto hace que el  cálculo de u para m grandes se lleve a cabo con rapidez.

[Graphics:Images/DifusMetExpli2D.nb_gr_16.gif]

Vamos a integrar usando  s=1/4. Esto  garantiza la estabilidad del método de resolución.

[Graphics:Images/DifusMetExpli2D.nb_gr_17.gif]

Veamos la solución en un instante tm cualquiera

[Graphics:Images/DifusMetExpli2D.nb_gr_18.gif]
[Graphics:Images/DifusMetExpli2D.nb_gr_19.gif]

Una gráfica tridimensional con la  solución en el instante t=20Δt:

[Graphics:Images/DifusMetExpli2D.nb_gr_20.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_21.gif]

Veamos la solución mediante un código de colores en el instante t=tm  (tm=0,5,10,20,30)

[Graphics:Images/DifusMetExpli2D.nb_gr_22.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_23.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_24.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_25.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_26.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_27.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_28.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_29.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_30.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_31.gif]

Veamos la solución mediante un mapa de líneas  isotérmicas

[Graphics:Images/DifusMetExpli2D.nb_gr_32.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_33.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_34.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_35.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_36.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_37.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_38.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_39.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_40.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_41.gif]

[Graphics:Images/DifusMetExpli2D.nb_gr_42.gif]
Ejercicios

(1) Integra el sistema usando s=1/2. ¿Qué sucede?
(2) Resuelve el problema empleando otras condiciones iniciales y/o de contorno (las que más te motiven/gusten/interesen).  Puedes, por ejemplo, hacer que estas condiciones dependan del tiempo.