“AÑO DEL BICENTENARIO DEL PERÚ:200 AÑOS DE

INDEPENDENCIA”

UNIVERSIDAD NACIONAL DE SAN MARTÍN

FACULTAD DE INGENIERÍA CIVIL Y ARQUITECTURA

ESCUELA PROFESIONAL DE INGENIERÍA CIVIL

Tema : Ecuaciones Diferenciales Parciales

Asignatura : Métodos Numéricos

Estudiante : Piero Sebastián Paredes Hidalgo - 72257252

Docente : Juan Carlos Ipanaqué Palacios

Ciclo : 2021-II

Tarapoto, San Martín-Perú

2021
DESARROLLO
1 .−¿Con el método explícito calcule la distribución de temperaturaenuna barra larga y¿
delgada que tiene una longitud de 10 cm y los siguientes valores :
cal
k =0.49 , ∆ x=2 cm y ∆ t=0.1 s . La temperatura
( s . cm .° C )
de los extremos de labarra se mantienen en 40 ° C y 60° C , y la inicial es de 100° C .
cal g
Considere que la barra es de aluminio con C=0.2174 y ρ=2.7 3 .
(g . °C ) cm

Solución
→ Ecuaciones parabólicas( Diferencias finitas)
2
∂ T ∂T
Ecuación de la conducción de calor :k = …(1)
∂x
2
∂t
La segunda derivada se representa mediante uan diferencia dividida finita centrada:
l l l
∂ T T −2T i +T i−1
2
k 2 = i +1 2
…(2)
∂x ∆x

Una diferencia dividida finitahacia adelante sirve para aproximar laderivada con
T l +1+ T li
respecto altiempo : ∂T = i …( 3 )
∂t ∆t

T li +1−2 T li +T li−1
T l+1 l
i +T i
Sustituyendo las ecuaciones : ( 2 ) y ( 3 ) en ( 1 )k =
∆ x2 ∆t
k ∆ t Datos
i =T i + λ ( T i+1 −2T i +T i −1 ) → donde : λ=
Que resulta:T l+1 l l l l
¿¿
T ( 0 , t )=40° C → Nodo T 0T ( 10 , t )=60 ° C → NodoT 5∆ t=0.1T 0=100 ° C ∆ x=2cm

k=
( 0.49
( s . cm .° C ) )
cal

=0.8348
cm
λ=
k∆t
2

( ( g .° C ) )( cm )
cal g s ¿¿
0.2174 2.7 3

Longitud de barra=10 cm→ con ∆ x=2 cm

∴ Existirá 4 nodos entre 0 cmhasta 10 cm→ 0 cm−2 cm−4 cm−6 cm−8 cm−10 cm
Cumplimiento de convergencia y estabilidad del ejercicio :
2
1 ∆x
∆t ≤ →0.1 ≤ 1.1978
2 k
i =T i + λ ( T i+1 −2T i +T i −1 ) → Dondel=0
Calcularemos laTempratura de cada nodo T l+1 l l l l

En t=0.1 sT 11=T 10+ λ ( T 02−2 T 10+T 00 ) =100+0.02087 ( 100−2 ( 100 ) + 40 )=98.7478° C

T 2=T 2+ λ ( T 3−2 T 2 +T 1 ) =100+0.02087 ( 100−2 ( 100 ) +100 ) =100 °C
1 0 0 0 0

T 13=T 03+ λ ( T 04 −2T 03 +T 02 )=100+ 0.02087 ( 100−2 ( 100 ) +100 )=100 ° C

T 14=T 04 + λ ( T 05−2 T 04 +T 30) =100+0.02087 ( 60−2 ( 100 ) +100 )=99.1652° C

En t=0.2 s
T 21=T 11+ λ ( T 12−2T 11 +T 10) =98.7478+ 0.02087 ( 100−2 ( 98.7478 ) + 40 )=97.5478° C
T 2=T 2+ λ ( T 3−2T 2 +T 1) =100+0.02087 ( 100−2 (100 )+ 98.7478 )=99.9738° C
2 1 1 1 1

