Ecuaciones diferenciales en derivadas parciales: la ecuación de Difusión

Ecuación de Continuidad Ley de Fick

κ es el coeficiente de difusión
⇥
+ ·J = 0 J J= r⇢
⇥t

(r, t)

⇤⇥
= r2 ⇥
⇤t

ρ representa densidad, calor (distribución de temperaturas), …

✓ 2 ◆
⇤⇥ ⇤ ⇥ ⇤2⇥ ⇤2⇥
= + 2+ 2
⇤t ⇤x2 ⇤y ⇤z

Caso unidimensional
⇤⇥ ⇤2⇥
=
⇤t ⇤x2
Métodos Computacionales
 Ecuación de Difusión unidimensional

Estabilidad: Los errores no crecen al avanzar en el tiempo.

Consistencia: Cuando los intervalos finitos h y τ tienden a cero la

ecuación de diferencias finitas utilizada debe tender
a la ecuación diferencial exacta.

Convergencia: Cuando los intervalos finitos h y τ tienden a cero la

solución obtenida numéricamente tiende a la solución
exacta.

Estabilidad y Consistencia = Convergencia

Métodos Computacionales 2
 Ecuación de Difusión unidimensional

Resolución numérica por el método explícito

⇥T ⇥2T
T(x,t) =
⇥t ⇥x2

xi = ih h, intervalo espacial
Discretización (diferencias finitas) Ti,n = T (xi , tn )
tn = n τ, intervalo temporal

Ti,n+1 Ti,n Ti+1,n 2Ti,n + Ti 1,n

=
⇥ h2
⇥
Ti,n+1 = Ti,n + {Ti+1,n 2Ti,n + Ti 1,n }
h2

⇥ 1
r⌘ (< ) coeficiente adimensional (dependiente de la discretización)
h2 2

Métodos Computacionales 3
 Ecuación de Difusión unidimensional

Resolución numérica por el método implícito: Crank Nicolson

⇥T ⇥2T
=
⇥t ⇥x2
⇢
Ti,n+1 Ti,n Ti+1,n 2Ti,n + Ti 1,n Ti+1,n+1 2Ti,n+1 + Ti 1,n+1
= +
⇥ 2 h2 h2

rTi 1,n+1 + 2(1 + r)Ti,n+1 rTi+1,n+1 = rTi 1,n + 2(1 r)Ti,n + rTi+1,n
⇥
r⌘ coeficiente adimensional (dependiente de la discretización)
h2
sin límite de estabilidad

Unidades de longitud y tiempo tales que L=1 y κ=1

T(0,t) = A(t) T(L,t) = B(t)

Condiciones de contorno de Dirichlet
Métodos Computacionales 4
 Ecuación de Difusión unidimensional

Resolución numérica por el método implícito: Crank Nicolson

rTi 1,n+1 + 2(1 + r)Ti,n+1 rTi+1,n+1 = rTi 1,n + 2(1 r)Ti,n + rTi+1,n

0 10 ..
1
.. .. .. .
B . . . CB C
B CB C
B CB C
B CB C
B r 2(1+r) r CB Ti 1,n+1 C
B CBB
C
C
B CB C
B CB C
B CB C
B r 2(1+r) r C B Ti,n+1 C
B CB C
B CB C
B CB C
B CB C
B CB T C
B r 2(1+r) r C B i+1,n+1 C
B CB C
B CB C
@ A@ A
.. .. .. ..
. . .
.

Métodos Computacionales 5
 Ecuación de Difusión unidimensional

Resolución numérica por el método implícito: Crank Nicolson

rTi 1,n+1 + 2(1 + r)Ti,n+1 rTi+1,n+1 = rTi 1,n + 2(1 r)Ti,n + rTi+1,n

0 10 ..
1 0
..
1
.. .. .. . .
B . . . CB C B C
B CB C B C
B CB C B C
B CB C B C
B r 2(1 r) r CB Ti 1,n C B bi 1 C
B CB C B C
B CB C B C
B CB C B C
B CB C B C
B r 2(1 r) r CB Ti,n
C B
C ⌘ B bi
C
B CBB C B
C
C
B CB C B C
B CB C B C
B CB C B C
B CB T C B bi+1 C
B r 2(1 r) r C B i+1,n C B C
B CB C B C
B CB C B C
@ A@ A @ A
.. .. .. .. ..
. . .
. .

Métodos Computacionales 6
 Ecuación de Difusión unidimensional

Resolución numérica por el método implícito: Crank Nicolson

rTi 1,n+1 + 2(1 + r)Ti,n+1 rTi+1,n+1 = rTi 1,n + 2(1 r)Ti,n + rTi+1,n

0 10 ..
1 0
..
1
.. .. .. . .
B . . . CB C B C
B CB C B C
B CB C B C
B CB C B C
B r 2(1+r) r CB Ti 1,n+1 C B bi 1 C
B CBB
C B
C B
C
C
B CB C B C
B CB C B C
B CB C B C
B r 2(1+r) r C B Ti,n+1 C = B bi C
B CB C B C
B CB C B C
B CB C B C
B CB C B C
B CB T C B bi+1 C
B r 2(1+r) r C B i+1,n+1 C B C
B CB C B C
B CB C B C
@ A@ A @ A
.. .. .. .. ..
. . .
. .

{Ti,n } {Ti,n+1 }

Métodos Computacionales 7
 Ecuación de Difusión unidimensional. Los extremos

Resolución numérica por el método implícito: Crank Nicolson

rTi 1,n+1 + 2(1 + r)Ti,n+1 rTi+1,n+1 = rTi 1,n + 2(1 r)Ti,n + rTi+1,n
T0,n+1 = T0,n

TN,n+1 = TN,n

rT0,n+1 + 2(1 + r)T1,n+1 rT2,n+1 = rT0,n + 2(1 r)T1,n + rT2,n

rTN 2,n+1 + 2(1 + r)TN 1,n+1 rTN,n+1 = rTN 2,n + 2(1 r)TN 1,n + rTN,n
0 1
b1 + 2rT0,n
B C
B .. C
B . C
B C
B C
B bi C
Renormalización de los términos independientes B 1 C
B C
B bi C
B C
B C
B bi+1 C
B C
B C
B .. C
B . C
@ A
Métodos Computacionales bN 1 + 2rTN,n 8