Procesos de combustion: modelos

webserver.dmt.upm.es

Procesos de combustion: modelos

Capitulo 16Procesos de combustion: modelosModelos te6ricos de combusti6nEn el capitulo anterior se ha presentado una descripci6n eminentemente fenomenol6gica de10s procesos de combusti6n. Se trata ahora de presentar algunos modelos sencillos que sirvenpara explicar y comprender mejor estos procesos, planteando las ecuaciones generales desdeun principio, siquiera para tener una visidn global del problema (hoy dia ya se han resueltoalgunos casos complejos con ayuda de ordenador).Correspondiendo a la caracterizaci6n de las llamas que se hizo antes, 10s modelos decombusti611 pueden tambiCn clasificarse atendiendo a diferentes aspectos:- Seg6n la dependencia espaciotemporal: uniforme, estacionario, unidimensional (ID),bidimensional (2D), tridimensional (3D), transitorio.- Segiin la condicidn de mezclado: premezclados (combusti6n homogCnea y ondas dedeflagraci6n y detonaci6n) y separados (llamas de difusi6n).- Seg6n la propagaci6n: homogCneo (en reposo, bien agitado), onda de deflagraci6n(subs6nica) y onda de detonaci6n (supers6nica).- Seg6n el tipo de flujo: laminar y turbulento.- Seg6n la fase del combustible: gaseosa, condensada (pero con evaporacibn), ycombusti6n heterogCnea propiamente dicha (en fase s6lida).- Seg6n la convecci6n: sin convecci6n, con convecci6n natural, con convecci6nforzada.- Seg6n 10s efectos de compresibilidad: incompresible, compresible.Sin embargo, 10s aspectos que se van a considerar en mis detalle a continuaci6n se hanclasificado de una manera mis sencilla y globalizadora, atendiendo a las disciplinas bisicasrequeridas para su comprensi6n, de la manera siguiente:


384 I. Martinez: TERMODINAMICA BASICA Y APLICADA- consecuencias de la Termodindmica del equilibrio:- composici6n de equilibrio- poder calorifico- temperatura de combusti6n adiabitica- consecuencias de la Termodinimica de la evoluci6n:- combustidn controlada por la cinCtica quimica:- combustor en reposo- combustor bien agitado- combusti6n controlada por difusi6n- llama bidimensional- llama axilsimCtrica- llama esf6rica- combusti6n controlada por onda de deflagraci6n o de detonaci6n- llama de premezclaLo m6s dificil de modelizar es la cinCtica quimica, y por eso las simplificaciones suelen serm6s Wsticas. En orden creciente de complejidad, se pueden citar 10s siguientes modelos:- Adici6n de calor a un gas inerte. Este es el modelo normalmente usado para el estudiode 10s ciclos termodindmicos correspondientes a 10s motores tCrmicos, per0 tambiCnes 6til para estudiar p.e. 10s campos de velocidades y temperatura en una llama dedifusi6n si se supone una deposici6n superficial de calor en la zona apropiada.- Reacci6n global irreversible monoetapa, que puede ser de primer orden (FjP), desegundo orden (F+O+P), o en general de orden arbitrario, correspondiendo a unareacci6n global F+v~O+V~N+V~P, donde F, 0 y N representan el combustible, eloxidante y una especie inerte, respectivamente, y P representa todos 10s productos(incluidos 10s inertes). El cambio de composici6n se analiza a partir de la ley develocidades de Arrhenius (ver mis adelante CinCtica quimica): DcF/Dt=-BacFacobc~exp(-T,lT), siendo B,, a, 6, c y Ta constantes dadas, y una leyconespondiente de deposici6n tCrmica asociada a la de producci6n de especies, dadapor DqlDt=-(h,lp)Dc~IDt. Muchas veces se supone que la reacci6n no genera mdscantidad de sustancias (F+vO+(l+v)P) y que ademds 10s efectos de dilataci6ntCrmica son tambiCn despreciables pese a1 gran salto de temperatura.- Otros modelos mis complicados se estudian con ayuda de hip6tesis de estadoestacionario para 10s radicales mis reactivos, hip6tesis de equilibrio parcial para lasreacciones mds rdpidas, desarrollos asint6ticos (entre 10s que cabe destacar 10s deenergia de activacidn alta), truncamientos en el ndmero de especies o en el ndmero dereacciones simples a contemplar, etc., y normalmente con ayuda de ordenador.N6tese tambiCn que debido a1 aumento exponencial de la velocidad de reacci6n con latemperatura, las reacciones quimicas en las zonas fn'as pueden suponerse congeladas,y en las zonas muy calientes pueden suponerse en equilibrio (y en ellas, a1 ser lasreacciones de combusti6n muy irreversibles, no podri quedar a la vez combustible yoxidante).


Cap. 16: PROCESOS DE COMBUSTION: MODELOSEcuaciones generalesSe desea conocer la evoluci6n de 10s campos de concentraciones de especies, de velocidades,de temperatura y de presi6n en un cierto dominio fluido (finito o infinito), es decir,determinar:yi (2, t), .'(X', t), T(2, t), ~(29 f )donde, como es traditional en la teoria de la combusti6n, se ha usado la fracci6n mhsica yipara caracterizar la composici6n, aunque se puede pasar a fracciones molares Xi, o adensidades pi, o a concentraciones volumktricas ci, fticilmente a travCs de las relaciones:Como ya se analiz6 en el Cap. 2, esta descripci6n como medio continuo estri basada en queexiste una escala de tamaiios (particula fluida) en 10s cuales se puede considerar que sealcanza el equilibrio termodinhmico en tiempos inferiores a 10s de inter&, y por tanto suestado queda univocamente descrito (en ejes m6viles) dados la energia U, el volumen V y lascantidades de sustancia ni.Aunque en el Cap. 2 se supuso que las ni eran conservativas y en 10s sistemas reactantes nolo son, el observador siempre puede escoger un tiempo de muestreo mucho mayor que eltiempo caracten'stico de reaccidn quimica y considerar el sistema en equilibriotermodinrimico total, como se hizo en el Cap. 9, o escoger un tiempo de muestreo muchomenor que el tiempo caracteristico de reacci6n quimica y considerar el sistema en equilibriotermodintimico quimicamente congelado, es decir, con 10s T, p y pi que corresponderian a1estado de equilibrio con las ni conservativas, pero que ahora irtin evolucionando lentamentecon el tiempo. Este punto de vista de equilibrio local parcial es el que se contempla en ladinhmica de 10s sistemas reactantes.Excepto la ecuaci6n de la cinktica quimica, todas las demris ecuaciones generales que sepresentan a continuaci6n han sido ya desarrolladas en capitulos precedentes, pero se hapreferido presentarlas aqui de nuevo.Partiendo de las ecuaciones (5.1) a (5.3) aplicadas a cada componente (ver (1 1.20) y (1 1.21)),y con ayuda del Teorema del Transporte de Reynolds (5.9, se obtienen las ecuaciones deconservaci6n para un volumen de control:


386 I. Martinez: TERMODINAMICA BASICA Y APLICADAdonde wi es la velocidad de producci6n misica de la especie i por unidad de volumen, y no seha tenido en cuenta el posible acoplamiento radiativo ya que, aparte de ser extremadamentecomplicado, en 10s sistemas gaseosos que se van a analizar es despreciable para sistemas detamaiio usual. Ahora, tomando un volumen unitario inmdvil respecto a una referencia dada(CA=O) y aplicando el Teorema de Gauss-Ostrogradski (5.20), se pasa a una forma diferencialequivalente a las ecuaciones (5.22) a (5.24):(vdase tambidn (1 1.20) y n6tese qUe pi = pyi, qUe ci = c + cdi , Cyi = 1, CyiCdi = 0 YC wi = 0).Las relaciones adicionales que ligan las variables intermedias introducidas serian:- ecuaciones constitutivas del equilibrio: p = p(yi ,T, p) y h = h(yi, T, p) , donde 10scoeficientes que aparezcan (p.e. Mi, c,; ....) deben ser datos.- ecuaciones constitutivas de la evoluci6n (transporte, y transformaci6n quimica):Jdi = Jdi(yi,~,~), += Y(Yi,~,p), a = a(Yi,~,P) y w = w(yi,T,~), donde 10scoeficientes que aparezcan (p.e. Di, p, k,..) deben ser datos.- condiciones iniciales y de contorno para las variables independientes (16.1).Rara vez se hace uso del conjunto de ecuaciones generales anteriores, incluso en lamodelizaci6n numCrica, ya que en 10s problemas de inter& no aparecen a la vez todos 10sfen6menos en la misma proporci6n y es necesario economizar esfuerzos.Aunque las ecuaciones anteriores serian tambiCn aplicables a flujos multifisicos yturbulentos, en estos casos resulta impracticable buscar una soluci6n determinista delcomportamiento detallado y se ha de recurrir a mCtodos estadisticos para analizar 10spromedios, estableciendo hip6tesis adecuadas sobre las fluctuaciones.


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 387La aproximaci6n mis utilizada (sobre todo cuando la difusi6n es dominante) es el llamadomodelo de Shvab-Zel'dovich (1950), cuyas hip6tesis son las siguientes:1. Se supone regimen estacionario: alat=O.2. Se supone que la presi6n es uniforme: Vp = 0. Esto s610 implica que se van adespreciar las fuerzas misicas (la flotabilidad) y las viscosas, y que 10s gradientes develocidad no son muy grandes (p.e., a1 pasar el fluido a travCs de una onda dedeflagraci6n la velocidad se multiplica por 5 aproximadamente y la presi6ndisminuye un 2%).3. Se supone mezcla de gases ideales, y por tanto las ecuaciones constitutivas de lamateria en el equilibrio son:p=-CCX,M,RUT i=1= P = -- pM con M I ~ X ~ M ,TRui=lR~TC(Y~IM~)C4. Se suponen despreciables 10s flujos cruzados de transporte, y que las difusividadesmisicas Di de cada especie respecto a1 resto de la mezcla o a la especie dominanteson todas iguales entre si (es bastante aproximado, except0 para el hidr6geno) y a suvez iguales a la difusividad terrnica de la mezcla a; es decir: Di=a, o lo que es lomismo, Le=1 (12.37).Entonces, las ecuaciones constitutivas de la materia en la evoluci6n (ecuaciones detransporte) serAn:4 * 4pGiyi = -pDiVyi (ley de Fick) con vdi = vi - v (16.1 1)' - j(vivj + v jvi) - (ip - p ,)~.JS,7, - (ley de Stokes)(16.12)4q = -kVT + E(pyicdih,) (ley de Fourier generalizada)i=lEn esta dltima se puede poner I$ en lugar de hi, ya que el transporte asociado deenergia termica sera mucho menor que el de energia quimica, si el sistema es muyreactivo.5. Se considera una dnica reacci6n quimica, y fuertemente irreversible, es decir, cuyavelocidad hacia adelante (consumo de reactivos) es muy superior a la velocidad haciaatris (consumo de productos para regenerar 10s reactivos), y para la ecuaci6nconstitutiva de la evoluci6n (cinitica quimica) se toma:


I. Martinez: TERMODINAMICA BASICA Y APLICADAreactivos "iwi = M i[{- vi) n ) Baexp- -Ea (ley de Arrhenius)j=1 RuTreactivosproductosSi esta simplificaci6n resulta inapropiada en un problema (porque la velocidad haciala izquierda sea importante o porque haya que considerar "arias reacciones), seaplicaria la ecuaci6n general:R reactivos a productosn,, reactivos .= ~ ~ x [ v ~ ~ - v l ) k / , C?r=l j=1reactivos;productqscon las R ecuaciones x vjr M, -+ xvjr~, (16.15)j=l j=1siendo kfi y kbr las constantes de velocidad hacia adelante (f=forward) y hacia atris(b=backward), cuya expresidn desarrollada en terminos de la temperatura es como en(16.14), y K:=kf,lkbr es la constante de equilibrio quimico de la reaccidn r expresadaen concentraciones, como se deduce de (9.21). Para las reacciones elementales esC=3 (ver mis adelante Cinetica quimica), per0 en total C es del orden de decenas y Rpuede llegar a la centena.Con las simplificaciones anteriores las ecuaciones de balance de especies (16.6) se reducen a:La ecuacidn de la energia, pasando a entalpias con ht=e+plp y haciendo uso de la hip6tesisde regimen estacionario, de la de presi6n uniforme (implica ht= h) y de la de gases ideales, sereduce a V - (pch) = -V . a, que combinada con (16.13), (16.11) y (16.10) toma la forrna:Si ahora se hace el cambio de variables:


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 389es fficil comprobar que tanto 10s balances de especies como el balance de energia se reducena la misma ecuaci6n diferencial:y restando dos a dos entre las C+1 ecuaciones (16.19) se puede generar un sistemamatemfiticamente equivalente en el que s610 una de ellas retuviese el tCrmino independiente,que es el iinico no lineal.N6tese que 10s dos tCrminos del primer miembro de (16.19) representan la convecci6n y ladifusi6n, y el segundo miembro la generaci6n (el tCrmino de acumulaci6n local no aparecepor haber supuesto regimen estacionario).A veces es importante la utilizaci6n de las ecuaciones de conservaci6n no para un elementode volumen sin0 para el salto a travCs de una discontinuidad (llama, superficie liquido-gascon evaporaci6r1, etc.), las cuales se obtienen sin dificultad a partir de las anteriores. Por otraparte, algunos problemas son tan sencillos que ni siquiera es necesario recunir a 10s balancesen forma diferencial y se usan ecuaciones de conservaci6n globales (algebraicas).A continuaci6n, y antes de entrar en el estudio de la evoluci611, se van a considerar lasconclusiones que la Termodinimica del equilibrio impone a todo proceso de combusti6n.Termodindmica de la combustionSe tratan aqui 10s condicionantes que la Termodinimica del equilibrio impone a 10s procesosde combusti6n, siendo todos ellos consecuencia directa del estudio general sobre sistemasreactantes, realizado en el Cap. 9.La Termodinirnica del equilibrio no es capaz de predecir si un sistema abandonado en unascondiciones iniciales (y con ciertas condiciones de contorno) reaccionarfi o no, ni quCproductos se formarian si reaccionase, ni a quC velocidad. Lo que si define es si una reacci6ndada es posible o imposible (en esas condiciones), y si por otras razones se sabe quCproductos van a aparecer, la Termodinimica del equilibrio define la composicidn que sealcanzaria en el estado de equilibrio, y el balance energCtico, esto es, la energiaintercambiada con el exterior (normalmente en forma de calor) en funci6n de 10s estadosinicial y final, o la temperatura final si se prescribe el valor del intercambio energCtico.Com~osici6n de equilibrioEl estado de equilibrio de un sistema reactante queda descrito por un conjunto de valoresT,p,ni, siendo ni la cantidad de componente i (i=1,2, ... C) de masa molar Mi y C el n6mero deposibles componentes considerados en el modelo (se elige empiricamente). LaTermodinimica no enseiia quC productos aparecen; p.e., a1 reaccionar CHL, y O2 podrfaformarse C02, o CH30H, o C6HI2o6, etc., y aunque parece obvio que a alta temperatura lasmolCculas habrfin de ser muy sencillas, el hollin existe a altas temperaturas y no lo es. LaTermodinfimica enseiia la direccidn quimica en la que puede avanzar la reacci6n (viabilidad),


