3/10/2024
Ecuaciones Diferenciales a
Derivadas Parciales
PARABÓLICAS
Trabajo Práctico N° 9
Fecha: 3 de Octubre 2024
Prof. Práctico: Coggiola, Vivian (3D1 y 3D3)
EDP PARABÓLICAS
𝜕𝑦 𝜕²𝑦
=𝑘
𝜕𝑡 𝜕𝑥²
∆x
∆t
La figura 1 muestra dos
cotas dada por las ∆x
condiciones de frontera.
1
3/10/2024
Condición de frontera
Condición de frontera
1
2
(Condiciones en la
frontera)
Acotada
espacialmente
Variación de T respecto al espacio x Pudiendo presentar condiciones de frontera de:
(Δx) ➢ Dirichlet: cuando las temperaturas en los
bordes son conocidas en los extremos.
➢ Neumann: no se conocen todas las
temperaturas en los extremos.
Ejemplo 1
Sea una barra delgada, conductora de calor, aislada en todo su largo (0,5 m de largo ) y sometida en el extremo
izquierdo a una temperatura de 120ºC y aislada en su extremo derecho.
Determine la distribución de temperaturas T(x,0.5)
120℃
𝜕2 𝑇 𝜕𝑇
k = k=0.025 ∆𝑥 = 0.1 ∆𝑡 = 0,1 𝑇 𝑥, 0.1 = 120 + 20 𝑠𝑒𝑛(𝑥)
𝜕𝑥 2 𝜕𝑡
¿Qué datos tengo para resolver este problema?
Distribución espacial ∆𝑥 = 0.1 0 < 𝑥 < 0,5
0,0 0,1 0,2 0,3 0,4 0,5
0,1
Distribución temporal 0,2
∆𝑡 = 0,1 0,1 < 𝑡 < 0,5 0,3
0,4
Sabemos que la barra 0,5
tiene un largo de 0,5m y
un ∆x= 0.1 4
2
3/10/2024
Ejemplo 1
Sea una barra delgada, conductora de calor, aislada en todo su largo (0,5 m de largo ) y sometida en el extremo
izquierdo a una temperatura de 120ºC y aislada en su extremo derecho. 120℃
Determine la distribución de temperaturas T(x,0.5)
𝜕2 𝑇 𝜕𝑇
k = k=0.025 ∆𝑥 = 0.1 ∆𝑡 = 0,1 𝑇 𝑥, 0.1 = 120 + 20 𝑠𝑒𝑛(𝑥) Condición inicial
𝜕𝑥 2 𝜕𝑡
0,0 0,1 0,2 0,3 0,4 0,5
0,1 120,0000 121,9967 123,9734 125,9104 127,7887 129,5885 CONDICION INICIAL
0,2
0,3
0,4
0,5
T(0.0,0.1) = 120 + 20.sen(0.0) = 120
T(0.1,0.1) = 120 + 20.sen(0.1) = 121,99666= 121,9967
Trabajamos con 4 cifras decimales con
T(0.2,0.1) = 120 + 20.sen(0.2) = 123,97338= 123,9734 redondeo simétrico (excepto que se
aclare lo contrario)
5
T(0.5,0.1) = 120 + 20.sen(0.5) = 129,5885
CALCULAMOS LAS CONDICIONES DE FRONTERA ¿Cuántas CONDICIONES?
Sabemos que la barra tiene un largo de 0,5m y un ∆x= 0.1
Extremo izquierdo a una temperatura de 120º
0,0 0,1 0,2 0,3 0,4 0,5
0,1 120,0000 121,9967 123,9734 125,9104 127,7887 129,5885
0,2 120,0000
0,3
0,4
120,0000
120,0000
?
0,5 120,0000
Condición de frontera 1 Condición de frontera 2
¿Qué sucede con la frontera del extremo derecho?
RESOLVER POR M. EXPLICITO
(DIF FINITAS DIVIDIDAS) 6
3
3/10/2024
VERIFICACION DE CONVERGENCIA
A partir del Siendo ∆t = 0,1 nos aseguramos que el problema, a lo
1 Δ𝑥 2 1 0,1 2
problema ∆𝑡 ≤ . = . = 0,2 largo del tiempo converge a una solución y se vuelve
2 k 2 0,025
sabemos que estable
∆x = 0.1
k= 0.025 ∆t debe ser ≤ 0,2
¿Podemos resolver un problema sin tener los valores de ∆x y ∆t ?
Fórmula explícita de cálculo para Δ𝑡 1 Podremos resolver un problema, si asignamos un valor a ∆x que
verificar que el problema es 𝝺= 𝑘. ≤ permita calcular el ∆t apropiado para que el problema sea
convergente y estable. Δ𝑥 ² 2 convergente .
7
Continuamos resolviendo
Sabemos que la barra tiene un largo de 0,5m y un ∆x= 0.1
Extremo izquierdo a una temperatura de 120º
0,1 120,0000 121,9967 123,9734 125,9104 127,7887 129,5885
0,2 120,0000
0,3 120,0000
0,4 120,0000 𝑻𝒊𝒋+𝟏 𝑻𝒊𝒋+𝟏
0,5 120,0000
Puntos sin aislante Puntos con aislante
Necesitamos plantear 2
ecuaciones
4
3/10/2024
𝜕2𝑇 𝜕𝑇
k = k=0.025 ∆𝑥 = 0.1 ∆𝑡 = 0,1 𝑇 𝑥, 0.1 = 120 + 20 𝑠𝑒𝑛(𝑥)
𝜕𝑥 2 𝜕𝑡
Diferencia finita
𝜕2 𝑇 𝑇𝑖+1,𝑗 −2𝑇𝑖,𝑗 +𝑇𝑖−1,𝑗 dividida
K =
𝜕𝑥 2 ∆𝑥 2 centrada
Diferencia finita
𝜕𝑇 𝑇𝑖,𝑗+1 −𝑇𝑖,𝑗 dividida HACIA
= ADELANTE
𝜕𝑡 ∆𝑡
∆𝑡
k (𝑇 −2𝑇𝑖𝑗 + 𝑇𝑖−1𝑗 ) = 𝑇𝑖𝑗+1 − 𝑇𝑖𝑗
∆𝑥 2 𝑖+1𝑗
∆𝑡 ∆𝑡 0,1 𝑇𝑖−1,𝑗 𝑇𝑖𝑗 𝑇𝑖+1,𝑗
𝑇𝑖𝑗+1 = k (𝑇 −2𝑇𝑖𝑗 + 𝑇𝑖−1𝑗 ) + 𝑇𝑖𝑗 reemplazando k = 0.025( ) =0,25
∆𝑥 2 𝑖+1𝑗 ∆𝑥 2 0,12
𝑇𝑖𝑗+1 = 0,25(𝑇𝑖+1𝑗 −2𝑇𝑖𝑗 + 𝑇𝑖−1𝑗 ) + 𝑇𝑖𝑗 𝑻𝒊,𝒋+𝟏
Predicción tiempo
𝑇𝑖𝑗+1 = 0,25𝑇𝑖+1𝑗 − 0,5𝑇𝑖𝑗 + 0,25𝑇𝑖−1𝑗 + 𝑇𝑖𝑗 futuro
𝑇𝑖𝑗+1 = 0,25 𝑇𝑖+1𝑗 + 0,5𝑇𝑖𝑗 + 0,25𝑇𝑖−1𝑗 Ecuación para 𝑻𝒊𝒋+𝟏 sin aislante (1) 9
Calculamos 𝑻𝒊𝒋+𝟏 = 𝟎, 𝟐𝟓 𝑻𝒊+𝟏𝒋 + 𝟎, 𝟓𝑻𝒊𝒋 + 𝟎, 𝟐𝟓𝑻𝒊−𝟏𝒋
0,0 0,1 0,2 0,3 0,4 0,5 0,0 0,1 0,2 0,3 0,4 0,5
0,1 120 121,9967 123,9734 125,9104 127,7884 129,5885 0,1 120 121,9967 123,9734 125,9104 127,7884 129,5885
0,2 120 121,9917 0,2 120 121,9917
0,3 120 0,3 120
0,4 120 0,4 120
0,5 120 0,5 120
𝑇𝑖𝑗+1 = 0,25(123,9734) + 0,5(121,9967) +0,25(120) 𝑇𝑖𝑗+1 = 0,25(125,9104) + 0,5(123,9734) +0,25(121,9967)
𝑻𝒊𝒋+𝟏 = 𝟏𝟐𝟏, 𝟗𝟗𝟏𝟕 𝑇𝑖𝑗+1 = 123,963475 = 123,963475
0,0 0,1 0,2 0,3 0,4 0,5
0,1 120 121,9967 123,9734 125,9104 127,7884 129,5885
0,2 120 121,9917 123,9635
0,3 120
0,4 120
0,5 120
𝑇𝑖𝑗+1 = 0,25(127,7884) + 0,5(125,9104) +0,25(123,9734)
𝑇𝑖𝑗+1 = 125,89565 = 125,8957 10
5
3/10/2024
Calculamos los puntos con aislante ¿Qué punto es conocido? 𝑇𝑖+1𝑗 ó 𝑇𝑖−1𝑗
𝜕𝑇 𝑇𝑖+1𝑗 −𝑇𝑖−1𝑗
=0 2∆𝑥
=0 por lo tanto 𝑇𝑖+1𝑗 = 𝑇𝑖−1𝑗
𝜕𝑡
0,0 0,1 0,2 0,3 0,4 0,5
0,1 120 121,9967 123,9734 125,9104 127,7884 129,5885
A partir de la Ecuación 1) 0,2 120 121,9917 123,9635 125,8957 127,7689
0,3 120
𝑇𝑖𝑗+1 = 0,25 𝑇𝑖+1𝑗 + 0,5𝑇𝑖𝑗 + 0,25𝑇𝑖−1𝑗 (1) 0,4 120
0,5 120
Tengo como punto conocido el 𝑇𝑖−1𝑗 entonces
𝑇𝑖𝑗+1 = 0,5(129,5885) + 0,5(127,7884) = 128,68845=128,6885
reemplazo 𝑇𝑖+1𝑗 𝑝𝑜𝑟 𝑇𝑖−1𝑗
0,0 0,1 0,2 0,3 0,4 0,5
0,1 120 121,9967 123,9734 125,9104 127,7884 129,5885
𝑇𝑖𝑗+1 = 0,25 𝑇𝑖−1𝑗 + 0,5𝑇𝑖𝑗 + 0,25𝑇𝑖−1𝑗 0,2 120 121,9917 123,9635 125,8957 127,7689 128,6885
0,3 120 121,9867 123,9536 125,8810 127,5305
𝑇𝑖𝑗+1 = 2 ∗ 0,25 𝑇𝑖−1𝑗 + 0,5𝑇𝑖𝑗 0,4 120
0,5 120
𝑇𝑖𝑗+1 = 0,5 𝑇𝑖𝑗 + 0,5𝑇𝑖−1𝑗 𝑇𝑖𝑗+1 = 0,5(128,6885) + 0,5(127,7689) = 128,2287
𝑇𝑖𝑗+1 = 0,5 𝑇𝑖𝑗 + 0,5𝑇𝑖−1𝑗 Ecuación para 𝑻𝒊𝒋+𝟏 CON aislante (2)
11
0,0 0,1 0,2 0,3 0,4 0,5
0,1 120,0000 121,9967 123,9734 125,9104 127,7884 129,5885
0,2 120,0000 121,9917 123,9635 125,8957 127,7689
0,3 120,0000
0,4 120,0000
0,5 120,0000
𝑇𝑖𝑗+1 = 0,25 𝑇𝑖+1𝑗 + 0,5𝑇𝑖𝑗 + 0,25𝑇𝑖−1𝑗 Ecuación para 𝑻𝒊𝒋+𝟏 sin aislante (1)
𝑇𝑖𝑗+1 = 0,25(123,9734) + 0,5(121,9967) + 0,25(120) = 121,9917
𝑇𝑖𝑗+1 = 0,25(125,9104) + 0,5(123,9734) + 0,25(121,9967) = 123,9635
𝑇𝑖𝑗+1 = 0,25(127,7884) + 0,5(125,9104) + 0,25(123,9734) = 125,8957
𝑇𝑖𝑗+1 = 0,25(129,5885) + 0,5(127,7884) + 0,25(125,9104) = 127,7689
12
6
3/10/2024
0,0 0,1 0,2 0,3 0,4 0,5
0,1 120,0000 121,9967 123,9734 125,9104 127,7884 129,5885
0,2 120,0000 121,9917 123,9635 125,8957 127,7689 128,6885
0,3 120,0000 121,9867
0,4 120,0000
0,5 120,0000
𝑇𝑖𝑗+1 = 0,5 𝑇𝑖𝑗 + 0,5𝑇𝑖−1𝑗 Ecuación para 𝑻𝒊𝒋+𝟏 𝒄𝒐𝒏 aislante (2)
𝑇𝑖𝑗+1 = 0,5 129,5885 + 0,5 127,7884 = 128,6885
𝑇𝑖𝑗+1 = 0,25 𝑇𝑖+1𝑗 + 0,5𝑇𝑖𝑗 + 0,25𝑇𝑖−1𝑗 Ecuación para 𝑻𝒊𝒋+𝟏 𝒔𝒊𝒏 aislante (1)
𝑇𝑖𝑗+1 = 0,25 123,9635 + 0,5 121,9917 + 0,25 120 = 121,9867
13
Resultado final
14
7
3/10/2024
Ejemplo 2
Un fenómeno físico de propagación de calor es modelado mediante la siguiente ecuación:
𝜕𝑇 1 𝜕2 𝑇
𝑥, 𝑡 - 𝑥, 𝑡 = 0 ∆𝑥 = 0,5 ∆𝑡 =0,5 0≤ t ≤ 4 1≤x ≤ 3
𝜕𝑡 5 𝜕𝑥 2
𝑇 1, 𝑡 = 1 𝑇 3, 𝑡 = 1 𝑇 𝑥, 0 = cos 𝑥² − 4𝑥 + 3
∆t
∆x
𝑻𝒊𝒋+𝟏
(Puntos sin aislante )
¿Condiciones de frontera?
¿Condición inicial?
𝑇 𝑥, 0 = cos 𝑥² − 4𝑥 + 3 𝑇 1,0 = cos 1² − 4(1) + 3 = cos 0 = 1,0000 𝑇 1, 𝑡 = 1
𝑇 3, 𝑡 = 1
𝑇 1.5,0 = cos 1.5² − 4(1.5) + 3 = 0,731688 = 𝟎, 𝟕𝟑𝟏𝟕
𝑇 2.0,0 = cos 2² − 4(2) + 3 = cos −1 = 0,54030 … = 𝟎, 𝟓𝟒𝟎𝟑
15
𝑇 2.5,0 = cos 2.5² − 4(2.5) + 3 = 𝟎, 𝟕𝟑𝟏𝟕
Resolvemos por el método explícito (D.F.D) y discretizamos
Un fenómeno físico de propagación de calor es modelado mediante la siguiente ecuación:
𝜕𝑇 1 𝜕2 𝑇
𝑥, 𝑡 - 𝑥, 𝑡 = 0 0≤ t ≤ 4 1≤x ≤ 3 𝑇 1, 𝑡 = 1 𝑇 3, 𝑡 = 1 𝑇 𝑥, 0 = cos 𝑥² − 4𝑥 + 3
𝜕𝑡 5 𝜕𝑥 2
∆𝑥 = 0,5 ∆𝑡 = 0,5
𝑇𝑖𝑗+1 − 𝑇𝑖𝑗 1 𝑇𝑖+1𝑗 − 2𝑇𝑖𝑗 + 𝑇𝑖−1𝑗
− =0
∆𝑡 5 ∆𝑥 2
𝑇𝑖𝑗+1 − 𝑇𝑖𝑗 1 𝑇𝑖+1𝑗 − 2𝑇𝑖𝑗 + 𝑇𝑖−1𝑗 1∆𝑡
= 𝑇𝑖𝑗+1 − 𝑇𝑖𝑗 = (𝑇 −2𝑇𝑖𝑗 + 𝑇𝑖−1𝑗 )
∆𝑡 5 ∆𝑥 2 5∆𝑥 2 𝑖+1𝑗
1 ∆𝑡 0,5
reemplazando =( ) =0,4
5 ∆𝑥 2 0,52
𝑇𝑖𝑗+1 − 𝑇𝑖𝑗 = 0,4(𝑇𝑖+1𝑗 −2𝑇𝑖𝑗 + 𝑇𝑖−1𝑗 )
𝑇𝑖𝑗+1 = 0,4𝑇𝑖+1𝑗 − 0,8𝑇𝑖𝑗 + 0,4𝑇𝑖−1𝑗 + 𝑇𝑖𝑗
𝑻𝒊𝒋+𝟏 = 𝟎, 𝟒𝑻𝒊+𝟏𝒋 + 𝟎, 𝟐𝑻𝒊𝒋 + 𝟎, 𝟒𝑻𝒊−𝟏𝒋 Ecuación para 𝑻𝒊𝒋+𝟏 sin aislante (1) 16
8
3/10/2024
0≤ t ≤ 4 1≤x ≤ 3 𝑇 1, 𝑡 = 1 𝑇 3, 𝑡 = 1 𝑇 𝑥, 0 = cos 𝑥² − 4𝑥 + 3
1,0 1,5 2,0 2,5 3,0 1,0 1,5 2,0 2,5 3,0
3,0,0 1,0000 0,7317 0,5403 0,7317 1,0000 3,0,0 1,0000 0,7317 0,5403 0,7317 1,0000
0,5 1,0000 0,7625 1,0000 0,5 1,0000 0,7625 0,6934 1,0000
1,0 1,0000 1,0000 1,0 1,0000 1,0000
1,5 1,0000 1,0000 1,5 1,0000 1,0000
2,0 1,0000 1,0000 2,0 1,0000 1,0000
2,5 1,0000 1,0000 2,5 1,0000 1,0000
3,0 1,0000 1,0000 3,0 1,0000 1,0000
3,5 1,0000 1,0000 3,5 1,0000 1,0000
4,0 1,0000 1,0000 4,0 1,0000 1,0000
𝑇𝑖𝑗+1 = 0,4𝑇𝑖+1𝑗 + 0,2𝑇𝑖𝑗 + 0,4𝑇𝑖−1𝑗
𝑇𝑖𝑗+1 = 0,4(0,5403) + 0,2(0,7317) +0,4(1,0000) = 0,76246…=0,7625
𝑇𝑖𝑗+1 = 0,4(0,7317) + 0,2(0,5403) +0,4(0,7317) = 0,69342…=0,6934
17
1,0 1,5 2,0 2,5 3,0
3,0,0 1,0000 0,7317 0,5403 0,7317 1,0000
0,5 1,0000 0,7625 0,6934 0,7625 1,0000
1,0 1,0000 0,8299 0,7487 0,8299 1,0000
1,5 1,0000 0,8655 𝑻𝒊𝒋+𝟏 0,8655 1,0000
2,0 1,0000 0,8986 1,0000
2,5 1,0000 1,0000
¿Se puede calcular el punto
𝑇𝑖𝑗+1 indicado en la Tabla ? 3,0 1,0000 1,0000
3,5 1,0000 1,0000
4,0 1,0000 1,0000
1,0 1,5 2,0 2,5 3,0
3,0,0 1,0000 0,7317 0,5403 0,7317 1,0000
0,5 1,0000 0,7625 0,6934 0,7625 1,0000
1,0 1,0000 0,8299 0,7487 0,8299 1,0000
1,5 1,0000 0,8655 𝑻𝒊𝒋+𝟏 0,8655 1,0000
𝑇𝑖𝑗+1 = 0,4𝑇𝑖+1𝑗 + 0,2𝑇𝑖𝑗 + 0,4𝑇𝑖−1𝑗
2,0 1,0000 0,8986 1,0000
2,5 1,0000 1,0000
𝑇𝑖𝑗+1 = 0,4(0,8299) + 0,2(0,7487) +0,4(0,8299)
3,0 1,0000 1,0000
𝑇𝑖𝑗+1 = 0,81366…=0,8137
3,5 1,0000 1,0000
4,0 1,0000 1,0000
18
9
3/10/2024
Resultado final
t/x 1,0 1,5 2,0 2,5 3,0
0,00 1,0000 0,7317 0,5403 0,7317 1,0000
0,50 1,0000 0,7625 0,6934 0,7625 1,0000
1,00 1,0000 0,8299 0,7487 0,8299 1,0000
1,50 1,0000 0,8655 0,8137 0,8655 1,0000
2,00 1,0000 0,8986 0,8551 0,8986 1,0000
2,50 1,0000 0,9218 0,8899 0,9218 1,0000
3,00 1,0000 0,9403 0,9154 0,9403 1,0000
3,50 1,0000 0,9542 0,9353 0,9542 1,0000
4,00 1,0000 0,9650 0,9504 0,9650 1,0000
19
Ejemplo 3
Un fenómeno físico de propagación de calor es modelado mediante la siguiente ecuación:
𝜕𝑇 1 𝜕2 𝑇
𝑥, 𝑡 - 𝑥, 𝑡 = 0 ∆𝑥 = 0,5 ∆𝑡 = 0,25 0≤ t ≤ 1,5 1≤x ≤ 3 𝑇 1, 𝑡 = 1 𝑇 3, 𝑡 = cos 𝑡
𝜕𝑡 4 𝜕𝑥 2
𝑇 𝑥, 0 = cos 𝑥² − 4𝑥 + 3
∆x
t/x 1,0 1,5 2,0 2,5 3,0
0,00 1,0000 0,7317 0,5403 0,7317 1,0000 cos 0,0 = 1,0000
0,25 1,0000 0,9689 cos 0,25 =0,9689
0,50 1,0000 0,8776 cos 0,5 =0,9689
∆t 0,75 1,0000 0,7317 cos 0,75 =0,7317
1,00 1,0000 0,5403 cos 1,0 =0,5403
1,25 1,0000 0,3153 cos 1,25 =0,3153
1,50 1,0000 0,0707 cos 1,5 =0,0707
¿Condición inicial?
¿Condiciones de frontera?
𝑇 𝑥, 0 = cos 𝑥² − 4𝑥 + 3
𝑇 1, 𝑡 = 1
𝑇 1,0 = cos 1² − 4(1) + 3 = 𝑐𝑜𝑠 0 = 1
𝑇 3, 𝑡 = cos 𝑡
𝑇 1.5,0 = cos 1.5² − 4(1.5) + 3 = 0,731688 = 0,7317
20
10
3/10/2024
Resolvemos por Método Explícito (DFD)
Un fenómeno físico de propagación de calor es modelado mediante la siguiente ecuación:
𝜕𝑇 1 𝜕2 𝑇
𝑥, 𝑡 - 𝑥, 𝑡 = 0 ∆𝑥 = 0,5 ∆𝑡 = 0,25 0≤ t ≤ 1,5 1≤x ≤ 3 𝑇 1, 𝑡 = 1 𝑇 3, 𝑡 = cos 𝑡
𝜕𝑡 4 𝜕𝑥 2
𝑇 𝑥, 0 = cos 𝑥² − 4𝑥 + 3
𝑇𝑖𝑗+1 − 𝑇𝑖𝑗 1 𝑇𝑖+1𝑗 − 2𝑇𝑖𝑗 + 𝑇𝑖−1𝑗
− =0
∆𝑡 4 ∆𝑥 2
𝑇𝑖𝑗+1 − 𝑇𝑖𝑗 1 𝑇𝑖+1𝑗 − 2𝑇𝑖𝑗 + 𝑇𝑖−1𝑗 1∆𝑡
𝑇𝑖𝑗+1 − 𝑇𝑖𝑗 = (𝑇 −2𝑇𝑖𝑗 + 𝑇𝑖−1𝑗 )
= 4∆𝑥 2 𝑖+1𝑗
∆𝑡 4 ∆𝑥 2
1 ∆𝑡 0,25
reemplazando =1/4( ) =0,25
4 ∆𝑥 2 0,52
𝑇𝑖𝑗+1 − 𝑇𝑖𝑗 = 0,25(𝑇𝑖+1𝑗 −2𝑇𝑖𝑗 + 𝑇𝑖−1𝑗 )
𝑇𝑖𝑗+1 = 0,25𝑇𝑖+1𝑗 − 0,5𝑇𝑖𝑗 + 0,25𝑇𝑖−1𝑗 + 𝑇𝑖𝑗
𝑻𝒊𝒋+𝟏 = 𝟎, 𝟐𝟓𝑻𝒊+𝟏𝒋 + 𝟎, 𝟓𝑻𝒊𝒋 + 𝟎, 𝟐𝟓𝑻𝒊−𝟏𝒋 Ecuación para 𝑇𝑖𝑗+1 sin aislante (1) 21
Resultado final
t/x 1,0 1,5 2,0 2,5 3,0
0,00 1,0000 0,7317 0,5403 0,7317 1,0000
0,25 1,0000 0,7509 0,6360 0,7509 0,9689
0,50 1,0000 0,7845 0,6935 0,7767 0,8776
0,75 1,0000 0,8156 0,7371 0,7811 0,7317
1,00 1,0000 0,8421 0,7677 0,7577 0,5403
1,25 1,0000 0,8630 0,7838 0,7059 0,3153
1,50 1,0000 0,8775 0,7841 0,6277 0,0707
𝑻𝒊𝒋+𝟏 = 𝟎, 𝟐𝟓𝑻𝒊+𝟏𝒋 + 𝟎, 𝟓𝑻𝒊𝒋 + 𝟎, 𝟐𝟓𝑻𝒊−𝟏𝒋
𝑻𝒊𝒋+𝟏 = 𝟎, 𝟐𝟓 𝟎, 𝟓𝟒𝟎𝟑 + 𝟎, 𝟓 𝟎, 𝟕𝟑𝟏𝟕 + 𝟎, 𝟐𝟓 𝟏, 𝟎𝟎𝟎𝟎 = 𝟎, 𝟕𝟓𝟎𝟗
𝑻𝒊𝒋+𝟏 = 𝟎, 𝟐𝟓 𝟎, 𝟕𝟑𝟏𝟕 + 𝟎, 𝟓 𝟎, 𝟓𝟒𝟎𝟑 + 𝟎, 𝟐𝟓 𝟎, 𝟕𝟑𝟏𝟕 = 𝟎, 𝟔𝟑𝟔𝟎
𝑻𝒊𝒋+𝟏 = 𝟎, 𝟐𝟓 𝟎, 𝟕𝟑𝟏𝟕 + 𝟎, 𝟓 𝟎, 𝟓𝟒𝟎𝟑 + 𝟎, 𝟐𝟓 𝟎, 𝟕𝟑𝟏𝟕 = 𝟎, 𝟔𝟑𝟔𝟎
22
11
3/10/2024
Ejemplo 4
1. Sea una barra delgada, conductora de calor, aislada en todo su largo (0,5 m de largo ) y sometida en el extremo
izquierdo a una temperatura de 150ºC y aislada en su extremo derecho.
𝜕2 𝑇 𝜕𝑇
k = k=0.025 ∆𝑥 = 0.25 ∆𝑡 = 0,1 𝑇 𝑥, 0.1 = 150 + 50 𝑠𝑒𝑛(𝑥)
𝜕𝑥 2 𝜕𝑡
a. Construya la tabla e indique las condiciones iniciales y de frontera.
b. Obtenga la/s fórmula/s de cálculo por discretización con DF.
c. Determine la distribución de temperaturas T(x,0,3).
RESULTADO…
t/x 0,00 0,25 0,50
0,10 150,0000 162,3702 173,9713
0,20 150,0000 162,3394 173,0432
0,30 150,0000 162,2740 172,1869
23
24
12
