Repositorio Institucional de Documentos
Abstract
Se tiene una instalación para estudiar la fluidodinámica de un proceso de captura de dióxido de carbono basado en un ciclo de carbonatación-calcinación. Dicha instalación, presenta una problemática importante que es poder medir de manera continua y no intrusiva el caudal de partículas que circula por ella. El objetivo principal de este proyecto es realizar un estudio experimental y teórico de los ciclones tangenciales de la instalación para lograr resolver este problema. Tras este estudio, se obtendrá una correlación de la caída de presión en los ciclones en función de la carga de partículas con la que trabaja y la velocidad del gas. Aragonés de Villa-Ceballos, Mario; Martínez Santamaría, Ana
Full text
Análisisexperimentaldelainfluenciadelacargadesólidosenlacaídadepresiónenun ciclóntangencial 1
Análisisexperimentaldelainfluenciadelacargadesólidosenlacaídadepresiónenun ciclóntangencial 2
Análisisexperimentaldelainfluenciadelacargadesólidosenlacaídadepresiónenun ciclóntangencial 3 Agradecimientos ALuismiyAna,porsudirección,colaboraciónyasesoramiento.Pordarmela oportunidaddeintroducirmeenelmundodelainvestigación.Graciasalosdospor vuestraayudaconstanteparapermitirmerealizarycompletaresteproyectoconmás facilidad. Atodamifamilia,especialmenteamimadre,porsuapoyoincondicionalalo largodetodamicarrera.Siemprecreyeronenmisposibilidadesymeanimarontanto enlosbuenoscomoenlosmalosmomentos.AmiabueloÁngelquehoyporhoyya noestáentrenosotrosypormuypoconohavistocompletadomisueño,elcualle hubieradadomuchaalegría. Amiscompañerosycompañerasdelauniversidad,aloscualesconsidero grandesamigosporeltiempocompartidoylasexperienciasvividasjuntos.Gracias porvuestracompañíayapoyoincondicionalentodomomentoyhacerquetodohaya sidomuchomásdivertido. Amisamigosyamigasdetodalavida,porsucompañía,suapoyoabsolutoy sucompresiónentodomomento.Graciasporestarahísiempre. ARosaqueaunquehoydíanoestáenmivida,hasidounapersonamuy importanteparalograrllegaramimeta.Siempreteestaréagradecidoportuayuday comprensibilidad. Muchasgraciasatodos,sinvosotrostodoloqueheconseguidonohubiera sidoposible.
Análisisexperimentaldelainfluenciadelacargadesólidosenlacaídadepresiónenun ciclóntangencial 4
Análisisexperimentaldelainfluenciadelacargadesólidosenlacaídadepresiónenun ciclóntangencial 5 TABLADECONTENIDOS i. Índice5 ii.Índicedefiguras7 iii.Índicedetablas9 iv.Nomenclatura10 i.Índice 1.INTRODUCCIÓN:ALCANCEYOBJETIVO15 2.INFLUENCIADELACANTIDADDEPARTÍCULASENLACAÍDADEPRESIÓN DEUNCICLÓN21 2.1Introducción21 2.2Principiodefuncionamientodeunciclón23 2.3Clasificacióndeciclones25 2.4Geometríadeunciclón29 2.5Patróndeflujo30 2.6Eficienciadeunciclón33 2.7Caídadepresiónenunciclón35 2.8Modelosparaelcoeficiente,ζg,decaídadepresiónenvacío 37 2.8.1ModelodeShepherdyLapple(1939)38 2.8.2ModelodeAlexander(1949)39 2.8.3ModelodeCasalyMartines(1983)40 2.9Modelosparaelcoeficiente,ζs,decaídadepresiónencarga 41 2.9.1ModelodeBriggs(1946)41 2.9.2ModelodeSmolik(1975)42 2.9.3ModelodeKaskakov(1990)42 2.10ModelodeEnricoGriecoyLucaMarmo43 3.ANÁLISISEXPERIMENTAL45 3.1Descripciónteóricadeinstalación45 3.1.1Equiposprincipaleseinstrumentacióndelainstalación 45 3.1.2Funcionamientodelainstalación54 3.1.3Aplicacionesyventajas/desventajadelainstalación 56 3.2Planteamientosdelosensayos57 3.3Resultadosexperimentales61 3.3.1Comparaciónentrelasdiferentesalternativasteórico‐experimentales 61 3.3.2Seleccióndelacorrelaciónadecuadayajustedelamisma 69
Análisisexperimentaldelainfluenciadelacargadesólidosenlacaídadepresiónenun ciclóntangencial 6 4.CONCLUSIONES75 5.REFERENCIASBIBLIOGRÁFICAS79
Análisisexperimentaldelainfluenciadelacargadesólidosenlacaídadepresiónenun ciclóntangencial 7 ii.Índicedefiguras Figura1.DiagramageneraldelosprocesosylossistemasdecaptacióndeCO2 adaptadade[1]. Figura2.Movimientodelaspartículasenunciclón. Figura3.Vórticesenunciclóntangencial. Figura4.Curvasdeeficienciadelosciclonesenfuncióndeltamaño(diámetro)de partícula.Figuraadaptadade[4]. Figura5.Diferentesdisposicionesgeométricasenciclones[5]. Figura6.Dimensionesdeunciclón. Figura7.Perfildelavelocidadtangencialenunciclóntangencial.Figuraadaptadade [6]. Figura8.Perfildelasdiferentescomponentesdelavelocidadenunciclóntangencial deflujoinverso:(a)tangencial,(b)radial.(c)axial[6]. Figura9.Vistadelanteradelventilador. Figura10.Vistatraseradelventiladormásmotoreléctrico. Figura11.Conjuntodeválvulasdelcircuitodecarbonatación. Figura12.Lechofluidocirculantedelcircuitodecarbonatación. Figura13.Ciclóntangencial. Figura14.Filtrosdemanga. Figura15.Loopsealdelcircuitodecarbonatación. Figura16.Conexionesdelasloopsealsdeloscircuitosdecarbonataciónycalcinación Figura17.Conexionesdelaloopsealdelcircuitodecarbonatación. Figura18.Calentadoreléctrico(vistafrontal). Figura19.Calentadoreléctrico(vistalateral). Figura20.Válvulademariposa(Tipo:cuartodegiro). Figura21.Instalacióndeexperimentación.
Análisisexperimentaldelainfluenciadelacargadesólidosenlacaídadepresiónenun ciclóntangencial 8 Figura22.Conexióndelaloopsealconlechofluidocirculantedelcircuitode carbonatación. Figura23Configuracióndelcircuitodecarbonatacióndelainstalación. Figura24.Representacióndelosresultadosdelacaídadepresiónenvacio. Figura25.Representacióndelosresultadosdelacaídadepresiónenvacioteórica respectodelaexperimental. Figura26.Representacióndelosresultadosteóricosdelacaídadepresiónencarga respectodelosdatosexperimentales. Figura27.Representacióndelosresultadosdelacaídadepresiónencargaenfunción delamisma. Figura28.Representacióndelosresultadosdelacaídadepresiónencargaen funcióndelavelocidad. Figura29.Representacióndelosresultadosteóricosdelacaídadepresiónencarga respectodelosdatosexperimentales.
Análisisexperimentaldelainfluenciadelacargadesólidosenlacaídadepresiónenun ciclóntangencial 9 iii.Índicedetablas Tabla1.Modelosteóricosdelacaídadepresióndelosciclones. Tabla2.Datosgeométricosdelosciclonesdelainstalación. Tabla3.Condicionesdeoperacióndelciclóndecarbonatación. Tabla4.Rangodeparámetrosdeoperacióndelosdostiposdeensayos. Tabla5.NúmerosdeEulerenvaciocalculadospordiferentesmodelosteóricos. Tabla6.Resultadosexperimentalesdelosensayosenvaciorealizadosenelciclónde lainstalación. Tabla7.ValoresdelafunciónR2delosdiferentesmodelosteóricosanalizados. Tabla8.Resultadosexperimentalesdelosensayosencargarealizadosenelciclónde carbonatacióndelainstalación. Tabla9.ValoresdelafunciónR2losdiferentesmodelosteóricosanalizados. Tabla10.ValoresdelafunciónR2delosmodelosteóricoseinvestigadosanalizados.
Introducción:Alcanceyobjetivo 16 Porestemotivo,unaalternativaquehoyporhoyestáteniendomucha consideracióneslacapturayalmacenamientodedióxidodecarbono.Existen diferentesprocesosqueseexplicanacontinuación,peroelobjetivodetodoselloses elmismoreducirlasemisionesdedióxidodecarbonoalaatmósfera.Peroesta alternativapresentatambiénunaproblemáticaimportantequeeslanecesidadde almacenarelCO2capturado,yaquesuaplicaciónesmuyreducida.Además,en algunoscasosestossistemasdecapturanoevitanquesuusoreduzcaelrendimiento globaldelaplanta,debidoaquedichossistemasparaoperarrequierenabsorberuna ciertaenergía. Enlaactualidad,elCO2separadoenplantasindustrialesseutilizapara satisfacerlasdemandasdesuspropiosprocesosynoparaalmacenarlo.Lacapturade CO2hastaelmomentohasidoaplicadaenpequeñasplantasdegeneraciónde energíayaqueparalasgrandesplantasseencuentraenfasededesarrollo.La capturaenprocesosindustrialescorrientessellevaacabodesdehaceunos80años aunquelamayoríadelCO2queescapturadoesvertidoalaatmósferaporquenohay incentivosorequerimientosaseralmacenado.Eldióxidodecarbonoescapturadoen procesoscomolapurificacióndelgasnaturalyproduccióndegassintéticocon contenidodehidrógenoparaelaboraramoniaco,alcoholycombustibleslíquidos sintéticos[1].LamayoríadelastécnicasempleadasparalacapturadeCO2enlos ejemplosmencionadossonsimilaresalausadaenlacapturaenprecombustión. Existentresprocesosprincipalesdecapturadedióxidodecarbonocomoson:en Postcombustión,enPrecombustiónyoxicombustión. Lossistemasdecapturaenprecombustiónalteranelcombustibleantesdesu combustiónenunreactorconvaporyaireuoxígeno,paradarlugarunamezclade monóxidodecarbonoehidrógenoquesellamagasdesíntesis[1].Esteprimerproceso escomplejoycaro.Elmonóxidodecarbonoreaccionaconelvaporenunsegundo reactor,enelcualseproduceCO2ymáshidrógeno.Estamezclapuedesepararseen unflujodeCO2yotrodehidrógeno,frecuentementeporunprocesodeabsorción físicaoquímica.Elresultadoesuncombustiblericoenhidrógenoelcualactualmente puedeserusadoenalgunasaplicaciones,talescomo:pilasdecombustible,turbinasde gas,fabricacióndefertilizantes,calderasyhornos.Estatécnicaseaplicaenciclos combinadosdegasificaciónintegrada(CCGI). Lossistemasdecapturaenoxicombustiónutilizanoxígenopuroenvezdel aire,elcualestácompuestomayormentepornitrógeno,paralacombustióndel combustible.Elobjetivoesproducirungasdecombustióncompuesto principalmenteporaguaydióxidodecarbono.Estegasdecombustiónqueseorigina tienealtasconcentracionesdeCO2yaquenoexistenitrógeno[1].Posteriormente,el vapordeaguaseextraefácilmenteporcondensación.
Introducción:Alcanceyobjetivo 17 Porotrolado,esteprocesorequieredelaseparaciónpreviadeloxígenodel aireparaobtenerungasconunapurezadel95%al99%[1].Loquesetraduceenun altocosteparalaaplicacióndeestatecnología.Eldesafíodecómoseparareloxígeno delrestodelaireselacriogenizacióndelaire.Esteprocesodecapturapuede emplearseenhornos,calderasysistemasdeturbinasdegas. Lacapturadedióxidodecarbonoenpostcombustiónconsisteensepararel CO2delosgasesproducidostraslacombustión.Elobjetivodeestatecnologíaes poderaplicarlaengrandescentralesdecarbón,cicloscombinadosdegasnaturaly biomasayaquesegúnlosestudiosrealizadoseseconómicamenteviableynose necesitahacergrandesmodificacionesenlacentraldondesevainstalar.Existen diferentesmétodosdecapturadeCO2enpost‐combustiónaunquelamayoríade ellossehallantodavíaenfasesdedemostraciónoinvestigación. Algunosdeestosmétodossebasanenunprocesodeabsorciónquímica,enel cualutilizanunsolventeparacapturarlafraccióndeCO2existenteenlosgasesde combustión.Unadeestastécnicasesellavadoquímicoqueempleaunsolvente líquidocomoporejemplo,lamonoetanolamina[1].Otraconsisteenunciclode carbonatación‐calcinaciónquepormediodeunsorbentesólido,comoeselóxidode calcio,capturaeldióxidodecarbono.Dichoprocesoseexplicaconmayor detenimientoposteriormente.Peroambastécnicasestántodavíaenfasede desarrollo.Otrosprocesosenpostcombustión,queactualmenteseencuentranenfase investigación,sonlacriogenizaciónquecondensaelCO2existenteenlosgasesde combustión,laeliminacióndeCO2conotrosadsorbentessólidos,opasarlacorriente degasesdelacombustiónatravésdeunamembranayatraparelCO2[2].A continuaciónenlafigura1semuestranunarepresentaciónesquemáticadelos procesosdecapturadedióxidodecarbono. Figura1.DiagramageneraldelosprocesosylossistemasdecaptacióndeCO2adaptadode[1].
Introducción:Alcanceyobjetivo 18 Enesteproyectoseestudiaunainstalación,quesebasaenunsistemade capturadedióxidodecarbonoenpostcombustión,conunciclodecarbonatacióncalcinación.Unciclodecarbonatación‐calcinaciónesenunprocesodeabsorción química,enelcualseutilizaóxidodecalcio(CaO),conocidocomocal,comosorbente regenerabledeCO2.Esteprocesosellevaacaboendosetapas(reacciones):la primeradeellasseconocecomocarbonataciónylasegundacomocalcinación. Lareaccióndecarbonatación(reac.1)serealizaenreactoratemperaturas cercanasa650ºC.ElóxidodecalcioabsorbeelCO2existenteenunacorriente gaseosa.Deestaformaseevitaqueeldióxidodecarbonoseaexpulsadoala atmosfera.Finalmente,enestaetapaseobtieneunsólidodecarbonatocálcico (CaCO3)paraqueposteriormenteserecupereeldióxidodecarbonoenelcalcinador. Porotrolado,seproducecaloryaqueesunareacciónexotérmica. Encuantoalareaccióndecalcinación(reac.2),lasegundaetapa,serealizaen unreactoratemperaturaspróximasa900ºC.Endichoreactorseproduceuna reacciónquímicaenlacualseliberaelCO2existenteenlacaliza(CaCO3).Esta reacciónesendotérmicaporloqueesnecesariounaportedecalor,elcualseráel únicoconsumoenergéticodeesteciclo.Lasreaccionesquímicasqueseproducen sonlassiguientes: Reaccióndecarbonatación: 1800 ](reac.1) Reaccióndecalcinación: 1800 (reac.2) Porotrolado,losreactoresenloscualesserealizanestasreaccionespueden serdecualquiertipo:calderasconvencionalesdecombustión,quemadoresdegaso lechosfluidos.Estaúltimaopcióneslamejordetodasporqueunlechofluidoesun equipoquepermiteunabuenamezcladegasysólidos.Ademáspuedecontrolarque latemperaturaseahomogénea.Estosdosfactoressonimportantesparatenerunas reaccionesfavorablesyconellounbuenrendimientodelciclo.Porlotanto,la instalaciónenlacualseexperimentó,comoseexplicaráunpocomásadelante, constadedoslechosfluidos.
Introducción:Alcanceyobjetivo 19 HayquemencionarqueestanuevatecnologíadecapturadeCO2pormedio deunciclodecarbonatación‐calcinaciónpuedeseraplicada,porejemplo,en centralestérmicasdecarbónocicloscombinadosexistentesynuevosyaqueel consumodeenergíanecesarioparalaregeneracióndelsorbente(calcinación)se podríarecuperarintegrandolosflujosdecaloraaltatemperaturadelciclodecaptura enunciclodevapor.Ademásconlautilizacióndeestatécnicatambiénesposible capturaróxidosdeazufre(SOx). Lainstalacióndelaquesedisponeesdeflujofrío,enellanoseproducen reaccionesquímicas,yaquesolosepretendeanalizarlafluidodinámicadeeste proceso.Laproblemáticaquepresentaestainstalacióndecapturadedióxidode carbonoespodermedirdemaneracontinuaynointrusivaelcaudaldepartículas. Actualmenteestamedidadelcaudaldepartículasesnocontinuadebidoaquenose puedemedirentodomomentoeintrusivaporqueafectaalfuncionamientodela instalaciónproduciendocambiosdepresionesypérdidasdepartículas,comose explicarámásadelanteenestamemoria.Elobjetivoprincipaldeesteproyectoes realizarunestudioexperimental(enlainstalación)yteóricodelosciclones tangencialesdelainstalaciónparalograrresolveresteproblema.Endichoestudio,se obtendráunacorrelacióndelacaídadepresiónenlosciclonesenfuncióndelacarga departículasconlaquetrabajaylavelocidaddelgasalaentrada.Deestamanera, comotantolacaídadepresióncomolavelocidadsepuedenmedirdeformacontinua ynointrusiva,sepodrácalcularelcaudaldepartículasdeformacontinuayno intrusiva,apartirdedichacorrelación,comoseexpondráconmásdetenimientoalo largodelamemoria.Porotrolado,comoambosciclonestieneslasmismas dimensionesgeométricasyoperaránbajocondicionesdetrabajoiguales,entoncesla correlacióndelacaídadepresiónseráválidaparalosdos.
Introducción:Alcanceyobjetivo 20 Lasetapasquesesiguenparalograrelobjetivodeesteproyectoseenumeran acontinuación: Enprimerlugar,serealizaunestudiodelprincipiodefuncionamiento, clasificación,patróndeflujoyeficienciadeseparacióndelosciclonestangenciales. Posteriormente,serealizaunestudiomásexhaustivoparalacompresióndelacaída depresiónqueseproduceenunciclón.Porotroladoycomoseveráensu correspondienteapartado,seelaboraunadocumentacióndelosdiferentesmodelos teóricosparaelcálculodelacaídadepresiónenunciclóntangencial.Paraello,seha buscadoeinvestigadoenbibliografíaespecíficaendichotema.Contodoesto,en estaprimeraetapasepretendeconocerlasprincipalescaracterísticasqueserigenen unciclóntangencial. Ensegundolugar,serealizanunaseriedeensayosexperimentalespara obtenerlacaídadepresióndelciclón.Estosensayosseefectúanbajodiferentes condicionesdetrabajo,variandolavelocidaddelgasalaentradadelciclónylacarga departículas.Semidelapresiónalaentradaysalidadelciclón,paraconocerla pérdidadecarga. Enlatercerayúltimaetapaseevalúanlascorrelacionesexistentesdelos diferentesmodelosdecaídadepresiónanalizados.Seguidamente,seestudiayajusta unacorrelacióndelacaídadepresiónpordiferentescriteriosdecálculoyseelige aquellaquemejorseaproximaalosdatosexperimentales.Paraconcluir,esta correlaciónsevalidacomparándolajuntoconlosmétodosteóricosexaminadosy resultadosexperimentalesparademostrarquelacorrelaciónobtenidaeslamás idóneaparadeterminarlacaídadepresiónlosciclonesdelainstalación.Deeste modo,conestacorrelaciónsepuededeterminardemaneracontinuaynointrusivael caudaldepartículasquecirculaporlainstalación,cumpliendoelobjetivodeeste proyecto.
Introducción 21 2.INFLUENCIADELACANTIDADDEPARTÍCULASENLACAÍDADEPRESIÓNDE UNCICLÓN 2.1Introducción Losciclonessondispositivoscentrífugosqueseutilizanparalaseparaciónde partículassólidassuspendidasenunacorrientegaseosa.Actualmente,estosequipos sonlosmásempleadospararealizardichatarea.Principalmente,debidoasubaja inversión,fácilmantenimiento(noposeenpiezasmóviles)yquepuedentrabajaren condicionesdefuncionamientocomplicadas.Básicamente,unciclóntangencialestá compuestoporuncilindroverticalconunfondocónico.Normalmentetienen entradatangencialydescargaaxialaunquetambiénexistenconentradaydescarga axial. Losciclonesencomparaciónconotrosdispositivosdeseparaciónde partículassólidascomolascámarasdesedimentacióngravitacional,filtro,lavadores yprecipitadoreselectrostáticopresentanlassiguientesventajas[3]: Elmaterialcapturadoquedasecoy,normalmente,aprovechable. Bajocostedeinversiónymantenimiento. Porlogenerallosciclonestieneneficienciasmayoresquelacámarade sedimentacióngravitacional. Puedenserusadosencondicionesextremas,enparticularaaltastemperaturas ypresionesyconsustanciasquímicamenteagresivas. Notienenpartesmóvilesysonmuyrobustos Lacaídadepresiónesconstantecuandosetrabajaenlasmismascondiciones deoperación. Puedeserconstruidoapartirdediferentesmateriales,comoporejemplode láminadeacero,metalesfundidos,aleaciones,aluminio,plásticos,cerámicas, etc. Puedeserequipadoconrevestimientosdematerialescapacesderesistirla erosiónycorrosión. Enalgunosprocesos,elciclónpuedeencargarsedesólidosviscososopegajosos conayudadeunriegolíquidoadecuado. Losciclonestambiénsepuedenutilizanparaseparargotasdelíquidoen suspensiónenungasypartículassólidasensuspensiónenunlíquido.Eneste últimocasosellamanhidrociclones.
Introducción 22 Porotrolado,losciclonestambiénpresentanunaseriededesventajas[3] frenteaotrostiposdedispositivosdeseparacióndepartículassólidascomolos anteriormentemencionados: Losciclonesnosonlosequiposmásadecuadosparasepararpartículascon diámetrosmenoresde5μm,aunqueenalgunoscasosconcretospartículas muchomáspequeñaspuedenllegaraserseparadasconunaceptable rendimiento. Tienerendimientosbajoscuandooperafueradesupuntodefuncionamiento óptimo,esdecir,contamañosdepartículasmáspequeñosqueparalosquefue diseñado. Espropensoaerosionarseproduciéndoseundesgasteyensuciamientodel mismo,sobretodosilossólidosqueestáprocesandosonabrasivosoviscosos. Tienealtoscostesdeexplotacióndebidoalaenergíarequeridaparasuperarla pérdidadepresión. Losciclonespresentaneficienciasmenoresquelosfiltrosdemangas,lavadores yprecipitadoreselectrostáticos. Sielciclónnoesdiseñado,mantenidoyoperadocorrectamenteentonces trabajarápordebajodelrendimientoesperado.
Principiodefuncionamientodeunciclón 23 2.2Principiodefuncionamientodeunciclón Elfuncionamientodelosciclonessebasaendaralflujodegasesunarápida rotación.Deestaforma,laspartículassondesplazadasporlaaccióndelafuerza centrífugahacialasparedesdelciclón,loquepermitequelaspartículasmáspesadas (mayordiámetroodensidad)seaglomerenyseseparendelacorrientegaseosa. Comosepuedeobservarenlafigura2,laspartículasmáspesadasseaproximanhacia lazonadelaparedmientrasquelasmenospesadassequedanenlazonacentral,no lleganasepararsedelacorrientegaseosaysalenporlasalidasuperiordelciclón.Por otrolado,conlaactuacióndelafuerzagravitacional,laspartículasmáspesadas desciendenporlasparedeshacialasalidainferiordelciclón. Paraexplicarlatrayectoriaquedescribeelgasenelinteriordeunciclón,hay quemencionarquedichatrayectoriadescribeundoblevórtice,endondeelgas desciendetrazandounaespiralenlaregiónexternadelinteriordelciclón. Posteriormenteelgasasciendelimpiodepartículastrazandootraespiralporla regióncentraldelinteriordelciclón,comomuestralafigura3. Figura2.Movimientodelaspartículas enunciclón. Inicialmente,elgasconpartículascirculaporlaentradadelciclónconuna trayectorialinealycuandoelgasentraenelcuerpodelciclónsumovimientocambia aunatrayectoriahelicoidal.Estemovimientoesocasionadoporeldiseñodela Figura3.Vórticesenunciclóntangencial.
Principiodefuncionamientodeunciclón 24 entradadelciclón,queproducequeelgasentreconcomponentetangencialhacialas paredesinternasdelcuerpodelciclón.Eltipodeseccióndeentradaquemásseusa paralograrestepropósitoeslasecciónrectangular,talycomopuedeverseenla figura3. Enlaprimerapartedelcuerpodelciclóndondesudiámetroesconstante,el gascargadoconpartículas,describeunatrayectoriaaxialmentedescendienteen espiral.Esenestazonadondeseproducelaseparacióndelaspartículasdemayor tamañodelgasdebidoalaactuacióndelafuerzacentrifuga,éstasseaglomeranyse deslizandescendiendoporlasparedesdelciclón.Lafuerzacentrífugavaríasu magnituddependiendodelavelocidadtangencialdelgasydelradiodelcuerpodel ciclón.Porloque,estafuerzacentrifugapuedesermuchomayorquelafuerza gravitacional. Lapartecónicadelciclónesdiseñadadetalformaqueelgasesforzadoa entrarenlaregióninternadelciclón,dondeelmovimientoenespiral(endirección axial)esdirigidohaciaarriba.Aldisminuirelradioseproduceunaumentodela velocidadtangencialdelgas,estoexplicaporquéenelconodelciclónesdondese presentalamayorrecoleccióndepartículas,especialmentelasdemenortamaño (diámetro).Laspartículassiguendescendiendodeslizándoseporlasparedesinternas hastallegaralasalidainferiordelciclóndondesonalmacenadasenunrecipiente determinadoorecirculadascomoeselcasodelainstalaciónanalizadaeneste proyecto. Finalmente,elgasasciendeconunatrayectoriaenespiralporlaregión centraldelciclón.Estacorrientedegasseencuentra(casiensutotalidad)limpiade partículasdetamañodepartículaparaelcualelciclónfuediseñadoqueseparase. Dichacorrientegaseosasedirigehacialasalidasuperiordelciclón.Estasalida,quese conocecomo“vortexfinder”,setratadeuntubodediámetroconstanteyde longitudsimilaralapartedelcuerpodelciclóndediámetroconstante,dependiendo deldiseñoquesequisierarealizar.Porlotanto,elgaslimpioatraviesael“vortex finder”ysaledelseparadordepartículasalsiguienteprocesodelimpieza, almacenamientooalaatmósferadependiendodelaaplicacióndequesetrate.
Clasificacióndeciclones 25 2.3Clasificacióndeciclones Hoyendía,losciclonessonundispositivodeseparadordepartículasyamuy estudiado,eldiseñodeéstossebasanormalmenteenfamiliasdeciclonesquetienen proporcionesdefinidas.Estosdispositivoscentrífugospuedenserclasificadossegún [3]: A. Laeficienciadelciclón. B. Laconfiguracióndelaentrada. C. Laformadesucuerpo. D. Ladireccióndelfluidoenlaentradaysalidadelosciclones. E. Ladisposicióngeométrica. Encuantoalosciclonesdeentradatangencialquesepuedenclasificarsegún sueficienciaexistentrestipos[4],comopuedenapreciarseenlasiguientefigura4: Ciclonesdealtaeficiencia. Ciclonesconvencionales. Ciclonesdealtacapacidad. Figura4.Curvasdeeficienciadelosciclonesenfuncióndeltamaño(diámetro)departícula. Figuraadaptadade[4]. 0 20 40 60 80 100 0 5 10 15 20 η% dp,µm Altaeficiencia Convencional Altacapacidad
Patróndeflujo 32 Ladistribucióndelflujoydelapresiónalolargodeunciclón,sepuede comprenderconlapresiónestáticaydinámica:py½ρv2,respectivamente,dondeρ esladensidaddelfluido.LaecuacióndeBernoulliparaunflujoestablesinfricción afirmaque[3]: 1 22 Eltérminodelaizquierdadeestaecuaciónesconocidocomo“Trinomiode Bernoulli”.Enelcasodeloscicloneselprimerytercertérminodedichotrinomioson losimportantes,yaquelasdiferenciasdealturasnosonmuygrandes. Estaecuacióndemuestraquelapresiónestáticaydinámicasepueden intercambiarenelcampodefluido.Enlasáreasdondelavelocidadesalta,lapresión estáticaserábajayviceversa.Porconsiguiente,examinandodenuevolafigura anteriorsepuedeafirmarquelapresiónestáticaenlazonainternadelcilindrode Stairmandesbajayfueradeél,esalta.Encuantoalapresióndinámica,es proporcionalalcuadradodelavelocidad,porloquesudistribuciónserásimilarala delavelocidad,cuandoéstaesaltalapresióntambiényviceversa. Figura8.Perfildelascomponentesdelavelocidadenunciclón:(a)tangencial,(b)radial.(c)axial[6].
Eficienciadeunciclón 3333 33 2.6Eficienciadeunciclón Laeficienciaorendimientodeseparaciónglobalsedefinecomoelporcentaje delamasadepartículasqueentraalciclónconlacorrientegaseosayquees separadadelamismaenelinteriordelciclón.Aunqueelrendimientodependedel tamañoydeladensidaddelaspartículassólidas,noessuficienteparacaracterizar adecuadamenteunciclón.Dichoesto,lasiguienteecuacióndeterminael rendimientodeseparaciónglobal,dondeMceslamasadepartículascapturadasyMa eslamasadepartículasqueentranenelciclón: 3 Existenbásicamentedosconceptosenlosquesebasanlosmodelosde eficienciadeunciclón:elconcepto“órbitadeequilibrio”(equilibrium‐orbit)yel concepto“tiempodevuelo”(time‐of‐flight)[3].Ambosdependendeltamaño(d)yla densidaddepartículas(ρp),incrementándosesudependenciaamedidaqueaumenta eltamañodedichaspartículas.Porconsiguiente,estosmodelossonapropiadospara determinarlaeficienciadeunciclón.Elprimerosebasaenunbalancedefuerzas sobrelapartículaqueestárotandoenunvolumendecontrol(cilindroidealde Stairmand,VC).Elsegundoconceptosefundamentaenlaideadesiunapartícula, inyectadaenalgunaposiciónlateralenlaentrada,tieneeltiemposuficientepara alcanzarlapareddeciclónyserrecogidaantesdellegaralaparteinferiordelciclón. Araízdeestosdosconceptoshanaparecidoalgunosmodeloshíbridosentreambos, considerandotantounintercambiodepartículasatravésdelvolumendecontrol comoatravésdelamigracióndepartículaalapareddeciclón. Enconclusión,sepuedecomentarquelosmodelosobtenidosaraízdel conceptode“tiempodevuelo”soncompletamentediferentesalosqueseobtienen conelconcepto“órbitaenequilibrio”,perodanvaloresdeeficienciacorrectosy similares,enunaampliagamadediseñosdeciclónycondicionesdefuncionamiento. Encuantoalaformadelafuncióndelaeficienciadelciclón,queesdeterminadacon estosmodelos,idealmenteesunescalóndevaloriguala1apartirdelpuntocríticoo decorte(d50valordetamañodelapartículaparaelcualelrendimientodelciclónes del50%).Sinembargo,tieneperfildecurva“S”,estecambiosedebeprincipalmente alassiguientescausas: Debidoalaturbulenciacreadaporlacorrientegaseosaalgunasdelas partículasquedeberíanseratrapadasenelciclónnolosony viceversa. Existenpartículaspequeñasquenodeberíandesersepararasperoal adherirseaotrasdemayortamaño,soncapturadasconellas.
Eficienciadeunciclón 34 Laspartículasacausadesuagitadomovimientosedancontinuos choquesconlasparedesdelciclónyentreellasmismas, produciéndoseundesgasteyreduccióndetamaño.Porloque, algunasdeestaspartículasquetendríanquesercapturadasnoloson ydanlugaraunamayorconcentracióndepartículasdemenor tamaño. Laposicióndelapartículaalaentradapuedeasucaptura.
Caídadepresiónenunciclón 35 2.7Caídadepresiónenunciclón Existeuncompromisoimportanteentrelacaídadepresiónylaeficienciade capturaenunciclóntangencial.Engeneral,cuantomayoreslaeficienciamayores sonlaspérdidasdepresiónyenconsecuenciasuconsumoenergéticotambién.Esto esdebidoengranparteaquesilavelocidaddelacorrientegaseosaenlaentradadel ciclónseincrementa,tantolaeficienciacomolapérdidadecargasonmayores.Enel casodelacaídadepresión,esteaumentoesproporcionaladichavelocidadal cuadrado. Elprocesoestándarquesellevaacaboenlaindustriaparaevaluarlapérdida decargaenunciclónesmedirlapresiónestáticaenlapareddesutuberíadeaguas arribaydeaguasabajo.Esteprocesoescomplicadorealizarloenelinteriordelciclón debidoalmovimientohelicoidalqueexisteenelflujo.Lacaídadepresiónenun ciclónsedividenormalmenteentrescontribucionesozonas.Porloquese consideraraquelacaídatotaldepresióneslasumadetalescontribuciones: 1. Laspérdidasenlaentradadelciclón. 2. Laspérdidasenelcuerpodelciclón. 3. Laspérdidasenel“vortexfinder”. Laspérdidasenlaentradaseoriginanporlaexpansiónocompresióndela corrientegaseosa,ladisminucióndelapresiónestáticayporlafricciónenlaentrada delciclón.Dichaspérdidas,sonfrecuentementedespreciablesencomparacióncon lasotrasdoscontribuciones.Lamenordeestaspérdidassedaenlosciclonescon entradatangencial.Encuantoalaspérdidasenelcuerpodelciclón,quesonmayores quelasanteriores,seproducenporlafricciónconlaspropiasparedesysugran trascendenciaradicaenquelimitanlaintensidadenelvortex.Laspérdidasenel “vortexfinder”sonlasmásimportantes.Enconclusióntodasestaspérdidasdanlugar enlacorrientegaseosaaunadisminucióndesupresióndinámicayestáticaentrela entradaylasalida. Lacaídadepresiónquesufreelflujodentrodeunciclón,Δp,esproporcional alcuadradodelavelocidaddelflujo,yaqueenunciclónelflujoesturbulento.Esta caídadepresiónescaracterizadaapartirdelnúmerodeEuler[3]: ∆ 1 2 4 Dondevzeeslavelocidadaxialdelacorrientedegasenlaentradadelciclón. Enelanálisisdelacaídadepresiónenunciclóntangencialinfluyenunaserie deparámetroscomoelnúmerodeFroude,Reynolds(operabajocondicionesdeflujo turbulento),larelaciónentreladensidaddepartículasyladelgas,etc…Perolos
Caídadepresiónenunciclón 36 factoresmásinfluyentessonlageometríadelciclónylacargadepartículas.Esto explicalaestructuraquetienenlamayoríadelosmodelosdecaídadepresiónenun ciclón.Enalgunosdeellos,losdosefectossonintroducidosindependientemente,yel productodeambosseráelnúmerodeEulerdelciclón,ζc,[6]: 5 DondeζgrepresentaelnúmerodeEulerquehabríaenausenciadepartículas (vacío),ζ gesunafuncióndelageometríadelciclón.Encuantoaζ s1estimala reduccióndelacaídadepresióncuandoelciclónoperaconunaaltacargade partículas.Enconclusiónatodoloexplicado,sepuedededucirquelacaídade presiónenunciclónesfuncióndesugeometría,delavelocidaddelflujoalaentrada asícomodeladensidadycaudaldepartículasquecirculanporél.
Modelosparaelcoeficiente,ζg,decaídadepresiónenvacío 37 2.8Modelosparaelcoeficiente,ζg,decaídadepresiónenvacío ComosehamencionadoanteriormentesesabequeelnúmerodeEulerpara lacaídadepresiónenvacíoenunciclóntangencialesfuncióndelasdimensionesdel ciclón,aunqueenalgunosmodelosestenúmerotambiéndependelaspropiedades delgasysuscondicionesdeoperación.Actualmenteexistenvariosmodelosquehan sidodesarrolladosparaestimarlacaídadepresiónenvacío,entreellosdestacanel de:ShepherdyLapple(1939),Alexander(1949),Stairmand(1949,1952),Barth (1956),Muschelknautz(1972),CasalYMartines(1983),etc… Sinembargo,existendiscrepanciasentreellossobrelacontribuciónmás relevanteenlacaídadepresión.Paraalgunosautores,comoShepherdyLapple (1939),Alexander(1949),Stairmand(1949,1952),CasalYMartines(1983),lamás importanteeslacorrespondientealazonadeentradayalcuerpodelciclón. MientrasqueotrosmodeloscomoBarth(1956)yMuschelknautz(1972),tratanlas pérdidasenlazonadel“vortexfinder”comolasmásimportantes. Algunosdeestosmodelossonempíricosyotrosteóricosperoningunode ellospredicelacaídadepresiónconprecisiónparaunampliorangodediseñosde ciclones.Además,sedebemencionarquealgunosmodelosdecálculocomoelde Stairmand,BarthyMuschelknautzsoncomplejos.Senecesitaconocercondiciones defuncionamiento(cinemática)quenosonfácilesdeobtenercomosonlasdel “vortexfinder”.Poresto,acausadesudificultadyaqueconotrosmodelosmás sencillossepuedenobteneryaestimacionesbuenasdelacaídadepresión,nosehan utilizadoparaelestudiodelciclóntangencialtratadoenesteproyecto.
Modelosparaelcoeficiente,ζg,decaídadepresiónenvacío 38 2.8.1ModelodeShepherdyLapple(1939) ElModelodeShepherdyLapple,queespuramenteempírico,determinauna simpleexpresióndelnúmerodeEuler,ζ g,enfuncióndelageometríadelciclón, puestoquesuponeunnúmerodeReynoldslosuficientealtocomoparadespreciar suinfluencia.Lasdimensionesqueutilizadichomodelosoneldiámetrodelcuerpo del“vortexfinder”ydelaentradadelciclón,enelcasodeunciclónconunaentrada tangenciallaanchuraylaaltura.Estemodeloesampliamenteutilizadodebidoasu buenarelacióndecálculoyaproximación.ShepherdyLappledieronlasiguiente expresiónparaelnúmerodeEuler[6]: 16 6 ExisteunmodelomuysimilaraesteparadeterminarelnúmerodeEuler,Eug, elcualseconocecomoelmodelodeCoker.Dichomodelotienelamismaestructura queelmodelodeShepherdyLapple.Paraunciclónconunaentradatangencial dependedelamismageometríaloúnicoquecambiaeselvalordelacorrelaciónde 16a9,47.Porlotantoestemodelotambiénserápuramenteempírico.Paraalgunos diseñosdeciclonesconestemodeloseobtienenmejoresaproximacionesdelacaída depresiónenvacíoqueconelmodelodeShepherdyLapple. Porotrolado,lavelocidadcaracterísticadelciclónqueseutilizaraen“ambos modelos”eselvalorpromediodevelocidadaxialenlaentradadelmismo,esdecir: ó 7 DondeQciclónenteselcaudaldelfluidoqueentraenelciclón.
Modelosparaelcoeficiente,ζg,decaídadepresiónenvacío 39 2.8.2ModelodeAlexander(1949) ElmodelodeAlexanderdeterminaelnúmerodeEuler,ζg,enfuncióndela geometríadelciclónperoadiferenciadelmodeloanterior,Alexanderintrodujola dependenciadelaviscosidaddelgasconlatemperaturaapartirdelexponentede vórtice“n1”.Porotroladoenestemodelodelageometríadelciclónseutilizanel diámetrodelcuerpodelciclónyeldel“vortexfinder”.Tambiénsehademencionar quesetratadeunmodelosemi‐empíricoyqueparasuestudiosenecesitanconocer máspropiedadesdelgasqueenelmodelodeShepherdyLapple.Peroaunasísu cálculonoresultadifícilysuaproximaciónesmejor.Alexanderdefiniólasiguiente expresióndelnúmerodeEuler[6]: 4.62 11 8 0.81 142 31 0.2211 1.529 10.67. 283.10 Encuantoalavelocidadcaracterísticadelciclónqueseutilizaenestemodelo eselvalormediodelavelocidadaxialalaentradadelciclón,vze,altratasedeun modeloquetambiénconsideralaspérdidasenlaentradaycuerpodelciclóncomo lasmásrelevantes.
Modelosparaelcoeficiente,ζg,decaídadepresiónenvacío 40 2.8.3ModelodeCasalyMartines(1983) ElmodelodeCasalyMartines,espuramenteempíricoydeterminaelnúmero deEuler,ζg,enfunciónsolodelageometríadelciclón,enconcretodeldiámetrodel “vortexfinder”ydelasdimensionesdelaentradaqueparaunciclóndeunaentrada tangencialsonlaanchuraylaaltura.Estemodelocomosepuedeobservarpartedel modelodeShepherdyLappleyapartedeobtenerotrovalordelacorrelación,como fueelde11,3ydeelevarlaparteenfuncióndelageometríadelaecuación, añadieronotracorrelaciónindependientedevalor3,33.Deestamanera,laexpresión delnúmerodeEulerqueseobtuvofuelasiguiente[6]: 3.3311.3 11 ConestastresmodificacionesCasalyMartinesmejoraronlaestimación presióndelcálculodelacaídadepresiónenvacíoenunrangodediseñosde ciclones,encomparaciónconelmétododeShepherdyLapple.Encuantoala velocidadcaracterísticadelciclónqueseutilizaenestemodeloeselvalormediode lavelocidadaxialalaentradadelciclón,vze,porquecomoocurreconlosmodelos explicadostratalaspérdidasenlaentradaycuerpodelciclóncomolasmás relevantes.
Modelosparaelcoeficiente,ζs,decaídadepresiónencarga 41 2.9Modelosparaelcoeficiente,ζs,decaídadepresiónencarga Existeconstanciadequeenunciclónlacaídadepresiónseveinfluenciada cuandoestáoperandoconunaconcentracióndesólidosalta,detalformaqueal aumentarlacarga,lacaídadepresióndisminuye.Esteefectoseatribuye principalmentealapropiainerciadelaspartículas.Porestemotivoseintroduceun coeficientedecorrección,ζs,alnúmerodeEulerdelciclón,ζc. Estareduccióndelacaídadepresióncuandoexisteunaaltaconcentraciónde partículashasidoestudiadaporvariosinvestigadoreseneltemacomo:Briggs(1946), Smolik(1975),Baskakov(1990),…Dichosautoresobtuvierondiferentesmodelosde cálculodeζs,ycomoocurrecuandosetrabajaenausenciadecarga(vacio)existen discrepanciasentreellos.Enestaocasióndichasdiscrepanciasradicanespecialmente enlacuantificacióndesuefecto. Porotraparte,aligualqueocurreconlosmodelosenvacíoningunodeestos modelospredicelacaídadepresiónencargaconunagranprecisiónparaunamplio rangodediseñosdeciclones. 2.9.1ModelodeBriggs(1946) ElmodelodeBriggsdeterminaunaexpresióndelnúmerodeEuler,ζ p,en funcióndelaconcentracióndepartículasyladensidaddelgas.Apesardeserla correlaciónmásantiguadetodas,conellaseobtienenunosresultadosmuy correctos.Estosedebeaqueesunmodelomuchomásconservadorqueelresto.Por todoesto,estacorrelaciónesampliamenteutilizadaporsubuenarelacióndecálculo yaproximación.Briggsobtuvolasiguienteexpresiónparaelcoeficientede corrección,ζs[6]: 1 10.0086 12 DondeCsieslaconcentracióndepartículasengsol/kggas.
Descripciónteóricadelainstalación 48 CICLÓNTANGENCIAL:Estainstalaciónconstadedosciclonestangencialesde igualescaracterísticasgeométricas,unoenelcircuitodecarbonatación(véase figura13)yotroeneldecalcinación.Encuantoasugeometría,seconoceenla siguientetabla: DATOSGEOMETRICOS H(mm)837 Hc(mm)439 Acentrada(m2)0,003 a(mm)43 b(mm)71 Dx(mm)93 Dc(mm)180 S(mm)112 Dd(mm)70 Tabla2.Datosgeométricosdelosciclonesdelainstalación. Portenerunconceptobásicodelfuncionamiento,aloperarenlos ensayossoloconelcircuitodecarbonatación,seanalizósoloelciclóndeeste circuitoconunacorrientedeaireycaliza.Porlotanto,lafunciónque desempeñadichociclónesladesepararlacalizadelacorrientedeaire. Figura13.Ciclóntangencial.
Descripciónteóricadelainstalación Descripciónteóricadelainstalación 49 FILTROSDEMANGA:Losfiltrosdemangasonequiposquetienencomofunción separarlaspartículassólidasquearrastraunacorrientegaseosamedianteun medioporoso.EnestainstalacióndecapturadeCO2sehadecididocolocarunoen lasalidadecadaciclón,conelfinde,recogerpartículasdecalizaquenosehan podidoatraparconsurespectivociclón.Losfiltrosdemangasseleccionadospara estainstalación,quepuedenobservarseenlafigura14,estánformadosporuna bolsadefibrasintéticaconformademanga.Dichabolsaescolocada(sincarcasa) enunossoportes. Figura14.Filtrosdemanga.
Descripciónteóricadelainstalación 50 VÁLVULANOMECÁNICAOLOOPSEAL:Dispositivodelainstalación,queestá construidodemetacrilato,funcionacomosellodepresionesentreelciclón(baja presión)ysucorrespondientelechofluidocirculante(amáspresión).Ademásde recircularlaspartículassólidasentreambosequipos.Estainstalacióntienedos loopseals,unoparacadacircuito.Enlafigura15puedeobservarselaloopsealdel circuitodecarbonatación. Estosequiposconstandedosentradas:unadepartículasqueesla salidadelciclónyotradeaireprovenientedelventilador,yunasalidade partículasdirigidaallechofluidocirculante.Dichoesto,unacaracterística importantedeestosequiposesquesegúnsedistribuyansusentradasysu salidacambianlaconfiguracióndelainstalacióndependiendodelestudioy ensayoquesevayaarealizarenella(véasefigura16).Encuantoasu funcionamientoenestainstalación,básicamenteconsisteenfavorecerla circulacióndelacalizaentreelciclónyellechofluidocirculantequeestán conectados. Enlosensayosrealizadosenestainstalación,comosehamencionado anteriormente,solosetrabajóconunodeloscircuitosenconcretoconeldel carbonatador.Porconsiguiente,lasalidadelaloopsealestádirigidaallecho fluidocirculantedelmismocircuito,véasefigura17.
Descripciónteóricadelainstalación Descripciónteóricadelainstalación 51 Figura15.Loopsealdelcircuitodecarbonatación.Figura16.Conexionesdelasloopsealsdeloscircuit o decarbonataciónycalcinación. Figura17.Conexionesdelaloopsealdelcircuitodecarbonatación.
Descripciónteóricadelainstalación 52 CALENTADORELÉCTRICO:Estedispositivodelainstalación,elcualsepuede observarenlasfiguras18y19,esunsencilloequipoeléctricoquetransforma laenergíaeléctricaenenergíacalorífica.Dichodispositivo,siseempleatienela funciónenlainstalacióndecalentarelairedefluidificacióndelmodelodel calcinador.Esdecir,elcalentadoreléctricoesúnicamenteunrequerimientodel escaladodelainstalaciónparamantenerlosnúmerosadimensionalesylos ratiosdedensidadentrelacorrientedeaireylaspartículasdecaliza. EQUIPOSDEMEDIDA(INSTRUMENTACIÓN):Paraelestudiodeestainstalación sehanempleadounaseriedeequiposdemedidadepresión,temperaturay velocidaddelacorrientedeaire.Estosequipossecolocanenposiciones estratégicasdelapropiainstalación,comoporejemploenloslechosfluidos circulantes,loopsealsasícomoenlosciclones.Deestamanera,sepueden obtenermedidas,quesonenviadasaunsistemadeadquisicióndedatos,que permitenconocerentodomomentocómoestáoperandolainstalación. Figura18.Calentadoreléctrico(vistafrontal).Figura19.Calentadoreléctrico(vistaperfil).
Descripciónteóricadelainstalación Descripciónteóricadelainstalación 53 Paramedirlapresiónseutilizantransductoresdepresióndiferencial.En lainstalaciónexistentransductoresdediferentesrangosdelectura/medida:de 0a10mbar,de0a50mbaryde0a100mbar. Paramedirlatemperaturaylavelocidaddelacorrientegaseosaseusa untermo‐anemómetrodesondadehilocaliente.Elrangodemedidade velocidaddeltermo‐anemómetroutilizadoesde0a20m/s. Enlafigura20sepuedeobservarlaválvulademariposadelcircuitode carbonatación,queseutilizaenlosensayosrealizadosconcargaparamedirla circulacióndepartículasconlaqueseestátrabajando. Figura20.Válvulademariposa(Tipo:cuartodegiro).
Descripciónteóricadelainstalación 54 3.1.2Funcionamientodelainstalación Enesteapartadoesimportantedetallarquesevaaexponerel funcionamientodeestainstalaciónenflujofrío(figura21),nodelfuncionamiento real,esdecir,conreaccionesquímicasquetendríaunainstalacióndeestetipo(ciclo decarbonatación‐calcinación).Porestemotivo,enestainstalacióndeflujofriose decidióoperarconunacorrientedeaireycalizacomosorbente,yaqueesbarata, disponibleyabundante. Elfuncionamientodelainstalación,dondesehanllevadoacabolosdistintos ensayos,sepuedeexplicarapoyándoseenelprocesoquerigeunciclode carbonatación‐calcinación.Enprimerlugar,enellechofluidocirculantede carbonataciónsemezclaunacorrientedeaire,provenientedelventilador,con carbonatocálcico(caliza).Estacorrientedeaireconpartículasdecalizaseconduce porunatuberíaalciclóndelcircuitodecarbonatación.Enesteciclónsecapturala calizacontenidaenlacorrientedeaire,ydichacorrientesaledelciclónlimpiade caliza.Antesdeabandonarlainstalaciónlacorrienteatraviesaunfiltrodemangas Figura21.Instalacióndeexperimentación.
Descripciónteóricadelainstalación Descripciónteóricadelainstalación 55 paracapturarpartículasdecalizaqueelciclónnohayasidocapazdeseparar.La caliza,queesatrapadaporelciclón,sedepositaenlaloopsealcircuitode carbonatación.Estedispositivohacecircularlacalizaallechofluidocirculantedel circuitodecalcinación.Porotrolado,sepuedepasarunacorrientedeaire provenientedelventiladorporlaresistenciaeléctricadelainstalacióndondese calientaantesdeintroducirseenellechofluidocirculante.Estecalentadoreléctrico seempleaúnicamenteparaconservarlosnúmerosadimensionalesylosratiosde densidadentrelacorrientedeaireylaspartículasdecaliza. Acontinuaciónenellechofluidocirculantedecalcinaciónsemezclala corrientedeairecalentadaconlacalizaparaluegodirigirsealciclóndelcircuitode calcinación.Endondelacalizaseseparadelacorrientedeaire.Finalmente,lacaliza seconducealaloopsealdelcircuitodecalcinaciónparareenviarlaalcircuitode carbonataciónyasíreciclarlaparavolverlaausar,encuantoalaire,seexpulsadel ciclónlimpiodecaliza.Perocomoocurreenelcircuitodecarbonatación,antesde salirdelainstalaciónlacorrienteatraviesaunfiltrodemangas,endonde,se capturanlaspartículasdecalizaqueelciclónnohasidocapazdeatrapar.
Descripciónteóricadelainstalación 56 3.1.3Aplicacionesyventajas/desventajadelainstalación EstesistemadecapturadeCO2conunciclodecarbonatación‐calcinaciónse aplicaráprincipalmenteenplantasdegeneracióneléctricatantodecarbóncomode gas,talescomo,centralestérmicasconvencionalesycicloscombinados, respectivamente.Yademásenprocesosindustrialesdondeexistanhornos,calderas, etc…Laproblemáticaquetieneestesistemaesqueaunestáenfasedeinvestigación ydemostración.Perosillegaaserválidosedaríaungrancambiobeneficiosoenel mundodelacapturadeCO2porquetienemuybuenasventajasencomparacióncon elrestodeprocesosbasadosprincipalmenteenpre‐combustiónyenoxicombustión. Porotrolado,estainstalacióndeflujofrioenconcretopresentalassiguientes ventajasydesventajas: Ventajas: Elsorbenteelegido,lacalizaesbarata,disponibleyabundante. Costedeinversióndelainstalaciónnoesalto. Lainstalaciónnopresentaningúnriesgoalahoradeoperarconella,ya que,esunprocesoquenoesnadapeligroso. Elmantenimientodeestainstalaciónesmuysimple,aunqueellono quitaquenoseaimportante. Permiterealizardistintasconfiguracionesdecircuitosdentrodela propiainstalaciónparallevaracabodiferentesensayosyestudiosen ella. Desventajas: Existeunconsumoenergéticoenlaresistenciaeléctrica(calcinador). Seproduceunadegradacióndelsorbente,esmuyimportantecalcularla paraconocerlacantidaddesorbenteconlaqueseestáoperando.
Planteamientosdelosensayos 57 3.2Planteamientosdelosensayos Laproblemáticaquepresentaestainstalacióndecapturadedióxidode carbonoespodermedirdemaneracontinuaelcaudaldepartículasquecirculaporel ciclón.Actualmente,estamedidanoescontinuayaquenosepuedeconoceren cualquierinstante.Además,esintrusivadebidoaquesumedidaafectaal funcionamientodelainstalación.Resolverestoesnecesarioparaconoceryanalizar lafluidodinámicadelsistema.Partiendodequesepuedemedirdemaneracontinua ynointrusivatantolacaídadepresióncomolavelocidaddelacorrientealaentrada delciclón,ycomosehaexplicadoanteriormente,sesabequelacaídadepresiónes funcióndelageometríadelciclón,ladensidaddelgas,lavelocidaddelacorrientea laentradayelcaudaldepartículasrespectivamente: ∆,,,,,16 Enconsecuencia,siseconocecómoserelacionan,sepodrácalculardemanera continuaynointrusivaelcaudaldepartículasquecirculaporelciclón.Elobjetivo principaldeesteproyectoesobtenerunacorrelacióndelacaídadepresiónenlos ciclonesenfuncióndelacargadepartículasconlaquetrabajaylavelocidaddelflujoa laentradayladensidaddelacorriente.Paraobtenerdichacorrelación,senecesita realizarunestudioexperimental(enlainstalación)yteórico(conlosmodelos)delos ciclonestangencialesdelainstalación.Esimportantecomentartambiénquecomo ambosciclonesoperaránbajolasmismascondicionesdetrabajoytienenlamisma geometría,estacorrelacióndelacaídadepresiónseráválidaparalosdosciclones.De ahí,queunaconfiguraciónadecuadadelainstalaciónpararealizarlosensayoshaya sidosóloelcircuitodecarbonatación,esdecir,ellechofluidocirculanteseconectóala loopsealdelmismocircuito,talycomomuestranlasfiguras22y23(siguientepágina). Encuantoalafunción,quehadesarrolladoelciclóndedichocircuitohasidola desepararlacalizadelacorrientedeaireycircularlahacialaloopseal,mientrasquela corrientedeairesaledelainstalacióncasilimpiadepartículas,porelloantesde abandonarlainstalaciónestacorrienteatraviesaunfiltrodemangas,endondese atrapapartículaspequeñasqueelciclónnohapodidocapturar.Principalmente,esto esnecesarioporseguridad,paraquelosoperariosdedichainstalaciónnorespiren partículasdecaliza,además,paraquelanavedondeestálainstalaciónnoseensucie conlacaliza.
Resultadosexperimentales 64 Porotrolado,enlafigura24sehanrepresentadográficamentelosvaloresde lacaídadepresiónobtenidosconlosmodelosteóricosexaminadosyvalores experimentalesenfuncióndelavelocidaddeentradaenelciclón. Figura24.Representacióndelosresultadosdelacaídadepresiónenvacio. Enlafigura25seexponenlosvaloresdelacaídadepresióncalculadaconlos diferentesmodelosteóricosenvacío,respectodelosvaloresexperimentales obtenidos.Estarepresentacióngráficainformavisualmentemuchomejorquela anteriordelaaproximacióndecadavalorteóricodelacaídadepresiónobtenidocon surespectivovalorexperimentaldelensayo. 0 200 400 600 800 1000 1200 0 5 10 15 20 ΔP(Pa) Vcy(m/s) ∆pcyclone[Pa] ∆pcyclone[Pa] ∆pteocyclone[Pa] ∆pteoSherpherd cyclone[Pa] ∆pteoAlexander cyclone[Pa] ∆pteoCasal cyclone[Pa] ∆pteoCoker cyclone[Pa]
Resultadosexperimentales 65 Figura25.Representacióndelosresultadosdelacaídadepresiónenvacioteóricarespectodela experimental. Porlotanto,examinandolatablaanteriorsepuedecomprobarquelos resultadosdelajusterealizadoenesteapartado,aligualquelosdelosdemás modelosteóricos,tieneunaaproximaciónmuybuenaalosresultados experimentales.Porotrolado,observandolafigura25sepuedeapreciarcómounos modelosseajustanmejorqueotrosalosvaloresexperimentales,esdecir,suspuntos estánmáspróximosalarectax=y(ΔPcexp=ΔPcteo).Enconclusión,elajustedelacaída depresiónenvacíoquemejorseasemejaalosresultadosexperimentales,yquepor elloeselmásidóneopararepresentarlacaídadepresiónenlosciclonesdela instalacióncuandooperanenausenciadepartículas,eselajusteteóricoestudiadoen esteapartado(ΔPteocyclone).Otrosajustesquetienenunabuenaaproximaciónsonel modelodeCasalyelmodelodeAlexander. 0 200 400 600 800 1000 1200 0 200 400 600 800 ΔPteo(Pa) ΔPexp(Pa) ∆pcyclone[Pa] ∆pteocyclone[Pa] ∆pteoSherpherd cyclone[Pa] ∆pteoAlexander cyclone[Pa] ∆pteoCasalcyclone [Pa] ∆pteoCokercyclone [Pa]
Resultadosexperimentales 66 3.3.1.2Pruebasencargaconlavariacióndevelocidad Losresultadosexperimentalesqueseobtuvieronenlaspruebasencarga variandolavelocidadycargadepartículasfueronlossiguientes: Prueba/ Ensayo vzlecho [m/s] vzi [m/s] ∆Pcexp [mbar] ∆Pcexp [Pa] Tª [ºC] ρ [kg/m3] mp [kgsol]t[s] Gs [kgsol/s] Csi [kgsol/kggas]ζsexp C10,5834,3440,42142,11531,509 1,1590,010600,00020,0110,962 C21,2099,0141,461146,146 35,352 1,1440,316600,00530,1680,785 C32,22316,5723,945394,534 39,145 1,1311,374200,06871,2040,635 C42,30517,1823,554355,390 33,517 1,1510,554100,05540,9200,522 C50,5954,4330,42542,45133,472 1,1510,010600,00020,0110,937 C60,9687,2161,011101,100 35,139 1,1450,028300,00090,0370,847 C71,37410,2431,463146,317 37,951 1,1350,234300,00780,2200,614 C81,87413,9663,078307,825 40,287 1,1260,652150,04350,9070,700 C92,30917,2154,574457,367 41,218 1,1230,978100,09781,6610,686 C100,5584,1620,32932,85336,708 1,1390,010600,00020,0120,832 C111,1248,3781,300130,045 39,671 1,1290,066300,00220,0760,820 C121,2469,2861,760175,976 40,497 1,1260,140200,00700,2200,906 C131,2829,5591,943194,281 41,100 1,1230,178200,00890,2720,945 C141,68712,5762,554255,428 42,356 1,1190,552150,03680,8590,721 C152,15016,0252,976297,571 43,735 1,1141,020160,06381,1730,520 C160,5614,1810,31931,92239,101 1,1310,010600,00020,0120,807 C170,8346,2171,194119,365 42,234 1,1190,038600,00060,0301,378 C181,1258,3851,478147,844 43,918 1,1140,160600,00270,0940,943 Tabla8.Resultadosexperimentalesdelosensayosencargarealizadosenelciclóndecarbonatacióndela instalación.
Resultadosexperimentales 67 Unavezconocidalacargadepartículasencadaensayosecalculólacaídade presiónenelciclónapartirdelosdistintosmodelosteóricosencarga,loscualesse expusieronenlosapartados2.8y2.9.Estosvaloresteóricosdelacaídadepresióny lafunciónR2decadamodeloseexpresanenlatabla9yfigura26. Tabla9.ValoresdelafunciónR2losdiferentesmodelosteóricos analizados. Figura26.Representacióndelosresultadosteóricosdelacaídadepresiónencargarespectode losdatosexperimentales. ‐800 ‐600 ‐400 ‐200 0 200 400 600 800 1000 0 100 200 300 400 500 ΔPTeo ΔPExp ΔPTeo/ΔPExp ∆pBriggs[Pa] ∆pSmolik[Pa] ∆pBaskakov[Pa] ∆pGrieco[Pa] Prueba/ Ensayo ∆Pc Briggs [Pa] ∆Pc Smolik [Pa] ∆Pc Baskakov [Pa] ∆Pc Grieco [Pa] R2=0,932190,735330,907690,94805
Resultadosexperimentales 68 EvaluandotantoelvalordelafuncióndeR2comolagráficamostradaenla figura26,querelacionanlosresultadosteóricosconlosexperimentales,sepuede deducirqueelmodeloteóricoencargaquemejorseasemejaalosdatos experimentales,sinduda,eselmodelodeGrieco.Existenotrascorrelacionescomo ΔPcBriggslacualobtieneresultadosmuyfavorables.Tambiénsepuedeobservarquela correlaciónΔPcSmolikeslaquepeoresvalorescalcula. Porlotanto,lacorrelacióndelacaídadepresiónquesevaaanalizartiene comobaseelmodelodeGrieco,comosemuestraenlaecuación(19): ∆ 1 2 22 Dondetanto“B”como“n”y“m”sonlasconstantesquesecalcularanbajoun procedimientodeajustecomoseexpondráenelsiguienteapartado3.3.2.Encuanto alnúmerodeEulerenvacíodelciclónsetrabajaráconelvalorobtenidoporel ajusterealizadoenelapartadoanterior,yaquecomosehademostradoeselmás adecuado.Comoseobservaysepretendíaestacorrelaciónesfunción(entreotros datos)delacargadepartículasylavelocidaddelacorrientegaseosaenlosciclones delainstalación.
Resultadosexperimentales 69 3.3.2Seleccióndelacorrelaciónadecuadayajustedelamisma Enesteapartadoseexplicaráelajustedelacorrelacióninvestigada.Se comentarátantoeltipodecorrelacióncomoelmétododeajusteporelcualfue calculada.TodoestotambiénpodráverseenelarchivodeExcelenelcualse realizarontodosloscálculosnecesarios.Seguidamente,losresultadosdeesta correlaciónsecomparanconlosresultadosobtenidosdelosmodelosteóricosy experimentales.Parafinalizarseelegiráaquelajustedelacorrelaciónquese aproximemejorquelosmodelosteóricosalosdatosexperimentales. Lacorrelaciónexaminadafue: ∆1 2 23 Estacorrelaciónescomplejayaqueconstadetresincógnitascomoson“B”, “m”y“n”.Elcálculodeestacorrelación,cómosevaexplicar,serealizóbajola suposicióndequeelciclónoperaavelocidadconstanteconcargavariable,utilizando losdatosexperimentales(tabla8).Endichocriterioelvalordeζg(vacío)conelquese trabajófue4,004(ajusteexperimental). Paraobtenerelvalorde“B”,“m”y“n”serealizaronlossiguientespasos: 1) Losensayosrepresentadosenlastabla5seagruparonenparejas,aquellos queteníanunavelocidad(vzi)muysimilar.Deestaclasificaciónse obtuvieronlassiguientesparejasdeensayos:C1conC10,C4conC9yC11 conC18. 2) Elvalorde“n”seobtuvoapartirdelasparejasdeensayoselegidasenel pasoanterior.Paraello,seaplicóelsiguienteprocesodecálculoylas siguientesecuaciones: ∆12 ⁄24 ∆12 25 Dondeelsubíndice“i”representalosvaloresexperimentalesdeunensayo yel“j”losdelotroensayodelapareja.Apartirdeestesistemadedos ecuacionesydosincógnitassecalculóelvalorde“nij”.Deunadeellas,por ejemplo,ΔPisedespejóvmparasustituirloenlaotrayaplicando propiedadeslogarítmicasseobtuvo:
Resultadosexperimentales 70 log ∆1/2 ∆1/226 Unavezresueltoestoparacadapareja,setomócomovalorde“n”parala correlaciónelvalormediodelasumadetodaslas“n”delasparejasde ensayosanalizadas. 3) SerealizóelsiguientecambiodevariableB´=B∙vzim,paraposteriormente comoseexplicará,calcular“B”y“m”delacorrelaciónpormediodeun ajustelogarítmico.Paracadaensayodelatabla8secalculósuvalorde “B´”,suponiendoelvalorde“n”obtenidoenelpasoanterior.Despejando “B´”delaecuación23setieneque: ´∆1 2 27 Dondeelsubíndice“i”representacadaensayo. 4) Enesteúltimopasosecalculóelvalorde“B”y“m”apartiruna aproximaciónlogarítmica.Peroantes,comolosvaloresdeB´delosensayos sonnegativosycomosesabequenoexisteellogaritmodeunnúmero negativo,esnecesariorealizarelsiguientecambiodevariable: B´=B∙vzim→ ‐B´=‐B∙vzim a´=a∙vzimsiendo,a´=‐B´ya=‐B Aplicandovariaspropiedadeslogarítmicasseobtieneque: a´=a∙vzim→ln(a´)=ln(a)+m∙ln(vzi) Porlotanto,elajustelogarítmicoconsistióenrepresentargráficamentelos valoresdelln(a´)frentealosdelln(vzi)paracadaensayodelatabla5. Dichosvaloresseaproximaronporlarectay=A+qxysecalcularonlos valoresdeAyqconelprogramainformáticoExcel2007.Para seguidamente,obtener“B”y“m”conesteprocedimiento: y=ln(a´) A=ln(a)→ a=eA→B=‐a q=m
Resultadosexperimentales 71 Finalmente,losvaloresde“B”,“m”y“n”queselograronfueron‐2,7813E‐05, 5,701y‐0,74754respectivamente.Detalmanera,quelacorrelacióncalculadase definecomo: ∆ 1 24,004 2,781310 ,28 Unavezdefinidaestacorrelación,secalcularonsusvaloresdelacaídade presiónparacompararlosconlosresultadosexperimentalesyteóricosquesehan obtenidoanteriormente.Paraelegiraquelajustequemejorseaproximealosdatos experimentalesseutilizólafunciónestadísticaR2,ademásdelarepresentación gráficacomosepuedeverenlasiguientetablayfiguras: Prueba/ Ensayo ∆PcBriggs[Pa] ∆PcSmolik[Pa] ∆PcBaskakov[Pa] ∆PcGrieco[Pa] ΔPcCal1[Pa] R20,932190,735330,907690,948050,95863 Tabla10.ValoresdelafunciónR2delosmodelosteóricoseinvestigadosanalizados. Figura27.Representacióndelosresultadosdelacaídadepresiónencargaenfuncióndelamisma. ‐800 ‐600 ‐400 ‐200 0 200 400 600 800 1000 0,0 0,5 1,0 1,5 2,0 ΔP(Pa) C(kgsol/kggas) ΔPciclon ∆pBriggs[Pa] ∆pSmolik[Pa] ∆pBaskakov[Pa] ∆pcyclone[Pa] ∆pGrieco[Pa] ΔpCal1[Pa]
Resultadosexperimentales 72 Figura28.Representacióndelosresultadosdelacaídadepresiónencargaenfunción delavelocidad. Figura29.Representacióndelosresultadosteóricosdelacaídadepresiónencargarespectodelosdatos experimentales. ‐800 ‐600 ‐400 ‐200 0 200 400 600 800 1000 0 5 10 15 20 ΔP(Pa) Vcy(m/s) ΔPciclon ∆pBriggs[Pa] ∆pSmolik[Pa] ∆pBaskakov[Pa] ∆pcyclone[Pa] ∆pGrieco[Pa] ΔpCal1[Pa] ‐800 ‐600 ‐400 ‐200 0 200 400 600 800 1000 0 100 200 300 400 500 ΔPTeo ΔPExp ΔPTeo/ΔPExp ∆pBriggs[Pa] ∆pSmolik[Pa] ∆pBaskakov[Pa] ∆pGrieco[Pa] ΔpCal1[Pa]
Resultadosexperimentales 73 Examinandolosresultadosdetodaslascorrelacionesanalizadassepuede apreciarquelacorrelación,ΔPcCal1,tieneelvalordelafunciónR2máspróximoala unidaddetodasyportantoelmejordetodos.Dichovaloressuperioralmodelode Grieco,comosehacomentadoanteriormente,eselquemejoresresultadosobtiene detodoslosmodelosteóricos.Porconsiguiente,sepuedededuciratravésdeesta funciónestadística(R2),quelacorrelaciónΔPcCal1eslaquemejorseajustadetodasa losdatosadquiridosexperimentalmente. Ademásanalizandolasrepresentacionesgráficasunaaunasepuedeobservar quelaúltimadeellaseslaquemayorymejorinformaciónproporcionaparaconocer cuáldetodaslascorrelacioneseslaquemejorseaproximaalosdatosconseguidos experimentalmente.Deestaforma,enestaúltimagráficasepuedeapreciarcómo paravaloresbajosdelacaídadepresióntodaslascorrelacionesseajustanmuy correctamente.Peroparavaloresaltosyengeneral,lasquemejorseaproximanalos resultadosexperimentalesson:ΔPcCal1,ΔPcBriggsyΔPcGrieco.Ydetodasellasla correlaciónquetienemáspuntoscercanosalarectax=y,esdecirΔPcexp=ΔPcexp,es ΔPcCal1.Porlotanto,sedemuestragráficamentequeestacorrelacióneslaque obtieneresultadosdelacaídadepresiónenelciclónmáspróximosalosdatos experimentales. Endefinitivaparaterminarlaselecciónyelajustedelacorrelaciónmás adecuadaparalosciclonesdelainstalación,porlorazonadoanteriormente,se deducequeesΔPcCal1,esdecir,lacorrelacióninvestigada.
Referenciasbibliográficas 80
Análisisexperimentaldelainfluenciadelacargadesólidosenlacaídadepresiónenun ciclóntangencial 81
Análisisexperimentaldelainfluenciadelacargadesólidosenlacaídadepresiónenun ciclóntangencial 82