390 I. Martinez: TERMODINAMICA BASICA Y APLICADAla composici6n en el equilibrio, la sensibilidad respecto a perturbaciones fisicoquimicas(estabilidad), la temperatura que se alcanzaria en la combusti6n adiabsitica o la energia aevacuar en la combusti6n isoterma, etc.El hecho de que la viabilidad de una reaccidn sea positiva, A>O (siendo A5Cvfii la afinidadquirnica), no significa que vaya a ocurrir en realidad; p.e. para la reacci6nCH4+202=C02+2H20 a P=298 K y p@=100 kPa es A@=0,82 MJ/molcH,, per0 tarnbiCnAe=O, 11 MJImolc~, para CH4+(1D)02=CH30H.Sistema reactante significa que ni puede variar por composici6n y descomposici6n de otrasentidades mis elementales (10s sitomos de las molCculas). Se supone conocida la composici6nde cada molCcula en funci6n de 10s itomos a partir de sus reacciones de formaci6n:siendo E el n6mero de entidades o elementos conservativos y C el n6mero de componentesreactantes. TambiCn se suponen conocidas las reacciones globales entre 10s componentes:siendo R el n6mero de reacciones independientes (a la relaci6n R=C-E se le llama regla deVolterra). Por ejemplo, si basta con considerar 10s sitomos de C y 0 como 6nicas entidadeselementales conservativas, y s610 hay molCculas de C, 02, CO y C02 el n6mero dereacciones independientes es R=C-E=4-2=2, que pueden elegirse de entre C+02=C02,C+(1/2)02=CO y CO+(1/2)O2=CO2).Una consecuencia de la regla de Volterra de gran importancia para la cinCtica quimica, es queen un sistema reactante s610 existen R velocidades de reacci6n independientes (que semodelizarsin p.e. con la ley de Arrhenius), y por tanto, s610 R de 10s wi en (16.6) sonindependientes: p.e. para el sistema [C,02,C0,C02], si se toman como independientes wco ywco2, se tendri que wc=2wco-wco2 y wo2=-2wco+2wco2. En general, las relaciones entrelas velocidades de producci6n de 10s demis compuestos a partir de las consideradas comoindependien tes sersin:llamadas relaciones de Shvab-Zel'dovich, y donde 10s Vi, son 10s coeficientesestequiomCtricos de las relaciones de composici6n de 10s demis productos i = E+l ... C apartir de 10s productos base r = l...R (Mi=C~i,rMr).La Termodindmica enseiia que la composici6n de equilibrio ni (ni = nio + vi{) de un sistemade C componentes, integrado por n, cantidades de elementos conservativos (ne es conocido apartir de la composici6n inicial, ni,), en presencia de un bafio a T y p dados, vienedeterrninada por las C=E+R ecuaciones:


