SISTEMAS TERMICOS I

SOFTWARE PARA EL CÁLCULO DE POTENCIA DE CALDERAS PIROTUBULARES DE 2 PASOS

JOSE DANIEL GARCÍA LOPEZ COD: 2070286

NELSON EDUARDO JAIMES ARDILA COD: 2070279

MANUEL FERNANDO JEREZ CARRIZALES COD. 2071379

CRISTIAN FELIPE PÉREZ BROKATE COD. 2070290

CAMILO ANTONIO SANCHEZ VANEGAS COD. 2070301

SERGIO ARMANDO RIVERO GEREDA COD. 2070299

UNIVERSIDAD INDUSTRIAL DE SANTANDER

FACULTAD DE INGENIERIAS FISICOMECÁNICAS

ESCUELA DE INGENIERÍA MECÁNICA

SISTEMAS TERMICOS I
 1. INTRODUCCION:

Se plantea un análisis sobre los catálogos que se presentan para la selección de las calderas,
principalmente de las calderas pirotubulares que son las más usadas para servicios, en hospitales,
hoteles, etc.
Mediante este documento y software se busca una comparación teórica de los datos
proporcionados por fabricantes y un estudio de la transferencia de calor que ocurre dentro de una
caldera, teniendo en cuenta su geometría, combustible, presión de trabajo, etc.

2. OBJETIVO

2.1. OBJETIVO GENERAL

• Determinar teóricamente el calor generado en una caldera pirotubular y compararla con

el valor ofrecido por los fabricantes.

2.2. OBJETIVOS ESPECIFICOS

• Determinar mediante los datos geométricos ofrecidos por los fabricantes la geometría de
tubos, longitud efectiva y demás datos necesarios para el análisis teórico de la
combustión en el interior de la caldera.
• Obtener la grafica de temperatura de gases vs longitud efectiva de la caldera.
• Determinar el comportamiento de la absortividad vs la longitud efectiva de la caldera.
• Obtener la grafica del calor ofrecido vs longitud efectiva de la caldera.
• Determinar la grafica del comportamiento de la emisividad a lo largo de la caldera.
• Observar el comportamiento de la temperatura de la superficie de los tubos vs la longitud
de la caldera.

3. MARCO TEORICO

3.1. CALDERA PIROTUBULARES

Se denominan pirotubulares por ser los gases calientes procedentes de la combustión, los que
circulan por el interior de tubos cuyo exterior esta bañado por el agua de la caldera.

El combustible se quema en un hogar, en donde tiene lugar la

transmisión de calor, de esta forma los gases resultantes, se les hace
circular a través de los tubos que constituyen el haz tubular de la
caldera, y donde tiene lugar el intercambio de calor por conducción y
convección. Según sea una o varias las veces que los gases pasan a
través del haz tubular, se tienen las calderas de uno o de varios pasos.
En el caso de calderas de varios pasos, en cada uno de ellos, los
humos solo atraviesan un determinado número de tubos, cosa que se
logra mediante las denominadas cámaras de humos. Una vez
realizado el intercambio térmico, los humos son expulsados al
exterior a través de la chimenea.
 4. SOFTWARE PARA EL CÁLCULO DE POTENCIA DE CALDERAS PIROTUBULARES DE DOS
PASOS
Para poder desarrollar el análisis teórico de la caldera se recurrió al software Matlab el cual ofrece
un ambiente de programación adecuado, este permite realizar cambios a los datos de análisis y de
esta forma establecer procesos iterativos.

4.1. DATOS REQUERIDOS

El software realizado permite el cálculo de la potencia de una caldera, utilizando una serie de
datos de entrada que se deben proporcionar. Estos son:

- Diámetro tubo principal[in] - Emisividad superficie de los tubos

- Diámetros de los tubos 𝑓𝑓𝑡𝑡 3
- Flujo de combustible[ ]
secundarios[in] ℎ𝑟𝑟
- Presión de funcionamiento[PSI]
- Numero de tubos secundarios
- Exceso de aire[%]
- Longitud Efectiva[in]
- Tipo de combustible[CH4 o C2H6]
- Diámetro exterior de la caldera[in]
- Longitud total de la caldera[in]

En este caso se dejo una variable que permite establecer el número de divisiones que se desean
hacer para desarrollar el cálculo del calor, esto se debe a que es necesario desarrollar un análisis
por sección de tubo (integral discreta). Esta variable se establece como “Resolución de cálculo”.

4.2. PLATAFORMA GRAFICA

Matlab dentro de su programación permite desarrollar un entorno amigable en el cual se puede

interactuar con las variables que hacen parte del análisis. Además, permite presentar los
resultados tanto gráficos como concretos. De esta forma la plataforma grafica utilizada en el
software es la siguiente:
Los resultados que entrega el software son:
• La Potencia de la caldera [BHP]
• El aporte de Potencia de cada paso
• Distintas graficas que permite un mejor análisis de la caldera:
 Grafica de Emisividad Vs longitud efectiva
 Grafica de Absortividad Vs longitud efectiva
 Grafica de Calor generado Vs longitud efectiva
 Grafica de Temperatura de gases Vs longitud efectiva
 Grafica de Temperatura de superficie Vs longitud efectiva

4.3. CODIGO USADO:

El código aplicado será explicado con detenimiento:

4.3.1. ESQUEMA:
 4.3.2. DEFINICION DE VARIABLES

En estos dos primeros casos se determinan todas las variables que son necesarias para llevar a
cabo el cálculo de la potencia teórica de la caldera. Estos datos deben ser suministrados por el
fabricante, el cual especifica en sus catálogos no solo la geometría sino también el combustible
usado y el consumo de este, en el caso en que no se cuenten directamente con todos estos datos,
se debe hacen un cálculo aproximado.

A continuación se presenta el código necesario para adquirir los datos de entrada, con los cuales
se realizara el análisis.
Dp=str2double(get(handles.edit1,'string')); % Diametro tubo principal
Ds=str2double(get(handles.edit3,'string')); %Diámetro tubos secundarios
Nts=str2double(get(handles.edit2,'string')); %Numero tubos secundarios
Le=str2double(get(handles.edit4,'string')); %Longitud Efectiva
Dec=str2double(get(handles.edit7,'string')); %Diametro exterior caldera
L_total=str2double(get(handles.edit9,'string')); %longitud total de la caldera
comb=get(handles.popupmenu1,'value'); %Tipo de combustible
N=str2double(get(handles.edit10,'string')); %Resolucion Calculo
ex=str2double(get(handles.edit12,'string')); %Exceso de aire
Es=str2double(get(handles.edit13,'string')); %Emisividad superficie
fcombft=str2double(get(handles.edit21,'string')); %flujocombustible
PSI=str2double(get(handles.edit22,'string')); %Presion funcionamiento
exc=ex/100; %Porcentaje de exceso

4.3.3. SELECCIÓN COMBUSTIBLE

