Nome: Matheus Silva

1 Equação do calor discretizada pelo método

Crunk-Nicolson
Equação do calor:

∂ 2T ∂ 2T
 
∂T
=κ + (1)
∂t ∂x2 ∂y 2

1.1 Discretização temporal

Apromimação do termo temporal:

z+1 z
∂T Ti,j − Ti,j
≈ (2)
∂t ∆t

1.2 Discretização espacial

Aproximação do termo espacial direção x:
" #
z z z z+1 z+1 z+1
∂ 2T 1 Ti+1,j − 2Ti,j + Ti−1,j Ti+1,j − 2Ti,j + Ti−1,j
≈ + (3)
∂x2 2 ∆x2 ∆x2

Aproximação do termo espacial direção y:

" #
z z z z+1 z+1 z+1
∂ 2T 1 Ti,j+1 − 2Ti,j + Ti,j−1 Ti,j+1 − 2Ti,j + Ti,j−1
≈ + (4)
∂y 2 2 ∆y 2 ∆y 2

" #
z+1 z z z z z+1 z+1 z+1
Ti,j − Ti,j κ Ti+1,j − 2Ti,j + Ti−1,j Ti+1,j − 2Ti,j + Ti−1,j
= +
∆t 2 ∆x2 ∆x2
" #
z z z z+1 z+1 z+1
κ Ti,j+1 − 2Ti,j + Ti,j−1 Ti,j+1 − 2Ti,j + Ti,j−1
+ + (5)
2 ∆y 2 ∆y 2

1
1.3 Equações de diferenças finitas
Passando ∆x2 e ∆y 2 como fator comum:
z+1 z
Ti,j − Ti,j κ  z z z z+1 z+1 z+1

= Ti+1,j − 2T i,j + Ti−1,j + Ti+1,j − 2T i,j + Ti−1,j
∆t 2∆x2
κ  z z z z+1 z+1 z+1

+ T − 2T + T + T − 2T + T (6)
2∆y 2 i,j+1 i,j i,j−1 i,j+1 i,j i,j−1

Multiplicando os dois lados por 2 e definindo λx = κ/∆2x e λy = κ/∆2y :

n+1 n n n n n+1 n+1 n+1

2Ti,j − 2Ti,j = λx Ti+1,j − 2Ti,j + Ti−1,j + Ti+1,j − 2Ti,j + Ti−1,j +
n n n n+1 n+1 n+1


λy Ti,j+1 − 2Ti,j + Ti,j−1 + Ti,j+1 − 2Ti,j + Ti,j−1 (7)

Expandindo os termos, mantendo de cada lado da igualdade os termos

do instante de tempo n + 1 e n:

n+1 n+1 n+1 n+1 n+1 n+1 n+1

2Ti,j − λx Ti+1,j + 2λx Ti,j − λx Ti−1,j − λy Ti,j+1 + 2λy Ti,j − λy Ti,j−1 =
n n n n n n n
2Ti,j + λx Ti+1,j − 2λx Ti,j + λx Ti−1,j + λy Ti,j+1 − 2λy Ti,j + λy Ti,j−1
(8)
Agrupando por fator comum:

n+1 n+1 n+1 n+1 n+1

2 (1 + λx + λy ) Ti,j − λx Ti−1,j − λx Ti+1,j − λy Ti,j−1 − λy Ti,j+1 =
n n n n n



2 (1 − λx − λy ) Ti,j + λx Ti+1,j + Ti−1,j + λy Ti,j+1 + Ti,j−1 (9)

1.4 Equações para cada ponto da malha 3x3

1.4.1 i=1, j=1
Atribuindo x = 1, y = 1 para obter a equação para o primeiro ponto da
malha:

n+1 n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T1,1 − λx T2,1 + T0,1 − λy T1,2 + T1,0 =
n n n n n



2 (1 − λx − λy ) T1,1 + λx T2,1 + T0,1 + λy T1,2 + T1,0 (10)

2
 Expandindo os termos com fator comum:

n+1 n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T1,1 − λx T2,1 − λx T0,1 − λy T1,2 − λy T1,0 =
n n n n n
2 (1 − λx − λy ) T1,1 + λx T2,1 + λx T0,1 + λy T1,2 + λy T1,0 (11)

Os termos com índice 0 são condições de contorno, portanto, valores

conhecidos. Posicionando para o lado direito da equação, obtêm-se a equação
para o ponto i = 1, j = 1:

n+1 n+1 n+1

2 (1 + λx + λy ) T1,1 − λx T2,1 − λy T1,2 =
n n n n n n+1 n+1
2 (1 − λx − λy ) T1,1 +λx T2,1 +λx T0,1 +λy T1,2 +λy T1,0 +λx T0,1 +λy T1,0
(12)

