Tema 9.- Convección forzada

Diapositiva 4 JM.Corberán, R. Royo (upv) 4 Tema 9: Convección forzada El desarrollo de la capa límite en función de la distancia viene caracterizado...

85 downloads 344 Views 411KB Size
Diapositiva 1 Tema 9: Convección forzada

CONVECCIÓN FORZADA

JM.Corberán, R. Royo (upv)

1

Diapositiva 2 Tema 9: Convección forzada

ÍNDICE

placas Flujo externo

incompresible compresible

tubo único normal a tubos haces no circulares circulares

laminar Flujo interno turbulento

JM.Corberán, R. Royo (upv)

2

Diapositiva 3 Tema 9: Convección forzada

1.FLUJO EXTERNO SOBRE PLACAS PLANAS 1.1.METODOLOGÍA DE CÁLCULO Evaluación de la temperatura de recuperación Cuando el flujo se frena, el fluido se calienta, siendo la temperatura representativa del fenómeno de intercambio de calor la temperatura de recuperación Tr, pudiéndose calcular mediante la siguiente expresión: Tr - T∞ Ec donde Ec es el número de Eckert:

Ts - T∞

2

2

Ec =

y r el factor de recuperación r= Pr (1/2) (flujo laminar) r = Pr (1/3) (flujo turbulento)

=r⋅

U∞ Cp ⋅ (Ts - T∞ )

en la mayoría de aplicaciones de ingeniería, para flujos de baja velocidad, Tr es similar a la temperatura del fluido T∞ , Tr = T∞ Propiedades a la temperatura correspondiente (Ts + T∞ ) 2 (T + T∞ ) + 0.22 ⋅ (Tr − T∞ ) T ′′ = s 2

flujo incompresible: Tm = flujo compresible: JM.Corberán, R. Royo (upv)

3

Diapositiva 4 Tema 9: Convección forzada

El desarrollo de la capa límite en función de la distancia viene caracterizado por el valor del número de Reynolds a lo largo de la placa

U∞ ⋅ x υ

< 5 ⋅105 → flujo laminar Re > 5 ⋅105 → flujo turbulento

Expresión del Número de Nusselt: es local, f(x), calculándose el coeficiente de convección según la expresión:

hx =

JM.Corberán, R. Royo (upv)

Nu x ⋅ k x 4

Diapositiva 5 Tema 9: Convección forzada

1.2. CORRELACIONES DE CÁLCULO

Q=

Li



hx ⋅ (Ts − Tr ) ⋅ w ⋅ dx

Li −1

ZONA FLUJO LAMINAR Comprobar que 0.6 ≤ Pr ≤ 50

Nusselt local Nux =0.332⋅(Re x )1/2 ⋅Pr 1/3 obteniendo h local, hx Nusselt promedio Nu 0 −L =0.664⋅(Rex =L )1/2 ⋅Pr 1/3 obteniendo hL (Nusselt obtenido a partir de h promediado)

L

=

1 h X • dx L ∫0

Q = hL ⋅ (Ts − Tr ) ⋅ A

ZONA FLUJO TURBULENTO Comprobar que 0.6 ≤ Pr ≤ 50

Si 5·105
Nu=0.0296·Re (0.8)·Pr1/3 (Hasta 108 con error inferior al 15%)

JM.Corberán, R. Royo (upv)

5

Diapositiva 6 Tema 9: Convección forzada

3.FLUJO EXTERNO SOBRE TUBOS

1.Determinación de las propiedades necesarias a 2.Determinación del

JM.Corberán, R. Royo (upv)

ReD =

Tm =

T∞ + TS 2

U∞ ⋅ D µ/ρ

6

Diapositiva 7 Tema 9: Convección forzada

2.1.TUBO ÚNICO SECCIÓN CIRCULAR Comprobar: 1
NuD = C ⋅ Re D

(m)

1

Pr ⋅ Pr ⋅ ( ∞ ) 4 Prs n

0.37 si Pr ≤ 10 0.36 si Pr > 10

Determinación de n = 

Determinación de C y m (tabla siguiente)

ReD 1-40 40-103 103-2*105 2*105-106 JM.Corberán, R. Royo (upv)

C 0.75 0.51 0.26 0.076

m 0.4 0.5 0.6 0.7 7

Diapositiva 8 Tema 9: Convección forzada

TUBOS DE SECCIÓN NO CIRCULAR

NuD = C ⋅ Re D

