ECUACIONES

DIFERENCIALES
PARCIALES

UMRB-2020
 ECUACIONES DIFERENCIALES PARCIALES A PARTIR DE LA
MODELACIÓN DE FENÓMENOS FÍSICOS
Ecuación de Calor
Una de las ecuaciones diferenciales parciales mas estudiadas es:

Qy

Qx
α = Coeficiente de difusividad térmica
k = coeficiente de conductividad térmica
p = densidad del material
Cp = Capacidad calorífica del material
T = temperatura
t= tiempo
Qz

Esta ecuación describe la conducción de calor a través de un material en

régimen transitorio en sus 3 dimensión.
Ecuación de difusión

Cy

Cx

D = Coeficiente de difusividad
C = concentración

Cz
Fluido

Representa la difusión molecular de masa entre fluidos en sus 3 dimensiones.

NOTA: Repasar el balance de donde se obtiene las ecuaciones

 Distribución en estado estacionario que pueden
caracterizarse por EDP elípticas. a) Distribución de
temperatura sobre una placa caliente.

a) Barra larga y delgada que está aislada, excepto en sus

extremos. La dinámica de la distribución unidimensional de
temperatura a lo largo de la barra puede describirse mediante
una EDP parabólica. b) La solución, que consiste en
distribuciones correspondientes al estado de la barra en
diferentes momentos.
DIFERENCIAS FINITAS:
ECUACIONES
PARABÓLICAS
• MÉTODO EXPLICITO
• MÉTODO IMPLÍCITO
• MÉTODO CRANK-NICOLSON
 • MÉTODO
EXPLICITO
Si se sustituyen las derivadas parciales por sus aproximaciones en
diferencias finitas:

Ecuación de
Calor

Diferencias finitas hacia adelante

Diferencias finitas centrales

Remplazando en la ecuación de calor tenemos la siguiente ecuación
algebraica:
Conducción de calor en una barra
aislada:
Consideremos una barra aislada longitudinalmente y de longitud L muy
grande en comparación con su grosor, la longitud de la barra se medirá con
la variable x, siendo en sus extremos x = 0 y x = L.
Sean además la siguientes condiciones:

a) T(x,0) = f(x) 0<x<L Condición Inicial (En t = 0)

b) T(0,t) = g1(t) t>0 Condición de frontera 1 (En x = 0)

c) T(L,t) = g2(t) t > 0 Condición de frontera 2 (En x = L)

 T(x1,t1) T(x2,t1) T(xn,t1)
T(x1,t2) T(x2,t2) T(xn,t2)

T(x1,tmax) T(x2,tmax) T(xn,tmax)

Resolver este problema numéricamente, significa encontrar los valores

de T en puntos seleccionados en la barra x1, x2 , … , xn a ciertos tiempos
escogidos t1 < t2 …. <tmax
La forma de analizar el problema es gráficamente, representando el dominio
de definición T como el rectángulo que se ilustra en el sistema coordenado
x-t de la figura inferior.

T(x2,t2)
Nodos:
Representando los nodos en una red construida en el dominio de definición
Ejercicio 1
:Resolvemos el problema planteado, con una temperatura a tiempo
cero de -10°C, y después con una temperatura 40°C:
Solución:
Primero: Construimos la malla en el dominio de definición, para ello
subdividimos el tiempo y la longitud en varios subintervalos.
L = 0,4 m, en 4 subintervalos de 0,1 m
T = 1 h, en 10 subintervalos de 0,1 h
t
j

CF2 T(0,4,t) = 40ºC

CF1 T(0,t) = 40ºC
3

T(i,j+1)
2

T(i,j) T(i+1,j)
1
i

0 x
0 1 2 3 4
CI T(x,0) = -10ºC
Valores de T en cada t
nodo j

40ºC 40ºC

40ºC 40ºC

CF2 T(0,4,t) = 40ºC

CF1 T(0,t) = 40ºC
T(i,j+1)
40ºC 40ºC

T(i,j) T(i+1,j)
40ºC 40ºC
T(i-1,j) i

40ºC x
Nodo Singular -10ºC -10ºC -10ºC -10ºC -10ºC
T(40+(-10))/2 = 15ºC CI T(x,0) = -10ºC Nodo Singular
T = 15ºC
Usando la ecuación hallada anteriormente:

t
j
40ºC

40ºC