T 23=T 13+ λ ( T 14−2T 13+ T 12 ) =100+0.02087 ( 99.1652−2 ( 100 ) +100 )=99.9825 ° C

T 4=T 4 + λ ( T 5−2 T 4+ T 3 ) =99.1652+0.02087 ( 60−2 ( 99.1652 )+100 )=98.3652 °C
2 1 1 1 1

Así localcularemos por variaciones de 0.1 s por cada intervalo

Nodos
Intervalos(l) ∆ t ( s)
0(0cm) 1(2cm) 2(4cm) 3(6cm) 4(8cm) 5(10cm)
0 0 40 100 100 100 100 60
1 0.1 40 98.7478 100 100 99.1652 60
2 0.2 40 97.5478 99.9738 99.9825 98.3652 60
3 0.3 40 96.3974 99.9233 99.9485 97.5982 60
4 0.4 40 95.2939 99.8502 99.8989 96.8625 60
5 0.5 40 94.2350 99.7561 99.8345 96.1565 60

2 .−¿Consideramos la ecuación, aproxime lassoluciones dela siguiente ecuación :¿

2 2
∂ u( ∂ u
2
x , y )+ 2 ( x , y ) =0 ;0< x <0,5 , 0< y <0,5
∂x ∂y
Conlas condiciones de frontera :
π
u ( 0 , y )=0 , u ( 0,5 , y )=200 y ,0 ≤ y ≤
2
u ( x , 0 ) =0 , u ( x , 0,5 )=200 x ,0 ≤ x ≤ π
1 1
Use h= y k=
8 8
Solución→ Ecuacioneselípticas ( Diferencias finitas ) .

∂2 u ∂2 u u −2u i. j +ui−1. j ui . j+1−2 ui . j +ui . j−1

Fórmula : 2
( x , y ) + 2
( x , y ) = i +1. j + =0
∂x ∂y ∆ x2 ∆ y2
ui+1. j−2 ui . j+ ui−1. j u i. j +1−2u i . j +u i. j−1
→ Como ∆ x=∆ y reemplazamos 2
+ 2
=0
∆x ∆x
ui +1. j +ui−1. j +ui . j+1 +ui . j−1−4 u i . j =0

Calculamos las variaciones de x e y : x i=a+ (i ) h=¿x 0=0+ ( 0 ) ( 0.125 )=0

x 1=0+ ( 1 ) ( 0.125 )=0.125x 2=0+ ( 2 ) ( 0.125 )=0.25x 3=0+ (3 )( 0.125 )=0.375
x 4 =0+ ( 4 ) ( 0.125 )=0.5

y i=c+ ( i ) k=¿ y 0=0+ ( 0 ) ( 0.125 )=0 y 1=0+ ( 1 )( 0.125 ) =0.125 y 2=0+ ( 2 )( 0.125 )=0.25
y 3=0+ ( 3 ) ( 0.125 )=0.375 y 4 =0+ ( 4 ) ( 0.125 )=0.5

Colocamos las fronteras :

 0.5 200x 200x 200x
0.375 0 u3 u6 u9 200y
0.250 0 u2 u5 u8 200y
0.125 0 u1 u4 u7 200y
0.000 0 0 0
0.000 0.125 0.250 0.375 0.5

Partimos de :→u i+1. j +u i−1. j +ui . j+1 +ui . j−1=4 ui . j4 u1=u 4 +0+u2 +0 → 4 u1 −u4 −u2=0
4 u2=u5 +0+ u3 +u1 → 4 u2−u5−u3−u1=04 u3=u6 + 0+200 ( 0.125 ) +u2 → 4 u3−u6 −u2=25
4 u4 =u7 +u1 +u5 +0 → 4 u 4−u7 −u1−u 5=04 u5=u 8+u 2+u 4 +u6 → 4 u5 −u8−u 2−u4 −u6=0
4 u6=u 9+ u3+ 200 ( 0.125 )+ u5 → 4 u6−u 9−u3−u 5=25
4 u7=200 ( 0.125 ) +u 4 +u8 +0 → 4 u 7−u4 −u8 =25
4 u8=200 ( 0.250 )+u 5+u 7 +u9 → 4 u8−u 5−u7 −u9=50
4 u9=200 ( 0.375 ) +u 6+ u8 +200 ( 0.375 ) → 4 u9−u6 −u8=150