Cap. 16: PROCESOS DE COMBUSTION: MODELOSCn, = c 1veini e = 1,2,. .., E (E ecuaciones conservativas)i=l(16.23)dG(Typ'ni) = 0 r = 1.2, ..., R(R ecuaciones de equilibria)T,,,siendo G(T,p,ni) el potencial de Gibbs.Las ligaduras de equilibrio se pueden poner de otra forma, en funci6n de la afinidad A y laconstante de equilibrio, a partir de la expresi6n general:CdG = - SO + Vdp + xpidni = -SdT + Vdp - C qdtri=l r=l 1Rde donde se deduce que las condiciones de equilibrio serfin Ar=O. Introduciendo la constantede equilibrio K(T,pB) como se hizo en (9.20):es decir, con A@(T,pB) = -Cvi~i'(T,p') 5 RTlnK(T,pB), se concluye que las R condiciones deequilibrio de (16.23) se pueden poner como:conforme se desarroll6 en (9.24) y (9.25). Normalmente se utiliza el modelo de gas ideal ypor tanto 10s coeficientes de actividad serfin la unidad, ~=1. La variable adimensionalKr(T,pe) se tabula parape=lOO kPa o se ajusta linealmente de la forma lnK=A+BIT, como seha hecho en (16.26). N6tese que a1 igual que la funci6n de Gibbs de un sistema con unareacci6n puede obtenerse a partir de las de formaci6n de cada compuesto por Gr=CviGfi, lasconstantes de equilibrio estfin relacionadas a travks de Kr=nKfiVi, siendo Vi 10s coeficientesestequiomCtricos de las reacciones de formaci6n.Como ejemplo, supdngase que se desea calcular la composici6n de equilibrio en lacombusti6n isobfirica y adiabfitica de acetileno en aire a 100 kPa. Para ello basta conocer lacomposici6n inicial, medida p.e. por la riqueza @ (supongamos que hay un 10% en volumen


392 I. Martinez: TERMODINAMICA BASICA Y APLlCADAde C2H2 en aire y por tanto @=1,32), y las constantes de equilibrio, que si se aproximan comolnK=A+BIT quedan especificadas a partir de las funciones de Gibbs de formaci6n y entalpiasde formaci6n de 10s compuestos que intervengan. Si se supone que s610 intervienen 10scomponentes del conjunto [C2H2,O2,N2,CO2,H20,CO,H2,H,O,0H,NO], se tienen C=l 1compuestos y E=4 entidades conservativas [C,H,O,Nl. La riqueza se define respecto a laestequiomem'a global:de tal manera que la ecuaci6n de la mezcla reactiva se puede poner como:siendo C la cantidad de gases presentes en 10s productos de la combusti6n de g moles decombustible. El problema queda planteado en funci6n de 13 incdgnitas (1 1 xi, C y latemperatura de combusti611 adiabfitica Tad) y las 13 ecuaciones siguientes:para H24 = C[2xcZH2 + 2xH20 + 2xH2 + XH + xOH]para o 2. $ = c[2xO2 + 2xC02 + XHZO + xCO + x0 + xNO]para c 2 ' )' 3,76 = C[2xN2 + xN0]- balance energetic0 (volumen de control en regimen estacionario con una entrada y unasalida):productosreactivosw+Q= Cli,hi(Tad>p)- znihi('earada,p Isiendo datos w (el trabajo aportado p.e. por un agitador externo) y 8 (el calor aportado siel proceso no es adiabfitico), y hi las entalpias de 10s compuestos quimicos separados.- equilibrio quimico; C-E=7 ecuaciones de acoplamiento quimico, del tip0 (9.20), p.e.:para C2H2 + $ 02 = 2C02 + H 20


Cap. 16: PROCESOS DE COMBUSTION: MODELOSpara H2 + f 0 2 = H20para CO + f o2 = C02para N2 + 0 2 = 2N0para H2 + 0 2 = 2H0para H2 = 2Hpara O2 = 20o cualquier combinaci6n'de ellas. Los valores de g@ y h@ que aparecen en estas ecuacionesde equilibrio son ficilmente accesibles en la bibliografia; s610 cabe sefialar que 10s valorespara el agua han de ser 10s del estado hipotdtico de vapor a pe=lOO kPa y p=298 K, esdecir h@(vapor)=-242 Wmol y g@(vapor)=-229 kJ/mol. La soluci6n de estas 13ecuaciones con 13 inc6gnitas en funcidn de 6 se representa en la Fig. 16.1.Poder calorifico de un combustibleSe denornina poder calorifico estiindar de un combustible (en su combusti6n con el oxigenodel aire), a1 calor que seria necesario evacuar, por unidad de cantidad de sustancia (o demasa, o de volumen) de combustible, en un reactor isobirico e isotermo en condicionesestandar (p@=100 kPa y p=298 K), entrando 10s reactivos separados y saliendo 10sproductos separados y en las mismas condiciones estiindar. Este poder calorifico asi definido,cambiado de signo, coincide pricticamente con la entalpia estandar de reaccidn h,@ definidaen el Cap. 9, aunque alli se requiri6 que la reacci6n se produjese en estado de equilibrio (parareacciones de combustidn las energias que entran en juego son tan grandes que esasprecisiones carecen de inter& priictico). Por tanto:


I. Martinez: TERMODINAMICA BASICA Y APLICADAFig. 16.1. Composici6n de equilibrio xi de las 11 especies consideradas, y temperatura adiabitica T d,en la combusti6n de acetileno y aire a 100 kPa, en funci6n de la riqueza 4. Para el rango deriquezas representado la fracci6n de acetileno en 10s productos es despreciable.uroductosreactivoscch?- ctiih? C c00PCS = -h, = - vihL= -c vir$ =-cncombustible i=l i=lsiendo hi@=hfi@ la entalpia estsindar de fonnaci6n de 10s compuestos tomando como nulas lasde 10s elementos y vi = tii / ticombustible 10s coeficientes que aparecen en la ecuacidn molar dela mezcla reactiva para la unidad de sustancia de combustible:uroductosque se llaman simplemente coeficientes de mezcla, y que son proporcionales a las fraccionesmolares xi ya que xi = vi I vi .Por ejemplo, el PCS del gas natural, supuesto metano puro, y sabiendo que la reaccidn esCH4+2O2=CO2+2H2O, sera: (valores tornados de la Tabla A3.8): PCS=-[(I).(-394)+(2).(-286)-(I).(-75)-(2).(0)]=891 kJ/molm,,,,=55,7 MJ/kgmetano.La contribuci6n a la entalpia de reacci6n de la entalpia de formaci6n de la molCcula delcombustible suele ser pequefia (un 8% en el ejemplo anterior), por lo que cuando no seconoce con precisi6n la estructura molecular del combustible (carbones, fuel6leos y otrosdestilados del petr6leo) se desprecia su efecto. Similarmente, si se ha determinado lacomposici6n de un fuel6leo como 86% de C y 12% de H en peso, su poder calonlfico seestima como el 86% del del carbon0 miis el 12% del del hidrbgeno, es decir,0,86.33+0,12.143=46 MJIkg. Por otra parte, el aire necesario para la combusti6nestequiomCtrica si que es independiente de la estructura molecular, y asi para este ejemplo se


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 395requeriria 0,86.(32/12)+0,12.(32/4)=3,25 kg de O2 por kilo de fuelbleo, y por tanto3,25 .(1 +3,3)= 14 kgairJkgfue161eo, para combustidn estequiomCtrica.Suele hablarse de poder calorifico superior e inferior (PCS y PCZ), refiriindose a si secontabiliza o no la entalpia esthndar de vaporizaci6n del agua producida (de toda, no s610 dela que estaria en fase liquida en equilibria con 10s gases de escape en condiciones esthndar), yla relaci6n entre ellos es: PC$=-hr8 y PCI~PCS-V~~~~~~~@,con hlvQ=44 kJ/mol=2,44 MJIkg.En realidad el PCI no es mhs que una variable auxiliar de gran utilidad para reducir cdculos,ya que normalmente interesa saber el poder calorifico en situaciones reales en las cuales latemperatura es superior a la del punto de rocio y por tanto toda el agua esth en forma devapor.. ,.Temperatura de co mbusti6n adiabaticaSe llama temperatura de combusti611 adiabhtica, Tad, a la que alcanzarian 10s productos de lacombusti611 en un sistema aislado tkrmicamente y mantenido a presi6n constante,inicialmente a T@ y p@. N6tese que aunque no fuese Cste el estado inicial, a1 ser el salto detemperatura grande, la contribuci6n de la energia sensible de 10s reactivos seria normalmentedespreciable. La temperatura alcanzada seria mayor si se mantuviera el volumen constante.El valor de Tad se obtiene del balance energetic0 para un volumen de control en regimenestacionario con una entrada y una salida (5.40) con Q = w = 0 y energias cineticas ypotenciales despreciables, esto es:reactivosproductosZ hi ('entrada 9 p) = Z 'ihi ('ad 7 p)reactivosreacf ivosreactivos8 88 8Z nihik ,P )+[ Z &hi(Tenfida.p) Z tiihik ,P )]=productosproductos-Q 8nihi(Tad,~)- 'ihik YP )](16-28)8 8donde hi(^ ,p ) = /$ es la entalpia estindar de formaci6n del compuesto i, y se hanfconsiderado separadas las sustancias tanto a la entrada como a la salida. Para el c6mputo de10s terrninos de variaci6n de entalpia de cada sustancia entre 10s dos estados termomecinicospuestos explicitos en (16.28) puede usarse el modelo de sustancia calorificamente perfects,aunque a1 ser el interval0 de temperaturas muy grande conviene tomar la variaci6n de cp conT. Normalmente la entalpia termomechnica aportada por 10s reactivos es despreciable (salvoque se use aire muy caliente), y para la entalpia termomecinica de 10s productos se suponengases ideales, per0 es importante destacar que como el agua esth condensada en el estadoQ 8(T ,p ), para que ello sea vdlido, es necesario tomar como referencia el estado hipotCtico devapor de agua a T8=298 K y p8=100 kPa.Para obtener Tad de (16.28), y despreciando la entalpia termomechnica de 10s reactivos,pueden seguirse varios caminos:


I. Martinez: TERMODINAMICA BASICA Y APLICADAPC1 =x vi [hi (Tad , pi=lproductosproductosi=lxproductosQ PC1PC1 = vicPi (Tad- T') Td =x+ productosi=lreactivos+ productossiendo PC1 = - vihy el poder calorifico de la reacci6n real (y no el del combustible ensu combusti6n completa con el aire) y vi = tii / ticombwtible 10s coeficientes de mezcla.Se usa (16.29) cuando se disponen de tablas o grfificos de hi-hi0, y el problema se resuelvepor aproximaciones sucesivas. Se usa (16.30) cuando en lugar de disponer de las entalpias sedispone de las capacidades tCrmicas (corno se ve, en este caso hay que hacer una integraci6nen cada iteracibn). Se usa (16.31) cuando se quiere aproximar con unos valores medios de cppara cada gas en todo el intervalo de temperaturas, que, aunque desconocido, no variarSidemasiado de un problema a otro en tCrrninos relativos (p.e. pueden tomarse como valoresmedios cp(C02)=54 J/mol, cp(H20)=47 Jlmol, cp(diat6micos)=34 Jlmol, cp(monoat6micos)=21 Jlmol). Se usa (16.32) cuando se quiere aproximar la capacidad tCrmica de 10s gases deescape por un cp medio de la mezcla en todo el intervalo de temperaturas (p.e. cp=34 Jlmol,ya que en la prfictica el nitr6geno es el componente mayoritario). En estos dos 6ltimos casosla temperatura de combusti611 adiabfitica es despejable explicitamente.,Aunque el balance energetic0 (16.28) puede interpretarse como que con el poder calorificode la reacci6n se calientan 10s gases de escape desde la temperatura estfindar hasta latemperatura adiabfitica, eso no tiene por quC corresponder a la realidad, donde lacomposici6n de equilibrio de 10s gases de escape iria variando con la temperatura de acuerdocon la constante de equilibrio (que liga la composici6n y la temperatura); per0 como laentalpia es funci6n de estado, si que se puede calcular por ese camino hipotktico.En cualquier caso, el problema se suele resolver por aproximaciones. Se empieza suponiendoque no hay disociaci6n (combusti6n completa) y se calcula la temperatura que se alcanzaria.Luego, si Csta resulta superior a unos 2200 K, se aiiaden las reacciones de disociacidncorrespondientes y se calcula la nueva composici6n de equilibrio a la temperatura anterior ycon ella un nuevo valor de la temperatura adiabfitica, y asi sucesivamente.En la Fig. 16.2 se esquematiza el resultado del cfilculo de la temperatura adiabfitica decombusti6n del metano en funci6n de la riqueza de la mezcla; la curva a corresponderia a lareacci6n ideal sin disociaci6n) con oxigeno puro, la curva b incluiria el efecto de un gasinerte (el nitr6geno del aire), y la curva c aiiadiria el efecto de la disociaci6n (aparici6n deCO Y H2).


Cap. 16: PROCESOS DE COMBUSTION: MODELOSFig. 16.2. Temperatura adiabitica de combusti6n del metano en funcidn de la riqueza de lamezcla; a, reacci6n ideal (sin disociaci6n) con oxigeno puro; b, efecto de un gasinerte (el nitr6geno del aire); c, efecto de la disociaci6n (aparicibn de CO y Hz).N6tese que la temperatura mixima se obtiene en la priictica con mezclas un poco msis ricasque la estequiomttrica; tambitn hay que tener en cuenta que un poder calorifico msis alto(mayor Ih,l) no implica una Tad mayor, pues el efecto de la disociaci6n puede ser dominante.Aunque pueda parecer parad6jic0, todavia hoy presenta m8s incertidumbre la medida directade la temperatura de combusti6n adiabitica que su c8lculo te6rico por el procedimientoanteriormente expuesto. La explicaci6n es la siguiente. Como se puede ver en el Cap. 21, latermometria por encima del punto de la plata (1235 K) se define por la ley de radiaci6n delcuerpo negro de Planck, y se rnide con aparatos que detectan la radiaci6n propia del objeto.Se pueden construir patrones de cuerpo negro a alta temperatura haciendo agujeros (para queno reflejen) en sMidos de alto punto de fusi6n (Tfp,ati,=2042 K, TfWOIframio=3700 K) y conellos calibrar la medida de altas temperaturas de superficies didas, per0 la medida de altastemperaturas en gases es todavia miis compleja porque a1 ser un sistema enrarecido presentamucha transmistancia y encima Csta depende fuertemente de la longitud de onda. Por ello, loque se hace es seleccionar un ancho de banda donde la transmitancia de la llama sea menor, yun tamafio de la muestra suficientemente grande para que la transmitancia sea realmentepequeiia, y considerar entonces que la emisi6n en esa banda (y la absorci6n) corresponde a lade cuerpo negro a la temperatura de la muestra.Como en la radiacidn propia de las llamas no siempre se encuentran bandas apropiadas, unmttodo muy utilizado es el de aiiadir alg6n otro elemento trazador, y el mejor es el sodio,cuyos vapores tienen una banda de absorci6n muy importante en el amarillo, cerca de 0,6 Fm(en realidad son dos lineas, una a 0,5890 pm y otra a 0,5896 pm). Se estima que laincertidumbre absoluta de este mttodo (el mejor) de medida experimental de la temperaturade combusti6n adiabfitica, denominado de la linea D del sodio, es de +20 K a 2000 K.Cinetica quimicaEl principal objetivo de la cinttica quimica es determinar la velocidad de reaccih i en unsistema de masa de control (para separar las variaciones de composici6n reactivas de lasconvectivas), la cual se define como (Cap. 9):


I. Martinez: TERMODlNAMlCA BASICA Y A PLICA DACpara una reacci6n 0 = C viMii=lsiendo ni la cantidad de sustancia i en la masa de control, [Mi] = ci su concentraci6n(moVms), pi su densidad (kg,m3), Wi la velocidad de produccidn m6sica (16.3) y Vi 10scoeficientes estequiomCtricos relativos a1 proceso reactivo que se considera, el cual debe serexplicitamente especificado, como se ha hecho en (16.33).Cuando se trataba de una reaccidn en equilibrio no tenia sentido distinguir entre reactivos yproductos, pero para 10s estudios cineticos lejos del equilibrio es muy importantedistinguirlos, por lo que las reacciones cineticas se escriben siempre asi:Por ejemplo, para 3H -+ H2+H se tendr6 C = 2 (1 = H, 2 = H2), vl' = 3, V$ = 0, vll' =1,v2" =1, vl = v1" - v1 ' = -2, v2 = v211 - v2' = 1. A veces en lugar de la flecha simple (+) sepone flecha doble (s), pero s610 por ahorrar escritura, pues lo que se quiere decir es que sevan a considerar dos reacciones cineticas, la que procede a derechas y su contraria, y no unareacci6n en equilibrio; de hecho, las velocidades de reacci6n serin muy distintas (s610iguales en el estado de equilibrio) y habri que distinguirlas adecuadamente (p.e. consubindices 1 y 2,6 f y b, del inglCs forward y backward, 6 + y -, 6 ' y ", etc.)La velocidad de reacci6n dependeri de la presi6n, la temperatura y la concentraci6n de cadaespecie (sea reactiva, no reactiva o catalizadora). Pensando en tertninos cinkticosmicrosc6picos, la velocidad de reacci6n seri proporcional a la frecuencia de colisiones de 10sreactivos, la cual seri proporcional a sus concentraciones, luego parece razonable que unprimer modelo para la dependencia de la velocidad de reacci6n con la concentraci6n apresi6n y temperatura constantes sea:que se denomina ley de acci6n de masas, llamindose a k constante de velocidad de reacci6n.La variaci6n de la composici6n de una especie reactante se obtiene combinando (16.33-16.35) para dar, segiin se mate de una sola reaccidn o de varias:CRCdc e dci-~'~,~)k,nc~"'j-? (16.36)dt Vj=1 dl r=l j=1= vi- = (viW-vi')kncjvj' o bien - = X(V"~,,Pero el modelo cinCtico basado exclusivamente en la probabilidad de choque de lasmolCculas estables que aparecen en (16.34) no ser6 muy realista porque no distingue entre


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 399choques elristicos e inelisticos, ni considera la aparici6n de radicales libres durante laevoluci6n. Para comprobarlo basta resumir algunos resultados experimentales, a una presi6ny temperatura dadas, en funci6n de la composici6n:Varias conclusiones se pueden deducir de estos y otros resultados:- La ley de acci6n de masas no es directamente aplicable a estas reacciones globales.- Unas veces la velocidad s610 depende de las concentraciones de 10s reactivos,como en (16.37), pero otras depende tambiCn de la concentraci6n de productos,como en (16.38).- A6n hay otras veces que la velocidad de reacci6n queda mejor descrita en funci6nde las concentraciones de otras especies transitorias que se detectanexperimentalmente, aunque no figuren ni como reactivos ni como productos en lareaccidn global, como pudiera haber sido la concentraci6n de molCculas OH en(16.37) o de hidr6geno at6mico en (16.38).- Si en lugar de considerar una reaccidn global como las anteriores, se considera unconjunto de reacciones elementales apropiadas, a veces es posible deducir lasvelocidades de reacci6n experimentales mediante aplicaci6n directa de la ley deaccidn de masas a las reacciones elementales.El 6ltimo punto anterior viene a resolver el problema de la cinCtica quimica, pues bastarri conencontrar el conjunto de reacciones elementales equivalentes a una dada, y aplicar la sencillaley (16.35) para determinar su velocidad de reaccidn, aunque la bdsqueda puede ser penosaporque no se sabe cdmo hacerla, ni si es 6nica, ni siquiera si existirri solucibn.Mecanismo de una reacci6nSe llama mecanismo de una reacci6n dada a1 conjunto de pasos o reacciones elementales que,aplicrindole la ley de acci6n de masas, es capaz de predecir la velocidad de reacci6n real.Seg6n eso, una reacci6n elemental es aquella cuya velocidad de reacci6n queda bien descritacon la ley de acci6n de masas, lo cual dependeri del rango de presiones, temperatura yconcentraciones existentes.Se llama orden de una reacci6n cuya velocidad dependa de las concentraciones en la formaempirica 6 = knciai , a la suma de 10s exponenetes, Ihi. Para reacciones elementales el


400 I. Martinez: TERMODINAMICA BASICA Y APLICADAorden coincide con la molecularidad o niimero de particulas que chocan. Asi, una mismareaccidn, p.e. A + B + C puede ser de primer o de segundo orden, como luego se veri, per0si es una reaccidn elemental, entonces es de primer orden.Las reacciones de primer. orden sirven para modelizar procesos de descomposicidn del tip0A + P en la que una molCcula A da lugar a uno o varios productos P. aunque tambiCn sirvepara reacciones bimoleculares A+B -+ P si las concentraciones son muy diferentes (cB >>cA). Se tiene entonces que:~ C A - k't- = -kcACB = -klcA -+ CA = C A,~dty la constante de reacci6n k' se mide experimentalmente a partir de:Pero las reacciones de primer orden no son posibles mis que si A era una molkculapreviamente excitada por una reacci6n de segundo orden o superior, ya que si A es unamolCcula estable no se descompone naturalmente. Fue Lindemannl en 1922 quien primer0propuso un mecanismo para explicar la ley de velocidades de descomposici6n de lasreacciones unimoleculares:A + M A * +M (excitaci6n tkrmica, ripida)A * +M 4 A + M (desactivacidn tkrmica, ripida)A* 3'~ (descomposici6n, lenta)Isiendo A* un estado energkticamente excitado de la molCcula A, que se llama compleipactivado (de tiempo de vida - 10-9 s), y M cualquier tip0 de molkcula, que normalmente seraA (aunque puede ser P o una especie inerte). Aplicando (16.36) a (16.41) se obtiene:y eliminando CA* a partir de (16.42) se llega a una ecuaci6n diferencial ordinaria de segundoorden que no es elemental, asi que, bien esa o bien el conjunto (16.42-43) habri de resolversenumCricamente (lo que es muy sencillo per0 poco ilustrativo). Pero si se hace la hipdtesis deestado estacionario, es decir que a1 ser A* muy activo no se acumula y, por tanto, se consumea la misma velocidad que se produce (lo cual, evidentemente no puede ser vilido a1 principioni a1 final de la reaccibn, per0 si durante la mayor pane del tiempo, entonces dcA*ldt = 0) y elsistema (16.42-43) se reduce a1. Lindemann, F.A., Trans. Faraday Soc. 17, p. 598, 1922.


Cap. 16: PROCESOS DE COMBUSTION: MODELOSque enseiia que la velocidad de una reacci6n monomolecular A + P, (6 / V = -dcAldt =kIIci), puede ser de 1% 22" grado; suelen ser de l9 grado a1 principio o si la presi6n es alta, yde 2"ado al final o si la presi6n es baja.Las reacciones de segundo orden sinen para modelizar procesos elementales de formaci6nde complejo activado, A+M+A*+M, de isomerizaci6n, A+M+B+M y de descomposici6n,A+M+B+C+M. Se comprende que 10s choques binarios implicados en las reacciones de 2"orden ser6n mucho m6s frecuentes que 10s ternarios, y estos que 10s cuaternarios, y por tantolas reacciones de 2"rden son las m6s corrientes. Estas tambiCn sirven para modelizarA+B+C + P cuando cc >> CA + CB, resultando entonces una ley de velocidades:que sine tambiCn para determinar k' a partir de medidas de concentraci6n en funci6n deltiempo.Las seacciones de tercer orden son menos frecuentes, pero a1 tener m6s grados de libertadpara disipar la energia sobrante de la reordenaci6n de enlaces, son m6s productivas, es decir,es mayor la fracci6n de choques ternarios que da lugar a cambio quimico que la de binarios,porque en Cstos la mayoria de las veces no se podr6 disipar esa energia y se volverin aromper 10s enlaces ya forrnados. Ejemplos tipicos de reacciones ternarias son:2H+M+ H2+M, 2N0+02+ 2N02, etc.En el estudio tedrico de la velocidad de una reaccidn global en funci6n de las velocidades delas reacciones elementales del modelo propuesto, conviene tener en cuenta las diferentesescalas de tiempo que normalmente aparecen. Por ejemplo, considCrese el esquema lineal(reacciones sucesivas):donde 10s ndmeros corresponden a las posibles reacciones a considerar (10s impares hacia laderecha y 10s pares a1 reds). Si un paso es muy lento (p.e. k5kl,ks, es decir, B es muy reactivo, entonces se puede considerar esecomponente en estado estacionario, es decir, dcB/dt ser6 despreciable en dcB/dt=klcA-k3cB,except0 en el instante inicial, en que dominarrin las condiciones realmente impuestas. Porotro lado, si un paso es muy ripido en 10s dos sentidos (p.e. k3,k4>>kl,k2,k5,k6) entonces sepuede considerar que esa reacci6n est6 en equilibrio auimico parcial y por tanto k3cg--k4cC=0.


I. Martinez: TERMODINAMICA BASICA Y APLICADAReacciones en cadenaComo ya se ha visto, para explicar la dependencia de la velocidad de reacci6n con lacomposici6n hay que postular un conjunto de reacciones elementales. De ellas, unas serin deiniciaci6n o generaci6n de comDuestos intermedios activos (10s que aparecen en elmecanismo y no en la reacci6n global), que en general serdn radicales libres (entidadesat6mico-moleculares con electrones no apareados y por tanto inestables a p y T ambientes,como el H, OH, CH3) aunque otras veces serin iones (sobre todo a alto T o en disolucidnacuosa). Otras reacciones elementales serin de propagaci6n de 10s compuestos intermedios(unos dan otros) y finalmente habri reacciones elementales de terminacidn que consumanesos compuestos intermedios.La propagacibn de la reacci6n puede ser por una secuencia o cadena lineal de pasoselementales, o por una secuencia o cadena ramificada donde la tasa de reproducci6n deradicales es superior a Ia tasa de consumo.Algunos ejemplos de mecanismos propuestos pueden ser ilustrativos. El mis sencillo es el dela descomposici6n del ozono en oxigeno, 203 + 302, que para concentraciones moderadas(cO3 = 0,25) a presi6n y temperatura ambiente, puede resumirse asi (10s valores de B, y To seexplican en el punto siguiente):B, [m3.m0l-1.s-11 , Ta [K]Iniciacibn: O3 + M + O2 + 0 + M 2.5 x lo8 11430Terminacibn: O3 + 0 + 202 5.2 x lo6 2090(16.47)aunque el modelo mejora aiiadiendo 10s pasos 20+M 2 02+M e incluso 10s contrarios a(16.47).El mecanismo de reacci6n mis clisico (con llama y todo) es el de la combusti611 delhidr6geno con bromo, H2(g) + Br2(g) + 2HBr(g), que es aplicable en un amplio rango depresiones, temperaturas y concentraciones:Iniciaci6n:Br2 + M -+ 2Br + M\Propagaci6n: Br +H2 + HBr +H 'H+Br2 + HBr +BrH+HBr + H2+BrTerminaci6n: 2Br + M + Br2 + M-,el cual, aplicando la hip6tesis de estado estacionario para 10s radicales (Br y H) proporcionaexactamente la ley de velocidad de reaccidn experimental expuesta en (16.38).Un mecanismo para la combustidn del hidr6geno en oxigeno H2 + (1/2)02 + H20, aunquecon un rango de aplicaci6n muy restringido (grandes vollimenes y bajas presiones), es:


Cap. 16: PROCESOS DE COMBUSTION: MODELOSm 3Iniciaci6n: H2 + O2 + 20H lo1* - 19500 Kmol sm 3Propagacibn: OH +H2 + H20 +H 1020 - 3000 Kmole sRarnificaci6n: H + O2 + OH + 0m32.2x1020 - 8300K } (16.49)mol - s.Para vol6menes pequefios, muchas reacciones de terminacibn ocurren en la pared a1desactivarse 10s compuestos intermedios; esta es la causa de que las mezclas H2+02 noexploten a unos 800 K si la presi6n es baja (


404 I. Martinez: TERMODINAMICA BASlCA Y APLICADAque se puede aplicar la ecuaci6n de van't Hoff, que en terminos de la constante de equilibrioK(T) toma la forma (9.21): alnK/a(l/T) = -h,lR, por lo que si la energia de activaci6n Ea secomporta como h, y varia poco con la temperatura, se podrci poner:a ln kl - -Ea, ---- + kl = B, exp- ~ ~ a ,1y si ambos procesos de activaci6n estcin en equilibrio y se les puede aplicar la ley de acci6nde masas (16.35)6, = 6, con 2~ = klcAcM 52 = k2cBcccM kl(T) - KC(T) (16.52)v ' v k2 (TIsiendo KC la constante de equilibrio expresada en concentraciones, seg6n se deduce de (9.20).energia1 AT coordenada de reaccidnFig. 16.3. Esquema de las energias de activacidn necesarias para la formaci6n de compuestosintermedios en una reaccidn del tipo A + B+C.Los valores de Ba y Ea se determinan experimentalmente mediante ajuste de medidas deconcentraci6n, como en (16.40) y (16.45), aunque las incertidumbres tipicas son de un ordende magnitud y se prefiere medir en una sola direcci6n y calcular la otra a partir de (16.52),incluso muy lejos del equilibrio.La medida experimental de la composicidn de productos intermedios en las reacciones decombusti6n es muy complicada porque s610 existen en zonas de espesor =lo4 m (aunquepuede aumentarse uno o dos drdenes de magnitud si se reduce mucho la presi6n) y sustiempos de vida son de 10-3 s a 10-lo s. Para ello se utiliza alguno de 10s mktodos tipicos deanalisis de reacciones rcipidas:- Llenado ripido. Se introduce la mezcla reactiva a Tamb en un dep6sito precalentado y seespera a ver si arde o no. El resultado depende mucho del estado superficial de la paredinterior (que van'a con el tiempo), y s610 es interesante para reacciones con tiempos deretardo grandes, y a6n asi hay que descontar el tiempo de la transmisi6n del calor.


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 405- Compresidn brusca. Se dispone la mezcla en un dispositivo cilindro-Cmbolo y secomprime, midiendo p, T y ni, y observando la aparici6n de llama (a veces aparecenllamas frias durante la etapa inicial). A1 ser la solicitaci6n creciente, aparecen efectoscruzados del tiempo de retardo y la temperatura de autoignici6n.- Tubo de choque. La reacci6n se provoca con una onda de choque que viaja a unos lo3rn/s, por lo que con la medida sincr6nica de detectores separados =1e3 m se puedeapreciar variaciones de composici6n correspondientes a s de intervalo.- Fot6lisis. S610 vale para reacciones que puedan iniciarse por destello. Se irradia lamuestra con el destello de un flash que a la vez dispara con un cierto retardo (=lo-6 s) unhaz de laser para el anilisis espectrosc6pico.- Ensanchamiento de las bandas de radiacidn producidas por NMR (resonancia magnCticanuclear) o ESR (resonancia de espin electr6nico).Combustor en reposoDentro de 10s procesos de combusti6n controlados por la cinCtica quimica, se handesarrollado tres modelos simples (correspondientes a otros tantos mCtodos experimentales):- Combustor en reposo: masa de control con campo de velocidades identicamente nulo,volumen fijo, mezcla reactante homogCnea en rCgimen transitorio.- Combustor bien agitado: volumen de control fijo y de masa fija (sale la misma queentra), movimiento turbulento para uniformizar el estado interior, mezcla reactantehomoginea en regimen transitorio o estacionario.- Combustor unidimensional: volumen de control fijo infinitamente largo, movimientoturbulento para uniformizar el estado interior, mezcla reactante homogknea en cadaseccibn, per0 evolucionando en rkgimen estacionario longitudinalmente.Es ficil ver que el ijltimo modelo coincidiri con el primer0 sin mis que cambiar las variablesespacial y temporal.El primer modelo que se propuso para analizar el combustor en reposo es el modelo deSemenov (1928), en el cual se considera una mezcla reactiva de gases ideales y de capacidadtCrmica constante (no dependiente de 10s cambios de temperatura y de composici6n), y que laaceleraci6n tCrmica es tan grande que el sistema alcanza altas temperaturas cuando todavia lacomposici6n y la presi6n son pricticamente las iniciales, por lo que la estequiometriaparticular de la reacci6n es irrelevante en esta fase. De todos modos, se supondrd que lareacci6n es F+V~O+V~N+V~P, donde las especies F, 0, N y P son el combustible, eloxidante, las especies inertes, y 10s productos, respectivamente; estos 6ltimos se tratan comouna sola especie. Como no se consideran variaciones apreciables de la composici6n, estemodelo se dice que es de tip0 termico.


406 I. Martinez: TERMODINAMICA BASICA Y APLICADACon este modelo, el consumo de combustible, seglin la ley de Arrhenius, y en funci6n de lacantidad de combustible nF, o de su concentraci6n CF, o de su fracci6n misica y ~, ser6 engeneral :6- wi a b c -dc, - -dco - -dcN - dcp ----=c Fc Oc N B a expA=--- -- ----- -m ~YFV viMi T dt vodt vNdt vpdt MFV dtdonde B,, a, b, c, Ta y B (definida por (16.53)) son constantes conocidas, m es la masaencerrada y Mi la masa molar de la especie i.El balance energttico seri en general dU I dt = Q = -U*A(T -Tab), siendo U la energiaintema y U* la transmitancia tCrmica por unidad de Area A. La energia intema tendrii unacomponente quimica (por avance de la reacci6n) y otra tQrnica (por el calentamientoresultante) tal que se podri poner dUldt=-murOdyFldt+mc,dTldt, siendo ur=urO la energiaintema estindar de reaccibn, que coincide con la entalpia estindar de reacci6n hrO si no hayincrement0 molar en la reacci6n (Cvi=O, aunque aqui da igual porque se desprecian 10sefectos de concentraci6n), y por tanto:CB ~ Y F d~dt dt-mb - + mc, - = -U*A(T - T ~ )Si se introducen ahora las variables adimensionales:se llega a1 sistema:do -0ad z 0- = a exp- - P(0 - 1) = Q~~~~- otkrmjSi se supone que el combustor es perfectamente adiabiitico (PO) la evoluci6n tCrmica sericomo la mostrada en la Fig. 16.4a7 apareciendo un tiempo caracteristico de inflamaci611, z,,,que se puede calcular con la aproximaci6n de Frank-Kamenetski?, como se resume acontinuaci6n. Se trata de desarrollar el tkrmino exp(0J0) en torno a1 punto el, lo cual sehace asi:con lo que la ecuaci6n (16.56) se integra fricilmente para obtener el tiempo critic0 zcr:2. Frank-Kamenetskii, D.A., "Diffusion and heat uansfcr in chemical kinetics", Plenum Press, 1969.


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 407o 'cr 7 1 8 oFig. 16.4. a) Evoluci6n 'explosiva' dc un sistema reactante aislado de gran energia de activaci6n. b) Calorliberado quimicament PqUh para varias presiones (masas) en el reactor y calor evacuado a travh delas paredes Q , en ~ func~on ~ ~ de la temperatura interior, 8. c) Evoluci6n de un sistema no aisladottrmicarnente (puede ser explosiva o no, seglin sea la presi6n).d8 -8, x a - dx I 8, -- = a exp- ) - - a e-eaex +d z 8 d z-+ jF-I/dr dx - a8, -+ zcr =- e 'a'a (16.58)En el caso normal en que exista transmisi6n de calor a1 ambiente, representando grificamenteel segundo miembro del balance energktico en (16.56), como se ha hecho en la Fig. 16.4b,puede apreciarse que el sistema dejaria de acelerarse tkrmicamente cuando la curva deenergia quimica liberada $ui,,, cortase a la de evacuacidn de calor QtLrrn, por lo que sdlo seproduciria la combusti611 si la presi6n es superior a un valor critic0 pcr, a1 cual corresponderiauna Tcr (punto de tangencia en la Fig. 16.4b). En la Fig. 16.4~ se han representado tresposibles evoluciones: la critica, una subcritica y otra supercritica. La relaci6n entre pCr y Tc-se obtiene resolviendo el par de ecuaciones Q ~ y dQqUh ~ I dB = dQlkr,,, ~ Id8 ~ que =delimitan el punto de tangencia en la Fig. 16.4b, obtenikndose como solucidn aproximadapara energias de activacidn grandes, 8, >> 1:OCr E 1 + - In Pamb - -a eaPcr1'a + creea 8,r(1+ 2 )=a+b+c %ra+b+cque se ha representado en la Fig. 16.5a. A la temperatura correspondiente a OCr se le llama deFrank-Kamenetskii.Sin embargo, este modelo tkrmico s610 predice un limite de presi6n para cada temperatura,mientras que en la prictica hay zonas en la que aparecen tres limites, con una especie depeninsula inestable como la mostrada en la Fig. 16.5b. Esta peninsula puede explicarse si seaiiade un esquema cinktico-quimico de embalamiento por reaccidn en cadena, como en(16.49). En efecto, Sean F, P, R y N el combustible, el producto, un radical y una especieinerte, respectivamente, y considkrese el siguiente mecanismo:


408 I. Martinez: TERMODINAMICA BASICA Y APLlCA DA- PcrPamb- PcrparnbcombustionFig. 16.5. Inestabilidad de un sistema reactante de alta energia de activaci6n (combusti6n): a) modelotkrmico, b) modelo completo termico y quimico.R k3 >P (terminaci6n en la pared)R+N+N k4 ,P (terminaci6n en la fase gaseosa)donde a mide el grado de ramificacibn. En el caso estacionario ser8:y para que la concentraci6n de radical sea finita y positiva, habri de ser:Pero, para una cierta temperatura, el 6ltimo tCrmino de (16.61) vm'a con la presi6n de unaforma parab6lica (crece indefinidamente tanto para presiones muy bajas como para presionesmuy altas), dando lugar a la especie de peninsula representada en la Fig. 16.5b.Combustor bien agitadoSe va a analizar aqui la evoluci6n de un sistema reactante en un recipiente de V=cte con unaennada por donde se aportan 10s reactivos en condiciones conocidas y un orificio de salida.La hip6tesis crucial va a ser la de suponer el recipiente (combustor, tambiCn llamado reactor)bien agitado y por tanto sin gradientes termodinimicos apreciables. Puede haber transmisi6nde calor, per0 el aporte de trabajo (agitador) se supone despreciable. La salida de productostendr6 lugar en las condiciones interiores.


Cap. 16: PROCESOS LIE COMBUSTION: MODELOS 409Las ecuaciones de balance de especies para el volumen de control (V=cte) en regimen noestacionario, que se deducen de mod0 anilogo a (5.35) con (16.3) en lugar de (5.1), seri:(suponiendo que haya R reacciones de interes).La del balance energetico, anilogamente a (5.37), seri:Se van a introducir ahora las hipdtesis siguientes: una sola reaccibn, de tipo Arrhenius, y condos 6nicas especies 1 y 2 (p.e. H2+H2=H2+H+H, o bien 0=2H-H2), siendo ambas gasesperfectos con propiedades sirnilares (M y cp, y por tanto coincidentes las funciones misicas ylas molares) y estando el reactor en flujo misico estacionario: riz, =me. Con estassimplificaciones, y llamando m a la masa (constante) de gases dentro del reactor e y a lafracci6n misica de la especie reactiva (p.e. H2), el balance de especies (16.62) pasa a ser:siendo para el ejemplo propuesto v= v"- 9=- 1 y v1=2. La concentraci6n de la otra especie (H)seri 1 -y.El balance energetico (5.37) sera como el desarrollado en (16.54), esto es,dUldt=mduldt=-mhr@dy/dt+mcvd7'/dt, y por otra parte he-h=-(ye-y)hr@+cp(Te-T); es decir:Ahora se introducen las variables adimensionales:donde 0, (temperatura de activaci6n adimensional) y a (poder calorifico adimensional) serinn6meros grandes3, 2 es la relaci6n entre el calor intercambiado con el exterior y la energiatCrmica sensible del flujo de aporte, siendo U* el coeficiente global de transmisi6n de calorQ = -U*A(T - Tomb), puede ser grande o pequefio, y Da es el llamado nlimero de Damkohler(1936), que mide la relaci6n entre el tiempo de residencia y el tiempo de reacci6n y puede ser3. A1 producto pea se le llama nlirnero de Zel'dovich, Ze, que en variables dimensionales es:


I. Martinez: TERMODINAMICA BASICA Y APLICADAgrande o pequeiio4.Con esta nomenclatura, y recordando que para el combustible es v=-1,(16.64-65) toman la forma:las ecuacionessiendo y la relaci6n de capacidades timicas (siempre de orden unidad). Ahora se va aanalizar el comportamiento del reactor en rCgimen estacionario y luego se van a ver algunosefectos no estacionarios.Caso estacionari~El problema (16.67) en el caso homogkneo (sin las derivadas temporales) admite unainterpretaci6n analitica sencilla que se ilustra en la Fig. 16.6. Se toma el niimero deDamkohler Da como variable independiente proporcional a1 tamaiio del reactor, y seobtienen 10s limites siguientes. Cuando Da4, entonces y+ye y 8+(@e+A)/(pA)~eequil,tirm, mientras que cuando Da+-, entonces y 4 y B+aye/(pA)=eequil,quh.TambiCn se puede determinar analiticamente el rigimen de inflamaci6n y el de extinci6n queaparecen en la Fig. 16.6a, mediante desarrollo en serie de la forma 8=8equil,tkrm+C6)rPan ey=ye-Cyaan para Da


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 41 1Es importante darse cuenta que el fen6meno de inflamaci6n viene gobernado por la relaci6nentre el calor generado en la reacci6n y el calor evacuado por las paredes, debido a lavariaci6n exponencial de la velocidad de reacci6n con la temperatura (ley de Arrhenius), y nocomo se creia antiguamente debido a que la velocidad de reacci6n tenia un salto brusco a latemperatura de inflamacibn.Si se considera ahora la evoluci6n temporal del sistema reactante dentro del reactor, y parasimplificar se supone me = m, = 0 y en las ecuaciones (16.64) y (16.65) se introducen lasvariables adimensionales:el sistema (16.67) pasa a ser:- dy - vV -0,d z --y exp- 8vl-- de - -ead T ay expg-P(e-l)que coincide practicamente con el (16.56) y que, a1 igual que aquC1, puede estudiarse por elmetodo general de 10s sistemas aut6nomos (en 10s que las variaciones no dependenexplicitamen te del tiempo).El comportamiento puede ser subcritico o supercritico, como se ilustra en la Fig. 16.7, seg6nque la transmisi6n de calor supere o no la deposici6n de energia quimica (puede verse que elcaso critic0 corresponde a una refrigeraci6n pale,).tiempotiernpoa)tiempob)tiempoFig. 16.7. Evolucidn de un sistema reactante encerrado en un reactor homogCneo:temperatura adimensional, 8, y fraccidn misica de reactive, y, en funci6n deltiempo. Esta evolucidn temporal es igual a la evolucidn espacial en un reactor deflujo estacionario. a) caso subcntico (mucha transmisidn de calor por lasparedes), b) caso supercritico (transmisicin de calor insuficiente).