(m)

⋅ Pr

1 ( ) 3

C y m se determinan a partir de la tabla siguiente, según la geometría y el número de Reynolds ReD: G E O M E T R IA

ReD

C

m

C ua d rad o

u

D

u

3

5

0.2 4 6

0.5 8 8

3

5

0.1 0 2

0.6 7 5

0.1 6 0

0.6 3 8

0.0 3 85

0.7 8 2

0.2 2 8

0.7 3 1

5x 10 -1 0 D

5x 10 -1 0

H ex á go no

u

D

u

3

5x 10 -1 .95x 10 4

D

1.9 5x1 0 -10

4

5

P la c a vertic al

u JM.Corberán, R. Royo (upv)

3

4x 10 -1 .5x1 0

4

8

Diapositiva 9 Tema 9: Convección forzada

2.2. BATERÍAS DE TUBOS: USO GENERALIZADO EN INTERCAMBIADORES

Haz de cuatro hileras y tres tubos por hilera

JM.Corberán, R. Royo (upv)

9

Diapositiva 10 Tema 9: Convección forzada

1º Determinación de UMÁX según la geometría ReDM =

Sr

D

A1

D

Sr U∞ T∞

U∞ T∞

SD

SL

SL

U MAX ⋅ D µ/ρ

A2

A1

A2

A1

Disposición de los tubos en un banco. (a) Alineados (b) Alternados

Aplicando la ecuación de continuidad en cada caso tendremos: S

TUBOS ALINEADOS U∞ · ST= Umáx·(ST-D) , de donde Umáx = U∞ ⋅ S −T D T TUBOS ALTERNADOS U∞ · ST= Umáx·AMin, de donde Umáx = U∞ ⋅

ST A MIN

siendo Amin = Mín(A1,2A2) JM.Corberán, R. Royo (upv)

10

Diapositiva 11 Tema 9: Convección forzada

Si el número de HILERAS ≥ 20: 1º Evaluar Pr∞ a T∞ , PrS a TS 2º Determinar el resto de propiedades a Tm =

Ts + T∞ 2

3º Determinar C y n según el valor de ReDmáx y la geometría. 4º Comprobar

1000 < ReDM < 2 * 106  0.7 < Pr < 500 1

5º NuD = C ⋅ RenDM ⋅ PrT

m

JM.Corberán, R. Royo (upv)

0.36

 Pr∞  4   Pr  S

11

Diapositiva 12 Tema 9: Convección forzada

Tabla 7.7: Constantes de la ecuación 7.67 para batería de tubos en flujo cruzado (16) CONFIGURACIÓN

ReD,MÁX

C

M

Alineados

10-102

0.80

0.40

Alternados

10-102

0.90

0.40

Alineados

102-103

Alternados

102-103

Alineados

Aproximar a un único (aislado) cilindro

103-2x105

0.27

0.63

103-2x105

0.35 (ST/SL)1/5

0.60

103-2x105

0.40

0.60

Alineados

2x105-2x106

0.021

0.84

Alternados

2x105-2x106

0.022

0.84

(ST/SL<0.7)a Alternados (ST/SL< 2) Alternados (ST/SL> 2)

a Para S /S >0.7 la transferencia de calor es ineficiente y los tubos alineados no deben usarse JM.Corberán, R. Royo (upv) T L

12

Diapositiva 13 Tema 9: Convección forzada

Si el número de HILERAS <20: se corrige el NuD calculado anteriormente con el factor:

(NuD )N hileras (NuD ) 20 hileras cuyos valores aparecen en la siguiente tabla (para ReD >1000):



1

2

3

4

5

7

10

13

16

alineados

0.7

0.8

0.86

0.9

0.92

0.95

0.97

0.98

0.99

alternados

0.64

0.76

0.84

0.89

0.92

0.95

0.97

0.98

0.99

JM.Corberán, R. Royo (upv)

13

Diapositiva 14 Tema 9: Convección forzada

3.FLUJO INTERNO 3.1. INTRODUCCIÓN

En flujo interno, el perfil de temperaturas depende de la posición en la sección transversal y a lo largo del conducto: T=T(x,r). TS R

U

r

T

Distribución de la temperatura

Tb

Distribución de temperaturas en el flujo de tuberías

JM.Corberán, R. Royo (upv)

14

Diapositiva 15 Tema 9: Convección forzada

