TEMA 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

1. INTRODUCCION
CONDICIONES PARA EL RÉGIMEN PERMANENTE:

- CONDICIONES DE CONTORNO INDEPENDIENTES DEL

TIEMPO

- SISTEMA HA ALCANZADO EL EQUILIBRIO TÉRMICO CON

DICHAS CONDICIONES DE CONTORNO

∂T
≅0
∂t
CASOS OBJETO DE ANÁLISIS:

- UNIDIMENSIONALES
- BI- Ó TRI- DIMENSIONALES

- SIMPLES, HOMOGÉNEOS
- COMPUESTOS

- CONDUCTIVIDAD CONSTANTE
- CONDUCTIVIDAD VARIABLE

- SIN GENERACIÓN INTERNA

- CON GENERACIÓN INTERNA

2. SISTEMAS SIMPLES UNIDIMENSIONALES SIN

GENERACIÓN
GRADIENTES DE TEMPERATURA SIGNIFICATIVOS SÓLO EN
UNA DIRECCIÓN.

EJEMPLOS:
• PAREDES DE HORNOS.
• PAREDES DE CÁMARAS FRIGORÍFICAS.
• FLUJO A TRAVÉS DE TUBERÍAS O CONDUCTOS.
TEMA 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

PLACA PLANA
T1 HIPOTESIS
T2
• G=0
• Medio homogéneo
• k = cte.

• ECUACIÓN DIFERENCIAL:

L d2 T
2
=0
dx

• CONDICIONES DE CONTORNO:

 x = 0 T = T1

 x = L T = T2

• CAMPO DE TEMPERATURAS
T1
T − T2 T2
T( x ) = T1 − 1 ⋅x
L

• FLUJO DE CALOR

dT k W
q( x ) = cte. q( x ) = − k q= ⋅ ( T1 − T2 )  m2 
dx L

k⋅A
Q = q⋅ A = ⋅ (T1 − T2 )
L
TEMA 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

CILINDRO HUECO INFINITO

r2
r1 r2

r1 T2
T1

• ECUACIÓN DIFERENCIAL

d dT
(r )=0
dr dr

• CONDICIONES DE CONTORNO:

r = r1 T = T1

r = r2 T = T2

T2 − T1
T ( r ) = T1 + ln ( r / r1 )
• TEMPERATURA: ln ( r2 / r1 )

2 πkL
Q( r ) = (T1 − T2 ) = cte
• FLUJO DE CALOR: ln ( r2 / r1 )

1 k (T1 − T2 )
q( r ) =
• DENSIDAD DE FLUJO: r ln ( r2 / r1 )
TEMA 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

ESFERA HUECA

r2 r2
r1
T2
r1 T1

• ECUACIÓN DIFERENCIAL:

d 2 dT
(r )=0
dr dr

• CONDICIONES DE CONTORNO:

r = r1 T = T1

r = r2 T = T2

1 1
−
r r1
• TEMPERATURA: T(r ) = T1 + ( T − T2 )
1 1 1
−
r1 r2

4 πkr1 ⋅ r2
• FLUJO DE CALOR: Q (r ) = ( T1 − T2 ) = cte
r2 − r1

1 kr1 ⋅ r2
• DENSIDAD DE FLUJO: q(r ) = ( T1 − T2 )
r r2 − r1
2
TEMA 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

3. CONCEPTO DE RESISTENCIA TÉRMICA

Transferencia de calor análoga al flujo de electricidad:
I
V1 − V2
I=
R
V1 R V2

INTENSIDAD DE CORRIENTE ≈ FLUJO DE CALOR

DIFERENCIA DE POTENCIAL ≈ DIFERENCIA DE TEMPERATURA

RESISTENCIA ELECTRICA ≈ RESISTENCIA TERMICA

Por tanto expresamos el flujo de calor como:

∆T
Q = q⋅ A = DIFERENCIA DE TEMPERATURA
R RESISTENCIA TERMICA

PLACA PLANA CILINDRO HUECO ESFERA HUECA

L ln(r2 r1 ) r2 − r1
R= R= R=
kA 2π ⋅ k ⋅ L 4π ⋅ k ⋅ r1 ⋅ r2

