TEMA 4

Ecuaciones de variación para sistemas no

isotérmicos

Ecuaciones de energía
La ecuación de energía en los distintos sistemas coordenados
Ecuaciones adaptadas para procesos de convección natural
Resumen de ecuaciones
Flujo tangencial con generación de calor de origen viscoso
Enfriamiento por transpiración
Ecuación de energía

Balance de energía:

 Velocidad de   Velocidad de   Velocidad de 

 acumulación  entrada de energía  salida de energía 
     
 = - 
 de energía   cinética e interna   cinética e interna 
cinética e interna   por convección   por convección 

 Velocidad neta   Velocidad neta 

 de adición de  de trabajo comunicado 
    vz
+ - 
calor por por el sistema
z  z vy
    y
 conducción   a los alrededores  z

vx x
vx x x

vy vz
y  y z

y x
Velocidad de acumulación de energía cinética e interna:
  ˆ 1 2
x y z  U  v 
t  2 

Velocidad neta de entrada de energía cinética e interna por convección:

  1   1  
y z v x  Uˆ  v 2   v x  Uˆ  v 2  
  2 x  2  x  x 
  1 2  ˆ 1 2 
ˆ
x z v y  U  v   v y  U  v  
  2  y  2  y  y 
  1   1  
x y v z  Uˆ  v 2   v z  Uˆ  v 2  
  2 z  2  z  z 

Velocidad neta de adición de calor por conducción:


y z q x x
 qx x  x   y z qy y  qy y  y   x y qz z  qz z  z 
Velocidad neta de trabajo comunicado por el sistema a los alrededores:

a. Fuerzas gravitacionales:

 x y z v x g x  v y g y  v z g z 
b. Fuerzas de presión:

y z  pv x  x  x   pv x  x 
x z pv y   y  y   pv y  y 
x y  pv z  z  z   pv z  z 
c. Fuerzas viscosas:


y z  xx v x   xy v y   xzv z  x  x   xxv x   xy v y   xzv z  x 

x z  yx v x   yy v y   yzv z  y  y   yxv x   yy v y   yzv z  y 
x y  zx v x   zy v y   zzv z  z  z   zxv x   zy v y   zzv z  z 
Transformación a un sistema de coordenadas móvil:



Dt

D ˆ 1 2
    
   
U  2 v    .q   v .g    .pv   .  .v  

Restando la ecuación de energía mecánica se obtiene la Ecuación de energía calorífica:

DUˆ   

Dt
    
   .q  p .v   :  v 
[1] [2] [3] [4]

1. Velocidad de ganancia de energía interna por unidad de volumen,

2. entrada de energía interna por conducción,
3. aumento reversible de energía interna debido a la compresión
4. aumento irreversible de energía interna debido a la disipación viscosa.
La ecuación de energía en función de la temperatura

 Uˆ   Uˆ    p   ˆ ˆ
ˆ
dU   ˆ    
 Vˆ 
dV   dT  p T    dV  CV dT
 T   T Vˆ    T Vˆ 

DT   p    Generación
ˆ
CV
Dt
 
   .q  T  
  .v   : v
 T Vˆ
  
  
 de energía 

Simplificando para el caso de fluido newtoniano y conductividad calorífica

constante:

DT  p   
ˆ
CV
Dt
2
 k T  T  
  .v  v
 T 

Casos particulares

(a) Gas ideal.

 p  p DT  
  
 T Vˆ T
 ˆ
CV
Dt
2

 k  T  p  .v 
(b) Proceso a presión constante.
DT
p  cons tan te  Cˆ p  k  2T
Dt

(c) Fluido incompresible.

  DT
  cons tan te  .v  0  Cˆ p  k  2T
Dt
(d) Sólido.
  DT
v 0  .v  0  Cˆ p  k  2T
Dt
La ecuación de energía en coordenadas rectangulares
(en función de las densidades de flujo)
 T T T T   q x q y q z 
 Cˆv   vx  vy  vz       
  t  x  y  z    x y z 
 p   v x v y v z   v x v y v z 
T         xx   yy   zz 
  T    x  y  z    x y z 
  v v y   v x v z    v y v z  
   xy  x     xz      yz   
  y x    z x    z y  

La ecuación de energía en coordenadas cilíndricas

(en función de las densidades de flujo)

 T T v  T T  1  1 q qz 
 Cˆv   vr   vz   (rq )  
z   z 
r
 t r r   r r r 
 p   1  1 v  v z   v r 1  v   v z 
T    ( rv r )    
  rr r      v r   zz 
  T   r  r r  z   r    z 
    v   1 v r   v z v r   1  v z v   
  r   r      rz      z   
  r  r  r    r z   r  z  
La ecuación de energía en coordenadas esféricas
(en función de las densidades de flujo)

 T T v  T v  T 
 Cˆv   vr   
  t  r r  r sen   
1  2 1  1 q  
 2 (r q r )  q sen    
 r r r sen   r sen   
 p   1  2 1  1 v  
T    ( r v r )   
v sen    
 T    r 2 r r sen   r sen   
 v r  1 v  v r   1 v  v r v  cot   
  rr             
  r  r  r   r sen   r r  
  v  1 v r v    v  1 v r v  
  r        r   r  r sen    r 
   r r  r   
 1 v  1 v  cot   
    v 
 r  r sen   r  
La ecuación de energía en coordenadas rectangulares
(en función de las propiedades de transporte, para ρ, μ y k constantes)

  T  T  T  T    2T  2T  2T 
Cˆv   vx  vy  vz k 2   2
  t  x  y  z    x  y 2
z 
 v 2  v 2  v  2    v  v 
2
y y
 2  x     
z
    
x
 
   x    y    z      y  x 

 v x v z   v y v z  
2 2
      
  z  x    z y  


La ecuación de energía en coordenadas cilíndricas

(en función de las propiedades de transporte, para ρ, μ y k constantes)

 T T v  T T  1   T  1  T  T 
2 2
Cˆv   vr   vz  k  r  2 2  2
 t r r  z   r r   r  r  z 
2
 v 2  1  v     v 
2
   v
 1  v 
2
2  r       v r     z         z

 r   r      z    z r  
2 2
 v v   1 v r   v  
 z  r    r 
 r z   r  r  r   
La ecuación de energía en coordenadas esféricas
(en función de las propiedades de transporte, para ρ, μ y k constantes)

 T T v  T v  T   1   2 T 
 Cˆv   vr     k  2 r  r r 
 t r r  r sen    r  
1   T  1  2T   v 2
 2  sen   2 2 2 
 2  r 
r sen      r sen     r 
2 2
 1 v  v r   1 v  v r v  cot   
      
 r  r   r sen   r r  
2
    v  1 v  2  1 v   v   
 r r
   r      r  
   r  r  r    r sen    r  r 
2
 sen    v   1 v   
     
 r   sen   r sen    
Componentes de la densidad de flujo de energía q (Ley de Fourier)

Coordenadas rectangulares:
T T T
qx  k qy  k qz  k
x y z

Coordenadas cilíndricas:
T 1 T T
qr   k q   k qz  k
r r  z

Coordenadas esféricas:
T 1 T 1 T
qr   k q   k qz  k
r r  r sen  
Ecuaciones adaptadas para procesos de convección natural

Limitaciones a la transformación de la ecuación de movimiento:

• Bajas velocidades de fluido
• Pequeñas variaciones de temperatura.

Fluido en reposo (ley de la hidrostática): 

  
p  ρg 

Desarrollo en serie de la densidad:  
 Dv  

     g T  T    
Dt
   
  . 
  g T T 


Ec. de movimiento: 
 
Dv  
   .    p  g 
Dt