1.4.2 i=2, j=1

Equações para o ponto i = 2, j = 1:

n+1 n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T2,1 − λx T3,1 − λx T1,1 − λy T2,2 − λy T2,0 =
n n n n n
2 (1 − λx − λy ) T2,1 + λx T3,1 + λx T1,1 + λy T2,2 + λy T2,0 (13)
n+1
O termo T2,0 é condição de contorno, portanto, valor conhecido. Pas-
sando para o lado direito obtêm-se a equação para o ponto i = 2, j = 1:

n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T2,1 − λx T3,1 − λx T1,1 − λy T2,2 =
n n n n n n+1
2 (1 − λx − λy ) T2,1 + λx T3,1 + λx T1,1 + λy T2,2 + λy T2,0 + λy T2,0 (14)

1.4.3 i=3, j=1

Equações para o ponto i = 3, j = 1:

n+1 n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T3,1 − λx T4,1 − λx T2,1 − λy T3,2 − λy T3,0 =
n n n n n
2 (1 − λx − λy ) T3,1 + λx T4,1 + λx T2,1 + λy T3,2 + λy T3,0 (15)

3
 n+1 n+1
Os termos T4,1 e T3,0 são condições de contorno, portanto, valores co-
nhecidos. Passando para o lado direito, obtêm-se as equações para o ponto
i = 3, j = 1:

n+1 n+1 n+1

2 (1 + λx + λy ) T3,1 − λx T2,1 − λy T3,2 =
n n n n n n+1 n+1
2 (1 − λx − λy ) T3,1 +λx T4,1 +λx T2,1 +λy T3,2 +λy T3,0 +λx T4,1 +λy T3,0
(16)

1.4.4 i=1, j=2

Atribuindo x = 1, y = 2:

n+1 n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T1,2 − λx T2,2 − λx T0,2 − λy T1,3 − λy T1,1 =
n n n n n
2 (1 − λx − λy ) T1,2 + λx T2,2 + λx T0,2 + λy T1,3 + λy T1,1 (17)
n+1
O termo T0,2 é condição de contorno, portanto, valor conhecido. Posici-
onando para o lado direito da equação, obtêm-se as equações para o ponto
i = 1, j = 2:

n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T1,2 − λx T2,2 − λy T1,3 − λy T1,1 =
n n n n n n+1
2 (1 − λx − λy ) T1,2 + λx T2,2 + λx T0,2 + λy T1,3 + λy T1,1 + λx T0,2 (18)

1.4.5 i=2, j=2

Atribuindo x = 2, y = 2:

n+1 n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T2,2 − λx T3,2 − λx T1,2 − λy T2,3 − λy T2,1 =
n n n n n
2 (1 − λx − λy ) T2,2 + λx T3,2 + λx T1,2 + λy T2,3 + λy T2,1 (19)

1.4.6 i=3, j=2

Equações para o ponto i = 3, j = 2:

4
 n+1 n+1 n+1 n+1 n+1
2 (1 + λx + λy ) T3,2 − λx T4,2 − λx T2,2 − λy T3,3 − λy T3,1 =
n n n n n
2 (1 − λx − λy ) T3,2 + λx T4,2 + λx T2,2 + λy T3,3 + λy T3,1 (20)
n+1
O termo T4,2 é condições de contorno, portanto, valor conhecido. Pas-
sando para o lado direito, obtêm-se a equação para o ponto i = 3, j = 2:

n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T3,2 − λx T2,2 − λy T3,3 − λy T3,1 =
n n n n n n+1
2 (1 − λx − λy ) T3,2 + λx T4,2 + λx T2,2 + λy T3,3 + λy T3,1 + λx T4,2 (21)

1.4.7 i=1, j=3

Atribuindo x = 1, y = 3:

n+1 n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T1,3 − λx T2,3 − λx T0,3 − λy T1,4 − λy T1,2 =
n n n n n
2 (1 − λx − λy ) T1,3 + λx T2,3 + λx T0,3 + λy T1,4 + λy T1,2 (22)
n+1 n+1
Os termos T0,3 e T1,4 são condições de contorno, portanto, valores co-
nhecidos. Posicionando para o lado direito da equação, obtêm-se a equação
para o ponto i = 1, j = 3:

n+1 n+1 n+1

2 (1 + λx + λy ) T1,3 − λx T2,3 − λy T1,2 =
n n n n n n+1 n+1
2 (1 − λx − λy ) T1,3 +λx T2,3 +λx T0,3 +λy T1,4 +λy T1,2 +λx T0,3 +λy T1,4
(23)