RESISTENCIA TÉRMICA PELÍCULA CONVECTIVA

TS − T∞
TS Q = h ⋅ A ( TS − T∞ ) =
T
∞ 1 (h ⋅ A )
h
R = 1 (h ⋅ A )
TEMA 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

4. SISTEMAS COMPUESTOS SIN GENERACIÓN

COEFICIENTE GLOBAL DE TRANSFERENCIA

PLACAS EN “SERIE”
L1 L2 Ln

T1

T2
k1 k2 Q kn
T1 T2

L1/k1A L2/k2A Ln/knA

Las n placas en “serie” pueden ser tratadas

análogamente a n resistencias en serie que se suman
R T = ∑ Ri

T1 − T2 ∆T
Q= =
L1
+
L2
k 1A k 2 A
+ ...+
Ln
k nA
∑ Re s

CAPAS CILÍNDRICAS CONCÉNTRICAS EN “SERIE”

Te

Ti hi
r4 he
r3
r1 r2
Ti − Te
Q=
1 ln(r2 r1 ) ln(r3 r2 ) ln(r4 r3 ) 1
+ + + +
hi Ai 2 πk 1L 2 πk 2 L 2 πk 3 L he A e
TEMA 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

SISTEMAS COMPLEJOS “SERIE” / “PARALELO”

La analogía eléctrica puede extenderse para resolver

problemas más complejos, con resistencias tanto en serie
como en paralelo:

La resolución es aproximada pues al ser distinta la

conductividad de los materiales en paralelo, el problema sería
bidimensional.

Esta metodología se utiliza como método aproximado en el

cálculo de resistencias térmicas en elementos de edificación:

Ver apartado 3 del Documento de Apoyo al Documento Básico

HE (Ahorro de Energía) del Código Técnico de la Edificación
(CTE), DA-DB-HE/1: “Cálculo de parámetros característicos de
la envolvente”
TEMA 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

COEFICIENTE GLOBAL DE TRANSFERENCIA (U)

Flujo de calor que transfiere un sistema compuesto sin
generación por unidad de área y por grado de
diferencia de temperatura entre sus extremos

∆T Q 1
Q = U A ∆T = U= =
∑R A ⋅ ∆T A ⋅ ∑ R

• PLACA PLANA

1
U=
1 L1 L 2 L 1
+ + + ...+ n +
hi k1 k 2 kn he

• CILINDRO HUECO

Q = U1 A 1 ∆T = U2 A 2 ∆T = ... = U4 A 4 ∆T

Se define un COEFICIENTE GLOBAL para cada

superficie. El siguiente está referido al área exterior:

1
U4 =
r4 r ln(r2 r1) r4 ln(r3 r2 ) r4 ln(r4 r3 ) 1
+ 4 + + +
hi r1 k1 k2 k3 he
TEMA 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

5. CONDUCTIVIDAD TÉRMICA VARIABLE

• b < 0: k ↓ si T ↑
k = k0 (1+ bT ) ≠ cte • b > 0: k ↑ si T ↑
• b = 0: k = k0 = cte

La ecuación queda:

d dT dT
(k ) = 0 ⇒ −k = q = cte
dx d x dx

T1

X=0 T = T1
X=L T = T2 T2

dT
Q = − Ak0 (1 + bT )
El flujo de calor es constante dx

  T + T2 
Si se define k m = k0 1+ b  1  = k0 [1+ b Tm ]
  2 

km A
El flujo de calor se expresa como Q= (T1 − T2 )
L
TEMA 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

6. SISTEMAS SIMPLES UNIDIMENSIONALES

CON GENERACIÓN

GENERACIÓN: Conversión de una forma de energía en

energía térmica.
Ejemplos:
• Efecto Joule en conductores eléctricos
• Reacciones químicas exotérmicas
• Reacciones nucleares
• Absorción de radiación electromagnética en
medios semitransparentes a dicha radiación
• Fraguado del hormigón
• Etc.

PLACA PLANA
d 2T G
Ecuación Diferencial:
2
+ =0
dx k
 x = − L T = T1
Condiciones de Contorno: 
 x = + L T = T2

GL2  x 2  T2 − T1 x T1 + T2
T ( x) = 1 − 2  + +
2k  L  2 L 2
C .Cont .
Gen .

