Intercambiadores de Calor
Intercambiadores de Calor
Intercambiadores de Calor
A mi tutor Don Rubén Ventas Garzón por haberme ayudado siempre que lo he
necesitado. Además de aconsejarme y guiarme a lo largo de todo el proyecto, con un
trato excelente que hacía las cosas más fáciles, sin duda una gran persona.
También agradecerles a todos mis compañeros de Universidad por estos años tan
buenos que nunca olvidaré. Cada compañero de alguna forma me ha aportado algo, no
doy ningún nombre porque los que han estado compartiendo clase y momentos conmigo
saben quienes son. Aunque a muchos no los volveré a ver, decir a todos que me
considero un afortunado de haber podido compartir todos los momentos de la carrera
con vosotros.
Por último dar las gracias a todos los profesores que he tenido a lo largo de la
carrera y que me han ayudado a finalizar el objetivo que hace cinco años comencé.
ÍNDICE
1. INTRODUCCIÓN ................................................................................... 1
TU UT
MECÁNICA ................................................................................................................ 6
UT
2.1.1.
TU CICLO DE COMPRESIÓN MECÁNICA ............................................... 6
UT TU UT
2.1.2.
TU CICLO DE ABSORCIÓN ........................................................................ 7
UT TU UT
2.2.
TU FLUIDOS DE TRABAJO .............................................................................. 8
UT TU UT
3.4. APLICACIONES............................................................................................... 19
TU UT
GENERADOR .......................................................................................................... 34
UT
I
5.1. CÁLCULO COEFICIENTE DE CONVECCIÓN ......................................... 37
TU UT
................................................................................................................... 117
................................................................................................................... 134
II
P.F.C.
1. INTRODUCCIÓN
El creciente consumo de energía demandado por la población mundial, convierte
en fundamental el estudio y la mejora de diseños de ahorro energético.
- La utilización del calor que ofrecen los paneles solares, u otros medios de
producción de energía renovables, con el consiguiente ahorro de dinero, debido
a la subida de los combustibles fósiles.
1
P.F.C.
2
P.F.C.
El Capítulo 7 recoge las conclusiones a las que se llega tras haber realizado el
proyecto y se plantean posibles trabajos futuros para completar este proyecto.
3
P.F.C.
1.2. OBJETIVOS
4
P.F.C.
2. MAQUINA DE ABSORCIÓN
Las máquinas de absorción son sistemas para la producción de frío mediante
producción de calor. Por ejemplo el calor producido en las plantas de cogeneración y el
calor residual de los motores térmicos, son fuentes de calor que pueden ser utilizados
para alimentar a las máquinas de absorción.
Las primeras noticias que se tienen de una máquina de absorción usada para la
refrigeración fue en el año 1774 construida por Nairne, posteriormente el diseño sería
perfeccionado por Carré en la segunda mitad del siglo XIX. Hasta la década de 1950
prácticamente el diseño no sufre variaciones.
Entre 1940 y 1950 las máquinas de absorción caen en desuso por el auge de los
refrigerantes halogenados y de la energía eléctrica barata que permitió el desarrollo de
los sistemas de refrigeración por compresión mecánica. Sería a partir de 1970 cuando
las máquinas de absorción vuelven a comercializarse, en el sector del aire
acondicionado de grandes edificios (sobre todo en Estados Unidos y Japón), debido al
aumento del precio de la energía eléctrica y por las limitaciones impuestas en el uso de
los refrigerantes tradicionales.
Los sistemas de absorción pueden ser de simple efecto, doble efecto, triple
efecto, etc. En este proyecto sólo se estudian los primeros ya que el banco de ensayos
trata de simular un sistema de simple efecto.
5
P.F.C.
6
P.F.C.
7
P.F.C.
Para poder resolver el ciclo de absorción debemos saber las propiedades de los
fluidos de trabajo, refrigerantes y absorbentes, con ello se determinará los pares
refrigerante-absorbente.
• Alto calor latente de cambio de fase: el agua y el amoniaco tienen un alto calor
latente de cambio de fase (en torno a 10 veces el de los refrigerantes
halogenados).
• Debe ser estable, no inflamable y que no sea tóxico, ni peligroso. En el caso del
NH3 está clasificado como de clase II (refrigerante de media seguridad), por lo
B B
• Precio y disponibilidad.
8
P.F.C.
2.2.3. DISOLUCIONES
El bromuro de litio es una sal de color blanca, con una temperatura de ebullición
(≈ 2.200 º C) muy superior a la del agua, por lo que la separación en el generador no es
difícil. Presenta el inconveniente de ser altamente corrosiva. Aunque es miscible con el
agua en concentraciones superiores al 60%, presenta un límite de solubilidad a
concentraciones superiores al 65%, dependiendo de la temperatura, que se manifiesta en
la formación de cristales.
9
P.F.C.
3.INTERCAMBIADOR DE PLACAS
El intercambiador de calor es un dispositivo diseñado para transferir calor de
manera eficiente de un fluido a otro, estando separados estos por una barrera sólida o
encontrándose en contacto. Según la forma de trabajo para conseguir transferir calor, se
tendrá la siguiente clasificación:
10
P.F.C.
Los fluidos frío y caliente se introducen por los orificios de conexión y circulan
por los canales que se forman entre placas de manera alternativa. Un fluido es
conducido por los canales impares mientras que el otro es conducido por los pares. Los
dos fluidos se encuentran así separados, sin poderse mezclar, por una delgada placa a
través de la cual se produce la transferencia de calor. La distribución de los fluidos por
sus canales correspondientes se hace mediante una serie de juntas en los canales
impares que no permiten la entrada del fluido que ha de circular por los pares, y en los
pares que no permite la entrada del de los impares. Generalmente, aunque existen
múltiples configuraciones, el flujo de ambos fluidos se hace en contracorriente (figura
3).
11
P.F.C.
3.1. COMPONENTES
Un intercambiador de placas consiste en una pila de placas de metal corrugadas,
que se mantienen unidas mediante presión en un bastidor y selladas por medio de una
junta, de manera que se forman una serie de pasillos interconectados a través de los
cuales se hacen circular los fluidos de trabajo, los componentes principales de un
intercambiador de placas se muestran en la figura 4.
12
P.F.C.
- Placas
Las corrugaciones son resaltes que forman los correspondientes canales entre
placas, cuyos bordes se apoyan por presión unos sobre otros o bien están soldadas entre
ellas. Cada placa se sitúa enfrente de otra girada 180º respecto a la primera, como se
puede apreciar en la figura 6, para placas de corrugaciones tipo “Chevron”.
13
P.F.C.
Figura 7. Apilamiento de las placas. En rojo se muestra los canales por los que circula el fluido
caliente y en azul el fluido frío
(referencia [3] de la bibliografía)
14
P.F.C.
• Mediante juntas (1): Las placas se unen ejerciendo presión entre ellas mediante
las barras del bastidor e interponiendo entre ellas juntas para garantizar la
estanqueidad. Se desmonta fácilmente desatornillando las barras de unión, con
lo que las operaciones de mantenimiento y limpieza se pueden efectuar
rápidamente.
• Soldado (2): Las placas se unen por medio de soldaduras de manera que no se
necesitan juntas para garantizar la estanqueidad. La ventaja principal de este tipo
es que pueden utilizarse en aplicaciones donde la presión es un obstáculo para
las juntas. Pueden trabajar hasta presiones de unos 30-32 bares. Por el contrario,
presentan el inconveniente principal de que no pueden ser desmontados para su
limpieza. Además ante un cambio en las necesidades del proceso no es posible
modificar el número de placas.
15
P.F.C.
• Un paso.
• Varios pasos (por lo general dos): pueden ser de dos pasos para cada fluido o tan
solo de un paso para uno de ellos y dos para el otro.
16
P.F.C.
• profundidad corrugación, b
• radio de curvatura, R
17
P.F.C.
18
P.F.C.
3.4. APLICACIONES
Actualmente los intercambiadores de placas están entrando en casi todos los
campos de la transferencia de calor. Por sus ventajas son ampliamente utilizados, salvo
en aquellas aplicaciones en las que la temperatura, la presión o los caudales son muy
elevados.
• Recuperación de calor.
19
P.F.C.
3.5.1. VENTAJAS
• Mayor coeficiente global de transferencia, hasta cinco veces superior a los que
se pueden conseguir para las unidades de tubo y carcasa.
20
P.F.C.
• Control de la temperatura. La forma de los canales por los que circulan los
fluidos reduce la posibilidad de zonas de retención o estancamiento y
sobrecalentamientos locales.
• Flexibilidad. Pueden adaptarse y ser utilizados para una diversa gama de fluidos
y condiciones, incluyendo fluidos viscosos.
21
P.F.C.
3.5.2. LIMITACIONES
• Los materiales de las juntas son de varios tipos de elastómeros que tienen un
límite máximo de temperatura de funcionamiento (para los materiales más
usados) es de 140 a 150 ºC. Los fluor-elastómeros pueden aumentar ese límite
hasta los 180 ºC, pero a cambio de un mayor coste. Existen algunos modelos que
usan juntas de fibras de amianto comprimido, para los cuales la temperatura
límite asciende hasta un máximo de 250 ºC. Por el hecho de no llevar estas
juntas, los intercambiadores termo-soldados son capaces de aguantar
temperaturas superiores a los 140-150 ºC de los de tipo junta.
• Para cada modelo de placa, la boca de entrada tiene un tamaño fijo que limita la
cantidad de fluidos de alto volumen específico (vapores y gases húmedos) que
pueden entrar al intercambiador, de modo que este tipo de intercambiadores casi
nunca se utilizan en servicios con gran condensación.
22
P.F.C.
4. INSTALACIÓN EXPERIMENTAL
El esquema típico de una instalación de absorción de simple efecto sería el
mostrado en la figura 2, en él capitulo de absorción, sin embargo el esquema de la
instalación experimental corresponde al expuesto en la figura 13.
23
P.F.C.
4.1. COMPONENTES
El equipo se trata de una máquina para investigación, automatizada y
computerizada, para ensayos con absorbedores de nitrato de litio-amoníaco. Está
construida con componentes de acero inoxidable 316, 304, cobre, hierro, vidrio pirex
teflón y EPDM. Todos los instrumentos que se encuentran en la zona bañada por la
disolución son de acero inoxidable 316. En la figura 14 puede verse la instalación
experimental con los componentes más destacados que se pasan a comentar a
continuación.
24
P.F.C.
25
P.F.C.
26
P.F.C.
Este dispositivo está comunicado con el ordenador del sistema SCADA, lo que
permite conocer el caudal volumétrico de gas en tiempo real.
27
P.F.C.
La línea de líquido que sale del generador (en nuestro caso sería todo el agua) se
introduce en un intercambiador de placas situado en posición vertical. Dicho
intercambiador posee las mismas características que los descritos anteriormente,
variando únicamente las dimensiones del mismo. Se encuentra aislado con algodón
cerámico, cinta de aluminio y poliuretano para reducir pérdidas. El líquido se introduce
por la conexión superior izquierda y es refrigerado mediante la disolución que proviene
del absorbedor a unos 40º C. Dicha solución se introduce por la conexión inferior
derecha. En el interior del intercambiador cada uno de los fluidos posee un canal
independiente de 0,20 (líquido) y 0,25 (disolución) litros de capacidad.
La presión y temperaturas máximas de diseño de este equipo son 22 bar a 75º C
y 19 bar a 225º C.
Este dispositivo esta comunicado con el ordenador del sistema SCADA, lo que
permite conocer el caudal de disolución concentrada y su densidad en tiempo real.
28
P.F.C.
- (10) Absorbedor
Las partes inferior y superior del depósito se encuentran selladas mediante dos
bridas que proporcionan, además, la sujeción para dos mirillas de vidrio pirex, a través
de las cuales se puede observar el proceso de absorción.
La presión y temperatura máximas de diseño de este equipo son 20 bar y 100 ºC,
respectivamente.
29
P.F.C.
30
P.F.C.
31
P.F.C.
- Circuito de disolución
32
P.F.C.
• Bomba de agua (11). Se controla el caudal de agua que circula por el circuito de
la disolución, mediante un medidor Coriolis.
• Bomba de aceite (1). Para esta bomba no se tendrá un Coriolis para poder medir
el caudal, por lo que se controlará la frecuencia.
Todos las variables a excepción de la potencia del calentador, son tomadas desde
el software ADKIR y se controlan desde el ordenador central.
33
P.F.C.
- Geometría
- Fluidos de trabajo
34
P.F.C.
35
P.F.C.
36
P.F.C.
• •
Q = mc ·Cpc ·(Tce − Tcs ) = m f ·Cp f ·(T fe − T fs ) (5.1)
donde:
• •
- mc y m f , es el gasto másico del fluido caliente y del fluido frío
respectivamente.
37
P.F.C.
Q
U= (5.2)
Atotal ·∆Tlm
donde:
(∆T1 − ∆T2 )
∆Tlm = (5.4)
⎛ ∆T ⎞
ln⎜⎜ 1 ⎟⎟
⎝ ∆T2 ⎠
donde:
∆T1 = Tce − T fs
(5.5)
∆T2 = Tcs − T fe
1 1 1
= + Ren _ c + Rcond + Ren _ f + (5.6)
U · A hc · A h f ·A
Debido a que las placas están construidas en acero inoxidable AISI 316, cuya
conductividad térmica es de 13,4 W/(m ºC), se tiene que para un espesor de placa de 0,4
mm, la resistencia de conducción de la placa vale:
e 0,0004 4 ⎡ºC ⎤
Rcond = = = (5.7)
k · A 13,4· A 134000· A ⎢⎣ W ⎥⎦
38
P.F.C.
1 1 1
= + (5.8)
U hc h f
los fluidos caliente y frío respectivamente. Los coeficientes de convección son función
de las propiedades físicas de los fluidos y de la geometría de las placas.
1
Nu = C ·Re ·Pr a 3
(5.9)
v·Lc ·ρ
Re = (5.10)
µ
µ ·c p
Pr = (5.11)
k
h·Lc
Nu = (5.12)
k
39
P.F.C.
donde:
1
h·Lc
Nu = = C ·Re a ·Pr 3 (5.14)
k
1 1 1
= 1
+ 1
(5.16)
U ka a
kf
C c ·Re ·Pr · c
c c
3
C f ·Re ·Pr ·
f
3
f
Lc Lc
40
P.F.C.
expone a continuación.
1 1 1
= 1
+ 1
(5.16)
U a k a
kf
C c ·Re ·Pr · c
c c
3
C f ·Re ·Pr ·f
3
f
Lc Lc
1
kc
Se multiplica en los dos lados de la ecuación por Re ca ·Prc3 · , teniendo:
Lc
⎛ ⎞
1
1
kc
1 ⎜ 1 1 ⎟
·Re c ·Prc · = Re c ·Prc ·k c ·⎜
a 3 a 3
1
+ 1
⎟ (5.17)
U Lc ⎜ a 3 a 3 ⎟
⎝ C c ·Re c ·Prc ·k c C f ·Re f ·Pr f ·k f ⎠
⎛ 1 ⎞
1 kc 1
1 ⎜ a
1 k c ·Re c ·Prc3 ⎟
· ·Re ca ·Pr = + c
3
·⎜ 1
⎟ (5.18)
U Lc Cc C f ⎜ a 3 ⎟
k
⎝ f ·Re f ·Pr f ⎠
Tomando un valor inicial para el exponente “a” es posible obtener los valores de
las constantes Cc y Cf , ya que la ecuación (5.18) presenta la forma de la ecuación de una
B B B B
recta:
1 1
y= + ·x (5.19)
Cc C f
41
P.F.C.
1 1 1
= 1
+ 1
(5.16)
U k kf
C c ·Re ca ·Pr · c
c
3
C f ·Re af ·Pr · 3
f
Lc Lc
Se pasa el sumando del coeficiente de convección del lado frío al otro lado de la
ecuación, junto al término de la transferencia de calor:
1 Lc Lc
− 1
= 1
(5.20)
U
C f ·Re af ·Pr ·k f
3
f C c ·Re ca ·Pr ·k c
c
3
⎛ ⎞ 1
⎜1 Lc ⎟ k ·Pr 3 1
⎜ − 1
⎟· c c = (5.21)
⎜U a 3 ⎟ Lc C c ·Re ca
⎝ C f ·Re f ·Pr f ·k f ⎠
1
= C c ·Re ca (5.22)
⎛ ⎞ 1
⎜1 Lc ⎟ k ·Pr 3
⎜ − 1
⎟· c c
⎜U a ⎟ Lc
⎝ C f ·Re f ·Pr f3 ·k f ⎠
⎛ ⎞
⎜ ⎟
⎜ ⎟
⎜ ⎟
⎜
Ln⎜
1 ⎟
⎟ (
= Ln C c ·Re ca ) (5.23)
⎛ ⎞ 1
⎜⎜ 1 Lc ⎟ k ·Pr 3 ⎟
⎜⎜ − ⎟· c c ⎟
⎜⎜ ⎜ U ⎟ Lc ⎟⎟
1
a 3
⎝⎝ C f ·Re f ·Pr f ·k f ⎠ ⎠
42
P.F.C.
⎛ ⎞
⎜ ⎟
⎜ ⎟
⎜ ⎟
⎜ ⎟
⎟ = Ln(C c ) + a·Ln(Re c )
1
Ln⎜ (5.24)
⎛ ⎞ 1
⎜⎜ 1 Lc ⎟ k ·Pr 3 ⎟
⎜⎜ − ⎟· c c ⎟
⎜⎜ ⎜ U ⎟ Lc ⎟⎟
1
a 3
⎝⎝ C f ·Re f ·Pr f ·k f ⎠ ⎠
introducen en la ecuación (5.24), de forma que se obtiene un nuevo valor de “a”. Este
nuevo valor de “a” se reintroduce de nuevo en la ecuación (5.18) hasta conseguir
converger al método de cálculo.
43
P.F.C.
Se tendrán cinco temperaturas distintas de salida del fluido frío, de todas ellas
solo para 95 ºC se toman datos puesto que para el resto ya estaban tomados (referencia
[8] de la bibliografía). Las temperaturas de salida del generador del fluido frío son las
siguientes:
- 55 ºC
- 65 ºC
- 75 ºC
- 85 ºC
- 95 ºC
- 2,5 kg/min
- 5 kg/min
- 7,5 kg/min
- 9 kg/min
- 2,5 kg/min
- 5 kg/min
- 7,5 kg/min
- 8,5 kg/min
44
P.F.C.
- 4,5 Hz
- 8 Hz
- 12 Hz
- 15 Hz
- 4,5 Hz
- 8 Hz
- 12 Hz
- 15 Hz
- 20 Hz
- 25 Hz
- 30 Hz
- 35 Hz
45
P.F.C.
5.1.1.3. IMPLEMENTACIÓN
(∆T1 − ∆T2 )
∆Tlm = (5.4)
⎛ ∆T ⎞
ln⎜⎜ 1 ⎟⎟
⎝ ∆T2 ⎠
46
P.F.C.
Se vuelve a exportar de nuevo datos de una hoja Excel a Mathcad (ANEXO B).
En este caso se exporta parte de una hoja de cálculo con las propiedades físicas de cada
fluido calculado para cada experimento a la temperatura media, donde:
Al igual que en el ANEXO A, la fila que esta en color rojo se elimina por las
razones que se comentaron anteriormente.
Para calcular el caudal del aceite [kg/s], no se tiene como en el caso del lado del
fluido frío un caudalímetro que muestre el caudal de cada experimento. Con lo que se
calculará a partir de la potencia de las resistencias del aceite (Qc) ofrecida por el
vatímetro, sustituyendo en la expresión (5.1):
• • Q
Q = mc ·Cp c ·(Tce − Tcs ) ⇒ mc = (5.1)
Cp c ·(Tce − Tcs )
Para poder calcular la velocidad en el canal del fluido frío y caliente y los
números de Reynods y Prandtl, es necesario definir la geometría del intercambiador
(figura 25) con los datos facilitados por el fabricante.
47
P.F.C.
Figura 25. Datos geometría del intercambiador necesarios para la realización del método de Wilson
Una vez definida la geometría del intercambiador ya se puede calcular todos los
datos que serán necesarios para el método de Wilson.
⎛ • ⎞
⎜ m ⎟
⎜⎜ nº ⎟
canales ⎟
v= ⎝ ⎠ (5.26)
(ρ ·b·w)
Acto seguido se calcula los números adimensionales de Reynolds y de Prandtl
mediante las expresiones (5.10) y (5.11) respectivamente.
v·Lc ·ρ
Re = (5.10)
µ
µ ·c p
Pr = (5.11)
k
48
P.F.C.
Q
U= (5.2)
Atotal ·∆Tlm
• •
Q = mc ·Cpc ·(Tce − Tcs ) = m f ·Cp f ·(T fe − T fs ) (5.1)
(∆T1 − ∆T2 )
∆Tlm = (5.4)
⎛ ∆T ⎞
ln⎜⎜ 1 ⎟⎟
⎝ ∆T2 ⎠
49
P.F.C.
1 1
y= + ·x (5.19)
Cc C f
Se vuelve a tratar a la ecuación como una recta, definiendo los valores de “x” e
“y”. En esta segunda ecuación se introducen los valores de Cf y Cc, calculados en la
B B B B
primera ecuación. Con ello tendremos definidos los valores de “x” e “y”, y mediante la
función “slope” de Mathcad, se tiene la pendiente de la recta que en este caso es nuestra
“a”, como ya se vio en la expresión (5.25).
50
P.F.C.
Como el valor de “a” es 0.85, el siguiente paso sería introducir ese nuevo valor
en la primera ecuación del método de Wilson y repetir todo el proceso hasta conseguir
que el valor de “a” converja. Otra forma de hacerlo sería mediante un programa que
fuese un bucle y no terminase hasta que converja, figura 28.
- a = 0.788
- Cf = 0.116
B B
- Cc = 0.164
B B
1 1 1
= 1
+ 1
(5.16)
U a k a
kf
C c ·Re ·Pr · c
c c
3
f
3
C f ·Re ·Pr ·
f
Lc Lc
51
P.F.C.
En Mathcad se denomina:
Como se puede observar los resultados deben de coincidir obteniendo una recta,
la mayoría de los datos se encuentran dentro de un rango de error del 10 %, no obstante
hay algún dato que parece no cumplir esto. Por ello se aplicara un criterio estadístico a
los datos en el siguiente apartado para poder eliminar algún dato que pueda ser erróneo.
52
P.F.C.
xi − x
DRi = (5.27)
σ
donde:
- x = la media
- σ = la desviación estándar
coeficiente de desviación
número de medidas (n)
(DR0)
B B
2 1,15
3 1,38
4 1,54
5 1,65
7 1,8
10 1,96
15 2,13
25 2,33
50 2,57
100 2,81
300 3,14
500 3,29
(referencia [9] de la bibliografía)
53
P.F.C.
Si DRi > DR0 se podría descartar el punto, en caso contrario se podrá retener. Si
B B B B
sea inadecuada o que el proceso que está siendo investigado es muy variable.
⎛ ⎞
⎜ ⎟
⎜ ⎟
⎜ ⎟
⎜ ⎟
= Ln(C c ) + a·Ln(Re c )
1
Ln⎜ (5.24)
⎛ ⎞ 1 ⎟
⎜⎜ 1 Dh ⎟ k ·Pr 3 ⎟
⎜⎜ − ⎟· c c ⎟
⎜⎜ ⎜ U ⎟ Lc ⎟⎟
1
a
⎝⎝ C f ·Re f ·Pr f3 ·k f ⎠ ⎠
Como ya se explicó la ecuación (5.24) se ajusta a una recta, por tanto se aplicará
el criterio de Chauvenet para los datos de esta recta (figura 30) para comprobar que no
existe ningún dato erróneo.
Figura 30. Gráfica de los datos experimentales de la ecuación (5.24) en su última iteración
54
P.F.C.
Para calcular DRi (5.27), las variables de la expresión para el caso que se estudia
B B
xi − x
DRi = (5.27)
σ
donde:
- xi = y2
B B
- x = f(x2)
∑ (x )
n
i −x
- σ= i =1
(5.29)
n
Por lo que se calcula DRi para cada punto de la recta y se compara con el valor
B B
Existiendo un punto que puede ser rechazado(señalado con un circulo rojo en la figura
32). Dicho punto corresponde con el punto señalado en azul en el ANEXO A y B.
Como puede observarse en el ANEXO A, se trata del experimento realizado de forma
consecutiva con un experimento que resultó erróneo. Por otro lado comparando con los
mismos experimentos para los cuales se mantiene la temperatura de salida del agua
constante a 95 ºC y el caudal de agua constante a 2,5 Kg/min, variando solo el caudal de
aceite, la potencia para este experimento se aleja bastante del resto, cuando debería ser
estable. Por ello se considera un error y se descarta este resultado.
Se vuelve a realizar de nuevo todo el proceso del criterio de Chauvenet para los
94 datos restantes. Rechazándose otros seis puntos más, la eliminación de estos puntos
aparecen señalados en color verde en el ANEXO A y B. En la tabla 5 se resume el
criterio de Chauvenet con los puntos descartados en el orden de descarte, así como la
desviación DR por la que son descartados.
55
P.F.C.
Los siete puntos descartados tienen como caudal de agua 2,5 Kg/min. esto es
lógico ya que para menor caudal es más difícil que el proceso sea estable,
produciéndose un error más alto. También destacar que de los siete puntos descartados,
cuatro son para los experimentos tomados con temperatura de salida del agua constante
a 95 ºC y caudal de agua constante a 2,5 Kg/min., que como ya se vio anteriormente se
tuvo un error de medición para estas variables y el caudal de aceite de 25 Hz. Por tanto
al aplicar Chauvenet se descartan más datos realizados manteniendo dichas variables
constantes y variando el caudal del aceite, posiblemente por no tener una potencia del
vatímetro constante.
x2 vs. y2
3,5
3
2,5
2
y2
1,5
1
0,5
0
3 3,5 4 4,5 5 5,5 6 6,5
x2
apartado 5.2.2. IMPLEMENTACIÓN, variando solo los datos que pasan de 95 a 88.
Los nuevos valores de “a”, “Cf” y “Cc” son: B B B B
- a = 0.761 B
- Cf = 0.145
B B B
- Cc = 0.182
56
P.F.C.
1 1 1
= 1
+ 1
(5.16)
U a k a
kf
C c ·Re ·Pr · c
c c
3
f
3
C f ·Re ·Pr ·
f
Lc Lc
Se recuerda que:
Figura 32. Resultados teóricos vs. Experimentales, para los datos una vez aplicado el criterio de
Chauvenet
57
P.F.C.
h = f ( x1 , x 2 ,..., x n ) (5.30)
2
⎛ ∂h
n
⎞
∆h = ∑ ⎜⎜ ·∆xi ⎟⎟ (5.31)
i =1 ⎝ ∂x i ⎠
Son dadas por el fabricante y serán de 0,1 ºC tanto para el fluido frío (agua)
como para el caliente (aceite).
∆T fe = ±0.248 %
∆T fs = ±0.182 %
∆Tce = ±0.181 %
∆Tcs = ±0.195 %
58
P.F.C.
Te + Ts
Tm = (5.32)
2
⎛ ∂T ⎞ ⎛ ∂T
2
⎞
2
(∆Te )2 + (∆Ts )2
∆Tm = ⎜⎜ m ·∆Te ⎟⎟ + ⎜⎜ m ·∆Ts ⎟⎟ = (5.33)
⎝ ∂Te ⎠ ⎝ ∂Ts ⎠ 2
∆Tmf = ±0.148 %
∆Tmc = ±0.13 %
Las propiedades físicas del agua y del aceite, se expresan como una función que
solo depende de la temperatura media (5.34).
Cada propiedad tendrá una ecuación diferente por lo que se tendrán distintas
incertidumbres para cada propiedad y para cada fluido.
- Incertidumbre densidad
2
⎛ ∂ρ ⎞
∆ρ = ⎜⎜ ·∆Tm ⎟⎟ (5.35)
⎝ ∂Tm ⎠
∆ρ f = ±0.005 %
∆ρ c = ±0.0046 %
59
P.F.C.
- Incertidumbre viscosidad
2
⎛ ∂µ ⎞
∆µ = ⎜⎜ ·∆Tm ⎟⎟ (5.36)
⎝ ∂Tm ⎠
∆µ f = ±0.65 %
∆µ c = ±0.217 %
- Incertidumbre Cp
2
⎛ ∂Cp ⎞
∆Cp = ⎜⎜ ·∆Tm ⎟⎟ (5.37)
⎝ ∂Tm ⎠
∆Cp f = ±0.0047 %
∆Cp c = ±0.0145 %
2
⎛ ∂k ⎞
∆k = ⎜⎜ ·∆Tm ⎟⎟ (5.38)
⎝ ∂Tm ⎠
∆k f = ±0.0487 %
∆k c = ±0.0037 %
• Incertidumbre caudal
•
∆ m f = ±4.805 %
• Q • Q
mc = = mc = (5.1)
Cp c ·(Tce − Tcs ) Cp c ·∆Tc
60
P.F.C.
2 2 2
⎛ • ⎞ ⎛ • ⎞ ⎛ • ⎞
•
⎜ ∂ mc ⎟ ⎜ ∂ mc ⎟ ⎜ ∂ mc ⎟
∆ mc = ⎜ ·∆Q ⎟ + ⎜ ·∆Cp c ⎟ + ⎜ ·∆Tc ⎟ (5.39)
⎜ ∂Q ⎟ ⎜ ∂Cp c ⎟ ⎜ ∂∆Tc ⎟
⎝ ⎠ ⎝ ⎠ ⎝ ⎠
donde:
•
∂ mc 1
= (5.40)
∂Q Cp c ·∆Tc
•
∂ mc −Q
= (5.41)
∂Cp c Cp c 2 ·∆Tc
•
∂ mc −Q
= (5.42)
∂∆Tc Cp c ·∆Tc 2
- Incertidumbre potencia
∆Q = ±0.5%
- Incertidumbre Cpc
∆Cp c = ±0.0145 %
61
P.F.C.
- Incertidumbre ∆Tc
2 2
⎛ ∂∆Tc ⎞ ⎛ ∂∆Tc ⎞
∆(∆Tc ) = ⎜⎜ ·∆Tce ⎟⎟ + ⎜⎜ ·∆Tcs ⎟⎟ = (∆Tce )2 + (∆Tcs )2 (5.44)
⎝ ∂Tce ⎠ ⎝ ∂Tcs ⎠
∆(∆Tc ) = ±8.115 %
•
∆ m c = ±8.13 %
• Incertidumbre velocidad
⎛ • ⎞
⎜ m ⎟
⎜⎜ nº ⎟
canales ⎟
v= ⎝ ⎠ (5.26)
(ρ ·b·w)
2
⎛ ∂v • ⎞ ⎛ ∂v ⎞
2
∆v = ⎜ ⎟
·∆ m + ⎜⎜ ·∆ρ ⎟⎟ (5.45)
⎜ • ⎟ ⎝ ∂ρ ⎠
⎝∂m ⎠
donde:
∂v 1
= (5.46)
∂m
•
ρ ·b·w·n º canales
•
∂v −m
= 2 (5.47)
∂ρ ρ ·b·w·n º canales
62
P.F.C.
∆v f = ±4.805 %
∆vc = ±8.13 %
v· L c · ρ
Re = (5.10)
µ
2 2
⎛ ∂ Re ⎞ ⎛ ∂ Re ⎞ ⎛ ∂ Re ⎞
2
donde:
∂ Re ρ ·Lc
= (5.49)
∂v µ
∂ Re v·Lc
= (5.50)
∂ρ µ
∂ Re − ρ ·v·Lc
= (5.51)
∂µ µ2
∆ Re f = ±4.84 %
∆ Re c = ±8.133 %
63
P.F.C.
µ ·c p
Pr = (5.11)
k
2 2
⎛ ∂ Pr ⎞ ⎛ ∂ Pr ⎞ ⎛ ∂ Pr
2
⎞
∆ Pr = ⎜⎜ ·∆µ ⎟⎟ + ⎜⎜ ·∆Cp ⎟⎟ + ⎜ ·∆k ⎟ (5.52)
⎝ ∂µ ⎠ ⎝ ∂Cp ⎠ ⎝ ∂k ⎠
donde:
∂ Pr Cp
= (5.53)
∂µ k
∂ Pr µ
= (5.54)
∂Cp k
∂ Pr − Cp·µ
= (5.55)
∂k k2
∆ Pr f = ±0.652 %
∆ Prc = ±0.218 %
⎛ ∂Nu ⎞ ⎛ ∂Nu
2 2
⎞
∆Nu = ⎜ ·∆ Re ⎟ + ⎜ ·∆ Pr ⎟ (5.56)
⎝ ∂ Re ⎠ ⎝ ∂ Pr ⎠
64
P.F.C.
donde:
∂Nu
1
∆Nu f = ±3.69 %
∆Nu c = ±6.19 %
Nu·k
h= (5.12)
Lc
⎛ ∂h ⎞ ⎛ ∂h
2 2
⎞
∆h = ⎜ ·∆Nu ⎟ + ⎜ ·∆k ⎟ (5.59)
⎝ ∂Nu ⎠ ⎝ ∂k ⎠
∂h k
= (5.60)
∂Nu Lc
∂h Nu
= (5.61)
∂k Lc
∆h f = ±3.691 %
∆hc = ±6.19 %
65
P.F.C.
aceite -3
2,8·10 – 8·10
P
-3
P
± 0,217 P P
v (m/s)
aceite 4,3·10-2 - 0,43 P ± 8,13P
66
P.F.C.
Para nuestro caso la pérdida de carga se medirá en el lado del aceite, colocando
un diferencial de presión conectado en los puertos de entrada y salida del aceite.
∆Pfric ·Lc
f = (5.62)
2·ρ ·v 2 ·H
donde:
Lc = 2·b (5.13)
67
P.F.C.
donde:
∆Pcol = ρ ·g ·H (5.64)
ρ ·v 2 i
∆Psin g = ·∑ K i (5.66)
2 1
68
P.F.C.
donde:
- K = coeficiente de pérdida
Se tendrá un codo de 90º normal y acoplado, cuyo diámetro interior es 22.4 mm.
(una pulgada son 25.4 mm). Se debe interpolar, pero el dato más pequeño de la tabla
mostrada en la figura 33 es para 1 pulgada. Por lo que el coeficiente de pérdida será 0.5
(Kc = 0.5).
69
P.F.C.
2
⎛ D2 ⎞
K eb = ⎜⎜1 − i2 ⎟⎟ (5.67)
⎝ D ⎠
donde:
⎛ D2 ⎞
K cb = 0.42·⎜⎜1 − i2 ⎟⎟ (5.68)
⎝ D ⎠
donde:
Con lo que la expresión final de las pérdidas singulares para nuestro caso será
según la expresión (5.66) la siguiente:
ρ ·v 2
∆Psin g = ·(K c + K eb + K cb ) (5.69)
2
ρ ·v 2 ⎛⎜ ⎛ D2 ⎞⎞
2
⎛ Di2 ⎞
∆Psin g = · K c + ⎜⎜1 − 2 ⎟⎟ + 0.42·⎜⎜1 − i2 ⎟⎟ ⎟ (5.70)
2 ⎜ ⎝ D ⎠ ⎝ D ⎠ ⎟⎠
⎝
⎛ ρ ·v 2 ⎛ ⎛ Di2 ⎞
2
⎛ D2 ⎞⎞⎞
∆Pfrie = ∆Pmed − ⎜ · K c + ⎜⎜1 − 2 ⎟⎟ + 0.42·⎜⎜1 − i2 ⎟⎟ ⎟ ⎟
⎜ (5.71)
⎜ 2 ⎜ ⎟
⎝ ⎝ ⎝ D ⎠ ⎝ D ⎠ ⎟⎠ ⎠
70
P.F.C.
Una vez se tenga calculado los resultados obtenidos se ajustan a una ecuación de
la forma:
f = C ·Re b (5.72)
Para este caso se tiene el problema que no existe un Coriolis que muestre con
exactitud el caudal de aceite, con lo que se desconoce dicho caudal. Tampoco puede
hacerse con la expresión de la potencia (5.1) ya que no existe transferencia de calor.
• •
Q = mc ·Cpc ·(Tce − Tcs ) = m f ·Cp f ·(T fe − T fs ) (5.1)
71
P.F.C.
72
P.F.C.
5.3.2. IMPLEMENTACIÓN
Al igual que para el cálculo del coeficiente de convección, el cálculo del factor
de fricción de Fanning se realizará en el programa matemático Mathcad. Se explicará
durante este apartado paso por paso la realización del programa, que se encuentra
íntegro en el ANEXO F.
73
P.F.C.
Figura 35. Datos geométricos, variables experimentales y propiedades físicas definidos en Mathcad
Se puede observar que para pasar los datos de perdida de carga medidos,
llamado “∆Pbar”, se multiplica por 105, pasando de bar a unidades del sistema
P P
⎛ • ⎞
⎜ m ⎟
⎜⎜ nº ⎟
canales ⎟
v= ⎝ ⎠ (5.26)
(ρ ·b·w)
74
P.F.C.
75
P.F.C.
∆Pfric ·Lc
f = (5.62)
2·ρ ·v 2 ·H
76
P.F.C.
v·Lc ·ρ
Re = (5.10)
µ
77
P.F.C.
No existe un método para hallar una correlación del factor de fricción como era
el caso del método de Wilson para el cálculo de la correlación del Nusselt.
La forma de proceder para calcular la correlación del factor de fricción, será por
zonas, es decir, como se ha podido observar en la figura 39 se tendrán tres zonas
diferenciadas por tanto cada zona se ajustará a un modelo de curva exponencial,
expresión (5.72), cada una con sus valores de “C” y “b”.
f = C ·Re b (5.72)
Por tanto se introducirán en Mathcad tres rangos distintos de valores, cada uno
correspondiente a la zona donde se encuentre según el Reynolds, y se calcularán para
cada rango el valor del coeficiente de fricción del modo explicado en el punto 5.3.2.
IMPLEMENTACIÓN. Una vez hallado el factor de fricción se ajustan los datos
mediante la función “genfit”.
Los valores de “C” y “b” para cada rango de Reynolds son los siguientes:
- C = 23,395
- b = -0,531
- C = 10,98
- b = -0,33
78
P.F.C.
- C = 3,251
- b = -0,107
Figura 40. Comparación correlaciones experimentales del factor de fricción vs. datos
experimentales del factor de fricción
79
P.F.C.
Se aplicará Chauvenet para cada zona, puesto que cada una presenta una
correlación diferente.
Observando la figura 40, se trata de la zona que presenta un ajuste algo menor
por lo que parece que pudiese existir algún dato erróneo.
Se calcula DRi (5.27) para cada dato experimental del factor de fricción
perteneciente a la zona de flujo laminar.
xi − x
DRi = (5.27)
σ
donde:
∑ (x )
n
i −x
- σ= i =1
(5.29)
n
Por lo que se calcula DRi para cada punto y se compara con el valor de DR0,
calculado mediante interpolación en la tabla 4 (para n = 51, DR0 = 2.5748). Siendo
todos los puntos válidos. El punto que presenta mayor DRi (señalado en la figura 41 con
un círculo negro) será de 2,227.
80
P.F.C.
Se realizará del mismo modo que en el caso de la zona laminar. Para la zona de
transición se tendrán solo 12 datos experimentales, por lo que el valor de DR0 para
n = 12, será de 2,028 (tabla 4).
El número de medidas experimentales del factor de fricción para Re > 200 son
las restantes, n = 77 (DR0 = 2,6996, tabla 4).
81
P.F.C.
∆(∆Pmed ) = ±3.3 %
Te + Ts
Tm = (5.32)
2
⎛ ∂T
2
⎞ ⎛ ∂T ⎞
2
(∆Te )2 + (∆Ts )2
∆Tm = ⎜⎜ m ·∆Te ⎟⎟ + ⎜⎜ m ·∆Ts ⎟⎟ = (5.33)
⎝ ∂Te ⎠ ⎝ ∂Ts ⎠ 2
∆Tm = ±0.559 %
82
P.F.C.
∆ρ = ±0.0047 %
∆µ = ±0.933 %
• Incertidumbre caudal
•
∆ m = ±3.011 %
2
⎛ ∂v • ⎞ ⎛ ∂vtubo ⎞
2
∆vtubo = ⎜ tubo ⎟
·∆ m + ⎜⎜ ·∆ρ ⎟⎟ (5.77)
⎜ • ⎟ ⎝ ∂ρ ⎠
⎝ ∂m ⎠
donde:
∂v 1
•
= (5.78)
∂m ⎛ Di2 ⎞
ρ ·⎜⎜ π · ⎟⎟
⎝ 4 ⎠
•
∂v −m
= (5.79)
∂ρ ⎛ D2 ⎞
ρ 2 ·⎜⎜ π · i ⎟
⎟
⎝ 4 ⎠
∆vtubo = ±3.011 %
83
P.F.C.
⎛ • ⎞
⎜ m ⎟
⎜⎜ nº ⎟
canales ⎟
v= ⎝ ⎠ (5.26)
(ρ ·b·w)
2
⎛ ∂v • ⎞ ⎛ ∂v ⎞
2
∆v = ⎜ ⎟
·∆ m + ⎜⎜ ·∆ρ ⎟⎟ (5.45)
⎜ • ⎟ ⎝ ∂ρ ⎠
⎝ ∂ m ⎠
donde:
∂v 1
= (5.46)
∂m
•
ρ ·b·w·n º canales
•
∂v −m
= 2 (5.47)
∂ρ ρ ·b·w·n º canales
∆v = ±3.011 %
ρ ·vtubo 2
∆Psin g = ·(K c + K eb + K cb ) (5.69)
2
⎛ ∂∆Psin g ⎞ ⎛ ∂∆Psin g
2 2
⎞
∆ (∆Psin g ) = ⎜⎜ ·∆ρ ⎟⎟ + ⎜⎜ ·∆vtubo ⎟⎟ (5.80)
⎝ ∂ρ ⎠ ⎝ ∂vtubo ⎠
donde:
∂∆Psin g 2
vtubo
= ·( K c + K eb + K cb ) (5.81)
∂ρ 2
84
P.F.C.
∂∆Psin g
= ρ ·vtubo ·( K c + K eb + K cb ) (5.82)
∂vtubo
∆ ( ∆Psin g ) = ±6.023 %
2
⎛ ∂∆Pfric ⎞ ⎛ ∂∆Pfric ⎞
2
∆ ( ∆Pfric ) = ±3.436 %
v·Lc ·ρ
Re = (5.10)
µ
2 2
⎛ ∂ Re ⎞ ⎛ ∂ Re ⎞ ⎛ ∂ Re ⎞
2
donde:
∂ Re ρ ·Lc
= (5.49)
∂v µ
85
P.F.C.
∂ Re v·Lc
= (5.50)
∂ρ µ
∂ Re − ρ ·v·Lc
= (5.51)
∂µ µ2
∆ Re = ±3,092 %
En último lugar una vez calculadas las incertidumbres de las variables que son
necesarias para calcular el factor de fricción, ya se esta en disposición de calcular la
incertidumbre de dicho factor.
∆Pfric ·Lc
f = (5.62)
2·ρ ·v 2 ·H
2
⎛ ∂f ⎞ ⎛ ∂f
2
⎞ ⎛ ∂f
2
donde:
∂f Lc
= (5.86)
∂∆Pfric 2·ρ ·v 2 ·H
∂f − ∆Pf ·Lc
= (5.87)
∂ρ 2·ρ 2 ·v 2 ·H
∂f − ∆Pf ·Lc
= (5.88)
∂v ρ ·v·H
∆f = ±3,437 %
86
P.F.C.
87
P.F.C.
5.4. RESULTADOS
En este apartado se realiza un breve resumen de los resultados hallados
anteriormente y los rangos de experimentación, para poder comparar con otros autores
en el siguiente capítulo.
- a = 0.761
- Cf = 0.145
- Cc = 0.182
Correlaciones:
Nu = 0.145·Re 0.761
f ·Pr 3
f 95 < Ref < 559 1,99 < Prf < 3,7 (5.89)
1
Nu = 0.182·Re 0c .761 ·Prc3 28 < Rec < 557 53 < Prc < 133 (5.90)
f = C ·Re b (5.72)
88
P.F.C.
89
P.F.C.
90
P.F.C.
Nu = 0.145·Re 0.761
f ·Pr 3
f 95 < Ref < 559 1.99 < Prf < 3.7 (5.89)
1
Nu = 0.182·Re 0.761
c ·Pr
c
3
28 < Rec < 557 53 < Prc < 133 (5.90)
donde:
91
P.F.C.
Troupe (agua)
35
30
25
20 Troupe
Nu
15 P.F.C.
10
5
0
0 200 400 600
Re
92
P.F.C.
donde:
Tentrada + Tsalida
Tm = (6.3)
2
Tmc + Tmf
Tw = (6.4)
2
93
P.F.C.
C1 = 0.348
m = 0.663
Kumar (agua)
35
30
25
20 Kumar
Nu
15 P.F.C.
10
5
0
0 200 400 600
Re
El error para los Reynolds más bajos del rango (Re = 95), el error es de
aproximadamente el 35 %. A medida que aumenta el Reynolds el error disminuye hasta
el 23 % para los valores más altos del rango experimental del número de Reynolds
(Re = 560).
94
P.F.C.
Nu = Nu + Nu 3 3 3
·Pr ⎜⎜
3 ⎟⎟ 1 ≤ Re ≤ 10 4 20 º ≤ β ≤ 62 º (6.5)
⎝ µw
1 t
⎠
donde:
Tentrada + Tsalida
Tm = (6.3)
2
Tmc + Tmf
Tw = (6.4)
2
donde:
β = ángulo de “chevron”
95
P.F.C.
donde:
Wanniarachchi (agua)
30
25
20
Wanniarachchi
Nu
15
P.F.C.
10
5
0
0 200 400 600
Re
96
P.F.C.
Wanniarachchi (aceite)
100
80
60 Wanniarachchi
Nu
40 P.F.C.
20
0
0 200 400 600
Re
Nu = C1 ·Re ·Pr
m 3
(6.9)
donde:
97
P.F.C.
C1 = 0.2946
m = 0.7
Thonon (agua)
35
30
25
20 Thonon
Nu
15 P.F.C.
10
5
0
0 200 400 600
Re
Para el caso del aceite el rango queda limitado para términos bajos de Reynolds
por la correlación de Thonon, no llegando a cubrir todo el rango experimental del aceite.
La comparación se muestra en la figura 47.
98
P.F.C.
Thonon (aceite)
100
80
60 Thonon
Nu
40 P.F.C.
20
0
0 200 400 600
Re
Como ya ocurrió con otras correlaciones, esta correlación está realizada teniendo
como fluido de trabajo el agua. Por tanto no se comparará con el aceite.
El exponente del Prandtl para el caso de Buonopane será mayor (0.4) que el
tomado experimentalmente (1/3). El coeficiente “C” (0.2536) será también mayor que el
obtenido experimentalmente para el fluido frío (0.145). Esto se equiparará con un
coeficiente “a” que será menor para la correlación de Buonopane (0.65) que para la
correlación experimental (0.761).
99
P.F.C.
Buonopane (agua)
25
20
15 Buonopane
Nu
10 P.F.C.
0
0 200 400 600
Re
Será una buena comparación, dentro del pequeño rango de Reynolds donde se
puede comparar. En esta comparación al contrario que en el resto hechas para el agua,
será la correlación hallada experimentalmente la que presente un mayor valor de Nusselt
que la correlación del autor (Buonopane en este caso). El error oscila entre el 5 % para
los valores del rango de Reynolds bajos (Re = 400) y el 8 % para Re = 557 el valor más
alto del rango de comparación.
Nu = 0.78·Re ·Pr 0 .5 3
50 ≤ Re ≤ 20000 (6.11)
100
P.F.C.
25
20
Maslov
15 and
Nu
Kovalenk
10 o
P.F.C.
5
0
0 200 400 600
Re
Figura 49. Comparación Nu vs. Re, para la correlación de Maslov y Kovalenko (agua)
Como para el resto de comparaciones con las otras correlaciones con el agua
(excepto Buonopane), la distribución es semejante, pero será de un rango mayor para la
correlación de Maslov y Kovalenko. El error varía entre un 40 % para valores pequeños
de Reynolds y un 4 % para valores grandes de Reynolds (dentro del rango
experimental).
100
80
Maslov
60 and
Nu
Kovalenko
40
P.F.C.
20
0
0 200 400 600
Re
Figura 50. Comparación Nu vs. Re, para la correlación de Maslov y Kovalenko (aceite)
101
P.F.C.
El rango superior queda limitado a Re > 400, puesto que es el máximo que
permite la correlación de Cooper (6.12).
Cooper (agua)
25
20
15 Cooper
Nu
10 P.F.C.
0
0 200 400 600
Re
102
P.F.C.
Cooper (aceite)
100
80
60 Cooper
Nu
40 P.F.C.
20
0
0 200 400 600
Re
Se resume en la tabla 12, los rangos para los cuales se compara y el rango de
error de las comparaciones realizadas para cada autor.
103
P.F.C.
104
P.F.C.
C2
f = (6.13)
(Re) p
donde:
105
P.F.C.
• Re < 10
donde:
- C2 = 50
- p =1
donde:
- C2 = 19.40
- p =0.589
• Re > 100
donde:
- C2 = 2.99
- p =0.183
106
P.F.C.
Figura 53. Comparación factor de fricción vs. Re, para la correlación de Kumar
[ ]
1
f = f13 + f t 3 3 1 ≤ Re ≤ 10 4 (6.14)
donde:
⎛ Acorrugada ⎞ −1
f1 = 1774·( β ) −1.026 ·⎜⎜ ⎟·Re (6.15)
⎝ A ⎟⎠
donde:
107
P.F.C.
donde:
Sustituyendo cada dato de las expresiones anteriores por los del intercambiador
que se caracteriza, se llega a comparar en la figura 54 los resultados obtenidos
experimentalmente con los ofrecidos por la correlación de Wanniarachchi para el factor
de fricción de Fanning.
Figura 54. Comparación factor de fricción vs. Re, para la correlación de Wanniarachchi
108
P.F.C.
f = C 2 ·Re − p (6.18)
donde:
• Re ≤ 160
donde:
- C2 = 45.57
109
P.F.C.
- p = 0.67
• Re > 160
donde:
- C2 = 3.7
- p = 0.172
Figura 55. Comparación factor de fricción vs. Re, para la correlación de Thonon
110
P.F.C.
Figura 56. Comparación del factor de fricción experimental en función del Reynolds con otros
autores
Como ya se vio la correlación que peor se ajusta a los datos experimentales será
la correlación de Wanniarachchi (de color negro en la figura 56). La que mejor se ajusta
a los datos obtenidos experimentalmente será la correlación de Thonon (de color marrón
en la figura 56), de la que se puede apreciar también como el cambio de zona laminar a
transición y turbulento se ajusta más a los datos del factor de fricción de Fanning
hallado experimentalmente que para las otras correlaciones. Esto se debe a que la
correlación de Thonon propone un cambio para Re = 160 y Kumar para Re = 100,
observándose que el cambio para el factor de fricción experimental se encuentra
aproximadamente alrededor de Re = 200, estando más próximo por tanto a la
correlación de Thonon.
111
P.F.C.
7.1. CONCLUSIONES
Como resultado de este proyecto se han extraído las conclusiones que se detallan
a continuación:
1 1 1
= + Ren _ c + Rcond + Ren _ f + (5.6)
U · A hc · A h f ·A
e 0,0004 4 ⎡ºC ⎤
Rcond = = = (5.7)
k · A 13,4· A 134000· A ⎢⎣ W ⎥⎦
112
P.F.C.
35
30
25
Troupe
20
Nu
P.F.C.
15
Buonopane
10
5
0
0 200 400 600
Re
113
P.F.C.
114
P.F.C.
8. BIBLIOGRAFÍA
[1] Mercedes de Vega Blázquez: “Refrigeración por absorción”. Departamento de
Ingeniería Térmica y de Fluidos. Universidad Carlos III de Madrid, 2005.
[4] Björn Palm y Joachim Claesson: “Plate Heat Exchangers: Calculation Methods for
Single and Two-Phase Flow”. Heat Transfer Engineering, 2006.
[9] Antonio Lecuona Neumann y Alejandro Calvo Carnero: “Cursos sobre aplicaciones
del ordenador a la medida y análisis de resultados experimentales y de ensayos”.
[11] F.S.K. Warnakulasuriya y W.M. Worek: “Heat transfer and pressure drop
properties of high viscous solutions in plate heat exchangers”. International Journal of
Heat and Mass Transfer, 2008.
[12] Zahid H. Ayub: “Plate Heat Exchanger Literature Survey and New Heat Transfer
and Pressure Drop Correlations for Refrigerant Evaporators”. Heat Transfer
Engineering, 2003.
115
P.F.C.
Fuentes de Internet:
[13] www.alfalaval.es
[14] www.delaval.es
[15] plantasquimicas.iespana.es
116
P.F.C.
117
P.F.C.
• •
•
Tagua [ºC] m agua m aceite tfs [ºC] tfe [ºC] tce [ºC] tcs [ºC] W [w] mf [kg/min]
[kg/min] [Hz]
4,5 75,31 58,25 78,32 68,72 3050,36 2,50
8 75,17 58,62 75,79 71,13 2962,18 2,50
2,5
12 75,03 58,71 75,34 72,39 2887,48 2,50
15 75,01 58,80 75,24 72,91 2842,65 2,50
4,5 75,50 57,28 90,50 67,15 5778,86 5,00
8 75,14 57,52 79,11 69,09 5887,76 5,00
5
12 75,14 57,88 77,14 70,96 5956,47 5,00
15 74,16 55,64 75,17 70,28 6334,04 5,00
75
4,5 74,87 54,99 99,98 65,60 10237,51 7,46
8 75,01 55,40 84,78 68,23 10077,58 7,49
7,5
12 75,06 55,75 80,69 70,19 9927,58 7,50
15 75,09 55,81 79,25 70,72 9893,94 7,50
4,5 74,92 56,49 104,92 66,70 11368,57 8,99
8 74,81 56,63 87,08 68,60 11194,77 9,01
9
12 74,96 56,84 82,35 70,26 11107,41 8,99
15 74,95 56,81 80,52 71,01 11121,15 9,00
4,5 85,02 63,93 88,08 75,92 3755,10 2,50
8 84,98 62,92 85,80 79,68 4019,30 2,50
2,5
12 85,14 62,87 85,56 81,46 3998,38 2,50
15 85,13 63,03 85,50 82,29 3962,22 2,50
4,5 85,18 62,49 99,61 73,59 7807,88 5,00
8 84,85 61,93 89,66 76,47 7894,27 5,00
5
12 85,09 62,36 87,52 79,18 7837,07 5,00
15 83,95 61,53 85,59 79,02 7782,70 5,00
85
4,5 85,21 66,63 108,12 75,42 9785,66 7,50
8 85,13 67,05 93,39 78,58 9280,03 7,50
7,5
12 84,52 65,94 89,44 79,45 9507,74 7,50
15 84,96 65,69 88,48 80,19 9860,80 7,50
4,5 84,06 65,60 113,52 74,70 11492,51 9,00
8 85,19 67,07 96,66 78,90 11110,74 9,00
9
12 84,87 66,80 91,71 79,95 11094,55 9,00
15 85,31 67,10 90,07 80,70 11052,31 9,00
118
P.F.C.
• • •
Tagua [ºC] m agua m aceite tfs [ºC] Tfe [ºC] tce [ºC] tcs [ºC] W [w] mf
[kg/min] [Hz] [kg/min]
4,5 95,02 74,15 96,39 86,81 3850,85 2,50
8 94,88 73,92 95,26 89,51 3808,08 2,50
12 95,00 74,02 95,10 91,03 3771,45 2,50
15 95,00 74,16 95,03 91,73 3767,14 2,50
2,5
20 94,83 73,97 94,87 92,30 3807,67 2,50
25 94,87 75,20 94,85 92,92 3466,45 2,50
30 94,99 75,76 94,99 93,40 3325,37 2,50
35 94,82 75,66 94,86 93,47 3304,86 2,50
4,5 95,31 77,87 102,09 86,74 6302,75 5,00
8 95,10 77,62 98,15 88,67 6367,59 5,00
12 95,13 77,64 96,61 89,92 6349,15 5,00
15 94,84 77,37 95,75 90,28 6339,60 5,00
5
20 95,25 77,65 95,74 91,53 6394,40 5,00
25 95,02 77,51 95,27 91,92 6213,70 5,00
30 95,19 77,59 95,38 92,54 6216,27 5,00
35 94,91 74,09 95,05 92,16 7154,09 5,00
95
4,5 95,10 79,86 107,27 87,04 8129,18 7,50
8 95,10 79,82 101,32 89,06 8171,04 7,50
12 95,08 79,75 98,47 89,91 8152,73 7,50
15 94,91 79,55 97,20 90,19 8170,57 7,50
7,5
20 94,98 79,55 96,31 90,93 8153,21 7,50
25 95,27 79,70 96,14 91,76 8196,51 7,50
30 95,20 79,49 95,85 92,09 8296,88 7,50
35 95,30 79,30 95,74 92,47 8213,41 7,50
4,5 95,11 78,06 111,79 86,45 10316,23 8,50
8 94,94 78,16 103,71 88,58 10164,73 8,50
12 95,40 78,70 100,43 89,90 10086,41 8,50
15 95,02 78,49 98,45 90,03 9891,02 8,50
8,5
20 94,88 77,44 97,18 90,30 10576,83 8,50
25 95,30 78,72 96,67 91,42 9896,08 8,50
30 95,13 79,24 96,04 91,81 9366,93 8,50
35 94,98 79,55 95,63 92,08 8998,05 8,50
119
P.F.C.
120
P.F.C.
ρ f [kg/m3]
P P ρ c [kg/m3]
P P µ f [Pa·s] µ c [m2/s]
P P µ c [Pa·s] Cp f [J/Kg·K] Cp c [J/Kg·K] K f [W/m·K] K c [W/m·K]
984,661 851,613 4,900E-04 7,031E-06 5,987E-03 4183,454 1975,778 0,651 0,115
984,663 851,549 4,900E-04 7,009E-06 5,969E-03 4183,452 1976,247 0,651 0,115
984,445 846,204 4,865E-04 5,521E-06 4,672E-03 4183,630 2015,509 0,651 0,114
984,501 850,038 4,874E-04 6,529E-06 5,550E-03 4183,584 1987,350 0,651 0,115
984,475 850,770 4,870E-04 6,755E-06 5,747E-03 4183,605 1981,969 0,651 0,115
984,436 850,931 4,864E-04 6,806E-06 5,792E-03 4183,637 1980,786 0,652 0,115
984,106 844,973 4,811E-04 5,250E-06 4,436E-03 4183,912 2024,552 0,652 0,114
984,345 849,467 4,849E-04 6,361E-06 5,403E-03 4183,712 1991,541 0,652 0,115
983,982 850,063 4,792E-04 6,536E-06 5,556E-03 4184,016 1987,164 0,652 0,115
983,878 850,457 4,776E-04 6,657E-06 5,661E-03 4184,105 1984,271 0,653 0,115
979,564 846,476 4,201E-04 5,584E-06 4,727E-03 4188,226 2013,510 0,661 0,114
979,500 846,511 4,194E-04 5,592E-06 4,734E-03 4188,293 2013,256 0,661 0,114
979,515 846,290 4,196E-04 5,541E-06 4,689E-03 4188,277 2014,878 0,661 0,114
979,496 846,175 4,193E-04 5,514E-06 4,666E-03 4188,297 2015,724 0,661 0,114
979,776 843,601 4,226E-04 4,971E-06 4,194E-03 4188,004 2034,628 0,660 0,114
979,807 846,164 4,229E-04 5,512E-06 4,664E-03 4187,972 2015,805 0,660 0,114
979,712 846,191 4,218E-04 5,518E-06 4,669E-03 4188,071 2015,604 0,660 0,114
980,581 846,910 4,322E-04 5,687E-06 4,816E-03 4187,182 2010,325 0,659 0,114
980,563 841,452 4,320E-04 4,580E-06 3,854E-03 4187,199 2050,413 0,659 0,114
980,418 844,858 4,302E-04 5,225E-06 4,415E-03 4187,346 2025,398 0,659 0,114
980,312 845,438 4,289E-04 5,350E-06 4,523E-03 4187,454 2021,138 0,659 0,114
980,285 845,684 4,286E-04 5,404E-06 4,570E-03 4187,481 2019,332 0,660 0,114
980,150 839,814 4,270E-04 4,314E-06 3,623E-03 4187,619 2062,446 0,660 0,114
980,141 844,134 4,269E-04 5,077E-06 4,285E-03 4187,628 2030,711 0,660 0,114
980,043 844,966 4,257E-04 5,248E-06 4,435E-03 4187,729 2024,599 0,660 0,114
121
P.F.C.
ρ f [kg/m3]
P P ρ c [kg/m3]
P P µ f [Pa·s] µ c [m2/s]
P P µ c [Pa·s] Cp f [J/Kg·K] Cp c [J/Kg·K] K f [W/m·K] K c [W/m·K]
980,053 845,259 4,258E-04 5,311E-06 4,489E-03 4187,719 2022,451 0,660 0,114
975,133 841,878 3,759E-04 4,654E-06 3,918E-03 4193,195 2047,288 0,667 0,114
975,447 841,478 3,787E-04 4,585E-06 3,858E-03 4192,823 2050,223 0,667 0,114
975,413 841,061 3,784E-04 4,514E-06 3,797E-03 4192,864 2053,284 0,667 0,114
975,369 840,852 3,780E-04 4,480E-06 3,767E-03 4192,915 2054,819 0,667 0,114
975,517 839,385 3,793E-04 4,248E-06 3,566E-03 4192,740 2065,600 0,666 0,114
975,781 841,301 3,816E-04 4,555E-06 3,832E-03 4192,430 2051,522 0,666 0,114
975,583 841,148 3,798E-04 4,529E-06 3,809E-03 4192,663 2052,648 0,666 0,114
976,168 841,713 3,851E-04 4,625E-06 3,893E-03 4191,979 2048,493 0,666 0,114
974,251 836,583 3,686E-04 3,857E-06 3,226E-03 4194,256 2086,175 0,668 0,113
974,148 839,721 3,677E-04 4,300E-06 3,611E-03 4194,382 2063,125 0,668 0,114
974,673 840,554 3,720E-04 4,431E-06 3,725E-03 4193,746 2057,007 0,667 0,114
974,617 840,614 3,716E-04 4,441E-06 3,733E-03 4193,814 2056,573 0,668 0,114
974,919 835,316 3,741E-04 3,698E-06 3,089E-03 4193,451 2095,487 0,667 0,113
974,126 838,746 3,676E-04 4,154E-06 3,484E-03 4194,408 2070,289 0,668 0,113
974,307 839,804 3,690E-04 4,313E-06 3,622E-03 4194,189 2062,522 0,668 0,114
974,080 840,045 3,672E-04 4,350E-06 3,654E-03 4194,464 2060,746 0,668 0,114
968,669 836,674 3,309E-04 3,869E-06 3,237E-03 4201,416 2085,506 0,673 0,113
968,794 836,249 3,316E-04 3,814E-06 3,189E-03 4201,249 2088,633 0,673 0,113
968,720 835,880 3,312E-04 3,767E-06 3,149E-03 4201,348 2091,341 0,673 0,113
968,673 835,709 3,310E-04 3,746E-06 3,131E-03 4201,412 2092,595 0,673 0,113
968,794 835,598 3,316E-04 3,732E-06 3,119E-03 4201,249 2093,412 0,673 0,113
968,366 835,438 3,293E-04 3,713E-06 3,102E-03 4201,825 2094,591 0,674 0,113
968,133 835,268 3,280E-04 3,692E-06 3,084E-03 4202,140 2095,837 0,674 0,113
968,225 835,284 3,285E-04 3,694E-06 3,086E-03 4202,015 2095,722 0,674 0,113
122
P.F.C.
ρ f [kg/m3]
P P ρ c [kg/m3]
P P µ f [Pa·s] µ c [m2/s]
P P µ c [Pa·s] Cp f [J/Kg·K] Cp c [J/Kg·K] K f [W/m·K] K c [W/m·K]
967,301 835,148 3,238E-04 3,678E-06 3,072E-03 4203,275 2096,717 0,674 0,113
967,459 835,693 3,246E-04 3,744E-06 3,129E-03 4203,058 2092,715 0,674 0,113
967,442 835,772 3,245E-04 3,754E-06 3,137E-03 4203,081 2092,137 0,674 0,113
967,634 835,907 3,255E-04 3,771E-06 3,152E-03 4202,818 2091,142 0,674 0,113
967,397 835,571 3,243E-04 3,729E-06 3,116E-03 4203,142 2093,611 0,674 0,113
967,525 835,593 3,249E-04 3,732E-06 3,118E-03 4202,968 2093,451 0,674 0,113
967,439 835,395 3,245E-04 3,708E-06 3,097E-03 4203,086 2094,905 0,674 0,113
968,727 835,587 3,313E-04 3,731E-06 3,118E-03 4201,339 2093,491 0,673 0,113
966,684 833,663 3,208E-04 3,506E-06 2,923E-03 4204,124 2107,630 0,675 0,113
966,698 834,728 3,209E-04 3,628E-06 3,028E-03 4204,105 2099,804 0,675 0,113
966,729 835,270 3,210E-04 3,693E-06 3,084E-03 4204,061 2095,821 0,675 0,113
966,858 835,539 3,216E-04 3,725E-06 3,113E-03 4203,884 2093,850 0,675 0,113
966,834 835,579 3,215E-04 3,730E-06 3,117E-03 4203,917 2093,551 0,675 0,113
966,681 835,400 3,208E-04 3,708E-06 3,098E-03 4204,129 2094,865 0,675 0,113
966,778 835,389 3,212E-04 3,707E-06 3,097E-03 4203,994 2094,945 0,675 0,113
966,809 835,316 3,214E-04 3,698E-06 3,089E-03 4203,951 2095,483 0,675 0,113
967,304 832,597 3,238E-04 3,391E-06 2,823E-03 4203,270 2115,455 0,674 0,113
967,328 834,210 3,239E-04 3,568E-06 2,976E-03 4203,237 2103,607 0,674 0,113
966,983 834,742 3,222E-04 3,630E-06 3,030E-03 4203,712 2099,704 0,675 0,113
967,187 835,243 3,232E-04 3,689E-06 3,081E-03 4203,431 2096,020 0,674 0,113
967,597 835,514 3,253E-04 3,722E-06 3,110E-03 4202,870 2094,029 0,674 0,113
967,010 835,349 3,224E-04 3,702E-06 3,093E-03 4203,674 2095,244 0,675 0,113
966,889 835,414 3,218E-04 3,710E-06 3,099E-03 4203,841 2094,766 0,675 0,113
966,834 835,452 3,215E-04 3,715E-06 3,103E-03 4203,917 2094,487 0,675 0,113
123
P.F.C.
124
P.F.C.
125
P.F.C.
126
P.F.C.
127
P.F.C.
128
P.F.C.
129
P.F.C.
130
P.F.C.
•
F [Hz] T [ºC] ∆Pmed [bar] ρ [kg/m3]
P P µ [m2/s]
P P µ [Pa·s] m [kg/min]
22,50 18,046 0,691 876,556 5,829E-05 0,051 51,009
25,00 18,596 0,792 876,257 5,543E-05 0,049 56,475
27,50 19,291 0,892 875,880 5,214E-05 0,046 61,891
27,50 20,114 0,874 875,434 4,862E-05 0,043 61,856
25,00 20,668 0,749 875,134 4,647E-05 0,041 56,403
22,50 20,999 0,637 874,954 4,525E-05 0,040 50,931
20,00 21,211 0,530 874,839 4,450E-05 0,039 45,432
17,50 21,340 0,433 874,769 4,405E-05 0,039 39,904
15,00 21,389 0,349 874,743 4,388E-05 0,038 34,343
12,50 21,383 0,270 874,746 4,390E-05 0,038 28,748
10,00 21,349 0,200 874,765 4,402E-05 0,039 23,118
7,50 21,286 0,140 874,799 4,424E-05 0,039 17,453
5,00 21,222 0,079 874,833 4,446E-05 0,039 11,752
16,00 18,582 0,487 876,265 5,550E-05 0,049 36,586
20,00 19,247 0,636 875,904 5,234E-05 0,046 45,469
24,00 20,215 0,796 875,380 4,822E-05 0,042 54,237
28,00 22,359 0,865 874,217 4,075E-05 0,036 62,837
12,00 93,520 0,097 835,634 3,737E-06 0,003 26,725
16,00 92,224 0,157 836,336 3,825E-06 0,003 34,785
20,00 90,369 0,229 837,342 3,957E-06 0,003 42,822
24,00 90,226 0,318 837,419 3,967E-06 0,003 50,700
28,00 89,831 0,422 837,634 3,997E-06 0,003 58,514
32,00 90,207 0,530 837,430 3,969E-06 0,003 66,177
36,00 90,832 0,649 837,091 3,923E-06 0,003 73,715
40,00 91,858 0,780 836,534 3,850E-06 0,003 81,098
44,00 92,075 0,923 836,417 3,835E-06 0,003 88,477
45,50 91,729 0,985 836,604 3,859E-06 0,003 91,280
20,00 94,926 0,228 834,871 3,645E-06 0,003 42,560
25,00 97,443 0,335 833,507 3,489E-06 0,003 52,116
40,00 98,787 0,751 832,778 3,410E-06 0,003 80,237
30,00 101,627 0,457 831,238 3,253E-06 0,003 61,290
16,00 105,683 0,157 829,039 3,047E-06 0,003 34,141
12,00 107,147 0,099 828,245 2,978E-06 0,002 26,241
20,00 111,380 0,216 825,950 2,791E-06 0,002 41,520
24,00 110,340 0,297 826,514 2,835E-06 0,002 49,197
28,00 110,868 0,393 826,227 2,813E-06 0,002 56,666
32,00 111,608 0,492 825,826 2,782E-06 0,002 64,010
36,00 112,420 0,601 825,386 2,748E-06 0,002 71,238
40,00 113,098 0,721 825,018 2,721E-06 0,002 78,375
44,00 113,773 0,852 824,652 2,694E-06 0,002 85,405
47,00 114,522 0,957 824,246 2,664E-06 0,002 90,570
131
P.F.C.
•
F [Hz] T [ºC] ∆Pmed [bar] ρ [kg/m3]
P P µ [m2/s]
P P µ [Pa·s] m [kg/min]
47,50 115,096 0,971 823,935 2,642E-06 0,002 91,353
16,00 117,920 0,148 822,404 2,537E-06 0,002 33,471
20,00 118,494 0,214 822,093 2,517E-06 0,002 41,026
24,00 119,132 0,292 821,746 2,494E-06 0,002 48,471
28,00 119,906 0,383 821,327 2,468E-06 0,002 55,799
32,00 120,792 0,476 820,847 2,437E-06 0,002 63,005
36,00 121,614 0,581 820,401 2,410E-06 0,002 70,108
40,00 122,526 0,697 819,906 2,380E-06 0,002 77,088
44,00 123,522 0,822 819,366 2,348E-06 0,002 83,943
48,00 124,584 0,959 818,790 2,315E-06 0,002 90,672
49,00 125,374 0,988 818,362 2,291E-06 0,002 92,247
16,00 129,730 0,142 816,000 2,164E-06 0,002 32,748
20,00 130,601 0,206 815,528 2,140E-06 0,002 40,121
24,00 128,535 0,285 816,648 2,197E-06 0,002 47,650
28,00 129,022 0,371 816,384 2,183E-06 0,002 54,880
32,00 130,023 0,463 815,842 2,155E-06 0,002 61,949
36,00 130,881 0,565 815,376 2,132E-06 0,002 68,922
40,00 131,865 0,677 814,843 2,105E-06 0,002 75,766
44,00 133,354 0,799 814,035 2,066E-06 0,002 82,416
48,00 135,375 0,925 812,940 2,015E-06 0,002 88,846
50,00 138,649 0,984 811,165 1,936E-06 0,002 91,601
48,00 141,410 0,907 809,667 1,874E-06 0,002 87,798
44,00 140,836 0,779 809,978 1,886E-06 0,002 81,219
40,00 138,277 0,662 811,366 1,945E-06 0,002 74,832
36,00 136,374 0,556 812,398 1,990E-06 0,002 68,197
32,00 135,465 0,456 812,891 2,013E-06 0,002 61,305
28,00 134,993 0,364 813,147 2,025E-06 0,002 54,254
24,00 134,625 0,279 813,347 2,034E-06 0,002 47,092
20,00 134,788 0,207 813,258 2,030E-06 0,002 39,790
16,00 53,896 0,220 857,118 9,379E-06 0,008 36,086
20,00 54,102 0,300 857,006 9,319E-06 0,008 44,512
24,00 54,345 0,394 856,874 9,250E-06 0,008 52,842
28,00 54,701 0,493 856,681 9,149E-06 0,008 61,070
32,00 55,267 0,618 856,374 8,994E-06 0,008 69,185
36,00 56,051 0,758 855,949 8,784E-06 0,008 77,177
40,00 57,055 0,902 855,405 8,528E-06 0,007 85,037
42,00 58,607 0,970 854,563 8,154E-06 0,007 88,808
40,00 59,668 0,894 853,988 7,913E-06 0,007 84,764
36,00 56,266 0,760 855,833 8,729E-06 0,007 77,158
40,00 50,906 0,932 858,739 1,032E-05 0,009 85,666
36,00 52,496 0,773 857,877 9,800E-06 0,008 77,494
132
P.F.C.
•
F [Hz] T [ºC] ∆Pmed [bar] ρ [kg/m3]
P P µ [m2/s]
P P µ [Pa·s] m [kg/min]
32,00 53,505 0,629 857,330 9,494E-06 0,008 69,319
28,00 53,996 0,500 857,064 9,350E-06 0,008 61,115
24,00 54,176 0,394 856,966 9,298E-06 0,008 52,850
20,00 54,131 0,303 856,990 9,311E-06 0,008 44,511
16,00 53,891 0,222 857,120 9,380E-06 0,008 36,087
12,00 84,341 0,108 840,610 4,440E-06 0,004 26,995
16,00 82,307 0,168 841,713 4,625E-06 0,004 35,198
20,00 77,778 0,246 844,169 5,083E-06 0,004 43,489
24,00 79,459 0,335 843,257 4,905E-06 0,004 51,415
28,00 80,791 0,441 842,535 4,771E-06 0,004 59,235
32,00 81,773 0,555 842,003 4,676E-06 0,004 66,963
36,00 82,149 0,679 841,799 4,640E-06 0,004 74,640
40,00 82,120 0,810 841,815 4,643E-06 0,004 82,266
44,00 81,819 0,944 841,978 4,671E-06 0,004 89,841
40,00 81,042 0,812 842,399 4,746E-06 0,004 82,392
36,00 80,512 0,684 842,686 4,798E-06 0,004 74,810
32,00 80,136 0,554 842,890 4,836E-06 0,004 67,111
28,00 80,135 0,445 842,891 4,836E-06 0,004 59,286
24,00 80,053 0,335 842,935 4,845E-06 0,004 51,378
20,00 80,193 0,245 842,860 4,830E-06 0,004 43,368
16,00 80,321 0,169 842,790 4,818E-06 0,004 35,274
12,00 80,265 0,111 842,820 4,823E-06 0,004 27,100
133
P.F.C.
134
P.F.C.
135
P.F.C.
136
P.F.C.
137
P.F.C.
138