1.4.8 i=2, j=3

Equações para o ponto i = 2, j = 3:

n+1 n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T2,3 − λx T3,3 − λx T1,3 − λy T2,4 − λy T2,2 =
n n n n n
2 (1 − λx − λy ) T2,3 + λx T3,3 + λx T1,3 + λy T2,4 + λy T2,2 (24)

5
 n+1
O termo T2,4 é condição de contorno, portanto, valor conhecido. Posi-
cionando para o lado direito da equação, obtêm-se a equação para o ponto
i = 2, j = 3:

n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T2,3 − λx T3,3 − λx T1,3 − λy T2,2 =
n n n n n n+1
2 (1 − λx − λy ) T2,3 + λx T3,3 + λx T1,3 + λy T2,4 + λy T2,2 + λy T2,4 (25)

1.4.9 i=3, j=3

Equações para o ponto i = 3, j = 3:

n+1 n+1 n+1 n+1 n+1

2 (1 + λx + λy ) T3,3 − λx T4,3 − λx T2,3 − λy T3,4 − λy T3,2 =
n n n n n
2 (1 − λx − λy ) T3,3 + λx T4,3 + λx T2,3 + λy T3,4 + λy T3,2 (26)
n+1 n+1
Os termos T4,3 e T3,4 são condições de contorno, portanto, valores co-
nhecidos. Passando para o lado direito, obtêm-se a equação para o ponto
i = 3, j = 3:

n+1 n+1 n+1

2 (1 + λx + λy ) T3,3 − λx T2,3 − λy T3,2 =
n n n n n n+1 n+1
2 (1 − λx − λy ) T3,3 +λx T4,3 +λx T2,3 +λy T3,4 +λy T3,2 +λx T4,3 +λy T3,4
(27)

6
1.5 Montando matriz com o sistema de equações gover-
nantes

2(1+λx +λy ) −λx 0 −λy 0 0 0 0 0
 −λx 2(1+λx +λy ) −λx 0 −λy 0 0 0 0

 0 −λx 2(1+λx +λy ) 0 0 −λy 0 0 0

 −λy 0 0 2(1+λx +λy ) −λx 0 −λy 0 0

 0 −λy 0 −λx 2(1+λx +λy ) −λx 0 −λy 0

0 0 −λy 0 −λx 2(1+λx +λy ) 0 0 −λy


0 0 0 −λy 0 0 2(1+λx +λy ) −λx 0


 0 0 0 0 −λy 0 −λx 2(1+λx +λy ) −λx

0 0 0 0 0 −λy 0 −λx 2(1+λx +λy

 n+1   n n n n n n+1 n+1 
T1,1 2 (1 − λx − λy ) T1,1 + λx T2,1 + λx T0,1 + λy T1,2 + λy T1,0 + λx T0,1 + λy T1,0
n+1  n n n n n n+1
T2,1  2 (1 − λx − λy ) T2,1 + λx T3,1 + λx T1,1 + λy T2,2 + λy T2,0 + λy T2,0 
n+1 n+1 
 n+1   n n n n n

T3,1  2 (1 − λx − λy ) T3,1 + λ x T 4,1 + λ x T 2,1 + λ y T 3,2 + λ y T 3,0 + λ x T 4,1 + λ y T3,0
n+1
 n+1   n n n n n

T1,2  
 n+1   2 (1 − λx − λy ) T1,2 + λx T2,2 + λx T0,2 + λy T1,3 + λy T1,1 + λx T0,2 

n n n n n
T2,2  = 
 n+1   2 (1 − λ x − λ y ) T 2,2 + λ x T 3,2 + λ x T 1,2 + λ y T 2,3 + λ y T 2,1


n n n n n n+1
T
3,2
 
 n+1   2 (1 − λx − λ y ) T3,2 + λ x T4,2 + λ x T2,2 + λ y T3,3 + λ y T3,1 + λ x T4,2


T  2 (1 − λx − λy ) T n + λx T n + λx T n + λy T n + λy T n + λx T n+1 + λy T n+1 
 1,3   1,3 2,3 0,3 1,4 1,2 0,3 1,4 
T n+1   2 (1 − λx − λy ) T2,3n
+ λx T3,3 n
+ λx T1,3 n
+ λy T2,4 n
+ λy T2,2 n
+ λy T2,4 n+1 
2,3
n+1 n n n n n n+1 n+1
T3,3 2 (1 − λx − λy ) T3,3 + λx T4,3 + λx T2,3 + λy T3,4 + λy T3,2 + λx T4,3 + λy T3,4
(28)