dT
Q = − kA =
dx
Q ( x ) = GAx +
kA
(T1 − T2 )
2L
Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

CILINDRO MACIZO
Ecuación Diferencial: 1 d  dT  G
r + = 0
r dr  dr  k

 dT
Condiciones de Contorno: r = 0 =0
 dr
r = R T = Ts

T (r ) =
G 2
4k
(
R − r 2 + Ts )
dT
Q ( r ) = − kA( r ) =
dr
Q ( r ) = Gπr 2 H
ESFERA MACIZA
1 d  2 dT  G
Ecuación Diferencial:
2 r + = 0
r dr  dr  k

 dT
Condiciones de Contorno: r = 0 =0
 dr
r = R T = Ts

T (r ) =
G 2
6k
( )
R − r 2 + Ts

dT 4
Q ( r ) = − kA( r ) = G πr 3
dr 3
Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

PLACA PLANA (GENERACIÓN Y k VARIABLE)

d  dT 
Ecuación Diferencial: k +G = 0
dx  dx 
k = k 0 (1 + bT )

 x = 0 T = T1
Condiciones de Contorno: 
 x = L T = T2

L
G
T − T1 k m (T2 ) x 2 x x
= + 1 − 
T2 − T1 k m (T ) L k m (T )
(T2 − T1 ) L  L 
L

 T +T 
donde: k m (T ) = k 0 1 + b 1 
 2 

(ver resolución numérica en hoja EXCEL)

(ver tres posibles resoluciones numéricas con EES)
Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

7. SISTEMAS MULTIDIMENSIONALES:
MÉTODOS DE RESOLUCIÓN

Conducción bi- ó tri- dimensional: Se establecen

gradientes de temperatura significativos en más de
una dirección espacial.

Ejemplos:
• Puentes térmicos en la edificación
Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE
Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE
Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE
Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

MÉTODOS DE RESOLUCIÓN

• Analíticos: Problemas simples académicos. Soluciones

exactas de referencia para validar otros métodos.
Ej.: Método de separación de variables + Principio de
Superposición.

• Gráficos: Actualmente en desuso por el auge de los

métodos numéricos debido a la potencia de cálculo de
los computadores. Se basa en la perpendicularidad de
isotermas y líneas de flujo

• Analógicos: La analogía con otros fenómenos

fácilmente medibles (magnitudes eléctricas) permite
cuantificar la conducción a partir de dichos fenómenos.

• Numéricos: No tienen limitaciones y permiten resolver

problemas complejos.
Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

8. SISTEMAS MULTIDIMENSIONALES:
MÉTODOS NUMÉRICOS

Posibilidades:
• Geometrías complicadas
• Sistemas compuestos (varios materiales)
• Propiedades variables (k variable, anisotropía,etc.)
• Generación interna variable
• Condiciones de contorno complejas
• Etc.

Bases de los métodos numéricos:

• Discretización del dominio (nodos)
• Planteamiento ecuaciones algebraicas en cada nodo

Tipos de métodos numéricos:

• Elementos Finitos: Buena adaptación a contornos
irregulares, elementos triangulares (tetraédricos en
3D), planteamiento de balances en forma integral en
elementos.
• Diferencias Finitas: Elementos adaptados al sistema
de coordenadas, peor adaptación a contornos
irregulares, cada elemento representa un nodo,
planteamiento de balances como en E.D.

Modos de plantear ecuaciones de nodo:

1. Aproximación Matemática: Sustitución de derivadas
por diferencias finitas
2. Aproximación Física: Balance de energía en cada
nodo
Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

APROXIMACIÓN MATEMÁTICA

Aproximación de la derivada primera:

• Diferencia hacia delante:

dF F ( xi + ∆x) − F ( xi ) F ( xi +1 ) − F ( xi )
F ′( xi ) = = lim ≅
dx xi ∆x → 0 ∆x ∆x
F ( xi +1 ) − F ( xi )
F ′( xi ) = + O{∆x}
∆x

• Diferencia hacia atrás:

dF F ( xi ) − F ( xi − ∆x) F ( xi ) − F ( xi −1 )
F ′( xi ) = = lim ≅
dx xi ∆x → 0 ∆x ∆x
F ( xi ) − F ( xi −1 )
F ′( xi ) = + O{∆x}
∆x