41 2 I. Martinez: TERMODINAMICA BASICA Y APLICADAEl caso de un reactor homogCneo a1 que se comunica trabajo, como p.e. en un dispositivocilindro-kmbolo de 10s usados para medir la temperatura de autoignici6n de mezclasreactivas, se estudiaria de una manera similar.Llamas de difusi6nSe consideran aqui las llamas que se establecen en las capas de mezcla de dos masasgaseosas de composici6n diferente, una contendri el combustible y tal vez gases diluyentes eincluso algo de comburente, mientras que la otra contendrsi el comburente, normalmentediluido con gases inertes como en el caso del aire. La masa gaseosa que contiene elcombustible puede provenir de la vaporizacidn in situ de combustibles condensados (p.e.combustidn de gotas, velas, candiles y fuegos naturales).Las llamas de difusi6n estin controladas por 10s fen6menos de transporte, pues mis bien setrata de dos procesos sucesivos, uno de mezcla no reactiva en proporci6n casiestequiomCtrica, seguido de un proceso de combusti6n. Estas llamas adolecen de lascaracteristicas de las ondas de deflagraci611, es decir, no tienen una velocidad ni un espesordefinidos. La combusti6n es menos eficiente (no se puede quemar deprisa) y da lugar aproductos inquemados (s610 parcialmente pirolizados) en 10s gases de escape, formsindosedep6sitos de hollin si la llama toca una superficie s6lida. Este deficiente quemado da lugar allamas muy luminosas (radian mucho mis que las de premezcla) por la radiaci6n tCrmicaamarillenta de las particulas s6lidas de hollin a alta temperatura.Como la iinica caracten'stica destacable de las llamas de difusi6n es su geometria, se van aestudiar tres casos tipicos: uno bidimensional (plano), otro axilsimCtrico (chorro) y otroesferico (gota).Llama bidimensionalConsidCrese un mechero bidimensional consistente en una ranura de separaci6n L con unaplaca adicional a una distancia intermedia aL (Fig. 16.8a), que acaba a una cierta altura z=0,y que separa dos corrientes unidimensionales de la misma velocidad, una de combustible Fposiblemente diluido con un gas inerte N (en Ocxx>L).Sean y ~, e yo, las fracciones misicas conocidas de combustible y oxidante en las corrientesincidentes (zIO) y sea yF(x,z) la fracci6n de combustible para z>0, la cual se desea calcularcon las hip6tesis introducidas por Burke y Schumann en 192g5:- El tiempo caracten'stico de mezclado por difusi6n es mucho mayor que el tiempocaracten'stico de reaccidn quimica (n6mero de Damkohler Da+=), lo que implicaque la llama es de espesor infinitesimal y que yFyo=O en todo el campo fluido, y enparticular en la llama, donde yF=yo=O.5. Burke, S.P. y Schumann, T.E.W., "Diffusion flames", Ind. Eng. Chcm. 20, p. 998-1004, 1928.


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 41 3Fig. 16.8. Llama de difusi6n (linea gruesa) en un mechero bidimensional: a) mezcla rica, b)mezcla estequiomCtrica, c) mezcla pobre. El par6metro b esd definido en (16.71) yes b=a-ll(l+~F,lyoJ.- La reacci6n es de una sola etapa irreversible y sin increment0 molar, F+tO=(l+v)P,por lo que las fracciones misicas coincidirin con las molares.- No se consideran variaciones de velocidad media (pese a que el calentamientoacelerari mucho la corriente) y se supondri que s610 hay transporte por difusi6n concoeficiente D; constante e igual para todas las especies (pese a que aumenta con latemperatura). En realidad, esta hip6tesis es equivalente a suponer nula la entalpia decombusti6n.La ecuaci6n diferencial serri la de un proceso estrictamente difusivo en ejes m6viles avelocidad vertical v, dyF / at = ~~d~~~ I ax2, o bien la de un proceso de difusi6n convectivaesta ionaria en ejes fijos, que serin 10s que aqui se adopten, DyF 1 Dt =vdyF 1 az =2 2Did yF / ax , que junto con las condiciones de contorno completan el problema:siendo xF(z) la posici6n de la llama, en la cual ha de ser yF=yo=O. La 6ltima condici6n de(16.70) expresa la igualdad de flujos estequiomCtricos hacia la llama por ambos lados. Lasoluci6n general para el campo de concentraciones del combustible y~ y del oxidante yo es6:con yF.yo =0, b = a--Yo,"YoYYo = YF,6. Strehlow, R.A., "Combustion fundamentals", McGraw-Hill, 1984.


41 4 I. Martinez: TERMODINAMICA BASlCA Y A PLICADA(siendo a el factor geomttrico especificado en la Fig. 16.8), de donde es fhcil obtener el perfilde la llama. En la Fig. 16.8 se han dibujado 10s tres casos que pueden aparecer en funci6n delparimetro b, que como se ve en (16.71) depende de la estequiometria, las proporcionesrelativas de combustible y oxidante, y la geometria.Llama axilsim4tricaLa caracten'stica mis relevante de las llamas de difusi6n axilsimdtricas es su longitud L, lacual se puede estimar para un chorro laminar de combustible gaseoso en aire (Fig. 16.9)calculando la distancia a la que la composici6n en el eje baja por difusi6n desde la unidadhasta la estequiomdtrica; es decir, como el recorrido de la corriente a velocidad V, durante untiempo t=(Dd2)2/Di, siendo Do el digmetro y Di el coeficiente de difusibn, por lo que lalongitud seri Llw~VJlo2/(4Di).Fig. 16.9. Estructura interna de una llama de difusi6n laminar axilsim6trica.Con ayuda de las ecuaciones del chorro de difusi6n obtenidas en el Cap. 12 se puede haceruna mejor estimaci6n de la longitud de la llama.En efecto, la disminuci6n de la fracci6n mdsica de combustible, y, a lo largo del eje es, seg6n(12.52):siendo V=ZD,~VJ~ el gasto volumttrico inyectado. Tambitn conviene recordar que para lamayoria de 10s gases el coeficiente de difusi6n Di es aproximadamente igual a la viscosidadcinemitica v (no confundir con 10s coeficientes de estequiometria). De (16.72) se desprendeque para un combustible dado la altura de la llama de difusi6n es proportional a1 gastovolumttrico, independientemente del digmetro o la velocidad del chorro.Se desea aqui ahora determinar el estado (composici6n y temperatura) en un punto genkricodel campo fluido. Para ello, sup6ngase que la reacci6n global en una llama de difusi6n sepuede modelizar con cuatro especies: el combustible F, el oxidante 0 (que normalmente serioxigeno, 02), 10s productos P (en realidad esta especie serd casi siempre una mezcla de C02y H20 de composicidn casi constante) y un gas inerte N (que normalmente sera N2).Sup6ngase que la reacci6n estequiomktrica es F+vm02=vpP siendo Vm el coeficiente


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 41 5estequiomCtrico del oxidante en fracciones misicas (vp no interviene en el razonamiento);p.e. en CH4+202=C02+2H20, vm1,,=2 y vdsic0=4 para el oxigeno.Pero antes de entrar en detalle en las ecuaciones del chorro laminar reactivo, se va ademostrar que basta una s61a variable f para determinar en cada punto el estado de la mezcla(T, P, ~i).ConsidCrese un punto del campo fluido, de composici6n yi (fracciones misicas). Aunque enrealidad se ha alcanzado esa composici6n en ese punto por difusi6n de especies en presenciade la reacci6n quimica, se puede pensar que se consigue por mezcla y reacci6n adiabitica ycompleta en un combustor como el de la Fig. 16.10, en el que se hace entrar una corriente 1de gasto misico f conteniendo combustible puro, otra corriente 2 de gasto misico 1-fconteniendo el oxidante diluido con el gas inerte, y saliendo una comente 3 de gasto misicounitario conteniendo productos de combustidn y el exceso de sustancias que no hanreaccionado.Fig. 16.10. Esquema de flujos de entrada y salida a un combuslor adiabatic0 en el que tienelugar una combusti611 completa.La suposici6n de que siempre hay reacci6n y que Csta es completa para cualquier relaci6n degastos implica que el reactor no puede ser homogCneo, sin0 que localmente se conseguirinlas condiciones de mezcla inflamable y luego se diluirin 10s productos con el resto de lacomente. Esta es una gran simplificaci6n del problema quimico, pues se considera p.e. que lacomposicidn de la especie 'productos' no varia con la temperatura, ni hay pir6lisis del excesode combustible, etc.En resumen, las condiciones impuestas a las concentraciones son:YF +YO +YP +YN =1(en cualquier punto)YI;, =I (en el conduct0 1)Yo, + YN, = 1 (en el conducto 2)Y F, YO, = 0 (en el conducto 3) Iy se supone que la entrada por 1 y 2 es a temperatura ambiente To.El balance de cualquier magnitud conservativa y aditiva, @,seri:y es ficil ver que son conservativas y ~, ypyoIv,,, y cpT-yFhr (entre otras), donde se hatomado el modelo de gas perfecto para la entalpia tCrmica y h, representa la entalpia de


41 6 I. Martinez: TERMODINAMICA BASICA Y APLICADAreaccidn por unidad de masa de combustible. Se supondri ademis que todas las especiestienen aproximadamente el rnismo cp.N6tese que aunque tanto y~ como yo irin disminuyendo a1 reaccionar, la diferencia ypyolv,se mantiene constante, debido a la relaci6n de estequiometria de la reacci6n.A partir de 10s valores de entrada y salida se puede despejar la fracci6n de gasto necesariapara que se den esas condiciones:y en funci6n de ella representar las variables de salida, como se muestra en la Fig. 16.11. Enparticular, puede comprobarse que:Fig. 16.1 1. Composici6n y temperatura dc salida dcl reactor dc la Fig. 16.9, en funci6n de lafracci6n de gasto f.En conclusi6n, la variable f, que mide la fraccidn en que se han mezclado la corriente decombustible y la corriente de oxidante en el reactor ideal, basta para determinar el estado dela mezcla en cada punto. Para calcular su distribuci6n espacial concreta, se procede de lamanera siguiente.Las ecuaciones que gobiernan la evoluci6n real de la llama de difusi6n estacionaria,axilsimktrica y alargada, (12.45-48), serhn:- balance de masa:


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 41 7- balance de especies (se detalla para el combustible F, y es similar para las otras):- balance de cantidad de movimiento axial:- balance de energia:siendo DF la difusividad mrisica de la especie F en la mezcla, que se suponeaproximadamente igual a la de las otras especies.Es fricil ver que, si ademas se introduce la simplificaci6n de considerar pcte yj.Up=Di=k/(pcp), es decir, la igualdad de 10s nlimeros adimensionales Pr=Sc=Le=l, se verificapara las tres variables~yr;-~~lv~, w y h=cpT-yFhr la misma ecuaci6n diferencial, como en(12.45--48), dando lugar a una soluci6n de semejanza similar a la de (12.49) con la mismavariable:que, en funci6n del n~mero de Reynolds Re=waolv, resulta ser:estando limitada su validez a zlDo>(3/32)Re, como cuando no habia llama.Este resultado indica que, con las simplificaciones asumidas, la estructura de una llama dedifusi6n laminar axilsimktrica es iddntica a la de un simple chorro de difusi6n laminar noreactivo, con el mismo mapa de velocidades y una correspondencia directa entre ypyolv,, yh=cpT-yFh, para la llama por una parte, y 10s campos de y y T del chorro no reactivo. Es fricilcomprender que esto no sera verdad en la practica por las siguientes discrepancias delmodelo:- Las grandes variaciones de temperatura ocasionan grandes variaciones de densidad, yse ha supuesto p=cte.- La dilataci6n tkrmica da origen a fuerzas de flotabilidad que tienden a estrechar lallama (lo mismo que se estrecha el chorro no reactivo si su densidad es inferior a la


41 8 I. Martinez: TERMODINAMICA BASICA Y APLICADAatmosfCrica). Sin embargo, es sorprendente el buen acuerdo entre la teoria y laprictica en cuanto a la longitud de la llama; en realidad, estudios rnis avanzados quetienen en cuenta la flotabilidad predicen la misma longitud de llama para el chorroaxilsimCmco (no asi para el caso bidimensional plano).- Si la combusti6n genera m8s cantidad de sustancia de la que consume (si Cvi>O),apareceri un flujo expansivo adicional a1 efecto tCrmico.- En realidad no seri /~/p=D~=k/(p~~) como se ha supuesto. Si la difusividad tCnnica esmenor que la misica (Le>l) la temperatura de la llama seri superior a la temperaturade combusti6n adiabitica, e inferior en el caso contrario.- Encubierto en el modelo de reactor de la Fig. 16.1 1 se ha supuesto que la cinCticaquimica era mucho rnis ripida que la difusi6n y por tanto que la reacci6n tenia lugaren una superficie de espesor nulo (Da+O).Con este modelo, la llama se sitlia donde se juntan y desaparecen el combustible y eloxidante, es decir, donde f=fesleq dada por (16.82), correspondikndole la superficie (ver(12.55)):siendo yooo la fracci6n misica de oxidante en la atm6sfera no perturbada. De (16.83) seobtiene una longitud de llama, Lum=z cuando rflam=O:(no confundir la v de la viscosidad con la v, del coeficiente de estequiometria misico) quepuede compararse con (16.72) y que enseiia una vez rnis que la altura es directamenteproporcional a1 gasto volumktrico y no a1 dirimetro o la velocidad por separado, per0 ahoratambiCn se ve que la longitud es menor si se utiliza oxigeno puro que si se usa aire. La llamapresenta en este modelo su mrixima anchura a una distancia (9116)L y vale:Podria parecer que la flotabilidad daria lugar a llamas rnis largas, per0 este efecto secompensa porque el estrechamiento debido a la flotabilidad aumenta 10s gradientes radiales,que aceleran el mezclado.Los valores tipicos de Re en llamas de difusi6n en rkgimen laminar estin entre 10 y 100, yaque si Re>200se genera turbulencia).