40ºC
Para el nodo i = 1, j = 0:
40ºC

i
40º
C -10ºC -10ºC -10ºC -10ºC -10ºC
x
 t
j
Para el nodo i = 2, j = 0: 40ºC

40ºC

40ºC

40ºC

i
40ºC
-10ºC -10ºC -10ºC -10ºC -10ºC x

t
j
40ºC
Para el nodo i = 3, j = 0:
40ºC

40ºC

40ºC

i
40ºC
-10ºC -10ºC -10ºC -10ºC -10ºC x
De este modo calculamost los valores de T de la segunda fila de
nodos. j

40ºC
40ºC

40ºC
40ºC

CF2 T(0,4,t) = 40ºC

CF1 T(0,t) = 40ºC

40ºC 40ºC

15ºC -10ºC 15ºC

40ºC 40ºC
i

x
15ºC -10ºC -10ºC -10ºC 15ºC
CI T(x,0) = -10ºC
Ahora deberíamos calcular los valores de T para la tercera fila y luego para
la cuarta y así hasta terminar las 10 filas en que subdividimos la hora.

Para el nodo i = 1, j = 1:
• MÉTODO
IMPLICITO
Ecuación de
Calor

Diferencias finitas centrales

Sustituyendo

donde
Conducción de calor en una barra aislada
Ejercicio 2
:
Resolvemos el problema planteado en la sección anterior:
 t
j

40ºC 40ºC

40ºC 40ºC

CF2 T(0,4,t) = 40ºC

CF1 T(0,t) = 40ºC 40ºC 40ºC

T(i-1,j) T(i,j) T(i+1,j)

40ºC 40ºC
i

40ºC
T(i,j-1)
x
Nodo Singular -10ºC -10ºC -10ºC -10ºC -10ºC
T(40+(-10))/2 = 15ºC CI T(x,0) = -10ºC Nodo Singular
T = 15ºC
 t
j
40ºC 40ºC

40ºC 40ºC

Para el nodo i=1, j=1 40ºC 40ºC

40ºC 40ºC

i
40ºC
40ºC

t -10ºC -10ºC -10ºC -10ºC -10ºC x

j
40ºC 40ºC

Para el nodo i=2, j=1 40ºC 40ºC

40ºC 40ºC

40ºC 40ºC

i
40ºC
40ºC
-10ºC -10ºC -10ºC -10ºC -10ºC x
 t
Para el nodo i=3, j=1 j
40ºC 40ºC

40ºC 40ºC

40ºC 40ºC

40ºC 40ºC

i
40ºC
40ºC
-10ºC -10ºC -10ºC -10ºC -10ºC x

El sistema formado será el siguiente:

Remplazando valores:

Resolviendo por cualquier método conocido:

• MÉTODO
CRANK-NICOLSON
El método de Crank-Nicolson ofrece un esquema implícito alternativo que tiene una exactitud de
segundo orden, tanto para el espacio como para el tiempo

Entonces, la primera derivada temporal se aproxima en tl+1/2 por

La segunda derivada en el espacio puede determinarse en el punto medio promediando las

aproximaciones por diferencias al principio (tl) y al final (tl+1) del incremento del tiempo

Sustituyendo las ecuaciones y en la ecuación y reagrupando términos, se obtiene

Ecuación de Calor en 2 dimensiones
Consideremos ahora una lamina delgada de espesor despreciable. La ecuación
que describe su comportamiento es la siguiente:

Esta ecuación permite hallar el valor de la Temperatura en cualquier posición

de la lamina y en cualquier tiempo en régimen transitorio.
 t x-t
t y -
j

y
CF2 T(x,Ly,t)
k

3 T(0 , y ,t)
x-y
CF t)
CI T(x,y,0)
x,y,
T(L
4
CF
i
x
CF1 T(x,0,t)
 t x-t
t y -

t2 y
j
T(i,j,k+1)

x-y
T(i,j+1,k)
t1
T(i-1,j,k)
T(i+1,j,k)
T(i,j,k)
to
x
T(i,j-1,k) i
Y luego remplazamos cada termino en la ecuación de calor bidimensional en
régimen transitorio
De donde la temperatura que no conocemos el la T(i,j,k+1)

Para hallar la temperatura T(i,j,k+1) tenemos: i=1, j=1, k=0

Se debe formar el sistema de ecuciones con cada i y j

FIN
Suerte para el siguiente parcial