Según el combustible seleccionado, se debe realizar el cálculo de la temperatura de llama inicial,
de las fracciones molares de los reactivos al igual que de los productos en función del exceso y de
la densidad del combustible de la siguiente forma:
if comb==1 %selección de el tipo de combustible, para este caso se
relaciona 1 con el CH4

TG=2329-(1485*exc)+(1124*exc^2)-(861.8*exc^3)+(638.6*exc^4)-(391.3*exc^5)+
(134.4*exc^6); Temperatura de llama del CH4 en función del exceso de
combustible

%determinacion de las fracciones molares de los productos en la combustión de CH4

en función del exceso

X_co2=1/(1+2+((1+exc)*7.52)+(2*exc)); %fraccion molar del bioxido de carbono

X_h2o=2/(1+2+((1+exc)*7.52)+(2*exc)); %fraccion molar del vapor de agua
X_n2=((1+exc)*7.52)/(1+2+((1+exc)*7.52)+(2*exc)); %fraccion molar del nitrogeno
X_o2=(2*exc)/(1+2+((1+exc)*7.52)+(2*exc)); %fraccion molar del oxigeno
denscomb=0.6472; %determinación de la densidad del CH4

else selección de el tipo de combustible, para este caso se relaciona con el

C2H6

TG=2382-1546.04*exc+1191*exc^2-929.5*exc^3+699.7*exc^4-433.3*exc^5+
149.5*exc^6; Temperatura de llama del C2H6 en función del exceso de
combustible

%determinación de las fracciones molares de los productos en la combustión de

C2H6

X_co2=2/(2+3+((1+exc)*13.16)+(3.5*exc)); %fracción molar del bióxido de

carbono
X_h2o=3/(2+3+((1+exc)*13.16)+(3.5*exc)); %fracción molar del vapor de agua
X_n2=((1+exc)*13.16)/(2+3+((1+exc)*13.16)+(3.5*exc)); %fracción molar del
nitrógeno
X_o2=(3.5*exc)/(2+3+((1+exc)*13.16)+(3.5*exc)); %fracción molar del oxigeno
denscomb=1.213; %determinación de la densidad del C2H6

end

Las formulas usadas en el código anterior se justifican mediante la siguiente metodología de

cálculo.

4.3.3.1. Temperatura de llama:

Determinación de las temperaturas de llama de los dos tipos de combustible implementados. Se

utilizo el programa EES y se muestra el cálculo tipo para el CH4. El archivo se encuentra adjunto y
se denomina ‘Temperatura de llama CH4.EES’

De esta forma se realiza el balance de combustión respectivo y se determina la Temperatura de

llama para un exceso dado, y se varía el exceso para obtener los demás valores de la temperatura.
“𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙 𝑑𝑑𝑑𝑑𝑑𝑑 𝐶𝐶𝐶𝐶4”

𝑄𝑄1 = 𝑁𝑁𝐶𝐶𝐶𝐶2 ∗ �𝐻𝐻𝑐𝑐𝑐𝑐 22 − 𝐻𝐻𝑐𝑐𝑐𝑐2 + 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑎𝑎𝑐𝑐𝑐𝑐2 � + 𝑁𝑁𝑂𝑂2𝑠𝑠[𝑗𝑗 ] ∗ �𝐻𝐻𝑜𝑜22 − 𝐻𝐻𝑜𝑜2 + 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑎𝑎𝑜𝑜2 � + 𝑁𝑁𝑁𝑁2[𝑗𝑗 ] ∗ �𝐻𝐻𝑛𝑛22 − 𝐻𝐻𝑛𝑛2 + 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑎𝑎𝑛𝑛2 � + 𝑁𝑁𝐻𝐻2𝑂𝑂
∗ �𝐻𝐻ℎ2𝑜𝑜2 − 𝐻𝐻ℎ2𝑜𝑜 + 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑎𝑎ℎ2𝑜𝑜 �
“𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸í𝑎𝑎 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝, 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 𝑎𝑎 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙(𝑖𝑖𝑖𝑖𝑖𝑖ó𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔)”

𝑄𝑄2 = 𝑁𝑁_𝐶𝐶𝐶𝐶4 ∗ (𝐻𝐻_𝑐𝑐ℎ4_1 − 𝐻𝐻_𝑐𝑐ℎ4 + 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑_𝑐𝑐ℎ4) + 𝑁𝑁_𝑂𝑂2𝑒𝑒[𝑗𝑗] ∗ (𝐻𝐻_𝑜𝑜2_1 − 𝐻𝐻_𝑜𝑜2 + 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑_𝑜𝑜2) + 𝑁𝑁_𝑁𝑁2[𝑗𝑗] ∗ (𝐻𝐻_𝑛𝑛2_𝑎𝑎𝑎𝑎𝑎𝑎 − 𝐻𝐻_𝑛𝑛2
+ 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑_𝑛𝑛2)
“𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸í𝑎𝑎 𝑑𝑑𝑑𝑑𝑑𝑑 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐, 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 𝑎𝑎 𝑙𝑙𝑙𝑙 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑑𝑑𝑑𝑑 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒”

𝑄𝑄1 = 𝑄𝑄2 “𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵 𝑑𝑑𝑑𝑑 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒í𝑎𝑎 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐ó𝑛𝑛”

“𝐷𝐷𝐷𝐷 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓 𝑒𝑒𝑒𝑒 𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 𝑑𝑑𝑑𝑑 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑢𝑢𝑢𝑢𝑢𝑢 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐ó𝑛𝑛, 𝑡𝑡𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 𝑑𝑑𝑑𝑑𝑑𝑑 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙
𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 ”

"𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸 𝑎𝑎 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠"

𝑇𝑇𝑒𝑒𝑒𝑒𝑒𝑒 = 25 + 273 “𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇 𝑑𝑑𝑑𝑑 𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝”

𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝐶𝐶𝐶𝐶2′; 𝑇𝑇_𝑒𝑒𝑒𝑒𝑒𝑒: 𝐶𝐶𝐶𝐶_𝑐𝑐𝑐𝑐2; 𝐻𝐻_𝑐𝑐𝑐𝑐2; 𝑆𝑆_𝑐𝑐𝑐𝑐2)
𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝑂𝑂2′; 𝑇𝑇_𝑒𝑒𝑒𝑒𝑒𝑒: 𝐶𝐶𝐶𝐶_𝑜𝑜2; 𝐻𝐻_𝑜𝑜2; 𝑆𝑆_𝑜𝑜2)
𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝑁𝑁2′; 𝑇𝑇_𝑒𝑒𝑒𝑒𝑒𝑒: 𝐶𝐶𝐶𝐶_𝑛𝑛2; 𝐻𝐻_𝑛𝑛2; 𝑆𝑆_𝑛𝑛2)
𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝐻𝐻2𝑂𝑂′; 𝑇𝑇_𝑒𝑒𝑒𝑒𝑒𝑒: 𝐶𝐶𝐶𝐶_ℎ2𝑜𝑜; 𝐻𝐻_ℎ2𝑜𝑜; 𝑆𝑆_ℎ2𝑜𝑜)
𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝐶𝐶𝐶𝐶4′; 𝑇𝑇_𝑒𝑒𝑒𝑒𝑒𝑒: 𝐶𝐶𝐶𝐶_𝑐𝑐ℎ4; 𝐻𝐻_𝑐𝑐ℎ4; 𝑆𝑆_𝑐𝑐ℎ4)

"𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒"

𝑇𝑇𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖 = 25 + 273 “𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇 𝑑𝑑𝑑𝑑 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑑𝑑𝑑𝑑𝑑𝑑 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐”

𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝐶𝐶𝐶𝐶4′; 𝑇𝑇_𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖: 𝐶𝐶𝐶𝐶_𝑐𝑐ℎ4_1; 𝐻𝐻_𝑐𝑐ℎ4_1; 𝑆𝑆_𝑐𝑐ℎ4_1)

𝑇𝑇𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖 = 25 + 273 “𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇 𝑑𝑑𝑑𝑑 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑑𝑑𝑑𝑑𝑑𝑑 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎”

𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝑂𝑂2′; 𝑇𝑇_𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖: 𝐶𝐶𝐶𝐶_𝑜𝑜2_1; 𝐻𝐻_𝑜𝑜2_1; 𝑆𝑆_𝑜𝑜2_1)
𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝑁𝑁2′; 𝑇𝑇_𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖: 𝐶𝐶𝐶𝐶_𝑁𝑁2_𝑎𝑎𝑎𝑎𝑎𝑎; 𝐻𝐻_𝑛𝑛2_𝑎𝑎𝑎𝑎𝑎𝑎; 𝑆𝑆_𝑛𝑛2_𝑎𝑎𝑎𝑎𝑎𝑎)

"𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸 𝑑𝑑𝑑𝑑 𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓ó𝑛𝑛 𝑑𝑑𝑑𝑑 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑢𝑢𝑢𝑢𝑢𝑢 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐"

𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑_ℎ2 = 0
𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑_𝑐𝑐ℎ4 = −74900
𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑_𝑐𝑐2ℎ6 = −84700
𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑_𝑛𝑛2 = 0
𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑_𝑜𝑜2 = 0
𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑_𝑐𝑐𝑐𝑐 = −110500
𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑_𝑐𝑐𝑐𝑐2 = −393500
𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑_ℎ2𝑜𝑜 = −241800

"𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁 𝑑𝑑𝑑𝑑 𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒"

𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑖𝑖 = 1; 40

𝑁𝑁𝑂𝑂2𝑒𝑒[𝑖𝑖] = 2 ∗ �1 + (𝑖𝑖 ∗ 0,02)� “𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁 𝑑𝑑𝑑𝑑 𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚 𝑑𝑑𝑑𝑑 𝑂𝑂2”

𝑁𝑁𝑁𝑁2[𝑖𝑖] = 2 ∗ �1 + (𝑖𝑖 ∗ 0,02)� ∗ 3,76 “𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁 𝑑𝑑𝑑𝑑 𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚 𝑑𝑑𝑑𝑑 𝑁𝑁2”
𝑁𝑁𝑂𝑂2𝑠𝑠[𝑖𝑖] = 2 ∗ (𝑖𝑖 ∗ 0,02) “𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁 𝑑𝑑𝑑𝑑 𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚 𝑑𝑑𝑑𝑑 𝑂𝑂2 𝑎𝑎 𝑙𝑙𝑙𝑙 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠”

𝐸𝐸𝐸𝐸𝐸𝐸

𝑁𝑁𝐶𝐶𝐶𝐶2 = 1 “𝑁𝑁𝑢𝑢𝑢𝑢𝑢𝑢𝑢𝑢𝑢𝑢 𝑑𝑑𝑑𝑑 𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚 𝑑𝑑𝑑𝑑 𝐶𝐶𝐶𝐶2”

𝑁𝑁𝐻𝐻2𝑂𝑂 = 2 “𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁 𝑑𝑑𝑑𝑑 𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚 𝑑𝑑𝑑𝑑 𝐻𝐻2𝑂𝑂”
𝑁𝑁𝐶𝐶𝐶𝐶4 = 1 “𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁 𝑑𝑑𝑑𝑑 𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑠𝑠 𝑑𝑑𝑑𝑑 𝐶𝐶𝐶𝐶4”

"𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 𝑎𝑎 𝑙𝑙𝑙𝑙 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙 (𝑇𝑇2), 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑠𝑠𝑠𝑠 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖 𝑦𝑦𝑦𝑦 𝑞𝑞𝑞𝑞𝑞𝑞 𝑒𝑒𝑒𝑒 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 𝐸𝐸𝐸𝐸𝐸𝐸 𝑒𝑒𝑒𝑒 𝑒𝑒𝑒𝑒 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒
𝑑𝑑𝑑𝑑 ℎ𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 𝑒𝑒𝑒𝑒 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐"

𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝐶𝐶𝐶𝐶2′; 𝑇𝑇2: 𝐶𝐶𝐶𝐶_𝑐𝑐𝑐𝑐2_2; 𝐻𝐻_𝑐𝑐𝑐𝑐2_2; 𝑆𝑆_𝑐𝑐𝑐𝑐2_2)

𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝑂𝑂2′; 𝑇𝑇2: 𝐶𝐶𝑝𝑝_𝑜𝑜2_2; 𝐻𝐻_𝑜𝑜2_2; 𝑆𝑆_𝑜𝑜2_2)
𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝐶𝐶𝐶𝐶′; 𝑇𝑇2: 𝐶𝐶𝐶𝐶_𝑐𝑐𝑐𝑐_2; 𝐻𝐻_𝑐𝑐𝑐𝑐_2; 𝑆𝑆_𝑐𝑐𝑐𝑐_2)
𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝑁𝑁2′; 𝑇𝑇2: 𝐶𝐶𝐶𝐶_𝑛𝑛2_2; 𝐻𝐻_𝑛𝑛2_2; 𝑆𝑆_𝑛𝑛2_2)
𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽(′𝐻𝐻2𝑂𝑂′; 𝑇𝑇2: 𝐶𝐶𝐶𝐶_ℎ2𝑜𝑜_2; 𝐻𝐻_ℎ2𝑜𝑜_2; 𝑆𝑆_ℎ2𝑜𝑜_2)
“𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑠𝑠𝑠𝑠 𝑢𝑢𝑢𝑢𝑢𝑢 𝑙𝑙𝑙𝑙 𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓ó𝑛𝑛 𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 𝑒𝑒𝑒𝑒 𝑐𝑐á𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐ó𝑛𝑛 𝑦𝑦𝑦𝑦 𝑞𝑞𝑞𝑞𝑞𝑞 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒
𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓ó𝑛𝑛 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑛𝑛𝑛𝑛 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 𝑒𝑒𝑒𝑒 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎”