25 25 75
u1 = u2= u3=
4 2 4
25 75 75 75 225
u 4= u5=25u6 = u7 = u8 = u9 =
2 2 4 2 4

3 .−¿ Determine ladistribución del calor en estado estable en una superficie metálica de ¿
dimensiones 0,5 m x 0,5 mutilizando una grilla de tamaño m=n=4 .
Las condiciones de borde en los ejes e son 0 y el calor enlos otros dos bordes crece
linealmente desde 0 a 100 grados .

Solución→ Ecuacioneselípticas ( Diferencias finitas ) .

ui +1. j +ui−1. j +ui . j+1 +ui . j−1−4 u i . j =0

Graficaremos la distribución de calor :→ Untamaño de grilla ointervalos :m=n=4

0.5 0 25 50 75 100
0.375 0 u3 u6 u9 75
0.250 0 u2 u5 u8 50
0.125 0 u1 u4 u7 25
0.000 0 0 0 0 0
0.000 0.125 0.250 0.375 0.5

Partimos de :→u i+1. j +u i−1. j +ui . j+1 +ui . j−1=4 ui . j4 u1=u 4 +0+u2 +0 → 4 u1 −u4 −u2=0
4 u2=u5 +0+ u3 +u1 → 4 u2−u5−u3−u1=04 u3=u6 + 0+25+u2 → 4 u3−u 6−u2=25
4 u4 =u7 +u1 +u5 +0 → 4 u 4−u7 −u1−u 5=04 u5=u 8+u 2+u 4 +u6 → 4 u5 −u8−u 2−u4 −u6=0
4 u6=u 9+ u3+ 50+u5 → 4 u6−u 9−u3 −u5=504 u7=25+ u4 +u8 + 0→ 4 u7−u 4−u 8=25
4 u8=50+u 5+ u7 +u9 → 4 u8−u 5−u7 −u9=504 u9=75+u 6 +u8 +75 → 4 u9−u 6−u8 =150

25 25 75
u1 = u2= u3=
4 2 4
25 75 75 75 225
u 4= u5=25u6 = u7 = u8 = u9 =
2 2 4 2 4

4 .−¿ Dadala función f (x )=2 x −5 x+3 ,determinar : ¿a ¿ El incremento de x en elintervalo x=−2 hasta x=2 .
2

b ¿ El incremento de la función y en el intervalo x =−2 hasta x=2 .

c ¿ El incremento de la función y en elintervalo x hasta x+ ∆ x

Solución:

a ¿Incremento de x en la función f ( x ) → ∆ x=x f −xi =2−(−2 )=4∴ ∆ x=4

b¿
∆ y =f ( x 2 )−f ( x 1) Por lo tanto calcularemos el incremento de y :
∆ y =f ( 2 )−f (−2 )∆ y =2 ( 2 )2 −5 ( 2 ) +3−( 2 (−2 )2−5 (−2 ) +3 )∆ y =−20

c¿
∆ y =f ( x +∆ x )−f ( x )∆ y =2 ( x + ∆ x )2−5 ( x +∆ x ) +3−( 2 ( x )2−5 ( x ) +3 )
∆ y =2 ( x +2 x . ∆ x +∆ x )−5 x−5 ∆ x +3−( 2 x −5 x +3 )
2 2 2