Cap. 16: PROCESOS DE COMBUSTION: MODELOSLlama esfhricaLa combusti611 de gotas es de gran importancia priictica (sirve para estudiar la combusti6n deliquidos pulverizados) y de enorme interks cientifico porque la geometria es muy sencilla (sise supone simetria esfkrica) y permite un andisis detallado. Desgraciadamente, en laexperimentaci6n no se puede conseguir la simetria esfkrica por causa de la gran flotabilidaddebida a 10s gradientes tkrmicos, la dilataci6n y el campo de fuerzas gravitatorio. Losmdtodos tradicionales de quemar gotas sueltas son tres. El miis sencillo es situar una gota enla punta de una fibra de cuarzo y encenderla con dos chispas simktricas para evitar fuerzasdescompensadas. Otro procedimiento es dejar caer una gota (y todo el equipo auxiliar) en untubo al vacio (para eliminar la resistencia aerodiniimica). El tercer mdtodo consiste ensuministrar combustible liquido por un tubito a una esfera de material cedmico poroso,ajustando el gasto para alcanzar un estado estacionario. En todos 10s casos es muy dificilaproximarse a la simetria esfdrica; la experimentaci6n en ingravidez parece ser la posibilidadmiis prometedora.Se supondrii que se trata de una gota de combustible en una atm6sfera oxidante (aunquetambikn podria tratarse de la llama entre una gota de oxigeno liquido y una atm6sfera demetano, p.e.). Si las gotas son coloidales m) se producirii una mezcla homogknea almenos en la zona de precalentamiento de la llama, que pasaria a ser del tipo de premezcla.Las gotas con tamaiios entre 10-5 m y 10-2 m arden con una llama de difusi6n mantenida porla vaporizaci6n que su calor de combusti611 genera. En las gotas mayores ya son muyimportantes 10s efectos de propagaci6n lateral de la llama y aparici6n de bolsas fluctuantes degases (como en la combusticin en cubeta). No se considera aqui la combusti611 de particulass6lidas esfkricas, que es muy compleja porque en unos casos hay combusti6n heterogknea(como en el carMn) y en otros casos 10s productos son particulas s6lidas (lo normal a1 ardermetales), ni la combustidn esfkrica de una mezcla gaseosa (p.e. con la esfera ceriimica antesdescrita), ni la combusti6n de liquidos mantenidos en estructuras s6lidas (por fuerzascapilares) como el pabilo de las velas o la mecha de un quinquet o un candil, que en ausenciade flotabilidad tambikn se aproximaria a la geometria esfkrica.El caso estrictamente esfkrico es sencillo y fue resuelto por Spalding (1953)7. Si se suponeque s610 intervienen cuatro especies (combustible F, oxidante 0, productos P y un gas inerteN), el campo de composici6n y temperatura esperados sera como el mostrado en la Fig.16.12.Se supone que ~610 hay una reaccibn, cuya estequiometria molar sera F+v002=vpP, que elproblema es casi-estacionario (el tamaiio de la gota varia lentamente en comparaci6n con lavelocidad de 10s gases si sus densidades son muy diferentes) y que las velocidades no sonmuy grandes (para poder considerar p=cte). Se suponen conocidos el tamaiio inicial de lagota rdt=O), las condiciones en el infinito, T, y yo, y se hace la hip6tesis de que la cinkticaquimica es mucho miis riipida que 10s fen6menos de transporte (como en el modelo deBurke-Shumann), por lo que la llama serii de espesor nulo y la temperatura en el radio de lallama serii Tad. Tambikn se consideran conocidas las propiedades de equilibrio y detransporte del fluido, y se usarii la aproximaci6n Pr=Sc=Le=l (compatible con la suposici6nde Tl&=Tad).7. D.B. Spalding, Fuel 32, p. 169,1953.