De esta forma se tabulan las temperaturas de llama para cada uno de los excesos y se realiza el
mismo cálculo para el combustible C2H6. El archivo correspondiente se denomina ‘Temperatura
de llama C2H6.EES’.

TEMPERATURA DE LLAMA
Exceso CH4 C2H6
0,02 2299 2351
0,04 2270 2321
0,06 2243 2293
0,08 2216 2265
0,1 2190 2238
0,12 2165 2212
0,14 2140 2186
0,16 2116 2161
0,18 2093 2137
0,2 2070 2114

Con los datos tabulados se realiza una regresión para dejar la Temperatura de llama en función del
exceso de aire.

Para CH4:
TG=2329-(1485*exc)+(1124*exc^2)-(861.8*exc^3)+(638.6*exc^4)-(391.3*exc^5)+
(134.4*exc^6)

Para C2H6:
TG=2382-1546.04*exc+1191*exc^2-929.5*exc^3+699.7*exc^4-433.3*exc^5+149.5*exc^6

4.3.3.2. Fracciones molares:

Por medio del balance molar de cada una de las ecuaciones químicas de combustión se pueden
determinar las fracciones molares de los gases de combustión en función del exceso de aire. El
análisis de la combustión se llevo a cabo de la siguiente forma:

 Análisis de combustión completa del CH4

1 ∗ 𝐶𝐶𝐶𝐶4 + 𝐴𝐴 ∗ (1 + 𝐸𝐸𝐸𝐸𝐸𝐸) ∗ (𝑂𝑂2 + 3,76 𝑁𝑁2 ) = 𝐵𝐵 ∗ 𝐶𝐶𝑂𝑂2 + 𝐶𝐶 ∗ 𝐻𝐻2 𝑂𝑂 + 𝐷𝐷 ∗ 𝑁𝑁2 + 𝐸𝐸 ∗ 𝑂𝑂2

De esta forma la combustión de un mol de CH4 arroja los siguientes resultados en función del
exceso:
𝐵𝐵 = 1
𝐶𝐶 = 2
𝐷𝐷 = 𝐴𝐴(1 + 𝐸𝐸𝐸𝐸𝐸𝐸) ∗ 3,76
𝐸𝐸 = 𝐴𝐴 ∗ 𝐸𝐸𝐸𝐸𝐸𝐸
 𝐴𝐴 = 2 𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣 𝑠𝑠𝑠𝑠 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 𝑎𝑎 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡é𝑠𝑠 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐ó𝑛𝑛 𝑠𝑠𝑠𝑠𝑠𝑠 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒

De esta forma se pueden determinar las fracciones molares de cada uno de los productos de la
combustión, para lo cual tenemos:
𝐵𝐵
𝑋𝑋𝐶𝐶𝑂𝑂2 =
𝐵𝐵 + 𝐶𝐶 + 𝐷𝐷 + 𝐸𝐸
𝐶𝐶
𝑋𝑋𝐻𝐻2𝑂𝑂 =
𝐵𝐵 + 𝐶𝐶 + 𝐷𝐷 + 𝐸𝐸
𝐷𝐷
𝑋𝑋𝑁𝑁2 =
𝐵𝐵 + 𝐶𝐶 + 𝐷𝐷 + 𝐸𝐸
𝐸𝐸
𝑋𝑋𝑂𝑂2 =
𝐵𝐵 + 𝐶𝐶 + 𝐷𝐷 + 𝐸𝐸

 Análisis de combustión completa del C2H6

Para este caso se realiza un análisis similar al realizado anteriormente:
1 ∗ 𝐶𝐶2𝐻𝐻6 + 𝐴𝐴 ∗ (1 + 𝐸𝐸𝐸𝐸𝐸𝐸) ∗ (𝑂𝑂2 + 3,76 𝑁𝑁2 ) = 𝐵𝐵 ∗ 𝐶𝐶𝑂𝑂2 + 𝐶𝐶 ∗ 𝐻𝐻2 𝑂𝑂 + 𝐷𝐷 ∗ 𝑁𝑁2 + 𝐸𝐸 ∗ 𝑂𝑂2
𝐵𝐵 = 2
𝐶𝐶 = 3
𝐷𝐷 = 𝐴𝐴(1 + 𝐸𝐸𝐸𝐸𝐸𝐸) ∗ 3,76
𝐸𝐸 = 𝐴𝐴 ∗ 𝐸𝐸𝐸𝐸𝐸𝐸
𝐴𝐴 = 3,5 𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣 𝑠𝑠𝑠𝑠 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 𝑎𝑎 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡é𝑠𝑠 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑠𝑠𝑠𝑠𝑠𝑠ó𝑛𝑛 𝑠𝑠𝑠𝑠𝑠𝑠 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒
Para este caso se determinan las fracciones molares de cada uno de la misma forma en
que se hizo para el CH4

FRACCIONES MOLARES PARA EL C2H6

Exceso N totales X_CO2 X_H2O X_N2
0,02 18,4932 0,108 0,162 0,726
0,04 18,8264 0,106 0,159 0,727
0,06 19,1596 0,104 0,157 0,728
0,08 19,4928 0,103 0,154 0,729
0,1 19,826 0,101 0,151 0,730
0,12 20,1592 0,099 0,149 0,731
0,14 20,4924 0,098 0,146 0,732
0,16 20,8256 0,096 0,144 0,733
0,18 21,1588 0,095 0,142 0,734
0,2 21,492 0,093 0,140 0,735
 FRACCIONES MOLARES PARA EL CH4
Exceso N totales X_CO2 X_H2O X_N2
0,02 10,5268 0,095 0,190 0,711
0,04 10,7136 0,093 0,187 0,713
0,06 10,9004 0,092 0,183 0,714
0,08 11,0872 0,090 0,180 0,715
0,1 11,274 0,089 0,177 0,716
0,12 11,4608 0,087 0,175 0,717
0,14 11,6476 0,086 0,172 0,718
0,16 11,8344 0,084 0,169 0,719
0,18 12,0212 0,083 0,166 0,720
0,2 12,208 0,082 0,164 0,721

4.3.4. CALOR ENTREGADO

Para poder determinar el calor entregado por la caldera se debe realizar el siguiente balance de
energía:
𝑄𝑄 = 𝑄𝑄𝑟𝑟𝑟𝑟𝑟𝑟 + 𝑄𝑄𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 = 𝑄𝑄𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒

4
𝜀𝜀𝜀𝜀 ∗ (𝜀𝜀𝜀𝜀 ∗ 𝜎𝜎 ∗ 𝑇𝑇𝐺𝐺𝑥𝑥 4 − 𝑎𝑎𝑔𝑔 ∗ 𝜎𝜎 ∗ �𝑇𝑇sup � )
𝑄𝑄 = 𝑑𝑑𝑑𝑑 ∗ + ℎ𝑖𝑖 ∗ 𝑑𝑑𝑑𝑑 ∗ (𝑇𝑇𝑇𝑇 − 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡)
𝜀𝜀𝜀𝜀 + 𝑎𝑎𝑔𝑔 + �𝜀𝜀𝜀𝜀 ∗ 𝑎𝑎𝑔𝑔 �
1 3
g ∗ (ρl − ρv ) Cpagua ∗ �Tsup − Tsat �
2
= �μ ∗ hfg ∗ � � ∗� n � � ∗ 𝑑𝑑𝑑𝑑
tensup Csf ∗ h ∗ Pr fg