• Diferencia central:
dF F ( xi + ∆x) − F ( xi − ∆x) F ( xi +1 ) − F ( xi −1 )
F ′( xi ) = = lim ≅
dx xi ∆x →0 2 ∆x 2 ∆x
F ( xi +1 ) − F ( xi −1 )
F ′( xi ) =
2 ∆x
+ O ( ∆x ) 2 { }

Aproximación de la derivada segunda:

dF dF
−
d 2F dx xi + 12 dx xi − 12 F ( xi +1 ) − 2 F ( xi ) + F ( xi−1 )
F ' ' ( xi ) = = lim ≅
dx 2 xi ∆x →0 ∆x ( ∆x ) 2
F ( xi +1 ) − 2 F ( xi ) + F ( xi −1 )
F ' ' ( xi ) =
( ∆x ) 2
+ O ( ∆x ) 2 { }
Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

Ecuación en diferencias finitas de un nodo (i,j):

∂ 2T ∂ 2T
+ 2 =0
∂x 2
∂y

∂ 2T Ti +1, j − 2Ti , j + Ti −1, j

≅
∂x 2 (∆x ) 2
∂ 2T Ti , j +1 − 2Ti , j + Ti , j −1
≅
∂y 2 (∆y ) 2

Sustituyendo las diferencias finitas de las derivadas

segundas en la ecuación diferencial y haciendo ∆x=∆y:

Ti +1, j + Ti −1, j + Ti , j +1 + Ti , j −1 − 4Ti , j = 0

Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

APROXIMACIÓN FÍSICA

Balance de energía en el volumen finito

correspondiente a un nodo interior (i,j):

Ti +1, j − Ti , j
qi +1 = k
∆x
T −T
qi −1 = k i −1, j i , j
∆x
T −T
q j +1 = k i , j +1 i , j
∆y
T −T
q j −1 = k i , j −1 i , j
∆y
(qi +1 + qi −1 )∆y + (q j +1 + q j −1 )∆x = 0

Sustituyendo las aproximaciones en diferencias finitas de

los flujos de calor en el balance y haciendo ∆x=∆y:

Ti +1, j + Ti −1, j + Ti , j +1 + Ti , j −1 − 4Ti , j = 0

Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

Balance de energía en el volumen finito

correspondiente a un nodo de un contorno convectivo
(i,j):

∆l
(q j+1 + q j−1 ) + (qi−1 + qc )∆l = 0
2

T i −1 , j − T i , j
q c = h (T o − T i , j ) q i −1 = k
∆l
Ti , j +1 − Ti , j T i , j −1 − T i , j
q j +1 = k q j −1 = k
∆l ∆l
Sustituyendo las aproximaciones en diferencias finitas de
los flujos de calor en el balance resulta la ecuación del
nodo del contorno:

 2 h ∆l  2 h∆l
Ti , j +1 + Ti , j −1 + 2Ti −1, j −  4 +  i, j
T = − To
 k  k
Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

RESOLUCIÓN DE LAS ECUACIONES EN

DIFERENCIAS FINITAS

Poniendo juntas las ecuaciones para todos los N nodos del

dominio se obtiene un sistema de N ecuaciones
algebraicas:
[A]⋅ {T } = {C}
[A]: Matriz (NxN) de coeficientes de las temperaturas en las
ecuaciones de nodos (matriz muy dispersa)
{T}: Vector (Nx1) de las temperaturas de los nodos
(incógnitas del problema)
{C}: Vector (Nx1) de términos independientes de las
ecuaciones de nodos (depende de las condiciones de
contorno)
Resolución del sistema de ecuaciones por inversión de la
matriz:

{T } = [A] ⋅ {C}
−1

A veces es más eficiente desde el punto de vista

computacional resolver el sistema mediante métodos
optimizados. Importante cuidar la numeración (ordenación)
de los nodos.

Ejemplos:

Ordenación buena Ordenación mala

 Tema 3: CONDUCCIÓN EN RÉGIMEN PERMANENTE

Matriz correspondiente a la numeración buena

 a11 a12 0 a14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
a a 22 a 23 0 a 25 0 0 0 0 0 0 0 0 0 0 0 0 0 
 21