I. Martinez: TERMODINAMICA BASICA Y APLICADAFig. 16.12. Perfiles de cornposicidn y y tcrnperatura Ten la cornbustidn de una gota con sirnetriaesfkrica. La llama de difusidn se sitiia en rl. Se supone que s610 intervienen cuatroespecies (cornbustiblc F, oxidante 0, productos P y un gas inerte N).Las inc6gnitas que se desea despejar son la posici6n rl de la llama y las condiciones en lasuperficie de la gota: To, y ~ y , droldt. Para ello se va a establecer el balance de oxidante en unradio r>rl y el balance de energia en un radio rorl debe haber un flujo de oxigeno por difusidn, dado por 47~~pD~dy~ldr,el cual debe compensar el flujo de oxigeno que desaparece por reaccibn, mds el flujoconvectivo hacia el exterior en el caso de que se generen en la reacci6n mis gases de 10s quese consumen (vp>vO). Si mF=4zro2pfiqdroldt es el consumo instantineo de combustible, elconsumo de oxidante serd mF voMo 1 MF y el flujo convectivo de oxidante hacia afuera serdmF(vP-vo)yoMmlMF, siendo Mm la masa molar de la mezcla. Por tanto, el balance deoxigeno seri:que integrada entre r=rl y r+.o a t = cte da:Por otra parte, el balance de flujos de energia en un radio r,


Cap. 16: PROCESOS DE COMBUSTION: MODELOSSumando (16.87) y (16.89) para eliminar rl, y multiplicando por mFro queda:A1 sustituir valores tipicos de hidrocarburos en (16.90) puede comprobarse que el segundosumando es mucho menor que el primero, por lo que se puede poner:e integrando desde el instante inicial t=O:que enseiia que el irea disrninuye linealmente con el tiempo. Pese a las grandessimplificaciones introducidas, la ecuaci6n (16.92) se ajusta muy bien a 10s resultadosexperimentales (*15%). A la constante P en ro2=r20,i~cial-~t se le llama constante deevaporaci6n, pues mucho antes de 10s experimentos de combusti6n ya se sabia que laevaporaci6n de gotas variaba linealmente con su Area, (11.35). TambiCn es ficil cornprobarque el consumo de combustible es directamente proporcional a1 radio, ya que2 2m = -p,, 4nr0 dro / dt = -p . 2nrodro / dt = p . 2nrJ, lo que permite predecir p.e. el gasto11qmasico necesario para el quemador de esfera cerimica porosa descrito anteriormente.Sustituyendo (16.91) en (16.87) se despeja el radio en el que se sit6a la llama, rl:f~qdonde se ha tenido en cuenta la igualdad de difusividades tkrmica y misica (Le=l).Sustituyendo valores tipicos en (16.93) puede verse que rllro=lO+lOO, y la soluci6n (16.92)podria haberse obtenido aproximadamente de una forma muy sencilla, sin tener en cuenta lacombusti6n, analizando simplemente la evaporaci6n de una gota en una atm6sfera inerte detemperatura T,=Td.Finalmente, conviene seiialar que todo este analisis se ha basado en que las propiedades delfluido no varian con la temperatura, habikndose despreciado incluso las variaciones dedensidad a1 integrar (16.88), por lo que no sera valid0 cuando la presi6n ambiente seacercana o superior a la critica, con10 ocurre en 10s motores cohete de combustible liquid0 yen 10s motores Diesel.


I. Martinez: TERMODINAMICA BASICA Y APLICADALlama de premezclaSe estudia aqui la estructura interna de una llama de premezcla plana en flujo laminar, esdecir, el detalle fino de la onda de deflagraci6n que se autopropaga a travCs de una mezclacombustible en reposo. Se tratari de calcular la velocidad de propagaci6n, el espesor de laonda y 10s perfiles de temperatura y composici6n.Experimentalmente, se pueden conseguir llamas planas estabilizrindolas sobre una rejilla o untap6n poroso, per0 s610 si la velocidad de la llama es pequeiia (


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 423Suponiendo que el espesor de la zona de reacci6n (C) es mucho menor que el de la zona deprecalentamiento (B), Cste se podri estimar mediante un balance energetic0 sencillo, sindifusi6n de especies ni reacci6n quimica. Este modelo tCrmico fue propuesto por Mallard yLe Chatelier en 18838.En ejes ligados a la onda y en rCgimen estacionario, el balance energetic0 ser6:donde se ha supuesto que no hay gradientes espaciales de presi6n (en realidad son s610 del2%), es decir, que el fluido no se acelera (vkase (5.31) y (5.32)) y por la misma raz6n, que nohay efectos viscosos. Llamando Vq a la velocidad de quemado (definida como la depropagaci6n normal de la llama respecto a la mezcla fresca) y LT a1 espesor de la capatCrmica de la llama, el anrilisis de 6rdenes de magnitud de (16.94) indica que VqILFaIL3 ypor tanto que V&=a. La difusividad tCrmica a se puede aproximar por la del aire, que sueleser el componente mayoritario; en condiciones ambientes, a=kl(pcp)=0,025/(1,2. 1000)=2.1V5 m2/s, per0 k crece con T1o y p decrece con T, por lo que a=pt2 y para unatemperatura media de unos 1500 K se llega finalmente a que VqLF0,2.10-3 m2/s, lo cualconcuerda razonablemente con 10s valores tipicos de Vq=0,5 m/s y LF1V3 m. Como lavelocidad de quemado serri proporcional a1 espesor LT, y a la velocidad de reacci6n e,tambiin se concluye que Vq = (a()ll2.La Ec. (16.94) puede integrarse en la forma:habiendo tomado el origen de distancias (x=O) en el punto que corresponderia a la separaci6nentre las zonas B y C (Fig. 16.13), que seri prdxima a la temperatura de autoinflamacidn dela mezcla fresca, aunque en este caso la mezcla en x=O contendri radicales y, pese a que suconcentraci6n es pequefiisima y no influyen en las propiedades tkrmicas, sus energias deactivaci6n son muy diferentes de las de 10s gases frescos y su influencia en las propiedadestermoquimicas puede ser apreciable.La relaci6n V&F~ tambiCn puede servir para estudiar c6mo vm'an estos parimetros con lapresi6n. Por una parte, la relaci6n entre la velocidad de quemado y el espesor de la llama seriproporcional a la velocidad de reacci6n: VqILT"{, la cual serri proporcional a pv--' y por tantoa pg-1. Por otra parte a=k/(pcp)-p-I, luego de Vq&p-1 y Vq/L~pv'-I se deduce finalmenteque Vq=p(v'-2)/2 y LTp-vM y como las reacciones elementales mris comunes son las de v'=2,para este caso la velocidad de la llama no varia apenas con la presi6n y el espesor de la llamaes proporcional a llp (esto es, a1 recorrido libre medio de las molCculas). En realidad, para10s hidrocarburos, si la velocidad es Vq=0,5 rn/s, su variaci6n con la presi6n es pequeiia (p.e.8. Mallard, E. y Le Cha~elicr, H.L., Ann. Mines 4, p. 379, 1883.


424 I. Martinez: TERMODINAMICA BASICA Y APLICADApara las mezclas metanolaire, la relacidn de velocidades a 1 MPa y a 0,l es 0,3), siendopositiva para velocidades mayores y negativa para velocidades menores.El espesor Lr de la zona de reacci6n (C) estri gobernado por el balance entre la generaci6n deradicales (cinCtica quimica) y su evacuaci6n (transporte por difusibn). En efecto, de (16.1 I),(16.14) y (16.16) se deduce:La zona C (Fig. 16.13), cuyo espesor es =0,1 mm a presi6n atmosfkrica, es la zona dereaccidn y no estsi en equilibrio termodinsimico (pese a que una molCcula sufre unas 105colisiones en esa zona), como es fdcil comprender a1 considerar la emisi6n lurninosa, que espor quimiluminiscencia (emisi6n debida a la recombinaci6n de radicales, que no correspondea la ley de Planck); la emisi6n en las mezclas ricas es debida a1 C2 y es verdosa, mientras queen las mezclas pobres es debida a1 CH y es violsicea.Los radicales aparecen a1 principio de la zona C (a la temperatura de autoinflamaci6n); lasmolCculas se descomponen en la zona C1 dando radicales de alta entalpia que se recombinanen la zona C2 para formar 10s productos con gran liberaci6n tCrmica (10s productos son debaja entalpia). Finalmente, cierta cantidad de radicales salen de la zona C a la D (10s de lacomposici6n de equilibrio a esa temperatura, aunque a veces salen mds, sobre todo en llamasa baja temperatura, donde 10s tiempos de relajaci6n son mayores).La falta de equilibrio termodinrimico puede comprenderse a1 considerar 10s diferentestiempos de relajaci6n de 10s diversos procesos que entran en juego. Asi, el equilibrio de 10smodos de traslaci6n (distribucibn maxwelliana de velocidades) se consigue en apenas 10colisiones tras una perturbaci6n inicial, mientras que el equilibrio de 10s modos de rotaci6nmolecular necesita unas 103 colisiones, el de vibraci6n unas 105 colisiones, el de disociaci6n([A][B]I[AB]=f(T)) unas lo7 colisiones, el de ionizaci6n ([A+][e-]l[A]=f(T)) unas 109colisiones, etc.En la zona D ya se puede considerar que se ha alcanzado el nuevo equilibrio termodinsimico,y, si la llama es grande y se evitan otras fuentes luminosas, se puede observar una tenueemisi6n rojiza correspondiente a la de cuerpo negro a la temperatura de la llama en lasbandas espectrales del H20 y C02.Ya se ha razonado por quC la velocidad de deflagracidn laminar depende poco de la presibn.Con la composici6n, la velocidad es msixima en la estequiometria o para mezclas un pocomris ricas, disminuyendo a la mitad o asi hacia 10s limites de ignici6n (limites depropagaci6n). Si en lugar de nitr6geno se usase He (en la misma proporci6n que el aire), lavelocidad mrixima pasa de ser 0,4 d s a 1,3 4s. Con la temperatura, la velocidad aumenta(debido a1 precalentamiento) casi cuadrriticaniente (p.e. para el metano en aire, Vq=O,l m/s a100 K, 0,4 m/s a 300 K, 1,6 m/s a 700 K).