De esta forma se presenta la ecuación de trasferencia de calor presente en el interior de la

caldera. Una parte asociada al intercambio de calor entre los gases de combustión y el tubo, y otra
parte asociada a la evaporación del agua.

𝜀𝜀𝜀𝜀 − 𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑒𝑒 𝑑𝑑𝑑𝑑𝑑𝑑 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡

𝜀𝜀𝜀𝜀 − 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑑𝑑𝑑𝑑𝑑𝑑 𝑔𝑔𝑔𝑔𝑔𝑔
𝑇𝑇𝐺𝐺𝑥𝑥 − 𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑑𝑑𝑑𝑑 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡
𝑎𝑎𝑔𝑔 − 𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴𝐴 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔
𝑇𝑇sup − 𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡
ℎ𝑖𝑖 − 𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝑑𝑑𝑑𝑑 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖
𝑑𝑑𝑑𝑑 − 𝐹𝐹𝐹𝐹𝐹𝐹𝐹𝐹𝐹𝐹𝐹𝐹𝐹𝐹𝐹𝐹 𝑑𝑑𝑑𝑑 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑎𝑎 𝑢𝑢𝑢𝑢𝑢𝑢 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑑𝑑𝑑𝑑 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡
μ − Viscocidad del liquido
g − Aceleracion de la gravedad
ρl − Densidad del liquido
ρv − Densidad del vapor
tensup − 𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖 𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣
 Cpagua − Calor especifico del agua
Tsat − Temperatura de saturacion
Pr − Prandlt del liquido
Csf − constante experimental dependiente de la combinacion superficie
− fluido tabla 10.3 Cengel
𝑛𝑛 − constante experimental dependiente de la combinacion superficie
− fluido tabla 10.3 Cengel
hfg − Entalpia de vaporizacion del agua

Teniendo de esta forma la ecuación relacionada con el calor transferido se acude a la solución
discreta de este tipo de ecuación, esto se logra dividiendo la caldera en un número discreto de
partes y se determina el calor entregado para cada una de estas divisiones. Como se puede
observar la ecuación presenta un dato que es desconocido, 𝑇𝑇sup , de esta forma se debe realizar
un proceso iterativo en el cual se varia la 𝑇𝑇sup hasta el punto en que la ecuación se cumpla (el
calor de radiación y convección interior debe ser igual al calor de ebullición).

De esta forma, con algunos de los datos iníciales entregados al programa se pueden determinar
algunos datos que se requieren para solucionar la ecuación presentada anteriormente:

dx=Le/N; %determinación de la división de la caldera

d=Dp*2.54/100; %conversión del diámetro de tubo principal en metros
ds=Ds*2.54/100; %conversión del diámetro de tubos secundarios en metros
dA=d*dx*pi; %área de cada una de las divisiones en el diámetro mayor
dAs=ds*dx*pi; %área de cada una de las divisiones en el diámetro menor
TGi(1)=TG; %Temperatura inicial del gas

4.3.4.1. Calculo del calor cedido en el tubo principal

Es así como en este punto se tienen todas las variables necesarias para desarrollar el análisis del
calor cedido por la caldera. Es importante destacar que se desarrollara el cálculo del calor del tubo
principal por aparte del calor suministrado en la segunda etapa (tubos secundarios) ya que es
necesario determinar la temperatura con la que el gas ingresa a los tubos secundarios.

%cálculo del calor entregado por la caldera en el primer paso (paso por el
%tubo principal)

contador=1; %se inicia un contador el cual permitirá construir los

vectores relacionados con la temperatura, emisividad, absortividad, calor, etc.

for i=1:N %En este caso se inicia el proceso de cálculo del calor
trasmitido en el paso por el tubo principal. Se hace desde 1 hasta N
ya que este es el número de veces que se debe desarrollar este
proceso

%Ahora es necesario determinar las propiedades de los gases de combustión y del

vapor con lo cual poder determinar el balance de energía presente en cada caso.

Para determinar las propiedades tanto de los gases de combustión como del vapor de
agua en el interior de la caldera se decidió elaborar diferentes funciones que permitieran
un cálculo adecuado y más organizado, de esta forma tenemos:
[Cp_agua,visco,Pr,densl,densv,tensup,Tsat,hfg]=PROPvapor(PSI); %Propiedades del
vapor de agua presente en el interior de la caldera

[PaLC,PaLW,EW,EC,Eg]=emisividad(TG,d,exc,comb); %Emisividad de
los gases a la temperatura de estos

[PM_gas]=PROPPM_gases(X_h2o,X_co2,X_n2,X_o2); %Determinación
del Peso molecular del gas, necesario para cálculos posteriores

[Cp_gas]=PROPCp_gases(TG,X_h2o,X_co2,X_n2,X_o2,PM_gas)*1000; %Determinación
del calor especifico del gas teniendo en cuenta el exceso, la temperatura y el
peso molecular.

[dens_gas]=PROPdens_gas(TG,X_h2o,X_co2,X_n2,X_o2,PM_gas); %se establece la

densidad de los gases a partir de propiedades previas de la combustión

[K_gases]=PROPk_gases(X_h2o,X_co2,X_n2,X_o2,TG); %Es necesario

además determinar la conductividad de los gases para el cálculo.

[visco_gases]=PROPvisco_gases(X_h2o,X_co2,X_n2,X_o2,TG); %Se
establece la viscosidad del gas
[Pr_gases]=PROPPr_gases(X_h2o,X_co2,X_n2,X_o2,TG); %Por ultimo
se determina el Prantl de los gases.

Las funciones y sus respectivas ecuaciones se determinaron de la siguiente forma:

4.3.4.2. Análisis función PROP_vapor

En el programa EES, se buscaron todas las propiedades del vapor, las cuales se dejaron en función
de la PRESION, teniendo como referencia la temperatura de saturación.

𝑃𝑃1 = 𝑃𝑃𝑃𝑃𝑃𝑃 ∗ 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐(𝑝𝑝𝑝𝑝𝑝𝑝; 𝐾𝐾𝐾𝐾𝐾𝐾)

𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 𝑑𝑑𝑑𝑑 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑎𝑎 𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑠𝑠𝑠𝑠 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 𝑒𝑒𝑒𝑒 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐

𝑚𝑚𝑚𝑚𝑢𝑢𝑤𝑤 = 𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉(𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠; 𝑇𝑇 = 𝑇𝑇1; 𝑃𝑃 = 𝑃𝑃1)

𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣 𝑑𝑑𝑑𝑑𝑑𝑑 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣 𝑎𝑎 𝑙𝑙𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒

𝑟𝑟ℎ𝑜𝑜𝑣𝑣 = 𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷(𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆; 𝑇𝑇 = 𝑇𝑇1; 𝑋𝑋 = 1)

𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑𝑑𝑑 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣 𝑎𝑎 𝑙𝑙𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒

𝑟𝑟ℎ𝑜𝑜𝑙𝑙 = 𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷(𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊; 𝑇𝑇 = 𝑇𝑇1; 𝑋𝑋 = 0)

𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 𝑒𝑒𝑒𝑒 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 𝑙𝑙𝑙𝑙 𝑙𝑙𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒

𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑎𝑎𝑙𝑙 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆(𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊; 𝑇𝑇 = 𝑇𝑇1)

𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑑𝑑𝑑𝑑𝑑𝑑 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎

𝐶𝐶𝐶𝐶 = 𝐶𝐶𝐶𝐶(𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊; 𝑇𝑇 = 𝑇𝑇1; 𝑃𝑃 = 𝑃𝑃1)

𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 𝑑𝑑𝑑𝑑𝑑𝑑 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑑𝑑𝑑𝑑𝑑𝑑 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣 𝑎𝑎 𝑙𝑙𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒

𝑇𝑇1 = 𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇(𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊; 𝑋𝑋 = 1; 𝑃𝑃 = 𝑃𝑃1)

𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑎𝑎 𝑎𝑎𝑎𝑎𝑎𝑎 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑑𝑑𝑑𝑑 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡

𝑃𝑃𝑃𝑃 = 𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃(𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊; 𝑇𝑇 = 𝑇𝑇1; 𝑃𝑃 = 𝑃𝑃1)

𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉 𝑑𝑑𝑑𝑑𝑑𝑑 𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 𝑎𝑎 𝑙𝑙𝑙𝑙𝑙𝑙 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒

ℎ𝑓𝑓𝑓𝑓 = ℎ𝑔𝑔 − ℎ𝑙𝑙

𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑑𝑑𝑑𝑑 𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓

ℎ𝑔𝑔 = 𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸(𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊; 𝑇𝑇 = 𝑇𝑇1; 𝑋𝑋 = 1)

 𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑑𝑑𝑑𝑑𝑑𝑑 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣

ℎ𝑙𝑙 = 𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸(𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊𝑊; 𝑇𝑇 = 𝑇𝑇1; 𝑋𝑋 = 0)

𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑑𝑑𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙

Como se observo anteriormente se dejaron todas las condición de las propiedades en

términos de la temperatura, con lo cual es permitido establecer una tabla con toda y cada una
de las propiedades en función de esta variable. Es importante destacar que algunas de estas
propiedades no solo son función de la temperatura sino también de la presión, sin embargo
solo se considero el estado de trabajo (Presión de trabajo) de la caldera la cual se establece en
el catalogo.

De esta forma con la tabla obtenida se realizó una regresión de 6to orden para cada
propiedad, y fueron estas las ecuaciones que se usaron en la función de MATLAB.

TABLA OBTENIDA EN EES

TABULACION DE LOS DATOS DETERMINACION DE LA REGRESION

De esta forma la estructura de la función en Matlab quedara de la siguiente forma:

function [Cp_agua,visco,Pr,densl,densv,tensup,Tsat,hfg]=PROPvapor(PSI)

hfg=2.26604478E+06-3.00016176E+03*PSI+1.50669335E+01*PSI^2-6.28828355E-
02*PSI^3+1.64955216E-04*PSI^4-2.41492422E-07*PSI^5+1.50133862E-10*PSI^6;

tensup=5.94730601E-02-2.17258820E-04*PSI+1.16275997E-06*PSI^2-4.82540042E-
09*PSI^3+1.26226981E-11*PSI^4-1.84478759E-14*PSI^5+1.14584475E-17*PSI^6

Cp_agua=4.20124862E+03+2.00288818E+00*PSI-7.91035729E-03*PSI^2+3.83573161E-
05*PSI^3-1.13381343E-07*PSI^4+1.86960730E-10*PSI^5-1.31006694E-13*PSI^6;

visco=2.63332344E-04-1.93487754E-06*PSI+1.53044586E-08*PSI^2-7.73074218E-
11*PSI^3+2.35613980E-13*PSI^4-3.95402693E-16*PSI^5+2.80476115E-19*PSI^6;

Pr=1.65395033E+00-1.16272925E-02*PSI+9.48159614E-05*PSI^2-4.81322741E-
07*PSI^3+1.47092009E-09*PSI^4-2.47246483E-12*PSI^5+1.75573717E-15*PSI^6;

densl=9.64097852E+02-9.44621307E-01*PSI+5.10104262E-03*PSI^2-2.30997637E-
05*PSI^3+6.66754407E-08*PSI^4-1.08343991E-10*PSI^5+7.52398472E-14*PSI^6;

densv=7.95570730E-02+3.71713008E-02*PSI-2.81883658E-05*PSI^2+1.25705793E-
07*PSI^3-3.38096852E-10*PSI^4+5.29827449E-13*PSI^5-3.59424410E-16*PSI^6;

Tsat=3.69592613E+02+1.13444974E+00*PSI-7.00644362E-03*PSI^2+3.26239232E-
05*PSI^3-9.54865174E-08*PSI^4+1.56429451E-10*PSI^5-1.09212918E-13*PSI^6;

4.3.4.3. Análisis función PROPCp_gases

En esta función se hicieron cálculos respectivos en EES de la misma forma que la anterior pero
esta vez, con la propiedad Cp para todos los gases de combustión en función de la
temperatura.

CALOR ESPECIFOC PARA LOS PRODUCTOS DE LA COMBUSTION

Cp_CO2 Cp_H2O Cp_N2 Cp_O2 Temperatura
61,19 53,17 36,44 38,59 2341
61,19 53,15 36,43 38,58 2326
61,18 53,12 36,43 38,56 2311
61,17 53,09 36,42 38,55 2296
61,16 53,07 36,41 38,54 2281
61,15 53,04 36,41 38,53 2266
61,14 53,01 36,4 38,52 2251
61,13 52,98 36,39 38,51 2236
61,12 52,96 36,39 38,5 2221
61,11 52,93 36,38 38,49 2206
61,1 52,9 36,37 38,48 2191
61,09 52,87 36,37 38,47 2176
61,08 52,84 36,36 38,46 2161
61,07 52,82 36,35 38,45 2146
61,06 52,79 36,35 38,44 2131
 61,05 52,76 36,34 38,43 2116
61,04 52,73 36,33 38,42 2101
61,03 52,7 36,32 38,41 2086

Nuevamente se realizaron regresiones de 6to grado para obtener las ecuaciones respectivas
de cada una de las propiedades, es así como la función en MATLAB se establece de la siguiente
forma:
function [Cp_gas]=PROPCp_gases(TG,X_h2o,X_co2,X_n2,X_o2,PM_gas)