Temperatura de masa TB: temperatura característica para cada sección. Corresponde a un promedio energético de la sección: R

TB =

∫ ρ ⋅ Cp ⋅ 2 ⋅π ⋅ r ⋅ T(r) ⋅ dr 0

R

∫ ρ ⋅ Cp ⋅ 2 ⋅π ⋅ r ⋅ dr 0

El calor transmitido por convección en cada posición se refiere a dicha temperatura de masa: dQ(x)=dA(x) h· (Ts(x) –TB(x))= P(x) dx h· (Ts(x) –TB(x)) Las propiedades necesarias para el cálculo del coeficiente de convección se evalúan a la temperatura de masa promediada entre la entrada y la salida del conducto:

JM.Corberán, R. Royo (upv)

TBi + TB0 2

15

Diapositiva 16 Tema 9: Convección forzada

Variación de temperatura a lo largo del conducto Suponiendo temperatura superficial interior y coeficiente de convección uniformes:

JM.Corberán, R. Royo (upv)

16

Diapositiva 17 Tema 9: Convección forzada

Siendo P el perímetro del tubo, tenemos: dQ = h ⋅ P ⋅ dx ⋅ (Ts − TB ) = m ⋅ Cp ⋅ dTB

dTB h⋅P =− ⋅ dx (TB − TS ) m ⋅ Cp h⋅ P ⋅ x

− TB − TS integrando = e m⋅Cp TBi − TS L

Q = ∫ dQ = h ⋅ P ⋅ L ⋅ DMLT 0

siendo DMLT la Diferencia Media Logarítmica de Temperatura: DMLT =

JM.Corberán, R. Royo (upv)

(TS - TBi ) - (TS - TB0 ) ≈ (TS − TB )promedio (TS - TBi ) ln (TS - TB0 ) 17

Diapositiva 18 Tema 9: Convección forzada

3.2. FLUJO LAMINAR

Re D =

U∞ ⋅ D < 2300 µ/ρ

-Todas las propiedades se evalúan a la temperatura de masa del fluido. Tubo corto (flujo no plenamente desarrollado). Comprobar:

(Re D ⋅ Pr/ L / D)1/ 3 > 2

0.48 < Pr < 16700 Ts = constante 1/3

  D  µ 0.14 Nu D = 1.86 ⋅    ⋅ Re D ⋅ Pr  ⋅ ( ) µ L    S Tubo largo (flujo desarrollado) (constantes válidas sólo para tubos cilíndricos)

 h⋅D

 = 3.66 Si TS es uniforme: Nu D =  k  f 

 h⋅D  Q  = 4.36  es uniforme: Nu D =  k A  f 

Si 

JM.Corberán, R. Royo (upv)

18

Diapositiva 19 Tema 9: Convección forzada

3.3. FLUJO TURBULENTO

ReD =

U∞ ⋅ D > 2300 µ/ρ

En flujo turbulento la región de entrada es mucho más corta, por lo que la siguiente correlación se considera válida para todo el conducto. Válida para:

L ≥ 10 D 0.7 < Pr < 160

NuD = 0.023 ⋅ Re D

0.8

⋅ Pr

n

0.3 si TS < TB (flujo se enfria) n 0.4 si TS > TB (flujo se calienta)

Nota: En conductos no circulares utilizar la definición de radio hidráulico:

Dh = JM.Corberán, R. Royo (upv)

4 ⋅ Area de paso Perímetro mojado 19

Diapositiva 20 Tema 9: Convección forzada

3.4.CORRECCIÓN POR EL EFECTO DE LA VARIACIÓN DE PROPIEDADES ENTRE LA TEMPERATURA DE LA SUPERFICIE Y EL FLUJO: LIQUIDOS GASES

Pr f = ( s )m PrB fB Pr Nu = ( s )n PrB Nu B Tipo de flujo Laminar

Turbulento

JM.Corberán, R. Royo (upv)

µs µB Ts TB µs µB Ts TB

µs m T ) ó ( s )m µB TB µ T ó ( s )n ó ( s )n µB TB ó (

Fluido

Sentido transferencia calor

m

n

Líquidos

Fluido se calienta Fluido se enfría

0.58 0.50

-0.11 -0.11

1

0

Gases

Calentándose o enfriándose

Líquidos

Fluido se calienta Fluido se enfría

0.25 0.25

-0.25 -0.11

Gases

Fluido se calienta Fluido se enfría

-0.2 -0.1

-0.55 0

20