Cap. 16: PROCESOS DE COMBUSTION: MODELOSMedida de la velocidad de deflagracion laminarSe trata de medir la velocidad normal (perpendicularmente a1 frente) de la onda respecto a 10sgases frescos Vq. Para ello se han propuesto varios mCtodos (Fig. 16.14) que puedenagruparse en:- llamas fijas (ancladas en mecheros: llama plana, llama bidimensional y llamaaxilsimCtrica)- llamas mbviles (en tubo, en pompa de jabbn, en bomba a volumen constante).En el mechero de llama plana (atribuido a Powlingg), Fig. 16.14a, se consigue estabilizar unallama de premezcla sobre un disco poroso (con orificios de menos de 1 mm) refrigerado,aunque hay que enfriarlo muy bien para que se estabilicen llamas de mis de 0,15 4 s devolocidad. Se enciende con gastos mayores, formindose una llama cbnica lejos del mechero,y luego se va reduciendo el gasto hasta que se hace plana. Representando la velocidadmedida en funcibn del grado de refrigeracibn y extrapolando a cero se calcula la velocidadnormal adiabitica.drejillab)C)disco porosoinerte -bamezcla 4combustible Imczclacombustiblelama iiinerte-wmezclacombustibletinerte-mmezclacombustible tFig. 16.14. Esquema de procedimientos para la medida de la velocidad de propagacidn de llamaslaminares: a) llama plana, b) llama bidimensional, c) llama axilsimktrica, d) llama haciaabajo en tubo, e) llama cn pompa dc jabbn, f) llama cn bomba a volumen constante.9. Powling, J., Fuel 28, p. 25, 1949.


426 I. Martinez: TERMODlNAMlCA BASlCA Y APLICADAEn el mechero de llama bidimensional (desarrollado por Lewis y von Elbelo), Fig. 16.14b, seancla la llama a una ranura de =1 cm de anchura y unos 10 cm de largo. Mediante unailuminaci6n adecuada se puede visualizar una zona de comportamiento bidimensional lejosde 10s bordes, y si se controla el perfil de velocidades de salida (parab6lico o mejor lineal)puede calcularse la velocidad normal (a partir del gasto volumCtrico y el ingulo del 'cono' otriingulo que forma la llama) simplemente como Vq=Vsena, siendo V la velocidad de salida(si el perfil no es constante conviene usar trazadores y medir V en cada punto) y a elsemiingulo en el virtice. Si se desea medir en llamas ricas conviene separar la llama dedifusi6n de la de premezcla, como se ha hecho en la Fig. 16.14b, aunque esto complica lainstalaci6n.En el mechero de llama axilsimCtrica (el tipico mechero de Bunsen), 16.14c, se procede igualque en el de llama bidimensional, pues en realidad Cste es una variante de aquC1, en la cual lainstalaci6n es mis sencilla (no hay que tomar precauciones en cuanto a la homogeneidad a lolargo de la ranura, ni hay problemas en 10s bordes), per0 el anilisis es algo mis complicadoporque la llama tiene curvatura (en 10s dos casos anteriores no tenia) y porque siempre semira a travCs de la propia llama, que actlia de lente. El diimetro tipico del tub0 es de 1 cm, yla incertidumbre en la velocidad suele ser de +20%.Para la medida en tubo suele usarse la propagaci6n hacia abajo, Fig. 16.14d, ya que lapropagaci6n hacia arriba viene muy distorsionada por la flotabilidad de 10s gases quemados,y en la propagaci6n en tubo horizontal la llama no es axilsimCtrica (en todos 10s casos, elfrente es parab6lico debido a1 efecto refrigerante de las paredes, y en el tub0 horizontal elvCrtice de la parabola esti por encima del eje del tubo). El diimetro del tubo suele ser deunos 5 cm, e influye algo en la velocidad medida. Por supuesto que la ignici6n ha de hacerseen el extremo abierto si se quiere medir V=Vq en la referencia ligada a1 laboratorio (enrealidad Vq=VAtubolANarna, siendo V la velocidad aparente, Atubo el Area de la secci6n del tuby All,, el kea de la llama), ya que si se produjese la ignici6n en el extremo cerrado lavelocidad aparente seria V=VqTd/Tah, debido a la expansi6n de 10s gases calientes. Encualquier caso, conviene practicar un pequeiio orificio en el extremo cerrado para amortiguarlas ondas aciisticas de presi6n que alli se reflejan.Para la medida en burbuja de jabbn, Fig. 16.4e, se hace una pompa de =1 cm de diimetroinsuflando la mezcla combustible en el liquid0 jabonoso, se introducen dos electrodos en elcentro y se provoca la ignici6n. La velocidad de deflagraci6n normal Vq se puede obtener dela velocidad aparente V como Vq=VTamDITad=V(rinicial/rfinal)3, siendo CS~OS 10s radios de lapompa de jab6n a1 principio y a1 final (conviene que la atm6sfera exterior sea inerte para queno siga ardiendo).La medida en bomba a volumen constante, Fig. 16.14f, es similar a la anterior,simplificindose enormemente las condiciones experimentales a cambio de una complicaci6nen el anilisis, pues ahora la expansi6n de 10s gases quemados esti restringida y tiende acomprimir la corona de gases frescos (y Cstos a aquCllos). La presi6n se puede consideraruniforme en cada instante, y su variaci6n con el tiempo (fricilmente medible) permite yacalcular la velocidad de deflagraci6n normal por medio dell:10. Lewis, B. y von Elbe, G., "Combustion, flamcs and explosions of gases", Acadcmic Press, 1961.11. Kuo, K.K., "Principles of combustion", John Wilcy & Sons, 1986.


Cap. 16: PROCESOS DE COMBUSTION: MODELOSaunque a veces se prefiere medir la variacidn de la presidn y el radio de la llama a la vez, porredundancia. En (16.97) R es el radio de la bomba, yla relaci6n de capacidades tdrmicas de10s gases (p1,4), pl la presidn inicial (normalmente la ambiente) y p2 la presi6n final de lacombusti6n adiabitica a volumen constante.La primera ecuacidn de (16.97) se obtiene igualando la expansi6n de la esfera de gasesquemados, 4xr2drldt, con el avance de la llama respecto a 10s gases frescos parados, 4m2Vqmis el avance debido a1 retroceso de 10s gases frescos por la compresi6n isentr6pica(deducida de dplp+ylV~IV~=O).La segunda ecuaci6n de (16.97) se obtiene planteando las ecuaciones de estado de gas idealpara la esfera quemada y la corona fresca en cada instante y a1 principio y a1 final,relacionando las presiones con la fracci6n de combustible quemado a travds del balanceenergdtico, donde se supone que todos son gases perfectos de igual capacidad termica. En laFig. 16.15 pueden verse algunos resultados experimentales12.Fig. 16.15. Resultados de un ensayo de cornbustidn en un recipiente esfkrico de radio R=80 mm. a)Variacidn de la presidn interior con el tiernpo; b) perfil de ternperaturas interiores en dosinstantes; c) variacidn con el ticrnpo de la posicidn de ues particulas trazadoras situadasinicialrnenle a 10%, 50% y 90% dcl radio; d) variacidn con el tiernpo de la velocidad normaldc qucrnado Vq y dc 10s dos terrninos que se restan en (16.97) para obtenerla.12. Takerno, T. y Iijirna, T., "Nonsteady flame propagation in closed vessels", Prog. Astronaut. Aeronaut. 76,p. 578-595, 1981.


I. Martinez: TERMODINAMICA BASICA Y APLICADARECAPITULACION1. En el capitulo anterior se puede encontrar una descripci6n fenomenol6gica de lascaracteristicas de 10s procesos de combusti6n y de algunos casos prricticos como laestabilizaci6n de llamas. Este capitulo es mis bisico porque trata de modelizarte6ricamente algunas de esas caracteristicas.2. En primer lugar se presenta una revisi6n de 10s tipos de anrilisis te6ricos que se handesarrollado, centrando el problema en la modelizaci6n de la cinktica quimica, que sueleser la hip6tesis determinante del tipo de aniilisis.3. Se presentan de manera muy sucinta las ecuaciones generales de la evolucidn local de 10ssistemas reactantes, pasando a una simplificaci6n muy utilizada en el aniilisis: el modelode Shvab-Zel'dovich.4. Las ligaduras que la Termodinrimica del equilibrio impone a 10s procesos de combusti6nse vuelven a tratar con mhs detalle de lo que cabia en un capitulo general sobreTermodiniimica de sistemas reactantes (Cap. 9). Los tres aspectos que se consideran son:la composici6n en el estado de equilibrio, el poder calorifico, y la temperatura decombusti6n adiabstica.5. Para la modelizacidn de la evoluci6n de 10s procesos de combusti6n7 se recurre a1 estudiode la cinCtica quimica, que enseiia a determinar la velocidad de reacci6n 5 a partir de unconjunto de reacciones elementales (mecanismo de la reacci6n) y de la ley de Arrhenius,en la forma:6. Se estudian en detalle dos problemas de combusti6n dominados por la cinCtica quimica,seglin el modelo de Semenov: el combustor en reposo, y el combustor bien agitado.7. De 10s problemas dominados por el transporte de calor y masa y donde 10s efectos decinktica quimica son de segundo orden, se presenta el mCtodo de Burke-Shumann paraanalizar las llamas de difusi6n bidimensional, axilsimktrica y esfkrica.8. Finalmente se estudia la llama de premezcla (es decir, la onda de deflagraci6n laminar),donde la cinktica quimica queda acoplada a1 trasporte de calor y masa para definir unavelocidad de propagaci6n caracteristica.PROBLEMASLos problemas de combusti6n completa se han agrupado en el Cap. 15, proponikndose aquis61o problemas en 10s que hay que determinar la composici6n, asi como 10s problemas miste6ricos de combusti6n controlada por la cinCtica quimica, por difusi611, y por propagaci6n deonda.


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 42916.1. En un combustor se inyecta a raz6n de un rnol de CO y un rnol de 02, a temperaturaambiente, y se provoca la reacci6n isobirica. Se pide:a) Determinar la constante de equilibrio en funci6n de la temperatura.b) Calcular la composici6n y temperatura de salida y el calor transferido, suponiendoque se ha producido 0,l rnol de C02.c) Calcular la composici6n y temperatura de salida y el calor transferido, suponiendoque se ha producido 0,9 rnol de C02.d) Calcular la composici6n y temperatura de salida en condiciones adiabriticas.e) Calcular la temperatura adiabritica correspondiente a la combustidn completa.Sol.: a) lnK=104-114,4(1-298lT) con T en kelvin; b) 5% C02, 46% CO, 49% 02,T=3980 K, Q=244 kJ; C) 58% C02, 6,5% CO, 353% 02, T=2600 K, Q=-87 kJ; d)43% C02, 18% CO, 39% 02, T=2750 K; e) T=3920 K.16.2. Se hace pasar 1 p kg/s de vapor por un lecho de carb6n a presidn ambiente,producikndose CO y H2. Se pide:a) Composici6n de salida a 100 OC y a 1000 OC.b) Calor intercambiado en ambos casos.c) Consumo de carb6n en ambos casos.Sol.: a) a 100 OC xHZo=0,99998, xco=O,OOOO1, xH2=0,00001, y a 1000 OC ~~~~=0,001,xco=0,499, xH2=0,499; b) ~=-0,03 W a 100 OC y ~=980 W a 1000 OC; c)m=0,25.10-6 kg/s a 100 OC y m=8,3.10-3 kg/s a 1000 OC.16.3. En una ciimara cerrada y rigida hay inicialmente 1 rnol de CO y 2 rnol de O2 a 25 OCy 200 kPa., y se inicia la reacci6n. Se pide:a) Calcular la presi6n y temperatura de combusti611 adiabiitica y la composici6n deequilibrio.b) Sabiendo que a1 cab0 de cierto tiempo la temperatura es de 1200 K, calcular lapresibn, la composici6n de equilibrio y el calor intercambiado.Sol.: a) Tad = 2870 K, p = 1,63 MPa, XO, = 0,606, xco, = 0,367, xco = 0,028; b) p = 0,67MPa, xo, = 0,6, xco, = 0,4, xco = 3. Q = -25,5 kJ1mol.16.4. Calcular la concentraci6n de equilibrio de NO en el aire a 300 K, 1000 K y 3000 Ksuponiendo que s61o hay N2, 02 y NO y que las variaciones de entropia sondespreciables, juzgando la bondad de las hip6tesis.Sol.: A 300 K xNo=10-15, a 1000 K ~ ~~=10-~, y a 3000 K xNo=lO-1, per0 a 3000 Kapareceriin otros compuestos, como NO2.16.5. Dentro de un recipiente esfkrico de acero, de 0,5 m de diiimetro, sumergido en unbafio a 500 K, existe un gas, inicialmente a 150 kPa, cuya velocidad de reacci6n semodeliza por la ley de Arrhenius con las constantes Ta=lO" K y Ba=106 m3.mol-1.s-l.El coeficiente global de transmisi6n de calor a travks de la pared es 10 W.m-2.K-1 yla entalpia de reacci6n -lo5 J.mol-1. Se pide:a) Presi6n y temperatura del estado critico de combusti6n (mixima temperatura deoxidaci6n lenta).b) Tiempo critico adiabiitico y variaci6n de la composici6n desde el estado inicial a1cn'tico.