Cp_h2o=47.52+0.003085*TG-3.0043e-7*TG^2+5.8055e-12*TG^3+1.1741727e-15*TG^4-
2.2818e-19*TG^5+1.617e-23*TG^6; %Calor especifico del H2O en función de la
temperatura (base molar).

Cp_co2=5.87747509E1+1.69378690E-3*TG-4.10871235E-7*TG^2+7.05593844E-11*TG^3-
8.32045444E-15*TG^4+8.31301123E-19*TG^5-4.62437735E-23*TG^6; % Calor
especifico del Co2 en función de la temperatura (base molar).

Cp_n2=3.49586547E+1+9.11541491E-4*TG-1.65241239E-7*TG^2+2.61350248E-11*TG^3-
3.45622253E-15*TG^4+3.34620828E-19*TG^5-1.71323646E-23*TG^6; % Calor
especifico del N2 en función de la temperatura (base molar).

Cp_o2=3.67749156E+1+9.03531759E-4*TG-9.03901171E-8*TG^2+2.13236048E-11*TG^3-
3.63578634E-15*TG^4+4.39677364E-19*TG^5-2.69753077E-23*TG^6; % Calor
especifico del O2 en función de la temperatura (base molar).

Cpmolar_gas=(X_h2o*Cp_h2o)+(X_co2*Cp_co2)+(X_n2*Cp_n2)+(X_o2*Cp_o2); %Calculo
del calor especifico de los gases en base molar.

Cp_gas=Cpmolar_gas/PM_gas; %Determinación del calor especifico en base másica.

En este caso se considera el cálculo de el calor especifico en base molar ya que se usaron las
fracciones molares para determinar el efecto de cada componente, de esta forma para determinar
el calor especifico másico es necesario dividir por la masa molar del gas.

NOTA: Todas las funciones adicionales de propiedades siguieron estas metodologías de cálculo.

4.3.4.4. Modificación a función EMISTOT

Se utilizo la función “emistot” pero esta se cambio de nombre, se usaron las modificaciones
hechas en un trabajo previo y se cambiaron las formula de las fracciones molares del CO2 y de
H2O:

function [PaLC,PaLW,EW,EC,Eg]=emistot(TG,d,ex,comb); %función que permite

el calculo de la absortividad de los gases la cual es función de la temperatura
del gas.

if comb==1%seleccion del combustible para el cálculo de la fracción molar del co2

y del H2O para el cálculo de la emisividad (CH4)

xc=1/(1+2+((1+ex)*7.52)+(2*ex));%x_C(:,2); %fracción del bióxido de carbono

xw=2/(1+2+((1+ex)*7.52)+(2*ex));%x_W(:,2); %fracción del vapor de agua
 else %selección del combustible para el cálculo de la fracción molar del co2 y
del H2O para el cálculo de la emisividad (C2H6)

xc=2/(2+3+((1+ex)*13.16)+(3.5*ex));%x_C(:,2); %fracción del bióxido de carbono

xw=3/(2+3+((1+ex)*13.16)+(3.5*ex));%x_W(:,2); %fracción del vapor de agua

end

4.3.4.5. Análisis Función AMISTOT

La función “amistot”, funciona de manera muy similar a “emistot” pero esta primera hace su
análisis a partir de la temperatura de superficie. Es asi como esta variable también entra dentro de
de la iteración del calor.
function [abs_g]=amistot(TG,tsup,d,ex,comb₧) %función que permite el calculo de
la absortividad de los gases la cual es función de la temperatura de la
superficie.

Una vez con las propiedades necesarias se procede a iniciar un proceso iterativo para igualar los
calores de ebullición del vapor y el calor proporcionado por los gases por convección y radiacion.
Egas1(i)=Eg; %se inicia el valor del vector de la emisividad del gas para
poder tabularlo para su posterior análisis grafico.

Seguidamente de la determinación de todas las

propiedades tanto de los gases como del vapor se procese
a hacer el análisis de intercambio de calor en cada uno de
los tramos del tubo primario, esto se realiza partiendo de
la suposición de la temperatura de superficie del tubo, la
cual es muy cercana a la temperatura del vapor.
for x=1:100 %se inicia un proceso secuencial en un valor de 1 el cual
determinara la variación de la temperatura de superficie hasta que los calores en
la ecuación del balance de energía sean iguales.

Suponemos una temperatura de superficie, con la cual se

calcula la absortividad y se encuentra el Q ebullición.
tsup=Tsat+x*0.1; %temperatura tentativa de la superficie de los tubos

[abs_g]=amistot(TG,tsup,d,exc,comb); %Calculo de la absortividad de los gases

con la temperatura de superficie supuesta

Ag1(contador)=abs_g; %Tabulación de los datos de la absortividad en un

vector para su posterior análisis grafico

%De esta forma se procede a calcular el calor de ebullición Q_ebullición (Q

DERECHA)

g=9.81; %gravedad

Csf=0.0130; % constante experimental dependiente de la combinación superficie

n=1; % constante experimental dependiente de la combinación superficie

qebu=visco*hfg*((g*(densl-densv))/tensup)^(1/2)*((Cp_agua*(tsup-Tsat))
/(Csf*hfg*Pr^n))^3; %calor de ebullición por unidad de área

Qder=qebu*dA;%calor de ebullición por sección

Una vez encontrado el Calor de ebullición, es

necesario determinar el calor emitido por los gases
por medios de radiación y convección en el interior
de los tubos.
%Caculo calor de radiación y convección de los gases (Q_izquierda)

%Es importante primero determinar el coeficiente de convección interior

H_interior

%Calculo del flujo de combustible y aire

fcomb=fcombft*7.86578333e-6; %Flujo de combustible en la caldera (se realiza una

conversión de unidades para trabajar en sistema internacional m^3/s)

mcomb=denscomb*fcomb; %Determinación del flujo másico de combustible que circula

por los tubos de la caldera

%Ahora es necesario determinar la masa de aire que reacciona con el combustible

en función del combustible seleccionado y el exceso.

if comb==1 %determinación de la masa de aire en la combustión de CH4

maire=(mcomb*(1+exc)*2*4.76*28.97)/16.043; %Flujo másico de aire del CH4

else %Determinación de la masa de aire en la combustión de CH4

maire=(mcomb*(1+exc)*3.5*4.76*28.97)/30; %Flujo másico de aire del C2H6

end