0 a32 a 33 0 0 a36 0 0 0 0 0 0 0 0 0 0 0 0 
 
a 41 0 0 a 44 a 45 0 a 47 0 0 0 0 0 0 0 0 0 0 0 
0 a52 0 a54 a 55 a56 0 a58 0 0 0 0 0 0 0 0 0 0 
 
0 0 a 63 0 a 65 a 66 0 0 a 69 0 0 0 0 0 0 0 0 0 
0 0 0 a 74 0 0 a 77 a 78 0 a 7,10 0 0 0 0 0 0 0 0 
 
0 0 0 0 a85 0 a87 a88 a89 0 a8,11 0 0 0 0 0 0 0 
0 0 0 0 0 a96 0 a98 a 99 0 0 a 9,12 0 0 0 0 0 0 
 
0 0 0 0 0 0 a10,7 0 0 a10,10 a10,11 0 a10,13 0 0 0 0 0 
 
0 0 0 0 0 0 0 a11,8 0 a11,10 a11,11 a11,12 0 a11,14 0 0 0 0 
0 0 0 0 0 0 0 0 a12,9 0 a12,11 a12,12 0 0 a12,15 0 0 0 
 
0 0 0 0 0 0 0 0 0 a13,10 0 0 a13,13 a13,14 0 a13,16 0 0 
0 0 0 0 0 0 0 0 0 0 a14,11 0 a14,13 a14,14 a14,15 0 a14,17 0 
 
0 0 0 0 0 0 0 0 0 0 0 a15,12 0 a15,14 a15,15 0 0 a15,18 
0 0 0 0 0 0 0 0 0 0 0 0 a16,13 0 0 a16,16 a16,17 0 
 
0 0 0 0 0 0 0 0 0 0 0 0 0 a17 ,14 0 a17 ,16 a17 ,17 a17,18 
 
0 0 0 0 0 0 0 0 0 0 0 0 0 0 a18,15 0 a18,17 a18,18 

Matriz correspondiente a la numeración mala

 a11 a12 0 0 0 0 a17 0 0 0 0 0 0 0 0 0 0 0 
a a 22 a 23 0 0 0 0 a 28 0 0 0 0 0 0 0 0 0 0 
 21
0 a32 a33 a34 0 0 0 0 a39 0 0 0 0 0 0 0 0 0 
 
0 0 a 43 a 44 a 45 0 0 0 0 a 4,10 0 0 0 0 0 0 0 0 
0 0 0 a54 a55 a56 0 0 0 0 a5,11 0 0 0 0 0 0 0 
 
0 0 0 0 a 65 a 66 0 0 0 0 0 a 6,12 0 0 0 0 0 0 
a 0 0 0 0 0 a77 a78 0 0 0 0 a 7,13 0 0 0 0 0 
 71 
0 a82 0 0 0 0 a87 a88 a89 0 0 0 0 a8,14 0 0 0 0 
0 0 a93 0 0 0 0 a98 a99 a9,10 0 0 0 0 a9,15 0 0 0 
 
0 0 0 a10, 4 0 0 0 0 a10,9 a10,10 a10,11 0 0 0 0 a10,16 0 0 
 
0 0 0 0 a11,5 0 0 0 0 a11,10 a11,11 a11,12 0 0 0 0 a11,17 0 
0 0 0 0 0 a12, 6 0 0 0 0 a12,11 a12,12 0 0 0 0 0 a12,18 
 
0 0 0 0 0 0 a13, 7 0 0 0 0 0 a13,13 a13,14 0 0 0 0 
0 0 0 0 0 0 0 a14,8 0 0 0 0 a14,13 a14,14 a14,15 0 0 0 
 
0 0 0 0 0 0 0 0 a15,9 0 0 0 0 a15,14 a15,15 a15,16 0 0 
0 0 0 0 0 0 0 0 0 a16,10 0 0 0 0 a16,15 a16,16 a16,17 0 
 
0 0 0 0 0 0 0 0 0 0 a17,11 0 0 0 0 a17 ,16 a17 ,17 a17 ,18 
 
0 0 0 0 0 0 0 0 0 0 0 a18,12 0 0 0 0 a18,17 a18,18 