Informe Práctica 1
Flujo interno
Autores:
Mikel Vicente Casanova
Jorge Moreno Santiago
David Lampurlanés Nadal
Oliver Cano Hopkins
Octubre 2024
ÍNDICE ÍNDICE
Índice
1. Introducción 2
2. Cálculos previos 2
3. Resultados 3
3.1. Perfiles de temperatura y velocidad axiales en las distintas ISO su-
perficies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
3.1.1. Velocidad axial . . . . . . . . . . . . . . . . . . . . . . . . . . 3
3.1.2. Temperatura axial . . . . . . . . . . . . . . . . . . . . . . . . 3
3.2. Perfiles de temperatura y de flujo de calor a través de la superficie a
lo largo del conducto . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
3.2.1. Temperatura . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
3.2.2. Flujo de calor . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
3.3. Evolución de la temperatura a lo largo de la coordenada axial . . . . 5
3.4. Coeficiente de convección a lo largo de la pared del conducto . . . . . 5
3.5. Otros cálculos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.5.1. Distancia donde el flujo se encuentra fluidodinámica desarrollado 7
3.5.2. Comparación perfil teórico con simulado . . . . . . . . . . . . 8
3.5.3. Valor tórico de la velocidad en el eje . . . . . . . . . . . . . . 9
1 UC3M
Departamento de Ingenierı́a Térmica y de Fluidos
2 CÁLCULOS PREVIOS
1. Introducción
El propósito de la práctica es analizar la transferencia de calor por flujo interno a
lo largo de un conducto. Para lo cual trabajaremos con el software Ansys fluent,
además el caso que estaremos trabajando es el C, tal que:
Fluido: Aceite
ρ = 757kg/m3
cp J/kgK = 1977J/kgK
K(W/mK) = 0, 086W/mK
µ = 0, 0001822kg/ms
D = 0, 05m
ṁ= 0, 015kg/s
Twall constant= 400o C
Tinlet = 280o C
2. Cálculos previos
Difusividad térmica: α = k
ρcp
= 5, 75 · 10−8 m2 /s
µ
Viscosidad cinemática : v = ρ
= 2, 407 · 10−7 m2 /s
v
Número de Prandtl: Pr = α
= 4, 188
ṁ ṁ
Velocidad media: um = ρAC
= 2 = 0, 01001m/s
ρ πD
4
ρum D
Número de Reynolds: Re = µ
= 2096, 44
Longitud a la que el flujo se desarrolla fluidodinámica:
z
f d,h
≈ 0′ 05 ReD → zf d,h ≈ 0′ 05 ReD D
D lam
zf d,h = 5, 24m
Longitud a la que el flujo se desarrolla térmicamente:
z
f d,t
≈ 0′ 05 ReD Pr → zf d,t ≈ 0′ 05 ReD Pr D
D lam
zf d,t = 21, 95m
2 UC3M
Departamento de Ingenierı́a Térmica y de Fluidos
3 RESULTADOS
3. Resultados
3.1. Perfiles de temperatura y velocidad axiales en las dis-
tintas ISO superficies
Las medidas de la temperatura y la velocidad se han realizado a lo largo de la
coordenada axial.
3.1.1. Velocidad axial
Figura 1: Velocidades axiales
Se observa que la mayor velocidad en todas las superficies se encuentra en el centro
del conducto, donde r = 0m. Ası́ como la menor velocidad mı́nima es, por continui-
dad, nula en el punto de contacto con la pared del tubo.
3.1.2. Temperatura axial
Figura 2: Temperaturas axiales
De igual forma, la temperatura en el punto de contacto con la pared de la tuberı́a
sera mayor. También se puede ver que la temperatura es mayor en todos los puntos
axiales conforme se avanza en las isosuperficies, además de que la temperatura en el
centro del del conducto se mantiene constante a lo largo de toda la tuberı́a, siendo
ası́ la temperatura de entrada y salida aproximadamente iguales para r < 0,02m.
3 UC3M
Departamento de Ingenierı́a Térmica y de Fluidos
3.2 Perfiles de temperatura y de flujo de calor a través de la superficie a lo largo del
conducto 3 RESULTADOS
3.2. Perfiles de temperatura y de flujo de calor a través de
la superficie a lo largo del conducto
3.2.1. Temperatura
Figura 3: Temperatura del fluido en contacto con la pared
Se observa como la temperatura aumenta rápidamente al inicio del conducto para
luego estabilizarse a la temperatura Twall .
3.2.2. Flujo de calor
Figura 4: Flujo de calor a través de la pared
En esta figura se ve como el flujo de calor a través de la pared disminuye a medida que
el fluido avanza, esto debido a la disminución del coeficiente convectivo producido
por la disminución del incremento de temperatura entre la pared y el fluido.
4 UC3M
Departamento de Ingenierı́a Térmica y de Fluidos
3.3 Evolución de la temperatura a lo largo de la coordenada axial 3 RESULTADOS
3.3. Evolución de la temperatura a lo largo de la coordenada
axial
Iso-superficie x (m) Temperatura media (K)
In 0 553
hd 5,24 570,572
S1 10 575,833
S2 20 582,865
td 21,95 583,887
S3 30 587,662
S4 40 591,346
S5 50 594,435
S6 60 597,023
Out 70 599,375
Cuadro 1: Temperaturas a lo largo del eje axial
Figura 5: Temperatura media a lo largo del eje axial
Como hemos visto en la figura 2 la temperatura al inicio será Tinlet y, progresiva-
mente, va aumentando debido al flujo de calor.
3.4. Coeficiente de convección a lo largo de la pared del
conducto
Para calcular el coeficiente de convección a lo largo de la pared usamos la siguiente
ecuación:
q ′′ (z)
h(z) = (1)
Tp (z) − Tm (z)
5 UC3M
Departamento de Ingenierı́a Térmica y de Fluidos
3.4 Coeficiente de convección a lo largo de la pared del conducto 3 RESULTADOS
Los valores de las tres incógnitas los podemos obtener de la simulación, siendo q ′′ (z)
el flujo de calor a través de la pared de la tuberı́a, Tp (z) la temperatura de la
pared (constante en nuestro caso), y Tm (z) la temperatura media en cada sección
(isosuperficie).
Igual que en el apartado anterior, calcularemos el coeficiente de convección en unas
secciones determinadas, con el fin de estudiar el comportamiento del fluido a lo
largo de toda la tuberı́a sin tener que hacer infinitos cortes. Como la tuberı́a mide 70
metros, hemos tomado la sección de entrada, la de salida y luego una cada 10 metros.
Si a estas les sumamos la xCr,t y la xCr,h nos quedan un total de 10 isosuperficies.
ISOSUPERFICIES POSICIÓN (m) FLUJO CALOR (W/m2 ) TM EDIA (K)
inlet 0 10872,5 553
hd 5,23894 583,51 559,13026
x10 10,1008 433,023 561,62638
x20 20,214 316,922 565,52441
td 22,1199 304,332 566,14867
x30 29,8554 266,03 568,59052
x40 40,274 232,586 571,19224
x50 49,6499 211,705 573,55124
x60 60,2915 193,954 575,6573
out 70 181,195 577,67522
Cuadro 2: Cálculo de temperatura media y flujo de calor
ISOSUPERFICIES POSICIÓN (m) hl oc(W/Km2 ) N uloc
inlet 0 90,60417 52,67684
hd 5,23894 5,124364 2,979281
x10 10,1008 3,888021 2,260477
x20 20,214 2,948781 1,714408
td 22,1199 2,848182 1,65592
x30 29,8554 2,547949 1,481366
x40 40,274 2,284561 1,328233
x50 49,6499 2,128785 1,237666
x60 60,2915 1,992486 1,158422
out 70 1,900817 1,105126
Cuadro 3: Calculo de coeficiente convectivo y Nusselt
Con los coeficientes de convección, usamos los datos de la segunda tabla (datos
iniciales del problema) para obtener el número de Nusselt, que se calcula:
hD
N uD = (2)
k
Comparando la última columna de la tabla 3 con el Nusselt teórico, que para una
Tpared =cte es N uD =3.66, vemos que, a partir del xCr,t han quedado demasiado pe-
6 UC3M
Departamento de Ingenierı́a Térmica y de Fluidos
3.5 Otros cálculos 3 RESULTADOS
queños, del orden de la mitad del valor teórico. Esto seguramente esté causado por
algún error al introducir las condiciones de contorno, que no hemos sido capaces de
detectar.
Figura 6: Variación del coeficiente convectivo con la posición
Podemos ver también, en el siguiente gráfico, que los valores del coeficiente convec-
tivo nos quedan demasiado pequeños y, consecuentemente, esto influye en el error
del número de Nusselt local.
3.5. Otros cálculos
3.5.1. Distancia donde el flujo se encuentra fluidodinámica desarrollado
Figura 7: Velocidad en el eje axial frente a la distancia
Como se puede comprobar, debido a algún error en la simulación el flujo no llega
a desarrollarse completamente dentro del tubo puesto que no llega atener un com-
portamiento asintótico cuando x=L. El resultado teórico de la distancia a la que el
flujo se desarrolla fluidodinamicamente es x=5,24m obtenido mediante la siguiente
formula
7 UC3M
Departamento de Ingenierı́a Térmica y de Fluidos
3.5 Otros cálculos 3 RESULTADOS
Figura 8: Perfil simulado
Figura 9: Perfil teórico
z
f d,h
≈ 0′ 05 ReD → zf d,h ≈ 0′ 05 ReD D (3)
D lam
3.5.2. Comparación perfil teórico con simulado
" 2 #
u(r) r
=2 1− (4)
um r0
El perfil de velocidades obtenido en la simulación para un punto cualquiera (en este
caso x=35m) es el siguiente:
Siendo el eje y la posición radial y el eje x la velocidad en cada una de sus coorde-
nadas.
El perfil teórico obtenido mediante cálculos teóricos es el siguiente:
Como podemos observar, los valores en los cortes con los ejes son bastante similares
a los obtenidos en la simulación. En la simulación el flujo tiene velocidades muy
8 UC3M
Departamento de Ingenierı́a Térmica y de Fluidos
3.5 Otros cálculos 3 RESULTADOS
parecidas entre si una vez se encuentra más cerca del eje, sin embargo, en los valores
teóricos va acercándose poco a poco al valor máximo de la velocidad.
3.5.3. Valor tórico de la velocidad en el eje
El valor teórico es u(r=0)= 0,020183454 m/s, mientras que en la simulación se puede
observar que el valor de la velocidad en esa coordenada es aproximadamente u(r =
0)=0,017 m/s por lo que podemos considerarlos lo suficientemente parecidos entre
sı́ como para darlos por correctos
9 UC3M
Departamento de Ingenierı́a Térmica y de Fluidos