mgases=mcomb+maire; %flujo másico de gases (masa de combustible sumado a

la masa de aire necesario para la combustión

vm=mgases/(dens_gas*pi*d^2/4); %Determinación de la velocidad media de los gases

Re=vm*d*dens_gas/visco_gases; %Calculo del Reynolds

Nu=4.36+(0.065*(d/Le)*Re*Pr_gases)/(1+0.04*((d/Le)*Re*Pr_gases)^(2/3));
%Determinación del Nusselt

hi=Nu*K_gases/d; %Calculo del coeficiente de convección interna

sigma=5.67e-8;

De esta forma ya se puede proceder a determinar el calor

cedido por el combustible de la siguiente forma:
Qizq=dA*(Es*Eg*sigma*TG^4-abs_g*sigma*tsup^4)/(Es+abs_g+(Es*abs_g))+hi*dA*(TG-
tsup); %Calor de radiación y convección interno cedido por los gases de
combustión a su paso por la caldera

Finalmente Se procede a comparar los dos calores

obtenidos, y se establece la condición para la cual el valor
 de estos calores es muy cercano, para nuestro caso se
considera que una diferencia de 20 Watts es insignificante
comparado con los valores del calor cedido total por
fracción de tubo “(delta_Q<10)|(delta_Q>-10)”. Cuando
esto sucede se procede a actualizar una serie de vectores y
en este punto termina el proceso iterativo, por lo tanto se
sale del FOR correspondiente a la iteración, y se procede a
pasar a la siguiente fracción del tubo.
delta_Q=Qizq-Qder; %diferencia del calor interno y externo

if (delta_Q<10)|(delta_Q>-10) %si la diferencia de calor es menor de

10 y mayor de -10 se procede a continuar con la actualización de los datos. Si no
se cumple la condicion, entonces se continua aumentando la temperatura de
superficie hasta que el delta_Q cumpla con la condición establecida

tsupi(i)=tsup; %Vector de temperatura de superficie para cada

una de las secciones

Qizqi(i)=Qizq; %Calor por ebullición real para la sección

analizada
Qderi(i)=Qder; %Calor por radiación y convección cedida por el
combustible en la sección analizada
Tsat;
Break %Se procede a salir del condicional, ya que la
condición fue cumplida
Break %Se procede a salir del “for” que aumenta la
temperatura de la superficie ya que se cumplió con la condición estipulada
end

end

Teniendo los datos para este primer tramo del tubo

principal de la caldera, se procede a actualizar los distintos
valores y continuar para el cálculo del siguiente diferencial.
TGi(i+1)=TG-Qizq/(mgases*Cp_gas); %Determinación de la temperatura del
combustible a la salida de la sección analizada anteriormente
TG=TG-Qizq/(mgases*Cp_gas);

Ag1(contador)=abs_g;

contador=contador+1; %Actualización del contador que permite tabular los

datos en vectores.
end

De esta forma se obtienen todos los calores asociados al tubo principal de la caldera, con lo cual se
procede a determinar el calor total producido en esta primera etapa.
Qtotal=0; %Se define la variable que llevara a su cargo el calor total del
tubo principal.

for i=1:N
 Qtotal=Qtotal+Qizqi(i); %Se procede a sumar todos los calores calculados
anteriormente y tabulados en el vector “Qizqui”.

end

4.3.5. Tubos secundarios:

El procedimiento para la determinación del calor cedido en el paso de los gases de combustión
por los tubos secundarios es el mismo que se hizo previamente para el tubo principal. Sin
embargo, es necesario tener en cuenta una serie de consideraciones:

• Se analiza un tubo por separado y luego se multiplica por el número de tubos

• El área de transferencia se ve afectada, esta depende del diámetro y numero de
tubos secundarios.
• El flujo másico de gases se divide proporcionalmente a la cantidad de tubos
presentes en la segunda etapa.

4.4. CATALOGO SELECCIONADO

Con la ayuda del catalogo de la compañía “HURST BOILER & WELDING COMPANY, INC”, se
obtuvieron los datos requeridos para hacer un análisis de resultados.

Este catalogo no tiene datos de numero de tubos y longitud efectiva, sin embargo proveen el valor
del área de transferencia de calor.

De esta forma se procedió a modificar el catalogo original, para poder obtener una ficha con las
especificaciones necesarias del software realizado, teniendo en cuenta consideraciones como la
invariancia del área del catalogo o dimensiones irregulares que puedan variar de forma
significativa los datos obtenidos.

Garantizando que el area de transferencia de calor sea la especificada por el catalogo, se calculo el
numero de tubos del segundo paso y se asumió una longitud efectiva como el 80% de la longitud
total de la casco.

𝐴𝐴𝐴𝐴𝐴𝐴𝑎𝑎𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 = 𝑝𝑝𝑝𝑝 ∗ 𝐷𝐷𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 ∗ 𝐿𝐿𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 + 𝑝𝑝𝑝𝑝 ∗ 𝑁𝑁𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 ∗ 𝐷𝐷𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 ∗ 𝐿𝐿𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒

En el catalogo se especifica, que para el consumo de combustible dado se utiliza GAS NATURAL,
con una composición de aproximada:
 COMPOSICION PORCENTAJE
Metano 75-95 %
Etano 0,3-12 %
Propano 0,1-4 %
Otros 5.

Basado: En calidad de gases – ECOPETROL

http://www.ecopetrol.com.co/documentos/CalidadGases2004.pdf
Por lo cual se obtienen resultados aproximados si se trabaja preferiblemente con el METANO
PURO.

CATALOGO ORIGINAL “Caldera Series 200 esp.pdf”

FICHA MODIFICADA “CALDERA FICHA.PDF”

 5. ANALISIS DE RESULTADO:

Los análisis de las calderas del catalogo se realizaron con un exceso de 0% de aire y con metano
puro.
• CALDERA 15 BHP

%𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸 = 12,87

• CALDERA 20 BHP
 %𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸 = 12,67

• CALDERA 50 BHP

%𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸 = 11,12
• CALDERA 150 BHP

%𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸 = 8,54

• CALDERA 800 BHP

 %𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸 = 5,08
5.1. TABLA DE RESULTADOS

POTENCIA TEORICA [BHP] POTENCIA SOFTWARE [BHP] % ERROR

15 13,06 12,87 %
20 17,46 12,67 %
50 44,44 11,12 %
150 137,19 8,54 %
800 759,40 5,08 %

6. CONCLUCIONES

• En el análisis realizado se obtuvieron datos que concuerdan con los datos del catalogo, en
este caso se presentan errores que varían del 1 al 13 %. Estos errores se presentan debido
a que el combustible usado es una aproximación del gas natural con que fue determinado
el valor del catalogo.
• El calor brindado por la caldera presenta un aumento en el momento en que los gases
pasan del tubo principal a los tubos secundarios. Esto se debe al aumento en área de
transferencia de calor (ver graficad de calor vs longitud efectiva).
• La temperatura de superficie de los tubos varía de 3 a 10°C respecto a la temperatura de
superficie. De esta forma el cálculo del calor de ebullición se llevo a cabo considerando
una ebullición nucleada adecuada para estos casos.
• Los catálogos de calderas no brindan toda la información necesaria para realizar el cálculo
teórico, por lo cual es necesario estimar algunos de los valores, teniendo en cuenta que no
se sobrepasen los datos entregados por el catalogo.

7. BIBLIOGRAFIA
• CENGEL, Yunus; Transferencia de calor, 3ed; Mc Graw Hill Interamericana:
Mexico:2005.
• BAUKAL, Charles; Industrial Burners Handbook: 2003.