430 I. Martinez: TERMODINAMICA BASlCA Y APLICADAc) Variacidn de la velocidad de reacci6n desde el estado inicial a1 cdtico.d) Estimar el n6mero de Rayleigh para ver si es importante la conveccidn interior (apartir de Ra=600 ya empieza a serlo).e) Variacidn de la temperatura y la composici6n con el tiempo, e influencia de latransmitancia tirrnica.Sol.: a) BCr = 1 + l/Oa = 1,05 + Tcr = 525 K, p = nRTIV + pcr = 158 kPa; b) T,, =exp(Ba)/(aBa) = 3. lo6 + tcr = 110 S; C) 4, / (, = ~X~(-BJB~-)/~X~(-B~) = 2,6; d) Ra =agATD3/(va) = lo8; e) en la Fig. P-16.5 se muestra la evolucidn temporal de latemperatura adimensional, 8 = TIT,, y la funcidn mbsica de productos combustibles,y, para cuatro valores de la trnasmitancia tCrmica: U* = 5 W.m-2.K-l, U*/2, U*/4 yU" = 0.Fig. P-16.5 (U* en W.m-2.K-I).16.6. En una atm6sfera inerte a 100 kPa y 25 "C, se tiene inicialmente una burbuja de jabdnde 1 cm de radio llena de una mezcla combustible de metano y aire en proporcidnestequiomCtrica. En un cierto instante se provoca la reacci6n con una chispa central.Se pide:a) Esquema de la configuraci6n en tres instantes del proceso (inicial, intermedio y final),indicando la nomenclatura que se va a usar para la geometna y las variablestermodinbmicas.b) Explicar por quC es de esperar que el radio de la llama varie linealmente con eltiempo, y hacer un esquema conjunto de la variaci6n del radio de la llama y del radiode la burbuja en funci6n del tiempo, deduciendo una expresi6n analitica para esta6ltima.c) Calcular la temperatura de la llama suponiendo combusti6n completa.d) Calcular el radio final.e) Calcular la velocidad observable de la llama y el tiempo que tarda en arder la burbuja,suponiendo que la velocidad de deflagraci6n es de 0,45 m/s.Sol.: a) drB/dt=Vq(l+TB/TF)=cte; C) TQd=2280 K; d) rB,2=rF,2=1,97 cm; e) t2=5 ms.


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 43116.7. En un reactor a presi6n constante se quema metano con oxigeno puro. Se pide:a) Suponiendo combusti6n estequiomktrica y completa, determinar la composici6n a lasalida, el poder calorifico y la temperatura de salida adiabritica.b) Sabiendo que en realidad se quema a raz6n de 2 moles de metano por cada 3 moles deoxigeno, y suponiendo que todo el hidr6geno se oxida completamente (y el carbon0parcialmente), determinar la nueva estequiomenia, la composici6n a la salida, elpoder calorifico y la temperatura de salida adiabhtica.c) Sabiendo que en 10s productos aparecen C02, CO, H20 y H2, y suponiendo que esthen equilibrio, determinar el n6mero de reacciones independientes necesarias pararnodelizar el proceso real, e indicarlas.d) Determinar las constantes de equilibrio de las reacciones anteriores en funci6n de latemperatura, aproximrindolas por lnK=A-BIT.e) Plantear el sistema de ecuaciones que permiten determinar la composici6n, el podercalorifico y la temperatura de equilibrio adiabritica en el proceso 6ltim0, indicando lasinc6gnitas y 10s datos.Sol.: a) (vol) 33% C02+67% H20, PCS=891 kJ/mol, Tad=5700 K; b) (vol) 33% C02+67%H20, PCS=608 kJ/mol, Td=4300 K; c) R=C-E=6-3=3, CH4+202=C02+2H20,C0+(1/2)02=C02, CO+H20=C02+H2; d)1 n2'~0,'H202 =-0,57+---X ~ ~ 4 X ~ 296650Te) hay 8 ecuaciones: las tres anteriores, y las cinco siguientes:donde C es la cantidad total de sustancia producida:y donde, por sencillez de escritura, se han sustituido 10s valores de entalpias deformaci6n y capacidades tirrnicas en la ecuaci6n del balance energktico.


432 I. Martinez: TERMODINAMICA BASICA Y APLICADA16.8. A una ciimara adiabitica de combusti6n llega una corriente estacionaria de 0,5 kg/minde octano liquido (C8H18) a 25 OC y 800 Pa. Se quema con 100% de exceso de aire(21% de 02, 79% de N2), de modo que en 10s productos se detecta la presencia deC02, 02, N2, NO y H20. A continuaci611, 10s productos pasan por una turbinaadiabitica de la que salen a 950 K y 200 Pa. Calciilese, sabiendo que la presi6n ytemperatura del ambiente son, respectivamente, 100 P a y 25 OC:a) La temperatura y la composici6n de equilibrio de 10s gases a la salida de la ciimara.b) La potencia extraida de la turbina.c) La irreversibilidad del procesod) Capacidad mhima de trabajo del Cf118 liquido, si la composici6n del ambiente es:XCO, = 0,0003, XO, = 0,2068, XN, = 0,7802 y XH,O = 0,0127.Sol.: a) 1500 K, C02: 6,47%, H20: 7,29%, 02: 10,08%, N2: 76,07%, NO: 0,09%; b)167 kW; c) 598 MJ; d) 5430 MJ.kmo1-1.16.9. Determinar la composici6n de equilibrio y la temperatura alcanzada en la combusti6nestequiometrica, adiabatica e isob~ca de propano en aire y en oxigeno puro, en 10ssiguientes supuestos:a) Suponiendo que s610 aparece C02 y H20.b) Suponiendo que ademis aparece CO, Hz, 02 y NO.c) Suponiendo que ademis aparece H, 0 y OH.Sol.:Tabla P- 16.9En aireEn oxigenoa b c a b cXCO, 0,1111 0,1003 0,1004 0,4286 0,0896 0,1396XOH 0,0027 0,1005T [K] 2324 2232 2219 6340 3420 308116.10. De medidas de la velocidad de reacci6n para 2N02 + 2N0 + 0 2 se ha obtenido laconstante de reaccidn, resultando ser k = 0,498.10-3 m3.mol-1.s-1 a 592 K yk = 4,743.10-3 m3.mol-l.s-l a 656 K. Se pide:a) Calcular la velocidad de reacci6n para c~o, = 3 mo1.m-3.b) Calcular las constantes E,, Ta y B, de la ley de Arrhenius.c) Representar la composici6n en funci6n del tiempo, sabiendo que inicialmente mo, =3 mob3, CNO = co, = 0 y que la reacci6n se hace a T = 656 K.


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 433Sol.:a) dm02/dt = -0,009 m~l.s-l.rn-~ a 592 K y dc~04dt = -0,085 rnol.~-l.m-~ a 656 K; b)E, = 0,11 MJ/mol, T, = EJR = 14000 K, B, = 4.103 m3.mol-l.s-l; c) Fig. P-16.10.Fig. P-16.1016.11. Se pretende estudiar la evoluci6n temporal en el proceso de combusti611 esfericadescrito en el problema P-15.27 (R = 0,4 m, pl = 100 Ha, T1 = 293 K, xl = 0,033). Sepide:a) Indicar quC relaci6n habri entre la velocidad de deflagraci6n laminar, Vq, y lavelocidad aparente de la llama a1 comienzo y a1 final de la combustibn.b) Relacionar Vq con la velocidad aparente de la llama y el perfil temporal de presiones.c) Relacionar la presi6n con la posici6n del frente de llama.d) Variaci6n de la presi6n con el tiempo cuando t + 0.e) Despreciando el efecto de la temperatura en Vq y suponiendo Vq = 0,5 d s representarla evoluci6n en funci6n del tiempo y en funci6n de la presi6n.Sol.: a) A1 principio la relaci6n entre el volumen de gases quemados y el del recipiente estan pequeiia que puede suponerse combusti611 a p = cte = 100 kPa y, por tanto, lallama se moverii con la velocidad Vq respecto a 10s gases frescos y Vq(l+T~l/7'~l)respecto a 10s gases quemados, siendo TF1 = T1 = 293 K y habiendo obtenido TBl =2310 K en el P-15.27 (p.e. si Vq = 0,5 ds, la velocidad aparente inicial seria0,5(1+2310/293) = 4,4 ds). A1 final de la combusti611 la velocidad aparente serApricticamente Vq (algo mayor porque ahora 10s gases frescos estin a 575 K (P-15.27)en vez de a 293 K), puesto que la velocidad de 10s gases frescos seri despreciable. Lavariaci6n de la velocidad aparente de la llama se ha esquematizado en la curva V1 dela Fig. 16-15d; b) Si r(t) es la posici6n instantinea de la llama, la relaci6n es:siendo dV~/dtl, la variaci6n temporal del volumen de gases frescos VF por compresi6nisentr6pica debida a la expansi6n de 10s gases quemados (pVi= cte); c)


434 I. Martinez: TERMODINAMICA BASICA Y APLICADAsiendo F=l+A = 30,3 10s moles de gases frescos consumidos por mol de combustible,ng la cantidad de combustible consumido, ng, la inicial, TF la temperatura instantheade 10s gases frescos, Tl la inicial, Vo el volumen del recipiente, f = ndnB1 la fracci6nmolar de combustible quemado, n = plpl y n2 r p dl la relaci6n entre la presi6n finaly la inicial (P-15.27), de donde se despeja:d) Si t + 0, entones drldt = Vq(l+TBl/TF,) = cte (no depende de t) y desarrollando elresultado anterior se obtiene:de donde se deduce que como r varia linealmente con t, Bfl tarnbiin varialinealmente con t y por tanto 6 sera proporcioanl a t3; e) Se integra numiricamente enp desde pl hasta p2, obtieniindose:TB(p) Se despeja de pVo = PnBRTB + F(nBl - nB)RTFsiendo P = 32,l 10s moles de productos generados por mol de combustible consumido(P-15.27). t(p) se obtiene por integraci6n dedonde dpldr puede obtenerse numkricamente como Apl[r(p+Ap)-r(p)]. En la Fig.P-16.1 la se representa la variaci6n con la presi6n y en la Fig. P-16.11 b la variaci6ncon el tiempo. Los vdores extremos son p2/p1 = 10,82, TF2fT1 = 1,97, TB1/T1 = 7,96,TBJT1 = 10,24 y t2 = 0,172 S.


Cap. 16: PROCESOS DE COMBUSTION: MODELOS 435Fig. P-16.11w16.12. Estimar el tamaiio critic0 de un dep6sito de fertilizante de nitrato am6nico para que seautoinflame en un ambiente a 50 OC, sabiendo que se descompone seg6n una ley deArrhenius de constantes Ta=1,6.104 K y Ba=6.1013 s-1. El coeficiente global detransmisi6n de calor a la atm6sfera es 10 W.m-2.K-1, la entalpia de reacci6n-380 kJ.mo1-I y la densidad 1750 kg/m3. Calcular tambikn el increment0 detemperatura antes de la combusti6n.Sol.: L = 11,2 m, Tcr = 329,5 K (en (16.56) a = 11,8 y P = 4,8.10-19).16.13. Considkrese la combusti611 cuasi-estacionaria y esferica de una gota de benceno 3 rnrnde diajnetro en aire. Se pide:a) Poder calorific0 y temperatura adiab8tica.b) Constante de quemado, P (r2-f-cia)/t (compirese con el valor experimental (&-dinicia! )/t=0,96mm2/s.c) Posici6n relativa de la llama.Sol.: a) PC1 = 3,136 MJ/molc,~,, Td = 2500 K; b)con k = ko,/m= 0,024- - 0,06 W.m-l.K-l, p1f4 = 879 kgIm3, cp -~,,(2500 K) = 1300 J.kgl.K-l, Teb = 353 K, hlv = 394 kJ/kg, resultando P =0,22. m2/s, o en tkrminos de diimetros, 0,88.10-6 m2/s (un 8% menor delexpeirmental); c) sustituyendo en (16.93) se obtiene rll&rgota = 70.

More magazines by this user
Similar magazines