Capítulo 4

CONDUCCIÓN BI-DIMENSIONAL,
REGIMEN PERMANENTE

T1 Ecuación de calor (k cte y sin

q T2<T1
generación, coordenadas cartesianas):
 2T  2T
2
0
x y 2

1
4.1- Método de separación de variables

Ej.: Placa plana retangular  2T  2T

 2 0
x y
2

Adimensionalizando la
y temperatura:
T  T1
T2,=1 
W T2  T1
 2  2
2
0
T1,=0 T1,=0 x y 2

CC :  (0, y)  0 ( x,0)  0
 (L, y)  0 ( x,W )  1
0 x Sep. variables - hip.:  (x, y) = X (x)Y (y)
0 T1,=0 L

2
Así,

CC :  (0, y)  0  C1  0
CC :  (x,0)  0  C2 sin x C3  C4   0
 C3  C4 (pois C2  0 pues  es función de x)
CC :  (L, y) = 0  C2C4 sin Le y  e  y  0
n
 sin L  0    (n  1,2,3,...)
L
nx ny / L  ny / L nx ny
  C2C4 sin
L
e  e  n L
 C sin sinh
L
3
Como el problema es lineal

nx ny
 (x, y)   C n sin sinh
n 1
L L

nx nW
C n :  (x,W)  1   C n sin sinh (1)
n 1
L L
- El cálculo de Cn se hace usando funciones ortogonales (cualquier
función f(x) puede ser escrito en terminos de una série infinita de
funciones ortogonales


f (x)   Angn (x) (2)

n1


 f (x)gn (x) dx   a gn (x) Angn (x) dx

b b
a (3)
n1

- Obs: g1(x),g2(x), …, gn(x) son ortogonales en el intervalo (a,b) si:

g
b
a m (x)g n (x)dx  0 m n
4
 El único término no nulo de la ec. (3) es cuando m=n:

 
b b
f (x)gn (x)dx An g n2 (x)dx
a a


b
f (x)g n (x)dx
 An  a

 a gn2 (x)dx
b 

f (x)   An gn (x)
n1
usando la ec. abajo:
Comparando las ecuaciones (1) y (2):

Entonces:
f (x)  1
L nx
nx  0
sin
L
dx 2 (1)n 1  1
gn (x)  sin  An  
2 nx
 0 sin L dx  n
L L

2(1) n 1  1
 Cn  n = 1,2,3,...
n sinh(nW / L) 5
Así,

2
 (x, y)   (1)n 1
 1 nx sinh(ny / L)
sin
 n1 n L sinh(nW / L)

=1

=0.75
=0.5
=0 =0.25 =0

=0

6
4.2 - Métodos aproximados
-Proporcionan soluciones aproximadas
-Las soluciones deben compararse con resultados
experimentales o de casos más simples
-Ej .: métodos gráficos y métodos numéricos

7
4.3 - Método de diferencias finitas
- La solución se obtiene en puntos discretos
- El primer paso es definir los puntos discretos, dividiendo el
dominio de interés en pequeños volúmenes. En cada volumen
tenemos 1 punto nodal  Malla computacional
Punto nodal

∆x m,n+1
∆y

m-1,n m,n m+1,n

y,n ∆y
m-1/2,n m+1/2,n
x,m m,n-1

∆x

8
- La ecuación de conducción de calor puede ser
aproximada por diferencias finitas:
-Bi-dimensional, regimen permanente

-Bi-dimensional, regimen permanente, condutividad térmica

constante

9
-Ecuación de calor 2D, regimen permanente y sin generación de
calor, obtenida a partir de diferencias finitas:
Sustituyendo (1) y (2) en la ecuación de conducción de calor, se tiene

Para una malla uniforme, (i.e., ∆x=∆y)

Para una malla uniforme, (i.e., ∆x=∆y) y sin generación de calor

Tm,n1  Tm,n1  Tm1,n  Tm1,n  4Tm,n  0 (4)

Ecuación algébrica aproximada - ec. calor obtenida por diferencias

finitas
 -La ecuación também puede ser obtenida a partir de un
balance de calor en el VC que contiene el punto nodal (m,n):

Asi, la ecuacion de energia 2-D, regimen permanente y con

generación de calor, para malla uniforme, queda:

12
-Las ecuaciones de acuerdo con la normativa vigente.
- Para puntos en el contorno, las ecuaciones deben ser obtenidas
separadamente, a partir de un balance de calor
- Ejemplo:

13
-Balance de calor en el punto nodal, considerando
convección en la superfície externa:

14
Ejemplo:
Usando el método de balance de energía, se produce la ecuación
de diferencias finitas para (m, n) en un plano situado en una
superficie aislada, con generación de calor interna

15
-Solución de las ecuaciones de diferencias finitas:
- La solución de los sistemas de ecuaciones puede ser
obtenida por métodos directos o iterativos
Métodos directos: implican un número fijo de operaciones
aritméticas. Son más utilizados para soluciones de
sistemas pequeños, pues usan mucha memoria y pasan
mucho tiempo computacional.

Métodos iterativos: más utilizados para las soluciones de

sitios grandes. Requisitos de computación reducidos

16
-Método de inversión de matrices

•Considere un sistema de N ecuaciones y N incógnitas

•A cada punto nodal (m,n) y asociado a un número I
•El sistema de ecuaciones algébricas a ser resuelto esta dado por:

Obs.: el cálculo de
matrices inversas
puede no ser
numericamente 17
eficiente
-Método iterativo de Gauss-Seidel
1. Las ecuaciones algébricas debem ser ordenadas tal que los
elementos diagonal sean mayores que los otros
(|a11|>|a12|,|a13|,…)sistema diagonal dominante (tasa de
convergencia maximizada)
2. Despues de esta reordenacion se debe explicitar la incógnita (T):
i1 N
a a ij k1 i:1,2,…,N
  Tj  
c
Ti 
k i ij k
Tj
a ii j 1 a ii a
j  i 1 ii
k:iteracción

3. Se debe proporcionar un valor inicial para T en cada punto nodal

4. Nuevos valores para T son calculados a cada iteracción k
5. Las iteracciones terminam cuando un critério de convergencia
es satisfecho ([Link]., i|T k-i T k-1| ≤ )

18
Comentários

Los resultados obtenidos numericamente debem ser siempre

verificados:
- Se pueden realizar balances de energia
-Comparaciones con resultados experimentales
-Soluciones de problemas con mismas características, sin
embargo con algunas simplificaciones ([Link]., geometrias mas
simples)
-Comparaciones con solucciones exactas conocidas
-Comparaciones con resultados numéricos de la literatura

A definición de la malla computacional debe llevar en cuenta el

hecho de que los resultados no pueden depender de la malla:
pruebas de malla deben ser realizados. La malla escojida debe
tener en cuenta la precisión de los resultados y el costo
computacional: 19
resultado mas precisomayor costo computacional