2 2 2 2
∆ y =2 x +4 x . ∆ x +2 ∆ x −2 x −5 ∆ x∆ y =4 x . ∆ x+ 2 ∆ x −5 ∆ x
5 .−¿ Aproxime lasolución de la ecuación de onda: ¿
2 2
∂ u ( x .t ) ∂ u( x.t)
2
−4 2
=0 0<x<1. t≥0
∂t ∂t
Con las condiciones de frontera

U (0, t) = u (1, t) =0 t>0,

Y las condiciones iniciales

∂u ( x .0 )
u ( x .0 ) = sin=( Πx ) . Y =0 0≤ x≤ 1.
∂x
Compare con la solución analítica a este problema:

U(x,t)= sin ( Πx ) cos( 2 Πt ).

Utilice h= 0.25 y k= 0.1. Y determine los valores de u (x, t) cuando t= 0.3

Solución
αk
→ Ecuaciones hiperbólicas ( Diferencias finitas ) λ= h=0.25k =0.1α 2=4 → α=2
h
2 ( 0.1 )
→ λ= =0.8→ Para la primera iteración :
0.25

u ( x i , t 1 )=f ( xi ) + k ( g ( xi ) ) λ=0.8 f ( x )=sen ( πx ) g ( x )=0→ Para la2 ° iteraciónen adelante :

[ ]
2 2 2 2
α k α k
u ( x i , t j +1) =2 1− 2
u ( x i , t i) + 2 [ u ( x i−1 , t j ) +u( x i +1 , t j ) ]−u(x i ,t j−1 )
h h

Calculamos :u ( x i , t 1 )=f ( xi ) + k ( g ( xi ) )u ( x i , t 1 )=sen ( πx ) +k ( 0 )=sen ( πx )

→ Para 2 °iteración en adelante :
u ( x i , t j +1) =2 [ 1−( 0.8 ) ] u ( x i , t i )+ ( 0.8 ) [ u ( xi−1 , t j )+ u ( xi +1 , t j ) ] −u ( x i , t j−1 )
2 2

u ( x i , t j +1) =0.72 .u ( x i , t i ) +0.64 [ u ( xi−1 , t j ) +u ( xi +1 ,t j ) ]−u ( x i , t j−1 )

u ( x i , t j +1) =0.64 u ( x i−1 ,t j ) + 0.72. u ( x i , t i) + 0.64 u ( x i+1 ,t j ) −u ( x i , t j−1 )

→ Para x , eltamaño de paso será cada h=0.25 x 0=0 x 1=0.25x 2=0.5x 3=0.75x 4 =1
→ Tamaño de paso en el tiempo :k =0.1Condiciones iniciales :→1 ° iteración
u ( x 1 , t 0 )=sen ( πx ) =sen ( π ( 0.25 ) ) =0.70710678u ( x 2 , t 0 )=sen ( πx )=sen ( π ( 0. 5 ) )=1
u ( x 2 , t 0 )=sen ( πx )=sen ( π ( 0. 75 ) )=0.70710678u ( x 3 ,t 0 )=sen ( πx )=sen ( π ( 1 ) )=0
2 ° Iteración:u ( x i , t j +1) =0.64 u ( x i−1 ,t j ) + 0.72. u ( x i , t i) + 0.64 u ( x i+1 ,t j ) −u ( x i , t j−1 )
u ( x i , t j +1) =0.64 ( 0 )+ 0.72 ( 0.7071067 8 )+ 0.64 ( 1 )−0.70710678=0.44201010

x
U(x,t) 0 0.25 0.50 0.75 1.0
0 0 0.70710678 1 0.70710678 0
0.1 0 0.70710678 1 0.70710678 0
t
0.2 0 0.44201010 0.6251 0.44201010 0
0.3 0 0.0112 0.0158 0.0112 0

De forma analítica=u ( x , t )=sen ( πx ) cos ( πx ) →t=0.3

x Analítica Aproximada e=¿ V real −V aprox ∨¿

0.25 -0.2185 0.0112 0.2299
0.5 -0.3090 0.0158 0.3248
0.75 -0.2185 0.0112 0.2297