UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE CENTRO DE TECNOLOGIA PROGRAMA DE PÓS-GRADUAÇÃO EM ENGENHARIA MECÂNICA DISSERTAÇÃO DE MESTRADO DESENVOLVIMENTO E CARACTERIZAÇÃO DE SOLUÇÕES TERNÁRIAS ÁGUA-GLICEROL-PROPILENOGLICOL COMO FLUIDO REFRIGERANTE SECUNDÁRIO Pedro Samuel Gomes Medeiros Orientador: Prof. Dr. Cleiton Rubens Formiga Barbosa Co-orientador: Prof. Dr. Francisco de Assis Oliveira Fontes Natal, Dezembro de 2012 DESENVOLVIMENTO E CARACTERIZAÇÃO DE SOLUÇÕES TERNÁRIAS ÁGUA-GLICEROL-PROPILENOGLICOL COMO FLUIDO REFRIGERANTE SECUNDÁRIO Pedro Samuel Gomes Medeiros Natal, Dezembro de 2012 Pedro Samuel Gomes Medeiros DESENVOLVIMENTO E CARACTERIZAÇÃO DE SOLUÇÕES TERNÁRIAS ÁGUA-GLICEROL-PROPILENOGLICOL COMO FLUIDO REFRIGERANTE SECUNDÁRIO Este trabalho corresponde à dissertação de mestrado apresentada ao Programa de Pós- Graduação em Engenharia Mecânica da Universidade Federal do Rio Grande do Norte, como parte dos requisitos necessários para obtenção do título de mestre em Engenharia Mecânica. Banca Examinadora: ___________________________________________________ Prof. Dr. Cleiton Rubens Formiga Barbosa Orientador – Presidente da banca examinadora – PPGEM/UFRN ___________________________________________________ Prof. Dr. Francisco de Assis Oliveira Fontes Co-orientador – PPGEM/UFRN ______________________________________________________ Prof. Dr. Enio Pedone Bandarra Filho Examinador externo à Instituição – UFU ______________________________________________________ Prof. Dr. Kleiber Lima de Bessa Examinador externo ao programa – UFRN ______________________________________________________ Prof. Dr. José Ubiragi de Lima Mendes Examinador interno – PPGEM/UFRN DEDICATÓRIA À minha querida e amada mãe Fátima Gomes por tudo o que me deu e proporciona em toda a minha vida... AGRADECIMENTOS Primeiramente agradeço ao Senhor Deus por tudo que foi proporcionado em minha vida, dando-me sabedoria e discernimento nas minhas atitudes. Agradeço aos meus amados pais Vanildo Medeiros e Fátima Gomes pelo total amparo na minha vida e na minha formação acadêmica, além da enorme compreensão dos momentos que estive ausente em suas vidas. Agradeço a minha querida Ana Claudia Bello pelo seu carinho e incentivo moral, principalmente na fase final deste trabalho. Ao senhor Manoel Alves (Melado) pelo apoio espiritual e otimismo na concretização dos meus trabalhos. Agradeço ao meu orientador, Prof. Cleiton Barbosa, por ter aceitado o desafio de orientar minhas ideias não convencionais, pela ampla paciência e sua contribuição na produção desta pesquisa. Ao CNPq pela concessão da bolsa de mestrado. Ao Prof. George Marinho pelo incentivo à minha iniciação a pesquisa e estímulo para realização de uma pós-graduação. Ao Prof. Francisco Fontes pela co-orientação deste trabalho. Ao Prof. Osvaldo Chiavone pelas suas contribuições sobre os conhecimentos da engenharia química essenciais ao desenvolvendo deste trabalho. Aos chefes dos laboratórios de Refrigeração, Metrologia, Mecânica dos Fluidos e NUPEG, respectivamente os professores Cleiton Barbosa, Luiz Pedro de Araújo, Ubiragi Mendes e Osvaldo Chiavone, pela disposição dos seus equipamentos na realização dos ensaios experimentais. Agradeço também a presença dos professores Ubiragi Mendes, Kleiber Bessa e Enio Bandarra pela formação da banca examinadora. Aos meus colegas de pós-graduação que estiveram comigo e acompanhando as dificuldades do desenvolvimento da minha pesquisa, dando-me apoio e força. A todos os meus sinceros agradecimentos e o sentimento de gratidão, pela contribuição de cada um para o desenrolar do meu mestrado. MEDEIROS, P. S. G. Desenvolvimento e Caracterização de Soluções Ternárias Água- Glicerol-Propilenoglicol como Fluido Refrigerante Secundário. 2012. 88 p. Dissertação (Mestrado em Engenharia Mecânica) – Programa de Pós-graduação em Engenharia Mecânica, Universidade Federal do Rio Grande do Norte, Natal-RN, 2012. RESUMO Cada vez mais o mundo está adotando uma matriz energética limpa e sustentável, com o uso da agricultura para produção de agroenergia e combustíveis verdes, como bioetanol e biodiesel. A produção do biodiesel gera um coproduto, a glicerina, em que as usinas produtoras têm dificuldades com o destino do seu excedente. Várias pesquisas estão sendo desenvolvidas para nortear diferentes usos do glicerol (glicerina pura). O glicerol possui total solubilidade com a água e pode ser usado como aditivo anticongelante aplicado como fluido refrigerante secundário, em sistemas de refrigeração indireta e com termoacumulação. Também, o glicerol é uma matéria-prima alternativa na produção de propilenoglicol, um álcool de grande aplicabilidade industrial inclusive como anticongelante. Porém, o melhor álcool anticongelante é o etilenoglicol, um álcool tóxico derivado do petróleo. As soluções água-glicerol (AG) e água-propilenoglicol (AP) possuem propriedades termofísicas de qualidade inferior e desequilibradas se comparadas às soluções água-etilenoglicol (AE). Desta forma, esta pesquisa inovadora teve como propósito o desenvolvimento e a caracterização de soluções ternárias água-glicerol-propilenoglicol (AGP) como fluidos secundários, com propriedades termofísicas desejáveis e competitivas com as soluções água-etilenoglicol. Equações preditivas simplificadas foram usadas para prever o comportamento das soluções AGP, onde as seguintes propriedades termofísicas foram avaliadas e estimadas teoricamente: ponto de congelamento, massa específica, calor específico e condutividade térmica. As concentrações para definir o ponto de congelamento das soluções AGP foram definidas a partir da Lei de Raoult das propriedades coligativas. A análise matemática inicial mostrou que as soluções AGP possuem propriedades mais equilibradas que as soluções AG e AP e competitiva com a solução AE. A comprovação experimental das soluções AGP foi feita a partir de ensaios para verificar suas propriedades (massa específica, condutividade térmica e viscosidade dinâmica), comparando com as soluções de referência AG e AP. Os resultados experimentais comprovaram as expectativas iniciais e viabilidade técnica do novo fluido secundário ternário. A grande vantagem dos fluidos AGP é que são atóxicos e derivam de fontes renováveis. Palavras-chave: Refrigerante Secundário Alternativo; Termoacumulação; Propriedades Termofísicas; Equações Preditivas Simplificadas; Glicerina do Biodiesel. MEDEIROS, P. S. G. Development and Characterization of Ternary Solutions Glycerol- Propylene Glycol-Water like Secondary Coolant Fluid. 2012. 88 s. Dissertation (Master’s degree in Mechanical Engineering) – Postgraduate Program in Mechanical Engineering, Federal University of Rio Grande do Norte, Natal-RN, 2012. ABSTRACT Increasingly the world is adopting a sustainable and clean energy matrix, with the use of agriculture to produce agro-energy and green fuels, such as bioethanol and biodiesel. The production of biodiesel generates a byproduct, glycerin, wherein the producing plants have difficulties with the fate of its surplus. Several studies are being developed to guide different uses of glycerol (glycerin pure). Glycerol has complete solubility in water and it can be used as antifreeze additive applied as a secondary refrigerant in indirect refrigeration systems and thermal energy storage. In addition, the glycerol is an alternative feedstock in the production of propylene glycol, an alcohol with large industrial applicability including as antifreeze. However, the best alcohol antifreeze is ethylene glycol, a toxic alcohol from petroleum derivate. The glycerol-water solutions (GW) and propylene glycol- water (PW) possess thermophysical properties of inferior quality and unbalanced compared to ethylene glycol-water solutions (EW). Thus, this innovative research had the purpose to development and characterization of ternary solutions glycerol-propylene glycol-water (GPW) as a secondary coolant with desirable thermophysical properties and competitive with the ethylene glycol-water solutions. Simplified predictive equations were used to predict the behavior of the solutions GPW, which the following thermophysical properties were evaluated and estimated theoretically: freezing point, density, heat capacity and thermal conductivity. The concentrations to set the freezing point of GPW solutions were defined from Raoult's Law of colligative properties. Initial mathematical analysis showed that the GPW solutions have properties more balanced than GW and PW solutions, and it is competitive with EW solution. The experimental confirmation of the GPW solutions was made from testing to check their properties (density, thermal conductivity and dynamic viscosity) compared to reference GW and PW solutions. The experimental results confirmed initial expectations and technical feasibility of new ternary secondary coolant. The more advantage is that GPW fluids are nontoxic and are derived from renewable sources. Keywords: Alternative Secondary Coolant; Thermal Energy Storage; Thermophysical Properties; Simplified Predictive Equations, Glycerin from Biodiesel. LISTA DE FIGURAS Capítulo 1 Figura 1.1 Glicerina gerada na produção de Biodiesel entre 2005-2011 16 Capítulo 2 Figura 2.1 Sistemas de refrigeração indireta e direta 22 Figura 2.2 Sistema de refrigeração indireta com termoacumulação 23 Figura 2.3 Tela do SecCool com os fluidos e suas propriedades 26 Figura 2.4 Diagrama de fase líquido-sólido do sistema água-glicerol 33 Figura 2.5 Diagrama de fase líquido-sólido ternário água-metanol-glicerol 35 Capítulo 3 Figura 3.1 Síntese do desenvolvimento das concentrações da solução AGP 48 Figura 3.2 Método de previsão das propriedades termofísicas da solução AGP 49 Figura 3.3 Procedimento geral da realização da fase de experimentos 50 Figura 3.4 Densímetro Anton Paar DMA 4500 51 Figura 3.5 Acoplamento da seringa com a amostra na entrada do densímetro 51 Figura 3.6 Analisador de Propriedades Térmicas KD2 Pro da Decagon Devices 53 Figura 3.7 Sensor tipo agulha KS-1 para medição em líquidos 53 Figura 3.8 Banho térmico de bloco seco Presys T-35N 54 Figura 3.9 Vista de topo do banho térmico mostrando o fosso interno 55 Figura 3.10 Bloco de prova para acondicionamento do fluido 55 Figura 3.11 Vista de topo do banho térmico com a proteção do fluido 56 Figura 3.12 Viscosímetro Quimis Q86M21 57 Figura 3.13 Spin 01 acoplado ao eixo de medição do viscosímetro 57 Figura 3.14 Controle de temperatura por meio de banho-maria entre béqueres 58 Capítulo 4 Figura 4.1 Comparativo entre os PC real e teórico das soluções AG e AP 61 Figura 4.2 Comparativo da Massa específica das soluções AGP com AP e AG 63 Figura 4.3 Comparativo do Calor Específico das soluções AGP com AG e AP 64 Figura 4.4 Comparativo da Condutividade Térmica AGP com AG e AP 65 Figura 4.5 Comparativo concentração de água nas soluções binárias e ternárias 67 Figura 4.6 Comparativo do CEV para os fluidos binários e ternários 68 Figura 4.7 Difusividade Térmica das soluções binárias e ternárias 69 Capítulo 5 Figura 5.1 Comparativo Experimental da Massa Específica AG, AP e AGP 72 Figura 5.2 Validação experimental de condutividade térmica da solução AP 74 Figura 5.3 Validação experimental de condutividade térmica da solução AG 74 Figura 5.4 Comparativo Exp. de Cond. Termica para os fluidos AP, AG e AGP 75 Figura 5.5 Validação da viscosidade dinâmica para a solução AP 76 Figura 5.6 Validação da viscosidade dinâmica para a solução AG 77 Figura 5.7 Comparativo da visc. dinâmica exp. dos fluidos binários e ternários 78 LISTA DE TABELAS Capítulo 2 Tabela 2.1 Propriedades Físico-Químicas do Etilenoglicol 27 Tabela 2.2 Propriedades Físico-Químicas do Propilenoglicol 28 Tabela 2.3 Propriedades Físico-Químicas do Glicerol 29 Tabela 2.4 Erro das Eqs (2.16) e (2.17) para a solução AP 39 Tabela 2.5 Erro das Eqs (2.16) e (2.17) para a solução AG 39 Tabela 2.6 Erro da Eq (2.19) para a solução AP 41 Tabela 2.7 Erro da Eq (2.19) para a solução AG 41 Tabela 2.8 Erro da Eq (2.21) para a solução AP 42 Tabela 2.9 Erro da Eq (2.21) para a solução AG 43 Tabela 2.10 Erro da Eq (2.23) para a solução AP 44 Tabela 2.11 Erro da Eq (2.23) para a solução AG 44 Capítulo 4 Tabela 4.1 Concentração dos Constituintes AGP para PC teórico de -11,91 °C 61 Tabela 4.2 Soluções AG, AP e AGP para Análise Teórica e suas Concentrações 62 Capítulo 5 Tabela 5.1 Validação da Massa Específica da solução AP 71 Tabela 5.2 Validação da Massa Específica da solução AG 72 Apêndice Tabela A1 Propriedades Água-Propilenoglicol (64,5/35,5) - ASHRAE 88 Tabela A2 Propriedades Água-Glicerol (60/40) - MELINDER 88 Tabela A3 Propriedades Termofísicas do Glicerol 88 LISTA DE ABREVIATURAS E SIGLAS ANP Agência Nacional de Petróleo, Gás e Biocombustível AE Solução Água-Etilenoglicol AG Solução Água-Glicerol AGP Solução Água-Glicerol-Propilenoglicol AP Solução Água-Propilenoglicol CEV Calor Específico Volumétrico [MJ/m³.K] EG Etilenoglicol GL Glicerol med Médio NUPEG Núcleo de Pesquisa em Petróleo e Gás da Univ. Fed. do Rio Grande do Norte PC Ponto de congelamento [°C] PG Propilenoglicol UFRN Universidade Federal do Rio Grande do Norte LISTA DE SÍMBOLOS C Calor Específico [kJ/kg.K] F Grau de liberdade k Condutividade Térmica [W/m.K] m Massa [kg] n Número de moles [mol] N Número de espécies químicas P Número de fases; Pressão R Constante universal dos gases T Temperatura [°C] u Potencial químico V Volume [m³] v Volume específico [m³/kg] Xi Fração mássica da substância i Yi Fração molar da substância i Z Propriedade termodinâmica extensiva ordinária Quantidade molar parcial Símbolos Gregos α Difusividade Térmica [mm²/s] µ Viscosidade Dinâmica [cP] ν Viscosidade Cinemática [cSt] ρ Massa Específica [kg/m³] SUMÁRIO CAPÍTULO I 1. INTRODUÇÃO ........................................................................................................... 15 1.1 Problemática e Motivação ........................................................................................... 15 1.1.1 Aplicação de Alcoóis como Aditivo Anticongelante ................................................ 16 1.1.2 Hipótese Investigativa: formulação de um novo fluido secundário .......................... 17 1.2 Objetivos ..................................................................................................................... 18 1.2.1 Geral .......................................................................................................................... 18 1.2.2 Específicos ................................................................................................................. 18 1.3 Estrutura da Dissertação .............................................................................................. 19 CAPÍTULO II 2. REVISÃO BIBLIOGRÁFICA ................................................................................... 21 2.1 Sistema de Refrigeração Indireta e Termoacumulação ............................................... 21 2.2 Propriedades e Critérios de Seleção dos Fluidos Refrigerantes Secundários.............. 23 2.3 Banco de Dados de Fluidos Secundários: SecCool ..................................................... 25 2.4 Características e Propriedades do Etilenoglicol, Proprilenoglicol e Glicerol ............. 26 2.4.1 Etilenoglicol .............................................................................................................. 26 2.4.2 Propilenoglicol .......................................................................................................... 27 2.4.3 Glicerol ...................................................................................................................... 28 2.5 Soluções e Propriedades Coligativas ........................................................................... 29 2.5.1 Desenvolvimento Matemático da Lei de Raoult ....................................................... 31 2.6 Introdução à Termodinâmica Química ........................................................................ 32 2.6.1 Diagramas de Fases de Sistemas Binários e Ternários ............................................. 32 2.6.2 Propriedades Parciais das Soluções ........................................................................... 35 2.7 Equações Preditivas das Propriedades Termofísicas ................................................... 37 2.7.1 Equações para Massa Específica ............................................................................... 38 2.7.2 Equação para Calor Específico .................................................................................. 40 2.7.3 Equação para Condutividade Térmica ...................................................................... 41 2.7.4 Equação para Viscosidade Dinâmica ........................................................................ 43 CAPÍTULO III 3. METODOLOGIA ........................................................................................................ 47 3.1 Método de Previsão das Propriedades Termofísicas da Solução Ternária AGP ......... 47 3.2 Métodos Experimentais das Soluções Ternárias AGP ................................................ 49 3.2.1 Medição da Massa Específica ................................................................................... 50 3.2.2 Medição da Condutividade Térmica ......................................................................... 52 3.2.3 Medição da Viscosidade Dinâmica ........................................................................... 56 CAPÍTULO IV 4. DESENVOLVIMENTO PREDITIVO DA SOLUÇÃO TERNÁRIA ÁGUA- GLICEROL-PROPILENOGLICOL ......................................................................... 60 4.1 Concentrações dos Constituintes a partir da Lei de Raoult ......................................... 60 4.2 Previsão da Massa Específica ...................................................................................... 63 4.3 Previsão do Calor Específico ...................................................................................... 65 4.4 Previsão da Condutividade Térmica ............................................................................ 66 4.5 Comparativo Teórico das Propriedades Termofísicas entre os Fluidos Binários e Ternários .................................................................................................................... 66 4.5.1 Concentração de água nas soluções binárias e ternárias ........................................... 66 4.5.2 Calor Específico Volumétrico ................................................................................... 67 4.5.3 Difusividade Térmica ................................................................................................ 68 CAPÍTULO V 5. VERIFICAÇÃO EXPERIMENTAL DAS PROPRIEDADES TERMOFÍSICAS DAS SOLUÇÕES AGP ............................................................................................. 71 5.1 Verificação Experimental da Massa Específica .......................................................... 71 5.2 Verificação Experimental da Condutividade Térmica ................................................ 73 5.3 Verificação Experimental da Viscosidade Dinâmica .................................................. 76 CAPÍTULO VI 6. CONCLUSÕES ............................................................................................................. 80 6.1 Sugestões para Trabalhos Futuros ............................................................................... 81 REFERÊNCIAS .......................................................................................................... 83 APÊNDICE .................................................................................................................. 88 Capítulo 1 Introdução 15 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Capítulo 1 INTRODUÇÃO 1.1. Problemática e Motivação O Brasil é um modelo no processo de conversão de energia não renovável para a renovável. Detém uma das matrizes energéticas mais limpas do mundo por expandir o uso do bioetanol e biodiesel não somente como substitutos do petróleo, mas transformando esses produtos em commodities mundiais, produzidos de uma forma ambiental e socialmente correta (VECCHIA, 2010). Em 2004, o Governo Federal lança o Programa Nacional de Produção do Biodiesel. Entrado em vigor em janeiro de 2005, o governo junto com a Agência Nacional de Petróleo, Gás e Biocombustível (ANP) introduziu o biodiesel na matriz energética nacional, mas não na forma pura, e sim a partir de um blend com o óleo diesel mineral (ANP, 2012a). O biodiesel é obtido a partir da transesterificação de triglicerídeos com um álcool de cadeia curta, resultando em éster (biodiesel) e glicerina bruta, em presença de um catalisador básico (DEMIRBAS, 2008). A glicerina bruta é um coproduto do processo formado por uma mistura de glicerol e impurezas (resíduos de catalisador, gorduras não reagidas, água, alcoóis). A glicerina purificada e destilada, com índice de pureza maior que 95%, tem diversos usos industriais como matéria-prima. No Brasil, desde 2005 o biodiesel é produzido e, consequentemente, tem-se glicerina como co-produto. A partir de 2008, a produção de glicerina no mesmo ano atingiu o valor de 124.415 m³ (ANP, 2012b). Em 2010 foi antecipado o uso do blend B5, cujo óleo diesel contém 5% de biodiesel, o que implicou na produção anual de 256.884 m³ de glicerina; ou seja, um aumento de mais de 100% em menos de dois anos. Em 2011, a produção anual de glicerina chegou à marca de 273.353 m³, um aumento de 6,41% em relação ao ano anterior. Na figura 1.1 é mostrado o gráfico da produção de glicerina anual no período de 2005-2011. Estima-se que cerca de 1 kg de glicerina bruta é produzido para cada 10 kg de biodiesel puro (THOMPSON, HE, 2006). Com o rápido crescimento da produção mundial de biodiesel, nos últimos anos foi gerado um grande excedente de glicerol. Isso tem provocado o fechamento de várias outras plantas industriais de glicerol. 16 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 1.1 – Glicerina gerada na produção de Biodiesel entre 2005-2011 Fonte: ANP, 2012b (adaptado) 1.1.1. Aplicação de Alcoóis como Aditivo Anticongelante Muito se tem pensado a respeito do excesso da glicerina. Várias pesquisas e estudos tentam buscar rotas alternativas para aumentar o grau de aplicabilidade desse produto. O glicerol é um álcool que possui solubilidade infinita com água e alcoóis e levemente solúvel em éteres e ésteres (SOLOMONS, FRYHLE, 2007). Em 1925, Leonard B. Lane escreveu um artigo sobre a total solubilidade do glicerol na água e que a mesma provoca alterações no ponto de congelamento da solução: surge então a aplicação do glicerol como agente anticongelante (LANE, 1925). Dossat (2004) explica que certas aplicações industriais de refrigeração necessitam do uso de fluidos refrigerantes secundários para transferência de calor. Esses fluidos são a base de alcoóis ou sais em solução com água. O pesquisador sueco Ake Melinder (MELINDER, 2007) defendeu em sua tese de doutorado a caracterização e utilização de soluções aquosas como fluidos refrigerantes secundários, usando diferentes alcoóis como agentes anticongelantes, dentre eles o glicerol. Melinder avaliou as propriedades termofísicas de soluções água-álcool(genérico), e percebeu que a adição de álcool como anticongelante altera as propriedades da água de forma negativa a custa da redução do ponto de congelamento, para abaixo de 0 °C. Dentre os alcoóis 0,07 9,5 36,7 124 172 257 273 0 50 100 150 200 250 300 M ilh ar es d e m ³ 2005 2006 2007 2008 2009 2010 2011 17 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN analisados, o etilenoglicol destaca o seu grande poder crioscópico e de alterações razoáveis das propriedades termofísica, justificando o seu largo uso industrial (MELINDER, 2007). Fink (2003) descreve também a aplicação do etilenoglicol como anticongelante devido a sua eficiência térmica, mas destaca que esse álcool tem seu uso limitado por causa da sua toxidade. Dossat (2004) e Fink (2003) falam do uso do propilenoglicol, um álcool atóxico, como uma alternativa ao etilenoglicol. A grande vantagem do propilenoglicol sobre o etilenoglicol é que o primeiro pode ser produzido a partir da glicerina, uma fonte renovável (MOTA, PESTANA, 2011). O etilenoglicol somente é produzido a partir do óxido de etileno, produto de origem petroquímica, ou seja, não renovável (DOW, 2012a). Na engenharia térmica moderna, fontes de armazenamento térmico têm sido soluções promissoras para o processo de conversão de energia de outras fontes em energia térmica, para uso posterior (ASHRAE, 2008). Para aplicações da refrigeração, a termoacumulação vem sendo utilizada em sistemas de água gelada que aplicam anticongelantes, permitindo o menor ponto de congelamento possível para aumentar o ganho térmico. As diversas soluções água-anticongelante têm como requesitos possuir um alto calor específico e alta condutividade térmica; no qual o etilenoglicol garante essas necessidades. O propilenoglicol em solução garante um alto calor específico, entretanto, tem baixa condutividade térmica. Já a solução água-glicerol se destaca pela boa condutividade térmica, mas peca por ter baixo calor específico. 1.1.2. Hipótese Investigativa: formulação de um novo fluido secundário A princípio, um método de compensar essas diferenças entre as soluções água- propilenoglicol e água-glicerol seria criar uma solução aquosa ternária contendo glicerol e propilenoglicol como anticongelantes. Partindo de uma hipótese que o uso concomitante de glicerol e propilenoglicol em solução possam garantir propriedades equilibradas, uma solução AGP (água-glicerol-propileniglicol) deverá possuir características e propriedades desejáveis e competitivas em relação à solução água-etilenoglicol. A hipótese investigativa consiste no princípio da adição das propriedades parciais de cada constituinte em solução para compor suas propriedades termofísicas finais. Não há estudos ou evidências sobre a aplicação das soluções AGP como fluido secundário em refrigeração. Também não há diagramas de equilíbrio termodinâmico líquido-sólido da solução ternária e nem informações e dados a respeito do seu comportamento como fluido 18 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN térmico. Somente uma pesquisa pioneira de caráter científico, tecnológico e de inovação é capaz de tentar viabilizar a solução AGP como uma alternativa à solução água-etilenoglicol. Portanto, sabendo-se que o etilenoglicol é tóxico e deriva de fontes não renováveis, que o glicerol é um co-produto do biodiesel e um excedente no mercado, e que o propilenoglicol pode ser produzido a partir do glicerol, nada mais viável que a solução ternária AGP para realçar a aplicação de fontes renováveis e dar uma rota alternativa ao uso do glicerol excedente do biodiesel. 1.2. Objetivos 1.2.1. Geral Desenvolver e caracterizar soluções ternárias de água-glicerol-propilenoglicol, como fluido secundário em sistemas e processos de refrigeração indireta com termoacumulação, com potencial de substituição à solução água-etilenoglicol. 1.2.2. Específicos • Prognosticar as concentrações dos constituintes do sistema ternário AGP para o ponto de congelamento arbitrado de -15 °C; • Estimar as propriedades termofísicas mais relevantes em um fluido térmico a partir de equações preditivas simplificadas para compreender o comportamento da solução AGP em função da temperatura; • Comparar as propriedades termofísicas estimadas de soluções AGP com os fluidos de referência água-glicerol e água-propilenoglicol, como também com a solução água-etilenoglicol; • Mensurar propriedades termodinâmicas e de transporte (massa específica, condutividade térmica e viscosidade dinâmica) do sistema AGP, tomando como fluidos de referências soluções reais de água-glicerol e água-propilenoglicol. 19 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN 1.3. Estrutura da Dissertação A organização do presente documento está estruturada em seis capítulos. Neste capítulo (capítulo 1) são apresentados a problemática e a motivação da pesquisa, definindo os objetivos geral e específicos da investigação. No capítulo 2 está a revisão bibliográfica, onde há um conjunto de informações referentes à aplicação de fluidos refrigerantes aquosos em sistemas de refrigeração indireta com termoacumulação, assim como um estudo físico-químico e termodinâmico de soluções aquosas e suas propriedades. Equações preditivas são apresentadas, como também sua aplicabilidade em fluidos secundários aquosos. O capítulo 3 mostra os materiais e métodos utilizados na pesquisa para definição das soluções AGP, de forma teórica e experimental. O desenvolvimento preditivo da solução ternária está no capítulo 4, onde se tem uma prévia do comportamento da solução água-glicerol-propilenoglicol em diferentes concentrações e em função da temperatura. Também nesse capítulo é realizado um comparativo com as soluções binárias a base de glicerol, propilenoglicol e etilenoglicol. A verificação experimental das principais propriedades termofísicas do fluido AGP proposto é apresentada no capítulo 5, em que as novas soluções são comparadas com as soluções binárias AP e AG. Finalmente no capítulo 6 estão as conclusões da pesquisa e sugestões para trabalhos futuros. Capítulo 2 Revisão Bibliográfica 21 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Capítulo 2 REVISÃO BIBLIOGRÁFICA 2.1. Sistema de Refrigeração Indireta e Termoacumulação Um sistema de refrigeração consiste numa máquina térmica que absorve trabalho para promover o abaixamento da temperatura do meio a ser refrigerado, por meio da transferência de calor deste ao fluido de trabalho da máquina. Quando o próprio fluido refrigerante do sistema é o responsável pela absorção direta do calor do meio a ser refrigerado, diz-se que o sistema é de expansão direta, ou simplesmente refrigeração direta. Neste caso, o evaporador está em contato direto com o espaço ou meio que está sendo refrigerado (DOSSAT, 2004). Porém, em algumas situações a aplicação de sistemas de refrigeração por expansão direta pode se tornar inconveniente ou de viabilidade técnico-econômica limitada ou inviável (DOSSAT, 2004). Exemplos são processos que requerem alta carga térmica com grandes percursos de tubulações e vários pontos de abastecimento de frio. Nesses casos, há uma grande distância entre a unidade condensadora (compressor e condensador) e a unidade evaporadora (dispositivo de expansão e evaporador), que inviabiliza o uso do sistema direto devido às perdas de eficiência em todos os dispositivos. A solução para os problemas relatados está na aplicação de um sistema secundário de refrigeração. Denominado de sistema de refrigeração por expansão indireta, o frio é gerado pelo sistema primário e então transferido, por meio de um trocador de calor, ao fluido refrigerante secundário que está em estado líquido e não muda de fase (PRUZAESKY, et al, 2008). O fluido secundário resfriado pelo fluido primário (sistema por expansão direta) é bombeado através de uma tubulação apropriada para os trocadores de calor finais, onde deverá absorver calor do espaço ou meio a ser refrigerado, promovendo assim o objetivo final da refrigeração (DOSSAT, 2004). Na figura 2.1 é mostrado um comparativo entre os sistemas direto e indireto. A temperatura de trabalho de um fluido refrigerante secundário pode chegar até -25 °C, ou de acordo com a necessidade de projeto (ASHRAE, 2006). Temperaturas usuais estão por volta de -10 °C a -20 °C, que são temperaturas encontradas em freezers e câmaras frigoríficas de armazenamento de alimentos (DOSSAT, 2004). 22 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 2.1 – Sistemas de refrigeração indireta e direta Fonte: Pruzaesky, et al, 2008 O confinamento do fluido secundário em tanques para o seu resfriamento cria a possibilidade de realizar um armazenamento térmico de energia em baixa temperatura para uso posterior. Partindo do princípio que a temperatura do fluido está sempre próxima ao seu ponto de congelamento, há uma possibilidade de que o frio gerado pelo sistema primário seja termoacumulado pelo fluido secundário e usado somente quando houver necessidade, flexibilizando os horários de funcionamento do sistema primário. O armazenamento térmico de energia, ou termoacumulação, é uma tecnologia empregada no acúmulo de energia térmica a partir de um meio condutor e absorvedor de calor (CEC, 1996). A grande vantagem da termoacumulação está na facilidade e prontidão do sistema ao absorver ou fornecer energia térmica quando requisitado (ASHRAE, 2008). Na figura 2.2 é mostrado um esquema do sistema de refrigeração indireta com termoacumulação. Para uma termoacumulação eficiente é necessário que o fluido secundário tenha boas propriedades termofísicas equilibradas de absorção e transmissão de calor simultaneamente. Como o armazenamento térmico e acúmulo de energia são processos volumétricos, elas dependem da capacidade da substância termoacumuladora reter energia interna em seu volume e da capacidade de transferir essa energia a outro meio (INCROPERA, DEWITT, 2002). Relações entre essas propriedades podem representar o quanto que uma substância consegue reter e transportar energia no tempo, ou seja, a velocidade com que o calor é armazenado e difundido pelo material (ÇENGEL, 2007). 23 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 2.2 – Sistema de refrigeração indireta com termoacumulação Fonte: Pruzaesky, et al, 2008 Muitas outras propriedades estão envolvidas para determinar a qualidade de um fluido em transferir e armazenar calor. Sendo assim, é necessário analisar qualitativamente e quantitativamente todas as propriedades relevantes de um fluido térmico para obter a máxima eficiência no armazenamento térmico e transporte de energia. 2.2. Propriedades e Critérios de Seleção dos Fluidos Refrigerantes Secundários Muitas são as propriedades intrínsecas aos fluidos. O princípio utilizado na seleção de um fluido térmico parte da premissa de que um fluido deve armazenar e transferir o máximo de energia térmica, ter fácil escoamento, quimicamente estável e não deteriorar materiais em contato. Autores como Melinder (2007), Dinçer e Kanoglu (2010), Stoecker (1998) e Torreira (2002) elencam e explicam diversas propriedades para seleção de fluidos térmicos para transferência de calor. As seguintes propriedades são as de natureza termofísica, que ainda podem ser divididas em propriedades termodinâmicas e propriedades de transporte: Propriedades Termofísicas – Termodinâmicas: • Ponto de Congelamento • Massa Específica • Calor Específico 24 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN • Calor Específico Volumétrico • Coeficiente de Expansão Térmica • Tensão Superficial • Pressão de Vapor Propriedades Termofísicas – Transporte: • Condutividade Térmica • Difusividade Térmica • Viscosidade Dinâmica • Viscosidade Cinemática De forma geral, um fluido refrigerante secundário ideal deve ter, dentre as propriedades termofísicas mais relevantes, baixo ponto de congelamento (ter as menores temperaturas possíveis em estado líquido para absorver o máximo de calor); alta massa específica (para ocupar o menor volume possível); alto calor específico (maior armazenamento térmico); alta condutividade térmica (otimização da taxa de transferência de calor); baixa viscosidade (facilidade ao escoar) e baixa pressão de vapor (baixa volatilidade). Avaliando as propriedades listadas, percebe-se que a água é a única substância que atende a esses requisitos, sendo assim um fluido ideal. Possui ótimas propriedades termofísicas para transporte e armazenamento de calor, boas propriedades químicas, além de ser de fácil aquisição e baixo custo, dentre outras propriedades já conhecidas. Entretanto, a água pura tem a grande desvantagem de solidificar a temperatura de 0 °C, e a refrigeração trabalha em temperaturas abaixo de 0 °C para obter a máxima transferência de calor. Diante deste dilema, a melhor decisão aplicada ao caso é tentar abaixar o ponto de congelamento da água por meio de adição de agentes anticongelantes. Segundo Fink (2003), anticongelante é definido como um aditivo que, quando posto em solução com a água, promove o efeito de redução do ponto de congelamento da solução. Trata-se de substâncias orgânicas ou inorgânicas totalmente solúveis em água. São usados em sistemas térmicos onde há a necessidade de trabalhar com fluidos a temperatura abaixo do ponto de congelamento da água. A adição do anticongelante reduz a qualidade das propriedades termofísicas. Os sais, como o cloreto de sódio, são de fácil obtenção e manipulação e não diminuem de forma significativa as propriedades termofísicas; porém não são aplicados atualmente devido ao grande poder de corrosão. 25 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Os alcoóis não trazem consigo o problema da corrosão, no entanto são caros e alguns deles são tóxicos. Os alcoóis listados por Melinder (2007) como anticongelantes são: metanol, etanol, etilenoglicol, propilenoglicol e glicerol. Os dois primeiros são voláteis e podem causar problemas de combustão; os alcoóis glicóis (etilenoglicol e propilenoglicol) e o glicerol não são voláteis e funcionam bem como agentes anticongelantes. 2.3. Banco de Dados de Fluidos Secundários: SecCool SecCool é um aplicativo computacional que compila informações de propriedades termofísicas de diversos fluidos refrigerantes secundários. Foi desenvolvido pelo projeto DESIK, uma parceria entre instituições de pesquisa da Dinamarca e empresas privadas de refrigeração para aplicação de simulação e otimização de processos de refrigeração. O aplicativo reúne dados experimentais consolidados de diversas bases de pesquisas, como ASHRAE e Dow Chemical, e do pesquisador Ake Melinder (SECCOOL, 2007). Esses dados têm toda credibilidade e formalidade experimental, servindo como parâmetros de seleção e projeto de sistemas de refrigeração indireta, além de comparar as diferentes soluções anticongelantes. Os dados compilados servem também como base comparativa para fluidos não listados. Soluções binárias aquosas de sais e alcoóis estão listadas na biblioteca do aplicativo, informando suas respectivas propriedades termofísicas: ponto de congelamento, massa específica, calor específico, condutividade térmica, viscosidades dinâmica e cinemática. As propriedades dos fluidos estão em função da concentração do anticongelante e da temperatura. Os dados experimentais foram inseridos de forma discreta na biblioteca e correlações matemáticas polinomiais são aplicadas para formar uma continuidade das informações (SKOVRUP, 2005). Na figura 2.3 é mostrada a tela do aplicativo: 26 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 2.3 – Tela do SecCool com os fluidos e suas propriedades Fonte: Seccool, 2007 2.4. Características e Propriedades do Etilenoglicol, Propilenoglicol e Glicerol 2.4.1. Etilenoglicol Etilenoglicol (EG) é o nome comum dado ao álcool 1,2-etanodiol (SOLOMONS, FRYHLE, 2007). O termo glicol deriva dos dois grupos hidroxilas (-OH) presentes na sua estrutura molecular. É um álcool não volátil, relativamente tóxico, higroscópico, incolor, e praticamente sem odor. Tem completa miscibilidade com a água e líquidos orgânicos polares (MEGLOBAL, 2008). A toxidez do etilenoglicol está relacionada à ingestão do produto, onde a rápida absorção pelo sistema digestivo causa danos ao sistema nervoso, tanto em animais quanto em humanos, levando a morte (DOW, 2009). Não há riscos quanto ao contato com a pele e nem por inalação. Por possuir boas propriedades físico-químicas, o etilenoglicol é considerado o melhor aditivo em fluidos térmicos a base de água, além de ser totalmente solúvel em qualquer 27 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN quantidade. O fato de não ser volátil evita a liberação de vapores, que contribui para o sua aplicação mesmo em altas temperaturas (DOW, 2012b). Na tabela 2.1 são mostrados os dados e propriedades físico-químicas do etilenoglicol: Tabela 2.1 – Propriedades Físico-Químicas do Etilenoglicol Fórmula C2H6O2 Massa Molecular 62 kg/kmol Ponto de Ebulição 197 °C Ponto de Fusão -13,4 °C Massa Específica a 25 °C 1110 kg/m³ Calor Específico a 25 °C 2,428 kJ/kg.K Condutividade Térmica a 25 °C 0,2579 W/m.K Viscosidade Dinâmica a 25°C 16,9 cP Solubilidade em Água Infinita Fonte: Dow, 2012b 2.4.2. Propilenoglicol Propilenoglicol (PG) é o nome ordinário do álcool glicol 1,2-propanodiol (SOLOMONS, FRYHLE, 2007), que, assim como o etilenoglicol, possui duas hidroxilas na sua estrutura molecular. É um líquido não volátil, incolor, praticamente inodoro, solúvel em água, higroscópico e atóxico (DOW, 2012c). Sua obtenção pode ser feita a partir da hidratação do óxido de propileno nas indústrias petroquímicas (fonte não renovável) ou pela hidrogenólise do glicerol (fonte renovável). A hidrogenólise consiste na remoção do terceiro radical hidroxila do glicerol e sua posterior hidrogenação para formação do propilenoglicol (CHUN, et al, 2007). Esse processo também produz o 1,3-propanodiol e em menor escala pode produzir etilenoglicol, apesar das dificuldades do processo de reação (CHUN, et al, 2007). A grande aplicabilidade do propilenoglicol o transforma numa matéria-prima valiosa. A sua utilização vai desde a indústria cosmética e de alimentos, até na fabricação de elastômeros, polímeros e solventes. Destaca-se também a aplicação do propilenoglicol como aditivo anticongelante em fluidos térmicos, principalmente na refrigeração alimentícia por ser uma substância atóxica (DOW, 2012c). Na tabela 2.2 são mostradas as propriedades físico-químicas do propilenoglicol: 28 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Tabela 2.2 – Propriedades Físico-Químicas do Propilenoglicol Fórmula C3H8O2 Massa Molecular 76 kg/kmol Ponto de Ebulição 187,4 °C Ponto de Fusão ≈ -60 °C Massa Específica a 25 °C 1033 kg/m³ Calor Específico a 25 °C 2,510 kJ/kg.K Condutividade Térmica a 25 °C 0,2061 W/m.K Viscosidade Dinâmica a 25°C 48,6 cP Solubilidade em Água Infinita Fonte: Dow, 2012c 2.4.3. Glicerol O glicerol (GL) é o nome do composto orgânico 1,2,3-propanotriol, que se trata de um álcool trihidroxilado, ou seja, com três radicais hidroxilas na sua molécula. É atóxico, não volátil, inodoro, incolor, de gosto doce e altamente higroscópico (KNOTHE, GERPEN, KRAHL, 2005). Também conhecido como glicerina, o glicerol, diferentemente dos alcoóis glicóis, é encontrado de forma combinado com ácidos graxos formando as moléculas de triglicerídeos. A obtenção de glicerol é feita a partir de reações químicas em triglicerídeos. Um dos métodos é o clássico processo de saponificação das gorduras, em que os triglicerídeos reagem com íons metálicos produzindo sabão e glicerina. Outro método é a produção de ésteres, onde os triglicerídeos reagem com alcoóis, resultando também em glicerina. Desta forma, o glicerol é oriundo de fontes renováveis, já que os triglicerídeos estão naturalmente presentes em animais e vegetais. Há diversas aplicações do glicerol, como matéria-prima ou na formulação de produtos. É de grande importância na indústria alimentícia e cosmética, como também na fabricação de espumas, lubrificantes, polímeros e borrachas (KNOTHE , GERPEN, KRAHL, 2005). A aplicação do glicerol como aditivo anticongelante vem desde o início do século XX. Em 1925, Leonard Lane fez pesquisa sobre o uso do glicerol em fluidos refrigerantes aquosos e estabeleceu o diagrama de fases água-glicerol (LANE, 1925). Na tabela 2.3 são descritas as propriedades do glicerol: 29 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Tabela 2.3 – Propriedades Físico-Químicas do Glicerol Fórmula C3H8O3 Massa Molecular 92 kg/kmol Ponto de Ebulição 290 °C Ponto de Fusão 17 °C Massa Específica a 25 °C 1261 kg/m³ Calor Específico a 25 °C 2,428 kJ/kg.K Condutividade Térmica a 25 °C 0,2579 W/m.K Viscosidade Dinâmica a 25°C 993,4 cP Solubilidade em Água Infinita Fonte: Çengel, 2007; Knothe, Gerpen, Krahl, 2005 Na tabela acima, o ponto de fusão do glicerol é dado numa temperatura próxima a do meio ambiente, e bem acima em relação aos alcoóis glicóis. Entretanto, o glicerol é raramente visto no estado sólido, somente em situações microcontroladas ou em temperaturas muito baixas. González, Longinotti e Corti (2011) realizaram um estudo sobre a viscosidade do glicerol em baixas temperaturas, onde por volta de -80 °C se apresenta em estado vítreo. Em situações normais, o glicerol se apresenta como um líquido super-resfriado num estado metaestável (KNOTHE, GERPEN, KRAHL, 2005). 2.5. Soluções e Propriedades Coligativas Segundo Castellan (1986), solução é uma composição homogênea a partir de duas espécies químicas totalmente dispersas a nível molecular. Soluções de dois componentes são chamadas de binárias, e as de três são ternárias. O componente em maior quantidade na solução é o solvente, enquanto que os demais constituintes, podendo ser mais de um, são chamados de solutos (CASTELLAN, 1986). A concentração dos constituintes (solutos e solvente) pode ser dada em base mássica ou base molecular (FELDER, ROUSSEAU, 2005). A fração mássica do componente i na solução é dada pela razão da massa de i pela massa total da solução, representada na seguinte equação: (2.1) 30 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN A fração molar do componente i na solução é a razão entre a quantidade de moles de i e a quantidade de moles totais da solução: (2.2) A representação percentual em massa de i é 100Xi, e a percentagem molar de i é 100Yi. Para o estudo dos fluidos secundários a água sempre será o solvente, enquanto que os aditivos anticongelantes são os solutos. Para o caso dos sais, nunca se deve extrapolar o limite de solubilidade do sal na água, isto é, a quantidade adicionada tem que se dissolver por completo formando apenas uma solução monofásica (DOSSAT, 2004). As propriedades termofísicas da solução água-anticongelante diferem das propriedades do solvente puro, a água. Há sempre mudanças em algumas propriedades termodinâmicas, tais como pressão de vapor, ponto de ebulição e ponto de congelamento. Essas mudanças são devido à ação do soluto (agente anticongelante) e dependem exclusivamente da sua massa molecular e da sua concentração no solvente (CASTELLAN, 1986; FELDER, ROUSSEAU, 2005). Essas mudanças nas propriedades são denominadas propriedades coligativas da solução e se referem à modificação das propriedades a partir das seguintes condições (FELDER, ROUSSEAU, 2005; ATKINS, PAULA, 2006): • O soluto é uma substância não volátil e de alto ponto de ebulição; • O soluto não reage com o solvente e nem se dissocia, ou seja, permanece molecularmente inalterado e apenas dissolvido no solvente. Segundo Atkins e Paula (2006), todas as propriedades coligativas estão relacionadas com o abaixamento do potencial químico do solvente provocado pela presença do soluto. As propriedades coligativas são aplicadas idealmente às soluções diluídas, ou seja, soluções em que a concentração do soluto é muito pequena e tende a zero (FELDER, ROUSSEAU, 2005). Porém, as alterações das propriedades coligativas se aplicam também quando a solução não é diluída, desde que a concentração do soluto não seja grande suficiente para alterar as condições químicas da solução. A Lei de Raoult é a explicação formal para o acontecimento dos fenômenos causado pelo soluto, em que o principal fator é a diminuição da pressão de vapor do solvente. O abaixamento da pressão de vapor tem duas consequências importantes: a redução da 31 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN temperatura de congelamento e o aumento da temperatura de ebulição da solução, quando comparadas ao solvente puro (BRADY, RUSSELL, HOLUM, 2002). O abaixamento do ponto de congelamento é denominado de efeito crioscópico. 2.5.1. Desenvolvimento Matemático da Lei de Raoult As considerações iniciais para a determinação matemática da Lei de Raoult é que a solução seja diluída. Tomando inicialmente que o solvente (água) esteja em estado sólido (gelo), o seu potencial químico em estado sólido deve ser igual ao potencial químico da solução diluída, já que a presença do soluto não irá alterar a energia livre de Gibbs do sistema (CASTELLAN, 1986). Desta forma, o equilíbrio pode ser definido matematicamente como: (2.3) onde u é o potencial químico em função somente da temperatura e pressão para o sólido, e também em função da fração molar do soluto na solução. Na equação (2.3), T é a temperatura de equilíbrio, ou seja, a temperatura de solidificação da solução. Como a pressão é constante na solução, T é apenas função da fração molar Y. A relação entre o potencial químico da solução e o potencial químico do solvente puro em estado líquido é uma função da temperatura, pressão e fração molar do soluto, onde é descrita como: (2.4) onde u é o potencial químico da solução e u° é o potencial químico do solvente líquido puro (em função da temperatura, pressão e fração molar), relacionados com a constante universal dos gases. Fazendo uma relação matemática da equação (2.4) com a equação (2.3) e por meio da equação de Gibbs-Helmholtz (não demonstrada), os compêndios de química geral e físico- química estabelecem uma relação final para o abaixamento crioscópico (BRADY, RUSSELL, HOLUM, 2002; CASTELLAN, 1986; ATKINS, PAULA, 2006): 32 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN [°C] (2.5) onde, ΔTf é o abaixamento crioscópico, Tf é o ponto de fusão do solvente puro, ΔHfus é o calor de fusão molar do solvente, R é a constante universal dos gases, e YA é a fração molar do soluto. Desta forma, a redução do ponto de congelamento depende exclusivamente das características do solvente e da sua concentração. Para soluções não diluídas, a aplicação da equação (2.5) não retrata o abaixamento crioscópico com exatidão, porém serve como referência do efeito de um soluto no solvente. 2.6. Introdução à Termodinâmica Química 2.6.1. Diagramas de Fases de Sistemas Binários e Ternários Diagrama de fases é um registro gráfico que relaciona os efeitos de temperatura, pressão e composição dos componentes presentes numa mistura, mostrando suas várias fases em equilíbrio termodinâmico (WALLAS, 1985). Fase é a forma homogênea como uma substância se apresenta a partir da sua composição e do estado que se encontra (ATKINS, PAULA, 2006). Os tipos de fases que podem existir em condições particulares são características da natureza química dos componentes (WALLAS, 1985). O número de fases presentes é dado a partir da regra de fases de Gibbs. Para sistemas com mais de uma fase e com um ou mais componente, o número de variáveis independentes que especifica o estado termodinâmico é dado pela seguinte relação: (2.6) onde P representa o número de fases presentes, N é o número de espécies químicas, e F é o grau de liberdade do sistema, ou seja, a quantidade de propriedades termodinâmicas intensivas e independentes necessárias para caracterizar o sistema (SMITH, VAN NESS, ABBOTT, 2007). 33 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Para determinação da quantidade de anticongelante para reduzir o ponto de congelamento da água, deve-se recorrer ao diagrama de fases do sistema água-anticongelante. O diagrama informa o percentual necessário em massa ou volume de anticongelante para que a solução adquira o ponto de congelamento desejado, para manter sempre no estado líquido. Se apenas um anticongelante for usado para o sistema, tem-se dois componentes (a água e o anticongelante); então o sistema terá três graus de liberdade, visto que se tem apenas uma fase presente (fase líquida) e o número de espécies químicas é dois; portanto, para definir o estado do sistema binário basta ter dados referentes à pressão, temperatura e concentração do anticongelante. Para uso da solução na refrigeração indireta com termoacumulação, a pressão atuante no tanque termoacumulador é a pressão ambiente; assim, um diagrama de fases a pressão constante é suficiente para determinar as fases presentes e a concentração do anticongelante. Na figura 2.4 é mostrado o diagrama de fases líquido-sólido da solução água-glicerol, desenvolvido em 1925 por Leonard Lane: Figura 2.4 – Diagrama de fases líquido-sólido do sistema água-glicerol Fonte: Lane, 1925 Para obter melhores resultados nas propriedades termofísicas das soluções aquosas, é possível o uso de dois alcoóis como agentes anticongelantes. Feldman e Dalhstrom (1936) fizeram estudos sobre a utilização concomitante de metanol e glicerol como aditivos %Glicerol em massa T em p er a tu ra ( °C ) Curva de Ponto de Congelamento da solução água-glicerol 34 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN anticongelantes. A explicação do uso dos dois alcoóis está na busca de um maior equilíbrio das propriedades termofísicas com a menor quantidade de soluto na solução. Um sistema ternário tem um grau de liberdade maior que um sistema binário. Fazendo uso da regra de fases de Gibbs (equação 2.6) para a fase líquida (P = 1) no sistema ternário (N = 3), F será igual a quatro. Tomando a pressão como constante no diagrama de fases, as variáveis temperatura e concentração de dois dos três componentes são necessárias para caracterizar o sistema (SANTOS, 2006). No sistema binário, para uma dada temperatura, a fração mássica de um dos constituintes é dada diretamente pela fração mássica do outro constituinte, em vez que: (2.7) Para um sistema ternário, mesmo mantendo a pressão constante, é preciso representar o diagrama de fases de forma tridimensional, uma vez que, para caracterizar a solução é necessário fixar a fração mássica de dois componentes, ficando o terceiro a ser determinado pela seguinte relação: (2.8) Diagrama de fases ternário tem aspecto diferente do diagrama binário, em que sua representação gráfica é a projeção planar das fases num triângulo equilátero, sendo que cada vértice representa um componente puro e os lados representam as variações das frações de cada um dos componentes (SANTOS, 2006). No interior do triângulo, as fases podem ser expressas por pontos, linhas e regiões, de acordo com a miscibilidade dos componentes (WALLAS, 1985). Assim, os pesquisadores construíram um diagrama de fase para os constituintes da solução água-metanol-glicerol, onde está mostrado na figura 2.5: 35 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 2.5 – Diagrama de fase líquido-sólido ternário água-metanol-glicerol Fonte: Feldman, Dalhstrom, 1936 2.6.2. Propriedades Parciais das Soluções As misturas em geral faz parte de nossas vidas, onde quase tudo é uma combinação de uma variedade de substâncias químicas. Muitas dessas misturas são soluções homogêneas em que há um equilíbrio de fases entre os seus constituintes. Um sistema homogêneo é um sistema com propriedades uniformes em toda sua extensão, ou seja, suas propriedades, tal como a densidade, têm o mesmo valor de um ponto para outro, no sentido macroscópico (PRAUSNITZ, LICHTENTHALER, AZEVEDO, 1999). Um sistema fechado é aquele que não troca matéria com o seu ambiente, embora possa trocar energia. Assim, num sistema fechado, não submetido à reação química, o número de moles de cada componente da solução é constante. Essa restrição pode ser expressa como (PRAUSNITZ, LICHTENTHALER, AZEVEDO, 1999): (2.9) onde ni é o número de moles do iésimo componente e N é o número de componentes presentes. Glicerol (%massa) Água (%massa) Metanol (%massa) 36 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN As soluções de água-álcool(genérico) não sofrem reações químicas e podem ser considerados como um sistema fechado no interior do tanque termoacumulador. Sistemas multicomponentes passam por mudanças em suas propriedades causadas pelo processo de mistura, onde as composições se tornam variáveis essenciais, assim como a temperatura e a pressão (SMITH, VAN NESS, ABBOTT, 2007). A química está na maioria das vezes preocupada com as reações provenientes sob condições em que a temperatura e a pressão são as variáveis que influenciam o processo. Portanto, é útil dispor de um conjunto de propriedades que descrevem o efeito de uma alteração na concentração de diversas quantidades termodinâmicas, para as condições de temperatura e pressão constantes. Essas propriedades são referidas como as quantidades parciais molares (OTT, BOERIO-GOATES, 2000). A quantidade molar parcial é definida como: (2.10) onde Z pode ser qualquer propriedade termodinâmica extensiva (energia interna, entalpia, entropia, volume) e que a derivada é para temperatura e pressão constantes. Para melhor compreensão da propriedade molar parcial, Z é uma função das propriedades termodinâmicas, ou seja, Z = f(P, T, n1, n2, ...), de modo que a diferenciação total de Z pode ser escrita como: (2.11) para um processo a T e P constantes. Substituindo a notação da equação (2.10) na equação (2.11) resulta em: (2.12) Espera-se que as propriedades parciais molares é função da composição, da temperatura e da pressão. Para um sistema a temperatura e pressão constantes, as propriedades parciais molares são apenas funções de composição (OTT, BOERIO-GOATES, 2000). A equação (2.12) está em base molar, podendo ser reescrita em base mássica. Para isso, deve haver uma relação de soma permitindo o cálculo das propriedades de mistura a 37 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN partir de propriedades parciais, desempenhando um papel oposto a equação (2.12) que permite o cálculo das propriedades parciais a partir das propriedades de mistura (SMITH, VAN NESS, ABBOTT, 2007). Assim, a equação (2.12) em base mássica é expressa da seguinte forma: (2.13) Segundo Smith, Van Ness e Abbott (2007), o conceito de propriedade parcial implica que a propriedade da solução representa um todo, ou seja, a soma de suas partes representadas como propriedades parciais das espécies que constituem a solução. Os constituintes de uma solução estão intimamente misturados, e devido às interações moleculares não podem ter propriedades privadas voltadas para si próprios. Entretanto, as propriedades parciais definidas nas equações (2.12) e (2.13) têm todas as características de propriedades das espécies individuais como elas existem na solução. Desta forma, elas podem ser especificadas como valores de propriedades para as espécies individuais (SMITH, VAN NESS, ABBOTT, 2007). Em síntese, propriedades parciais, como propriedades de solução, são funções da composição e das propriedades individuais de seus constituintes. Outras propriedades pertinentes às soluções, como propriedades de transporte, são derivadas do princípio aplicado da equação (2.12), em que a soma da influência dos constituintes na solução em função de suas concentrações determinam o comportamento final da solução. 2.7. Equações Preditivas das Propriedades Termofísicas A partir da termodinâmica química, as equações preditivas de propriedades termofísicas derivam do princípio básico das propriedades parciais de soluções. As equações descritas a seguir seguem a equação (2.12) para base molar ou a equação (2.13) para base mássica. 38 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN 2.7.1. Equações para Massa Específica Massa específica é uma propriedade termodinâmica intensiva de um sistema e representa a massa ocupada por uma unidade de volume: [kg/m³] (2.14) O inverso da massa específica é o volume específico, propriedade largamente aplicada na termodinâmica: [m³/kg] (2.15) Para determinar a massa específica de uma solução, experimentos em laboratório com equipamentos adequados para tal fim devem ser usados. Muitas soluções de uso geral foram caracterizadas e publicadas em manuais e livros especializados (FELDER, ROUSSEAU, 2005). Porém, para o estudo de uma nova solução ou de uma solução comum ainda não caracterizada, pode-se estimar o um valor médio de sua massa específica a partir de uma relação entre as frações mássicas dos componentes, Xi, e as massas específicas dos componentes puros, ρi, (FELDER, ROUSSEAU, 2005), por meio de dois métodos que seguem o modelo da equação (2.13). O primeiro método consiste em admitir a aditividade dos volumes, ou seja, o volume final da solução será igual ao somatório dos volumes dos componentes. Isto implica que o volume específico da solução é igual à média ponderada dos volumes específicos de cada componente relacionados à sua fração mássica na solução: [m³/kg] (2.16) onde a massa específica é dada pela equação (2.15) aplicada à equação (2.16). O segundo método é uma média ponderada das massas específicas dos componentes puros, multiplicando cada uma pela sua fração mássica: [kg/m³] (2.17) 39 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN A diferença entre os métodos está na forma em que a solução se comporta de acordo com a natureza química de seus componentes. Aconselha-se utilizar ambos os métodos como possíveis estimativas para a massa específica da solução, pois, como são equações gerais, não há uma margem de erro definida. Recomenda-se verificar experimentalmente qual método se encaixa melhor para prever o comportamento da solução (FELDER, ROUSSEAU, 2005). Para as soluções água-propilenoglicol (AP) e água-glicerol (AG), os dados das tabelas 2.4 e 2.5 mostram o comportamento dessas soluções à luz das equações (2.16) e (2.17), comparando os valores calculados com a base de dados do Seccool (2007). O cálculo foi feito variando a concentração do soluto numa temperatura constante de 25 °C. Os dados da massa específica da água, propilenoglicol e glicerol são, respectivamente nesta temperatura, 997,2 kg/m³, 1033 kg/m³, e 1261 kg/m³: Tabela 2.4 – Erro das Eqs (2.16) e (2.17) para a solução AP Conc% em Massa ρ Seccool (ASHRAE) [kg/m³] ρ Equação (2.16) [kg/m³] Erro% Eq (2.16) ρ Equação (2.17) [kg/m³] Erro% Eq (2.17) 20% 1016,4 1004,2 -1,2% 1004,4 -1,2% 30% 1025,6 1007,7 -1,8% 1007,9 -1,8% 40% 1033,4 1011,2 -2,2% 1011,5 -2,2% Tabela 2.5 – Erro das Eqs (2.16) e (2.17) para a solução AG Conc% em Massa ρ Seccool (Melinder) [kg/m³] ρ Equação (2.16) [kg/m³] Erro% Eq (2.16) ρ Equação (2.17) [kg/m³] Erro% Eq (2.17) 20% 1044,0 1040,7 -0,3% 1050,0 0,6% 30% 1069,5 1064,0 -0,5% 1076,3 0,6% 40% 1097,0 1088,3 -0,8% 1102,7 0,5% Percebe-se que as equações (2.16) e (2.17) para a previsão da massa específica da solução água-propilenoglicol têm valores semelhantes e inferiores aos do Seccool; enquanto que para a solução água-glicerol, a equação (2.16) fornece valores menores e a equação (2.17) valores maiores que os dados do Seccool. Em ambos os casos, o erro não foi superior a ±2,5%, sendo assim as equações preditivas têm capacidade de prever com boa exatidão a massa específica de soluções aquosas. 40 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN 2.7.2. Equação para Calor Específico Calor específico é uma propriedade termodinâmica intensiva que se refere à absorção e inércia térmica, como também é relacionado à energia interna de uma substância. Para sistemas gasosos, por serem compressíveis e dependentes do processo, o calor específico assume as formas de calor específico a volume constante ou de calor específico a pressão constante. Líquidos e sólidos são substâncias incompressíveis e há apenas um calor específico a ser aplicado. Para misturas de gases ideais, Himmelblau (1989) afirma que o calor específico molar a pressão constante é o somatório dos calores específicos molares dos componentes multiplicados pela fração molar: [kJ/kmol.K] (2.18) Felder e Rousseau (2005) afirmam que a entalpia total de um sistema é a soma das entalpias dos componentes puros. Desta forma, a entalpia associada ao processo de mistura dos componentes pode ser desprezada e se torna uma aproximação razoável para mistura de gases e líquidos similares, como também para soluções diluídas. Portanto, a equação para estimar o calor específico de uma solução líquida assume o formato da equação (2.18) aplicada a uma base mássica: [kJ/kg.K] (2.19) onde Xi é a fração mássica do componente i, Ci é o calor específico do componente i, e Cmed é o calor específico médio da solução. Percebe-se que tanto a equação (2.18) quanto a (2.19) possuem a mesma caracterização da equação das propriedades parciais de soluções. Segundo Teja (1983), há diversas formulações para estimar o calor específico de uma solução, sendo as equações mais convencionais a (2.18) e a (2.19). A formulação da equação (2.19) para um sistema binário foi testada para mais de 215 soluções, num total de 1083 pontos extraídos experimentalmente. A equação (2.19) na forma binária apresenta um erro máximo de 12,5%, onde esse nível de desvio é atingindo quando a equação é aplicada a soluções não aquosas. Como teste da equação (2.19), as tabelas 2.6 e 2.7 trazem o erro da equação preditiva em relação aos dados do Seccool para as soluções água-propilenoglicol e água-glicerol. O 41 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN calor específico da água é 4,179 kJ/kg.K a 25 °C. Os calores específicos do propilenoglicol e glicerol estão respectivamente nas tabelas 2.2 e 2.3. Tabela 2.6 – Erro da Eq (2.19) para a solução AP Conc% em Massa C Seccool (ASHRAE) [kJ/kg.K] C Equação (2.19) [kJ/kg.K] Erro% Eq (2.19) 20% 3,990 3,845 -3,8% 30% 3,869 3,678 -5,2% 40% 3,725 3,511 -6,1% Tabela 2.7 – Erro da Eq (2.19) para a solução AG Conc% em Massa C Seccool (Melinder) [kJ/kg.K] C Equação (2.19) [kJ/kg.K] Erro% Eq (2.19) 20% 3,856 3,826 -0,8% 30% 3,670 3,650 -0,6% 40% 3,471 3,474 0,1% A equação (2.19), para ambas as soluções, previu valores inferiores aos dados reais. Para a solução AP o erro da equação chegou a -6%; enquanto para a solução AG o erro não ultrapassou a casa do 1%. Portanto, mesmo com o erro de -6% para a solução água- propilenoglicol, o erro geral da equação (2.19) para as soluções aquosas foi inferior aos ±12,5% permitidos. 2.7.3. Equação para Condutividade Térmica Condutividade térmica é uma propriedade termofísica de transporte de calor que quantifica a habilidade de uma substância transferir energia térmica de um ponto a outro. Modelos para condutividade térmica têm sido determinados como função da composição por vários pesquisadores. As equações preditivas têm vários formatos matemáticos devido à complexidade do fenômeno envolvido e não há um consenso a respeito da utilização da composição em base molar ou mássica (GAITONDE, DESHPANDE, SUKHATME, 1978). As correlações matemáticas, em geral, se aplicam a apenas dois 42 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN componentes devido aos pequenos desvios das situações reais (POLING, PRAUSNITZ, O’CONNELL, 2004). Diferentes relações empíricas têm sido propostos para o cálculo da condutividade térmica de soluções líquidas binárias à pressão atmosférica. Provavelmente, o mais comum é a regra aditiva para a base mássica (LOSENICKY, 1968): [W/m.K] (2.20) O método da regra aditiva foi comparado com dados experimentais de soluções líquidas binárias, em que a equação (2.20) não se mostrou confiável para prever o comportamento das soluções. Em muitos casos foi necessária a inclusão de uma constante empírica para aumentar a confiabilidade do modelo preditivo (LOSENICKY, 1968). O modelo empírico simplificado proposto mais conhecido é a equação de Filippov: [W/m.K] (2.21) A grande vantagem da equação acima está na possibilidade de ajustar a constante 0,72 de acordo com os dados em sua validação no tipo de solução a ser estudada, além de ser composta em base mássica (POLING, PRAUSNITZ, O’CONNELL, 2004). A equação (2.21) aplica-se somente a soluções binárias líquidas, com uma incerteza esperada entre 4 e 8% (PERRY, GREEN, 2008). Aplicando a equação (2.21) as soluções água-propilenoglicol e água-glicerol, têm-se os seguintes erros preditivos mostrados nas tabelas 2.8 e 2.9. A condutividade térmica da água a 25 °C é 0,6076 W/m.K; e os dados do propilenoglicol e glicerol estão nas tabelas 2.2 e 2.3, respectivamente. Tabela 2.8 – Erro da Eq (2.21) para a solução AP Conc% em Massa k Seccool (ASHRAE) [W/m.K] k Equação (2.21) [W/m.K] Erro% Eq (2.21) 20% 0,5007 0,4810 -3,9% 30% 0,4523 0,4264 -5,7% 40% 0,4074 0,3776 -7,3% 43 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Tabela 2.9 – Erro da Eq (2.21) para a solução AG Conc% em Massa k Seccool (Melinder) [W/m.K] k Equação (2.21) [W/m.K] Erro% Eq (2.21) 20% 0,5291 0,5062 -4,3% 30% 0,4928 0,4625 -6,2% 40% 0,4549 0,4234 -6,9% A equação de Filippov previu valores inferiores aos dados das soluções presentes no banco de dados Seccool. O erro da equação (2.21) para as soluções AP e AG foi de até -7%, onde está dentro do limite aceitável de ±8%. 2.7.4. Equação para Viscosidade Dinâmica Viscosidade é uma propriedade física dos fluidos que caracteriza sua resistência ao escoamento. Pode-se definir viscosidade também como uma propriedade de transporte, onde a viscosidade dinâmica é a transferência molecular de quantidade de movimento linear entre as moléculas, a partir da aplicação de uma tensão de cisalhamento no fluido para compor o seu movimento, ou seja, o escoamento (BIRD, STEWART, LIGHTFOOT, 2002). Essa transferência de momento pode ser entendida também como o momento transferido pelo movimento macroscópico do fluido e que esse transporte convectivo de momento é proporcional à massa específica do fluido, formando a viscosidade cinemática (BIRD, STEWART, LIGHTFOOT, 2002). Para o estudo preditivo da viscosidade de soluções, apenas a viscosidade dinâmica é analisada matematicamente, pois a viscosidade cinemática pode ser definida a partir da viscosidade dinâmica pela seguinte relação geral: (2.22) onde ν é a viscosidade cinemática, µ a viscosidade dinâmica e ρ é a massa específica. A grande maioria dos métodos preditivos de viscosidade dinâmica assume valores das viscosidades dos componentes puros em relação à concentração destes componentes. Os vários modelos e métodos analíticos são para uma mistura de apenas dois componentes, e que 44 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN os mesmos devem apresentar semelhança química e viscosidades numericamente comparáveis (POLING, PRAUSNITZ, O’CONNELL, 2004). Por esse motivo, muitas das correlações se aplicam a óleos e hidrocarbonetos ou a sistemas puramente orgânicos. Para estimar a viscosidade de líquidos que não são hidrocarbonetos, a seguinte relação simplificada pode ser aplicada (PERRY, GREEN, 2008): (2.23) onde µmed é a viscosidade média da mistura, Yi é a fração molar do componente i na mistura, e µi é a viscosidade dinâmica do componente i. Essa equação é representativa para previsão geral da viscosidade dinâmica e erros por volta de 15% são admissíveis (PERRY, GREEN, 2008). Devido a grande diferença numérica de viscosidade entre água, propilenoglicol e glicerol, é provável que a equação (2.23) forneça valores no limite ou acima do erro admissível de 15%. Para isso, as tabelas 2.10 e 2.11 mostram o erro dos valores obtidos pela equação (2.23) frente aos dados do Seccool para as soluções água-propilenoglicol e água- glicerol. A viscosidade dinâmica da água a 25 °C é 0,89 cP; enquanto que do propilenoglicol e glicerol a mesma temperatura é de 48,6 cP e 993,4 cP, respectivamente. Tabela 2.10 – Erro da Eq (2.23) para a solução AP Conc% em Massa µ Seccool (ASHRAE) [cP] µ Equação (2.23) [cP] Erro% Eq (2.23) 20% 1,708 1,114 -34,8% 30% 2,508 1,288 -48,6% 40% 3,672 1,537 -58,1% Tabela 2.11 – Erro da Eq (2.23) para a solução AG Conc% em Massa µ Seccool (Melinder) [cP] µ Equação (2.23) [cP] Erro% Eq (2.23) 20% 1,521 1,235 -18,8% 30% 2,172 1,534 -29,4% 40% 3,174 2,001 -36,9% 45 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Sabendo que a equação (2.23) pode dar resultados com erros na ordem de 15% admissíveis, percebem-se pelas tabelas que a equação simplificada não se aplica as soluções analisadas. Os valores resultantes da equação preditiva são bem menores que os do Seccool e com erros maiores que os 15% até mesmo para baixas concentrações; para as altas concentrações, o erro atinge valores maiores que 50%. Como já foi dito, a equação (2.23) deve ser usada para componentes similares e sem grandes diferenças de viscosidade. Somente a viscosidade do propilenoglicol é 54 vezes maior que a da água e a do glicerol é mais de 1000 vezes maior que a da água. Essas discrepâncias de viscosidade dos alcoóis com a água também são motivos para o erro fora do valor admissível. Isto implica que a equação preditiva da viscosidade, que segue a tendência da equação (2.13) das propriedades parciais, não se aplica às soluções aquosas a base de propilenoglicol e glicerol. Capítulo 3 Metodologia 47 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Capítulo 3 METODOLOGIA Neste capítulo, definem-se os materiais e métodos utilizados no desenvolvido do novo fluido refrigerante ternário a base de água, glicerol e propilenoglicol. Duas fases da pesquisa foram estabelecidas: 1. Desenvolvimento do fluido ternário AGP aplicando equações preditivas para definir sua composição e propriedades termofísicas de forma teórica e preditiva. 2. Realização de experimentos para comprovar a viabilidade técnica do novo fluido ternário e analisar o comportamento de suas propriedades em comparação com fluidos binários. Todas as fases de análise dos dados, cálculo dos valores das propriedades e construção de gráficos e tabelas foram realizados em planilha eletrônica. 3.1. Método de Previsão das Propriedades Termofísicas da Solução Ternária AGP Definidas as equações preditivas válidas para o sistema água-álcool, a segunda fase inicia-se estimando a concentração mássica de cada constituinte da solução AGP a partir da Lei de Raoult, para o ponto de congelamento arbitrado de -15 °C. Essa temperatura está entre os valores adotados em sistemas de refrigeração indireta que variam entre -10 °C e -20 °C, conforme explicado na revisão bibliográfica. O valor médio de -15 °C implica inicialmente em menores erros de previsão nos cálculos das concentrações dos constituintes em relação à temperatura de -20 °C, como também, permite um uso prático mais flexível no sistema que a temperatura de -10 °C. Por se tratar de um sistema ternário, a determinação das frações mássicas do glicerol e proprilenoglicol serão definidas mantendo a fração mássica da água constante. Assim, é possível variar infinitos valores de cada aditivo anticongelante. Essa determinação será feita partindo do princípio que o abaixamento crioscópico depende apenas da fração molar do soluto e da sua massa molecular. Tratando-se da solução ternária AGP, o soluto será constituído por dois alcoóis, somando a participação de cada 48 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN álcool na solução, definindo assim a fração molar do soluto para o ponto de congelamento arbitrado. A representação esquemática desta etapa está descrito a seguir, na figura 3.1: Figura 3.1 – Síntese do desenvolvimento das concentrações da solução AGP Fonte: Autor A faixa de concentrações dos constituintes do fluido AGP parte inicialmente da solução água-glicerol para PG = -15 °C, em que a concentração em massa do glicerol é de 40%. O ponto de congelamento para essa mesma solução segundo a lei de Raoult é de -11,91 °C, para os mesmos 40% de glicerol. O propilenoglicol é adicionado à solução AG, mantendo o ponto de congelamento teórico constante. Assim é possível equilibrar a fração molar do soluto, onde a quantidade de glicerol presente na solução AGP é inversamente proporcional à quantidade de propilenoglicol, e vice-versa. Desta forma, o valor máximo de glicerol é de 40% em massa para 0% de propilenoglicol; e para 0% de glicerol, o valor máximo de propilenoglicol é de 35,5% para um PG teórico pela lei de Raoult de -11,91 °C. Definido a faixa de possíveis concentrações da solução AGP, foi priorizado as soluções com a maior concentração de glicerol, pois o mesmo é o produto resultante da purificação e destilação da glicerina bruta do biodiesel. As equações preditivas apresentadas têm grande confiabilidade apenas para sistemas binários. A previsão da condutividade térmica somente se aplica a dois componentes. Para estimar as propriedades do fluido AGP, foi considerada uma formulação pseudo binária, onde, admite-se inicialmente que uma solução binária de água-propilenoglicol seja representada como um único componente, pois esta solução é inteiramente homogênea e uniforme. O segundo componente a ser aplicado nas equações é o glicerol. A decisão de por o glicerol como o segundo componente deve-se ao fato de que os seus dados de propriedades termofísicas estão disponíveis na literatura em função da temperatura, como pode ser visto na tabela A3 do Apêndice. As propriedades termofísicas da solução água-propilenoglicol estão disponíveis em função da temperatura no aplicativo SecCool (2007), transcritas na tabela A1 do Apêndice. Ponto de Congelamento arbitrado de -15 °C Fração molar única do soluto composto por PG e GL Aplicação da Lei de Raoult Determinação das concentrações da solução AGP 49 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN A previsão da solução AGP deve ser feita em função da temperatura, pois somente assim é possível ter uma base comparativa com os fluidos binários, compreendendo o comportamento térmico das soluções no processo de refrigeração. Definidas as soluções AGP dentro da possível faixa de composição, foram selecionadas quatro formulações AGP para uma análise comparativa gráfica das propriedades massa específica, calor específico e condutividade térmica do sistema ternário com as soluções água-glicerol e água-propilenoglicol. Os dados tabulados das soluções binárias AP e AG estão no Apêndice, tabelas A1 e A2, respectivamente. Com base nos critérios expostos anteriormente, a concentração mássica mínima do propilenoglicol no SecCool para fornecer os dados das propriedades termofísicas da solução AP é de 11%, sendo essa a concentração mínima aplicada do proprilenoglicol na solução AGP. Para o cálculo preditivo da massa específica é empregado uma média aritmética das equações (2.16) e (2.17). Na análise dessas equações para as soluções aquosas com glicerol (seção 2.7.1) foi observado que uma possível média traga valores bem aproximados aos dados reais das soluções AG. Por esse motivo, a média será aplicada para previsão das soluções AGP. As propriedades combinadas calor específico volumétrico (CEV) e difusividade térmica dos fluidos AGP serão comparadas com as soluções água-glicerol, água- propilenoglicol e água-etilenoglicol, na temperatura constante de 25 °C. Na figura 3.2 é mostrado um diagrama de fluxo desta metodologia: Figura 3.2 – Método de previsão das propriedades termofísicas da solução AGP Fonte: Autor 3.2. Métodos Experimentais das Soluções Ternárias AGP Após o desenvolvimento da solução AGP e previsão de suas propriedades termofísicas, técnicas experimentais foram aplicadas para comprovar a eficiência real do novo Formulação de quatro soluções AGP Aplicação da solução AP como um único componente Previsão das Propriedades Termofísicas das soluções AGP em função da temperatura Comparativo AGP com as soluções AG, AP e AE 50 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN fluido e comparar nas mesmas condições com as soluções reais de referência água-glicerol e água-propilenoglicol. As propriedades analisadas experimentalmente foram massa específica, condutividade térmica e viscosidade dinâmica. Todos os experimentos foram realizados na Universidade Federal do Rio Grande do Norte (UFRN), no Núcleo de Tecnologia Industrial. As quatro soluções foram formuladas com amostras de 500 gramas. As massas dos constituintes foram calculadas a partir de cada concentração do fluido AGP. A água utilizada é do tipo destilada por processo de osmose reversa. O glicerol e o propilenoglicol utilizados apresentam teor de pureza de 99%. A montagem e mistura dos constituintes foram realizadas no Laboratório de Metrologia utilizando uma balança com resolução de 0,1 g, onde as quatro formulações foram pesadas e acondicionadas em garrafas do tipo pet. De maneira simplificada, a fase experimental pode ser representada da seguinte forma (figura 3.3): Figura 3.3 – Procedimento geral da realização da fase de experimentos Fonte: Autor 3.2.1. Medição da Massa Específica O processo de medição da massa específica foi efetuado no Núcleo de Pesquisa em Petróleo e Gás (NUPEG). A máquina utilizada é um densímetro modelo DMA 4500 da Anton Paar (figura 3.4) de medição instantânea da massa específica, onde possui controle de temperatura do fluido com exatidão de 0,03 °C, na faixa de medição de 0 a 90 °C. As Montagem das soluções AGP em balança Medição da massa específica em densímetro térmico Uso do banho térmico Medição da condutividade térmica Condicionamento em freezer Controle de temperatura por banho-maria Medição da viscosidade dinâmica 51 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN medições foram feitas em função da temperatura, no range de 5 a 25 °C, variação de 5 °C entre as medições. O tempo médio para medição é de aproximadamente 30 segundos. Figura 3.4 – Densímetro Anton Paar DMA 4500 Fonte: Autor A medição foi realizada com amostra de aproximadamente de 3 mL de cada fluido a partir de uma seringa conectada a entrada de medição. A seringa permanece acoplada durante todo o processo de medição (figura 3.5). Figura 3.5 – Acoplamento da seringa com a amostra na entrada do densímetro 52 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Para obter dados confiáveis, três medições foram feitas para cada temperatura, obtendo assim média aritmética e desvio-padrão dos valores. Os resultados experimentais foram validados a partir dos fluidos de referências AG e AP, confrontando com os seus dados do Seccool, disponível no Apêndice. Um comparativo dos resultados foi feito a partir de gráficos e tabelas que mostram o comportamento da massa específica dos fluidos AGP, AG e AP em função da temperatura. 3.2.2. Medição da Condutividade Térmica A medição da condutividade térmica foi realizada com o instrumento KD2 Pro da Decagon Devices (figura 3.6), disponível no Laboratório de Refrigeração. Trata-se de um sensor tipo agulha que deve ser imerso no fluido, medindo a condutividade térmica por diferenças de temperatura durante a medição. O sensor utilizado foi o KS-1, indicado para medição em líquido (figura 3.7). A faixa de medição do instrumento é de 0,02 a 2,00 W/m.K, com uma exatidão de 5%. O tempo de medição é de aproximadamente 3 minutos. A condutividade térmica dos fluidos AGP, AG e AP foi medida num range de temperatura entre -15 e 25 °C, variação de 5 °C. A cada ponto de temperatura três medições automáticas a cada 15 minutos foram realizadas. Para varrer essa faixa de temperatura, um banho térmico tipo bloco seco, localizado no Laboratório de Metrologia, foi usado. Trata-se do banho T-35N da Presys Instrumentos (figura 3.8), onde este equipamento possui uma faixa de operação de -35 a 125 °C, com exatidão no controle de temperatura de ± 0,1 °C. No interior do banho térmico há um fosso de 25,4 mm de diâmetro e 124 mm de profundidade (figura 3.9), onde é instalado um bloco de prova de ¾ pol. de diâmetro interno para acondicionamento do fluido e posterior resfriamento (figura 3.10). Na parte superior do banho, por ficar exposta, foram colocadas cinco camadas de fitas crepe que funcionou como isolante térmico e absorvedor de umidade do ar para não contaminar o fluido térmico (figura 3.11). O sensor KS-1 foi posto no centro do bloco de prova sobre camada protetora de fita. Para garantir total homogeneidade e estabilidade térmica no interior do fluido, para cada ponto de temperatura foi aguardado um tempo inicial de 15 minutos antes de realizar a primeira medição. 53 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 3.6 – Analisador de Propriedades Térmicas KD2 Pro da Decagon Devices Fonte: Autor Figura 3.7 – Sensor tipo agulha KS-1 para medição em líquidos Fonte: Autor 54 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Após as três medições em cada temperatura, cálculos de média aritmética e de desvio- padrão foram aplicados para ter o valor da condutividade térmica na respectiva temperatura. Com esses resultados, gráficos foram gerados primeiramente para validação do experimento, e depois para descrever o comportamento das soluções AGP frente às soluções de referências AG e AP. Figura 3.8 – Banho térmico de bloco seco Presys T-35N Fonte: Autor 55 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 3.9 – Vista de topo do banho térmico mostrando o fosso interno Fonte: Autor Figura 3.10 – Bloco de prova para acondicionamento do fluido Fonte: Autor 56 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 3.11 – Vista de topo do banho com a proteção do fluido Fonte: Autor 3.3.3. Medição da Viscosidade Dinâmica A medição da viscosidade dinâmica ocorreu a partir de um viscosímetro rotativo modelo Q86M21 da Quimis, localizado no Laboratório de Mecânica dos Fluidos (figura 3.12). A decisão da utilização deste equipamento foi apenas por questões estratégicas de tempo e disponibilidade. Este instrumento mede viscosidades na faixa de 100 a 600.000 cP com exatidão de ± 2%, num range de temperatura de 0 a 35 °C. O viscosímetro possui spins que são acoplados ao rotor para medição da viscosidade, de acordo com a faixa de medição da viscosidade adequada a cada spin. Para os menores valores de viscosidade, o spin 01 é utilizado na rotação de 60 rpm (figura 3.13). O ensaio de viscosidade dos fluidos secundários foi realizado em função da temperatura. Como não há banho térmico disponível para o viscosímetro, as garrafas que contêm os fluidos AGP, AG e AGP foram acondicionadas em freezer para executar o ensaio a temperatura controlada, aquecendo a partir de um banho-maria. O controle de temperatura foi feito colocando um béquer menor, que continha o fluido em análise, no interior de um béquer maior, onde água morna em pequenas quantidades preenchia o vazio formado (figura 3.14). 57 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 3.12 – Viscosímetro Quimis Q86M21 Fonte: Autor Figura 3.13 – Spin 01 acoplado ao eixo de medição do viscosímetro Fonte: Autor 58 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 3.14 – Controle de temperatura por meio de banho-maria entre béqueres Fonte: Autor Conforme a temperatura do ensaio aumenta, coloca-se água em temperaturas maiores para aquecer o fluido térmico dentro do béquer menor. As medidas de viscosidade foram feitas na faixa de temperatura de 5 a 25 °C, variação de 5 °C. Os valores de viscosidade medidos foram dentro de uma variação de ± 0,2 °C em cada ponto de temperatura, a fim de obter uma nuvem de dados sobre a viscosidade na temperatura desejada e menores erros de medição. Para cada fluido foram realizados cinco ensaios. Os valores finais de viscosidade foram obtidos fazendo uma média de cada nuvem de pontos para cada temperatura a partir de cada ensaio. Os valores de maior e menor média aritmética dentro dos cinco ensaios foram descartados, aproveitando apenas os três valores de média e realizando uma nova média aritmética para o valor final da viscosidade. Os dados da viscosidade dinâmica foram plotados em gráficos para validar o comportamento dos dados experimentais para as soluções de referências AG e AP, e depois comparados com os fluidos AGP para ter uma visão geral do comportamento das novas soluções ternárias. Capítulo 4 Desenvolvimento AGP 60 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM Capítulo 4 DESENVOLVIMENTO PREDITIVO DA SOLUÇÃO TERNÁRIA ÁGUA-GLICEROL-PROPILENOGLICOL 4.1. Concentrações dos Constituintes a partir da Lei de Raoult O sistema ternário água-glicerol-propilenoglicol não possui diagrama de fases do equilíbrio líquido-sólido, recorrendo a aproximações pela equação (2.5) da lei de Raoult, onde a fração molar do soluto pode ser um somatório das frações molares dos alcoóis que constitui o próprio soluto como um todo. As soluções AGP devem ter ponto de congelamento arbitrado em -15 °C, e como se trata de um sistema ternário, há infinitas combinações de concentrações de glicerol e propilenoglicol que resultam nesse ponto de congelamento. Para um ponto de congelamento de -15 °C estabelecido pela equação (2.5), as soluções água-glicerol e água-propilenoglicol requerem, respectivamente, uma concentração mássica de 45% e 42%; enquanto que as soluções reais requerem apenas 40% e 34% de concentração para os mesmos -15 °C de ponto de congelamento. Necessário se faz destacar uma análise do valor do ponto de congelamento para as soluções AG e AP quando se aplica suas concentrações para PC = -15 °C real na equação (2.5). Para ambas as soluções, o ponto de congelamento teórico resulta em valores na ordem de -11 °C. Na figura 4.1 é mostrado um comparativo entre os valores de ponto de congelamento real e pela lei de Raoult das soluções AG e AP, para concentrações respectivas de 40% e 34%. O gráfico mostra que a equação (2.5) fornece valores aproximados dos pontos de congelamento para as soluções AG e AP quando se aplica suas concentrações de 40% e 34%, as mesmas para um ponto de congelamento real de -15 °C. Isto ocorre porque a fração molar do propilenoglicol e glicerol são aproximados. O valor do PC teórico para água-glicerol é de -11,91 °C, para uma fração molar de 0,1154; enquanto que para água-propilenoglicol é de -11,22 °C, para uma fração molar de 0,1087. A diferença entre o PC da solução AG e AP é de 6,14% e a diferença do valor da fração molar entre AG e AP é de 6,16%, ou seja, não há diferenças percentuais entre as frações molares e as temperaturas do ponto de congelamento. Desta forma, para a equação de Raoult determinar as concentrações de glicerol e propilenoglicol das soluções AGP, considera-se PC = -11,91 °C para efeito de cálculo, visto que a solução AGP deve ter o maior teor possível de glicerol, ou seja, a fração molar definida 61 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM para a equação (2.5) é 0,1154. Também, essa fração molar da solução AG serve como um fator de segurança para que a solução AGP não congele em temperaturas maiores que -15 °C. Figura 4.1 – Comparativo entre os PC real e teórico das soluções AG e AP Fonte: autor Utilizando uma planilha eletrônica para efetuar os cálculos preditivos das concentrações, é possível determinar uma gama de concentrações dos constituintes para a solução AGP, tendo em vista um possível ponto de congelamento real de -15 °C. Na tabela 4.1 é apresentados algumas possíveis combinações de concentração mássica de água, glicerol e propilenoglicol: Tabela 4.1 – Concentração dos Constituintes AGP para PC teórico de -11,91 °C Conc% H2O Conc% GL Conc% PG 62,2 20,7 17,1 62,1 21,7 16,2 61,9 22,7 15,3 61,8 23,7 14,5 61,7 24,7 13,6 61,6 25,7 12,7 61,5 26,7 11,8 61,4 27,6 11,0 61,3 28,6 10,1 61,2 29,6 9,3 61,1 30,5 8,4 61,0 31,5 7,5 -15,62 -11,22 -15,39 -11,91 -18 -16 -14 -12 -10 -8 -6 -4 -2 0 PC Real PC Eq (2.5) P o n to d e C o n ge la m e n to [ °C ] AP AG 62 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM Todas as 16 combinações têm uma fração mássica de 0,1154 para os solutos glicerol e propilenoglicol, que implica num ponto de congelamento teórico de -11,91 °C. É possível que qualquer combinação resulte num ponto de congelamento real em torno de -15 °C. Para o cálculo preditivo das propriedades termofísicas das soluções AGP, a solução água-propilenoglicol é considerada como um único constituinte. Os dados da solução AP no SecCool são para uma concentração mínima de 11% de propilenoglicol em massa. Desta forma, as soluções AGP, que deve ter o máximo possível de glicerol, requerem um teor mínimo de 11% de propilenoglicol para sua análise preditiva, estabelecendo assim quatro combinações destacadas em negrito na tabela 4.1, que doravante serão os objetos de estudo. A referência para indicação das soluções AGP será a concentração mássica do propilenoglicol, onde a solução AGP 01 tem 13,6%; AGP 02 tem 12,7%; AGP 03 tem 11,8%; e AGP 04 tem 11%. Estabelecida as concentrações e suas respectivas denominações, pode ser feita a análise preditiva das propriedades termofísicas das quatro soluções AGP comparando com os dados reais das soluções água-glicerol e água-propilenoglicol. O ponto de congelamento teórico de -11,91 °C para as soluções AP requer uma fração molar de 0,1154 que implica numa concentração em massa de 35,5% de propilenoglicol. Segundo o SecCool, para esta concentração, o ponto de congelamento real da solução AP é - 16,8 °C. Apesar de ser maior que os -15 °C iniciais, essa concentração substitui a concentração anterior de 34% para efeito de análise comparativa entre os fluidos AGP e AG. Portanto, uma padronização na fração molar de 0,1154 é efetuada para todas as soluções binárias e ternárias, em que traduz num ponto de congelamento teórico de -11,91 °C pela lei de Raoult. Na tabela 4.2 são apresentadas as soluções a serem analisadas e suas respectivas concentrações mássicas percentuais: Tabela 4.2 – Soluções AG, AP e AGP para Análise Teórica e suas concentrações Sol/Comp H2O GL PG AG 60,0% 40,0% 0,0% AP 64,5% 0,0% 35,5% AGP01 61,7% 24,7% 13,6% AGP02 61,6% 25,7% 12,7% AGP03 61,5% 26,7% 11,8% AGP04 61,4% 27,6% 11,0% 63 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM 4.2. Previsão da Massa Específica Com as concentrações definidas das soluções AGP, AG e AP, a realização da previsão da massa específica torna-se uma tarefa simplificada, uma vez que as equações preditivas já foram validadas e têm-se os dados a respeito do glicerol e da solução água-propilenoglicol no Apêndice, respectivamente as tabelas A3 e A1. A previsão se dá com o uso dos dados da massa específica da solução AP para a concentração de propilenoglicol previsto nas soluções AGP, juntamente com os dados do glicerol. A análise é feita variando a temperatura entre 0 e 25 °C. Na figura 4.2 são mostradas as curvas comparativas da massa específica. O método utilizado é a média aritmética das equações (2.16) e (2.17), já que essa média irá equilibrar o valor final preditivo da solução AGP. Figura 4.2 – Comparativo da Massa específica das soluções AGP com AP e AG Fonte: Autor Avaliando o gráfico, percebe-se que as soluções AGP têm massa específica maior que a solução AP e menor que a AG, mesmo com a variação da temperatura. Um fluido secundário deve ter a maior massa específica possível, e neste requesito as soluções AGP estão em vantagem quando comparadas teoricamente com os dados reais da solução água- propilenoglicol. 1025,0 1035,0 1045,0 1055,0 1065,0 1075,0 1085,0 1095,0 1105,0 0 5 10 15 20 25 M as sa E sp e cí fi ca [ kg /m ³] Temperatura [°C] AG AP AGP01 AGP02 AGP03 AGP04 64 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM A solução AGP 04 possui maior massa específica que as outras soluções ternárias. Isso porque na AGP 04 há uma concentração maior de glicerol e o mínimo de propilenoglicol de 11%. Para as soluções AGP de forma geral, quanto maior o teor de glicerol, maior será sua massa específica em relação à solução água-propilenoglicol. 4.3. Previsão do Calor Específico A equação (2.19) para previsão do calor específico foi testada e aprovada para tal fim, conforme seção 2.7.2. Na figura 4.3 são representadas as curvas do calor específico teórico das quatro soluções AGP em comparação com os dados reais das soluções AP e AG: Figura 4.3 – Comparativo do Calor Específico das soluções AGP com dados reais AG e AP Fonte: Autor Percebe-se pela figura, que a solução AGP 01 possui calor específico maior que a AGP 04. Isso se deve ao fato da solução AGP 01 ter maior concentração de proprilenoglicol, já que este constituinte tem calor específico maior que o do glicerol. A solução AGP 04 tem o menor calor específico por ter o mínimo de propilenoglicol na sua formulação. Pode ser visto 3,300 3,400 3,500 3,600 3,700 3,800 0 5 10 15 20 25 C al o r Es p e cí fi co [ kJ /k g. K ] Temperatura [°C] AG AP AGP01 AGP02 AGP03 AGP04 65 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM também que os fluidos AGP trazem um equilíbrio nos valores do calor específico comparando com os dados das soluções AP e AG. 4.4. Previsão da Condutividade Térmica A previsão da condutividade térmica dos fluidos secundários é baseada na equação (2.21), onde a equação proposta por Filippov foi avaliada e validada para previsão em soluções binárias. Desta forma, na figura 4.4 pode ser visto o gráfico comparativo da condutividade térmica teórica dos fluidos AGP com dados reais, extraídos do SecCool, das soluções AP e AG: Figura 4.4 – Comparativo de Condutividade Térmica das soluções AGP com AG e AP Fonte: Autor Observando as curvas das soluções no gráfico, percebe que a condutividade térmica da solução água-propilenoglicol é inferior às demais, pois sabe-se que o propilenoglicol possui 0,4000 0,4100 0,4200 0,4300 0,4400 0,4500 0,4600 0 5 10 15 20 25 C o n d u ti vi d ad e T é rm ic a [W /m .K ] Temperatura [°C] AG AP AGP01 AGP02 AGP03 AGP04 66 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM baixa condutividade térmica quando comparado ao glicerol. O comportamento das soluções AGP mostra ser satisfatório, pois os resultados preditivos de condutividade térmica são maiores que os dados reais da solução AP. Dentre as soluções AGP analisadas, a solução AGP 04 possui a menor condutividade térmica prevista, mesmo tendo a menor concentração de propilenoglicol em solução. O cálculo preditivo foi de encontro à expectativa da condutividade térmica dos fluidos ternários, já que quanto maior a concentração de glicerol, maior a condutividade térmica. Observa-se também que não há diferenças entre as condutividades térmicas das soluções AGP. 4.5. Comparativo Teórico das Propriedades Termofísicas entre os Fluidos Binários e Ternários Definido os valores preditivos das quatro soluções AGP, é cabível um comparativo com as soluções água-etilenoglicol, água-glicerol, água-propilenoglicol, e um valor mediano das propriedades de AG e AP. Essa média entre as propriedades de AG e AP serve como uma meta a ser cumprida pelas soluções AGP no comparativo geral. Com uma temperatura padronizada de 25 °C, as propriedades combinadas calor específico volumétrico e difusividade térmica serão comparadas, permitindo um entendimento melhor das propriedades de armazenamento térmico (termoacumulação), transporte e difusão de energia térmica. 4.5.1. Concentração de água nas soluções binárias e ternárias Para o ponto de congelamento real de -15 °C, a solução água-etilenoglicol (AE) requer uma concentração em massa de etilenoglicol de 32% (SECCOOL, 2007). A concentração mássica de água nas soluções a base álcool indica basicamente a qualidade de suas propriedades termofísicas, onde quanto maior a concentração de álcool para o mesmo ponto de congelamento, maior é a redução das propriedades termofísicas. Comparando a concentração de água entre as soluções para PC = -15 °C, a solução AE se destaca como a que possui menor quantidade de soluto, como pode ser visto na figura 4.5: 67 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM Figura 4.5 – Comparativo concentração de água nas soluções binárias e ternárias Fonte: Autor Observa-se que as soluções AGP possuem baixo teor de água quando comparadas com as soluções AE e AP, sendo apenas um pouco maior que a solução AG. Porém, mesmo com uma concentração mais baixa, as soluções AGP exibem propriedades termofísicas satisfatórias e com valores equilibrados em relação às soluções AG e AP. 4.5.2. Calor Específico Volumétrico A análise da propriedade termodinâmica calor específico volumétrico (CEV) indica se um fluido é adequado para armazenar energia térmica, partindo da relação entre o seu calor específico e massa específica, onde CEV = ρC. Na termoacumulação é importante que um fluido secundário possua o maior valor numérico do CEV, pois assim irá acumular maior quantidade de energia em um volume menor de fluido. Sendo assim, um comparativo do CEV para as soluções AE, AG, AP e AGP deve ser efetuado a fim de verificar qual o fluido possui o maior efeito de acumulação térmica, dada a temperatura de 25 °C, como mostrado na figura 4.6. A solução água-propilenoglicol (AP) possui o maior CEV dentre os fluidos analisado. A explicação está no alto calor específico que o propilenoglicol possui em relação ao etilenoglicol e glicerol, onde este último apresenta o menor valor para o calor específico, que teve como consequência o menor CEV para a solução binária AG. 68% 60% 64,5% 61,7% 61,6% 61,5% 61,4% 50% 55% 60% 65% 70% AE AG AP AGP01 AGP02 AGP03 AGP04 C o n ce n tr aç ão M ás si ca % H 2 O 68 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM Figura 4.6 – Comparativo do CEV para os fluidos binários e ternários Fonte: Autor A perspectiva inicial das soluções AGP era de que o seu CEV pudesse ter um valor intermediário, numa média entre os valores do CEV das soluções AP e AG. No entanto, observa-se que a análise preditiva dos fluidos AGP traz um CEV maior que a média, onde a solução AGP 01 possui um CEV maior que o restante das soluções AGP. Neste caso, quanto maior for a concentração de propilenoglicol, mais a solução AGP irá armazenar energia térmica por unidade de volume. Porém, é importante destacar a diferença relativa entre o CEV maior e o menor entre as soluções. A diferença é de 2,6% entre as soluções AP e AG, que implica que na prática não há uma diferença significativa para a termoacumulação. 4.5.3. Difusividade Térmica A difusividade térmica avalia a capacidade de um fluido armazenar energia e depois transferi-la. Altos valores de difusividade representam altas taxas de condução de calor sem armazenamento térmico, enquanto que baixos valores significam retenção de energia e pouca transferência de energia. 3,823 3,808 3,907 3,857 3,894 3,889 3,884 3,878 3,700 3,750 3,800 3,850 3,900 C EV [ M J/ m ³. K ] 69 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM Na refrigeração indireta com termoacumulação, os fluidos secundários devem ter um equilíbrio na difusividade, já que é matematicamente expressa pela razão entre a condutividade térmica e o calor específico volumétrico. Significa dizer que um fluido secundário deve cumprir a missão de armazenar e depois transferir o melhor possível sua energia termoacumulada. Atualmente, a solução água-etilenoglicol é a mais utilizada por cumprir bem essa dupla missão. A solução água-propilenoglicol possui baixa difusividade, enquanto que a água- glicerol não fica tão atrás da solução a base de etilenoglicol. Para cobrir essa discrepância de valores entre as soluções AG e AP, os fluidos AGP dever ter um equilíbrio em relação às soluções binárias, como pode ser visto na figura 4.7: Figura 4.7 – Difusividade Térmica das soluções binárias e ternárias a base de EG, PG e GL Fonte: Autor A meta para a difusividade térmica das soluções AGP era ter valores medianos entre as soluções AG e AP. Observa-se que a difusividade estimada para os fluidos AGP são um pouco superior à solução AP, entretanto menor que meta da média AG-AP. Essa análise final permite identificar e comprovar que as propriedades termofísicas das soluções AGP são um pouco superior a solução AP e próximas a solução AG. A diferença relativa entre as soluções AE e AGP está na casa dos 6,7%; enquanto que entre as soluções AE e AP é de 10,1%. Assim os fluidos AGPs possuem difusividade competitiva em relação aos fluidos binários. 120,4 119,5 109,3 114,4 112,8 112,8 112,8 112,8 105,0 107,5 110,0 112,5 115,0 117,5 120,0 122,5 125,0 D if u si vi d ad e T é rm ic a [m m ²/ s] Capítulo 5 Verificação Experimental 71 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Capítulo 5 VERIFICAÇÃO EXPERIMENTAL DAS PROPRIEDADES TERMOFÍSICAS DAS SOLUÇÕES AGP A etapa experimental consiste em avaliar e comparar o comportamento real dos fluidos AGP com as soluções reais de água-glicerol e água-propilenoglicol nas mesmas condições. Os experimentos realizados foram das propriedades termofísicas termodinâmica (massa específica) e de transporte (condutividade térmica e viscosidade dinâmica). A análise é feita primeiramente validando cada experimento a partir dos dados das soluções AP e AG presentes no SecCool e disponíveis no Apêndice (tabelas A1 e A2), e depois comparando de forma conjunta as seis soluções AG, AP, AGP 01, AGP 02, AGP 03, e AGP 04 numa análise quantitativa-qualitativa. A concentração de glicerol na solução AG é de 40%, e a concentração de propilenoglicol na solução AP é de 35,5%. Esses dados de concentração das soluções binárias são os mesmos praticados na análise comparativa teórica com os fluidos AGP. 5.1. Verificação Experimental da Massa Específica O experimento compõe-se na medição da massa específica dos fluidos AGP e dos fluidos de referência AG e AP. A validação foi feita para verificar a exatidão do instrumento e da montagem das soluções. Nas tabelas 5.1 e 5.2 são mostrados os dados da massa específica experimental e do SecCool (tabelas do Apêndice), em função da temperatura, para as soluções água-propilenoglicol e água-glicerol respectivamente: Tabela 5.1 – Validação da Massa Específica da solução AP Temperatura [°C] ρ médio Exp [kg/m³] ρ desvpad ρ Tab A1 [kg/m³] erro% Exp/Tab 25 1025,7 0,3 1030,0 -0,42% 20 1028,3 0,3 1032,5 -0,41% 15 1030,8 0,3 1034,7 -0,38% 10 1033,1 0,3 1036,9 -0,36% 5 1035,4 0,3 1038,9 -0,34% 72 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Tabela 5.2 – Validação da Massa Específica da solução AG Temperatura [°C] ρ médio Exp [kg/m³] ρ desvpad ρ Tab A2 [kg/m³] erro% Exp/Tab 25 1097,2 0,4 1097,0 0,02% 20 1099,5 0,4 1099,3 0,02% 15 1101,7 0,4 1101,4 0,03% 10 1103,8 0,4 1103,5 0,02% 5 1105,7 0,4 1105,4 0,03% Os resultados experimentais da massa específica da solução água-propilenoglicol possuem valores inferiores aos dados da tabela A1. Porém, não há diferenças significativas entre os dados, pois a diferença é em média de 0,4% para menos dos resultados experimentais em relação aos dados reais (tabela A1). Para a solução água-glicerol, também não houve erros quando comparados os resultados experimentais com os dados da tabela A2. Praticamente os valores coincidem, tendo uma diferença mínima de apenas 0,02%. Assim, o instrumento e os procedimentos experimentais estão validados. Portanto, a massa específica das soluções AGP reais podem ser comparados com as soluções reais AG e AP. O gráfico da figura 5.1 mostra o comparativo: Figura 5.1 – Comparativo Experimental da Massa Específica AG, AP e AGP Fonte: Autor 1020,0 1030,0 1040,0 1050,0 1060,0 1070,0 1080,0 1090,0 1100,0 1110,0 5 10 15 20 25 M as sa E sp e cí fi ca [ kg /m ³] Temperatura [°C] AP AG AGP01 AGP02 AGP03 AGP04 73 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN As curvas do gráfico mostram que a massa específica das soluções AGP tendem a ter valores mais próximos da solução AG, sendo assim um ponto positivo. Dentre as soluções AGP, a AGP 04 possui a maior massa específica visto que possui maior teor de glicerol com relação às outras soluções AGP. Isto quer dizer que quanto maior a concentração do glicerol na solução AGP, maior será sua massa específica. Em vista disso do equilíbrio que uma solução ternária possa trazer em suas propriedades, todas as soluções AGP tiveram comportamento da massa específica condizente com a previsão. 5.2. Verificação Experimental da Condutividade Térmica Na previsão da condutividade térmica das soluções AGP, houve uma inversão quantitativa que fugiu da expectativa da condutividade térmica das soluções. A comprovação experimental tem por finalidade verificar se a condutividade térmica real das soluções AGP segue a mesma tendência qualitativa e quantitativa mostrada nos cálculos preditivos. A análise inicial consiste na validação do instrumento e dos procedimentos de medição. O experimento coletou dados de condutividade térmica numa faixa de temperatura de -15 a 25 °C, onde o valor de -15 °C é a temperatura limiar prevista do ponto de congelamento de todas as soluções. Nas figuras 5.2 e 5.3 são apresentadas as curvas dos resultados experimentais da condutividade térmica dos fluidos AP e AG, respectivamente, mostrando a faixa de incerteza do instrumento e confrontando com os dados das tabelas A1 e A2. Fazendo uma comparação conjunta do aspecto dos gráficos das soluções AP e AG, percebe-se algo no mínimo curioso. O instrumento de medição de condutividade térmica, o KD2 Pro, efetuou leituras de condutividade que não condiz com o comportamento real dos fluidos, nem no ponto de vista qualitativo. A condutividade térmica de qualquer solução aquosa é diretamente proporcional à temperatura, ou seja, quanto maior a temperatura, maior a condutividade térmica. Essa relação direta somente aconteceu na faixa de temperatura de -15 a 10 °C, onde, após esse valor de temperatura, os valores ao invés de aumentarem tiveram um decréscimo. Na faixa de -15 a 15 °C percebe-se que os valores experimentais estão dentro da incerteza do instrumento de ±5%. Vale lembrar ainda que o experimento foi realizado sob condições controladas, onde o fluido estava contido num tubo em um banho térmico, de temperatura 74 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN constante e estabilizada. Sendo assim, a validação só ocorre na faixa de temperatura de -15 a 15 °C, em relação aos dados extraídos do SecCool (tabelas A1 e A2). Figura 5.2 – Validação experimental de condutividade térmica da solução AP Fonte: Autor Figura 5.3 – Validação experimental de condutividade térmica da solução AG Fonte: Autor 0,340 0,360 0,380 0,400 0,420 0,440 0,460 -15 -10 -5 0 5 10 15 20 25 C o n d u ti vi d ad e T é rm ic a [W /m .K ] Temperatura [°C] K Exp K Tab A1 0,340 0,360 0,380 0,400 0,420 0,440 0,460 -15 -10 -5 0 5 10 15 20 25 C o n d u ti vi d ad e T é rm ic a [W /m .K ] Temperatura [°C] K Exp K Tab A2 75 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Entretanto, como os fluidos AP, AG e AGP foram testados sob as mesmas condições, é possível avaliar de forma qualitativa o desempenho da condutividade térmica em função da temperatura, na faixa aceitável e validada de -15 a 15 °C. O comportamento gráfico da nuvem de pontos, apesar de errôneo, mostra uma tendência particular para o experimento e passível de ser comparado, como pode ser visto na figura 5.4: Figura 5.4 – Comparativo Experimental de Cond. Térmica para os fluidos AP, AG e AGP Fonte: Autor Analisando qualitativamente os dados no gráfico, não é possível definir o comportamento das soluções AGP em relação às curvas preditivas da figura 4.4, onde há grandes flutuações nas medidas de condutividade térmica dos próprios fluidos AGP, tendo ainda uma sobreposição dos seus valores com os resultados da condutividade térmica da solução AG. Esse perfil não foi observado nos cálculos preditivos nem nos pontos da figura 4.4. Mesmo com validação parcial do experimento e com o comportamento anômalo experimental, as curvas mostram que as soluções AGP são viáveis e com valores consideráveis de condutividade térmica para um fluido refrigerante secundário, com o mesmo nível de competitividade da solução AG para a condutividade térmica. 0,38 0,39 0,40 0,41 0,42 0,43 0,44 0,45 0,46 -15 -10 -5 0 5 10 15 C o n d u ti vi d ad e T é rm ic a [W /m .K ] Temperatura [°C] AP AG AGP1 AGP2 AGP3 AGP4 76 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN 5.3. Verificação Experimental da Viscosidade Dinâmica A viscosidade dinâmica foi a única propriedade termofísica que não passou na validação de sua equação preditiva simplificada. Por se tratar de um modelo geral, a equação não foi capaz de descrever o comportamento das soluções aquosas a base de alcoóis. Porém, para efeito de comparação da viscosidade dos fluidos AGP com os fluidos de referência AG e AP, um estudo experimental foi realizado. O viscosímetro disponível, o Quimis Q86M21, não é adequado para medições de fluido com viscosidade abaixo de 50 cP. Entretanto, nas análises preliminares foi verificado que é possível representar uma curva qualitativa de viscosidade dinâmica em função da temperatura. Para efeito de validação, nas figuras 5.5 e 5.6 são mostradas as curvas da viscosidade experimental, com a incerteza do instrumento de ±2%, e os dados do SecCool para as soluções AP e AG, respectivamente os dados das tabelas A1 e A2: Figura 5.5 – Validação da viscosidade dinâmica para a solução AP Fonte: Autor 2 3 4 5 6 7 8 9 10 5 10 15 20 25 V is co si d ad e D in âm ic a [c P ] Temperatura [°C] µ Exp µ Tab A1 77 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 5.6 – Validação da viscosidade dinâmica para a solução AG Fonte: Autor Observando as curvas é perceptível que os dados experimentais são maiores que os do SecCool, mas mostram a mesma tendência de curva, sendo assim o experimento validado apenas para uma análise qualitativa. O glicerol tem viscosidade dinâmica mais de 1000 vezes maior que a da água a 25 °C. O propilenoglicol possui viscosidade 54 vezes maior que a da água, para a mesma temperatura. No entanto, as soluções AP e AG, para um ponto de congelamento de -15 °C possuem viscosidade dinâmica de valores bem aproximados ou até semelhantes. Como a soluções AGP contêm glicerol e propilenoglicol em proporções semelhantes aos fluidos binários, é possível afirmar que não deve haver diferenças significativas entre as viscosidades dos fluidos ternários e binários. Isto pode ser comprovado observando os pontos de viscosidade dinâmica em função da temperatura na figura 5.7: 2 3 4 5 6 7 8 9 10 5 10 15 20 25 V is co si d ad e D in âm ic a [c P ] Temperatura [°C] µ Exp µ Tab A2 78 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – PPGEM/UFRN Figura 5.7 – Comparativo da viscosidade dinâmica experimental dos fluidos binários e ternários Fonte: Autor Os pontos de viscosidade dos seis fluidos são praticamente coincidentes em toda a faixa de temperatura. A única restrição é para a solução AP, onde na temperatura de 5 °C apresenta o maior valor de viscosidade dentre os fluidos analisados. Porém, as pequenas diferenças numéricas existentes da viscosidade dinâmica dos fluidos AGP, AP e AG são desprezíveis, fazendo que as soluções ternárias tenham o mesmo nível de resistência ao escoamento que as soluções binárias a base de glicerol e propilenoglicol. 4,5 5 5,5 6 6,5 7 7,5 8 8,5 9 5 10 15 20 25 V is co si d ad e D in âm ic a [c P ] Temperatura [°C] AP AG AGP01 AGP02 AGP03 AGP04 Capítulo 6 Conclusões 80 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM Capítulo 6 CONCLUSÕES Com base na fundamentação teórica, nos cálculos preditivos desenvolvidos e nas propriedades termofísicas mensuradas experimentalmente sobre os fluidos AGP investigados, conclui-se que:  Foi confirmada a hipótese investigativa inicial de que soluções ternárias de água- glicerol-propilenoglicol (Fluidos AGP) possam exibir características e propriedades equilibradas em relação às propriedades das soluções água- propilenoglicol e água-glicerol;  Também foi confirmado que as soluções AGP são competitivas e capazes de ser uma alternativa às soluções água-etilenoglicol;  Os aditivos anticongelantes que compõe as soluções AGPs podem ser obtidos a partir da glicerina oriunda do biodiesel, constituindo numa rota alternativa ao uso deste coproduto excedente. Desta forma, as soluções ternárias contribuem para o desenvolvimento de uma matriz energética mais sustentável e potencialmente atrativa;  Os fluidos ternários desenvolvidos nesta pesquisa (AGP 01, AGP 02, AGP 03 e AGP 04) são tão eficientes quanto a solução água-propilenoglicol para a termoacumulação, e também melhores que a solução AP em difusão térmica. Assim, os fluidos AGP têm um perfil mais adequado para serem aplicados a sistemas de refrigeração indireta com termoacumulação (aplicação em chillers industriais);  O artifício de considerar a água-propilenoglicol como um único componente da solução ternária transformou a solução ternária num sistema pseudo binário. Este artifício possibilitou o uso de equações preditivas binárias ao sistema ternário AGP sem perdas na qualidade da previsão de suas propriedades;  Não existe uma única concentração ideal para um dado ponto de congelamento da solução água-glicerol-propilenoglicol. Por se tratar de um sistema com um grau de liberdade maior que os sistemas binários, existem outras combinações possíveis de 81 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM concentração dos solutos, cujos valores dependem das características desejáveis da solução e de sua respectiva aplicação;  As viscosidades dinâmicas dos fluidos AGP têm comportamento semelhante às soluções AP e AG, ou seja, a substituição de um fluido binário a base de glicerol ou propilenoglicol pelo AGP não acarreta em mudanças no sistema hidráulico de refrigeração.  As quatro soluções AGP desenvolvidas e selecionadas na pesquisa possuem propriedades termofísicas equilibradas e versáteis que garantem sua aplicabilidade tanto como fluido de transporte quanto para termoacumulação. 6.1. Sugestões para Trabalhos Futuros O aprofundamento na análise experimental do comportamento das soluções AGP é fundamental para assegurar a sua viabilidade e como uma possível alternativa à solução água- etilenoglicol. Determinados aspectos necessitam de uma investigação mais a fundo e mais confiável e, desta forma, como proposta de continuidade deste trabalho, sugere-se: • Fazer um estudo termodinâmico das soluções AGP para determinação exata do equilíbrio líquido-sólido, pois assim será possível indicar com exatidão as concentrações necessárias dos alcoóis para dado ponto de congelamento; • Realizar os ensaios experimentais em equipamentos de maior confiabilidade; • Comprovar experimentalmente a previsão do calor específico e da difusividade térmica das soluções AGP; • Estudar diferentes equações preditivas de viscosidade dinâmica viáveis para soluções aquosas; • Obter dados das propriedades termofísicas em toda a faixa de aplicação de temperatura dos fluidos secundários; • Gerar um banco de dados com os resultados experimentais do fluido AGP para facilitar o processo de seleção de fluidos refrigerantes secundários. Referências 83 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM Referências ANP – AGÊNCIA NACIONAL DE PETRÓLEO, GÁS E BIOCOMBUSTÍVEL. Anuário Estatístico 2011. Agência Nacional de Petróleo, Gás e Biocombustíveis – ANP. Disponível em: . Acesso em: 18 de setembro de 2012 (b). ANP – AGÊNCIA NACIONAL DE PETRÓLEO, GÁS E BIOCOMBUSTÍVEL. Biodiesel: introdução. Agência Nacional de Petróleo, Gás e Biocombustíveis – ANP. Disponível em: . Acesso em: 17 de setembro de 2012 (a). ASHRAE – AMERICAN SOCIETY OF HEATING, REFRIGERATING AND AIR- CONDITIONING ENGINEERS. 2006 ASHRAE Handbook Refrigeration. Atlanta-GA, USA: American Society of Heating, Refrigerating and Air-conditioning Engineers, 2006. ASHRAE – AMERICAN SOCIETY OF HEATING, REFRIGERATING AND AIR- CONDITIONING ENGINEERS. 2008 ASHRAE Handbook HVAC Systems and Equipment. Atlanta-GA, USA: American Society of Heating, Refrigerating and Air- conditioning Engineers, 2008. ATKINS, P.; PAULA, J. de. Atkins’ Physical Chemistry. 8th edition. New York-NY, USA: W. H. Freeman and Company, 2006. BIRD, R. B.; STEWART, W. E.; LIGHTFOOT, E. N. Transport Phenomena. 2nd edition. New York-NY, USA: John Wiley and Sons, 2002. BRADY, J. E.; RUSSELL, J. W.; HOLUM, J. R. Química: a matéria e suas transformações. 3ª edição. Volume 1. Rio de Janeiro, RJ: LTC, 2002. CASTELLAN, G. Fundamentos de Físico-Química. Rio de Janeiro, RJ: LTC, 1986. CEC – CALIFORNIA ENERGY COMMISSION. Source Energy and Environmental Impacts of Thermal Energy Storage. Sacramento-CA, USA: California Energy Commission, 1996. ÇENGEL, Y. A. Heat and Mass Transfer: a practical approach. 3rd edition. New York-NY, USA: McGraw-Hill, 2007. CHUN, H. Z.; et al. Chemoselective catalytic conversion of glycerol as a biorenewable source to valuable commodity chemicals. Chemical Society Reviews, v. 37, n. 1, p. 527- 549, nov. 2007. 84 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM DEMIRBAS, Ayhan. Biodiesel: a realistic fuel alternative for diesel engines. London, UK: Springer, 2008. DINÇER, I.; KANOGLU, M. Refrigeration Systems and Applications. 2nd edition. Chichester-WS, UK: John Wiley and Sons, 2010. DOSSAT, R. J. Princípios de Refrigeração. Curitiba, PR: Hemus, 2004. DOW. Dow Ethylene Oxide. Dow Chemical Company. Disponível em: . Acesso em: 15 de setembro de 2012 (a). DOW. Dow Ethylene Glycol: physical properties. Disponível em: . Acesso em: 03 de outubro de 2012 (b). DOW. Etilenoglicol: ficha de informações de segurança do produto químico. Guarujá, SP: Dow Chemical do Brasil, 2009. DOW. Propilenoglicol Dow Grau Industrial: ficha de dados técnicos. Disponível em: . Acesso em: 03 de outubro de 2012 (c). FELDER, R. M.; ROUSSEAU, R. W. Elementary Principles of Chemical Processes. 3rd edition. Hoboken-NJ, USA: John Wiley and Sons, 2005. FELDMAN, H. B.; DAHLSTROM, W. G. Freezing Points of the Ternary System Glycerol-Methanol-Water. Industrial and Engineering Chemical, november 1936, pp. 1316- 7, 1936. FINK, J. K. Oil Field Chemicals. Burlington-MA, USA: Gulf Professional Publishing, 2003. GAITONDE, U. N.; DESHPANDE, D. D.; SUKHATME, S. P. The Thermal Conductivity of Liquid Mixtures. Ind. Eng. Chem. Fundam., vol. 17, no. 4, 321-325, 1978. GONZÁLEZ, J. A. T.; LONGINOTTI, M. P.; CORTI, H. R. The Viscosity of Glycerol- Water Mixtures Including the Supercooled Region. J. Chem. Eng. Data, 56, 1397-1406, 2011. HIMMELBLAU, D. M. Basic Principles and Calculations in Chemical Engineering. 5th edition. Englewood Cliffs-NJ, USA: Prentice Hall, 1989. INCROPERA, F. P.; DEWITT, D. P. Fundamentos de Transferência de Calor e de Massa. 4ª edição. São Paulo, SP: LTC, 2002. 85 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM KNOTHE, G.; GERPEN, J. V.; KRAHL, J. (editors). The Biodiesel Handbook. Champaign, EUA: AOCS Press, 2005. LANE, L. B. Freezing points of glycerol and its aqueous solutions. Industrial and Engineering Chemistry. Vol. 19, no. 9, 1925. LOSENICKY, Z. Thermal Conductivity of Binary Liquid Solutions. The Journal of Physical Chemistry. Vol. 72, No. 12, pp. 4308-4310, November, 1968. MEGLOBAL. Ethylene Glycol: product guide. London, UK: Meglobal, 2008. MELINDER, A. Thermophysical Properties of Aqueous Solution Used as Secondary Working Fluids [doctoral thesis]. Stockholm, Sweden: Royal Institute of Technology, Dept. of Energy Technology; 2007. MOTA, C. J. A.; PESTANA, C. F. M. Co-produtos da Produção de Biodiesel. Rev. Virtual Quim, 3 (5), 416-425, 2011. OTT, J. B.; BOERIO-GOATES, J. Chemical Thermodynamics: principles and applications. Volume 1. San Diego-CA, USA: Elsevier, 2000. PERRY, R. H.; GREEN, D. W. (editors). Perry’s Chemical Engineers’ Handbook. 8th edition. New York-NY, USA: McGraw-Hill, 2008. POLING, B. E.; PRAUSNITZ, J. M.; J. P. O’CONNELL. The Properties of Gases and Liquids. 5th edition. New York-NY, USA: McGraw-Hill, 2004. PRAUSNITZ, J. M.; LICHTENTHALER, R. N.; AZEVEDO, E. G. Molecular Thermodynamics of Fluid-Phase Equilibria. 3rd edition. Upper Saddle River-NJ, USA: Prentice Hall, 1999. PRUZAESKY, F. C.; et al. Novas soluções para fluidos secundários em refrigeração. In: MMA. Uso de Fluidos Naturais em Sistemas de Refrigeração e Ar-condicionado. Publicação Técnica. Brasília-DF: Ministério do Meio Ambiente – MMA, 2008. p. 63-78. SANTOS, R. G. dos. Transformações de Fases em Materiais Metálicos. Campinas, SP: Editora da Unicamp, 2006. SECCOOL. Properties. Software Version 1.33 [S.I.]: Desik Project Group, 2007. SMITH, J. M.; VAN NESS, H. C.; ABBOTT, M. M. Introdução à Termodinâmica da Engenharia Química. 7ª edição. Rio de Janeiro, RJ: LTC, 2007. SOLOMONS, T. W. G.; FRYHLE, C. B. Organic Chemistry. 9th edition. Hoboken-NJ, USA: John Wiley and Sons, 2007. 86 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM STOECKER, W. F. Industrial Refrigeration Handbook. New York-NY, USA: McGraw- Hill, 1998. SKOVRUP, M. J. SecCool Datafit: users manual. Lyngby-KGS, Denmark: IPU, 2005. TEJA, A. S. Simple Method for the Calculation of Heat Capacities of Liquid Mixtures. Journal of Chemical and Engineering Data, Vol. 28, No. 1, pp. 83-85, 1983. THOMPSON, J. C.; HE, B. B. Characterization of Crude Glycerol from Biodiesel Production from Multiple Feedstocks. Appl. Eng. Agric. Vol. 22, 261-265, 2006. TORREIRA, R. P. Fluidos Térmicos: água, vapor, óleos térmicos. Curitiba, PR: Hemus, 2002. VECCHIA, R. O meio ambiente e as energias renováveis: instrumentos de liderança visionária para a sociedade sustentável. Barueri, SP: Manole, 2010. WALLAS, S. M. Phase Equilibria in Chemical Engineering. Newton-MA, USA: Butterworth-Heinemann, 1985. Apêndice 88 Pedro Samuel Gomes Medeiros – Dissertação de Mestrado – UFRN/PPGEM Apêndice Propriedades termofísicas da solução água-propilenoglicol, com 35,5 % de concentração mássica de propilenoglicol (PC real = -16,8 °C), em função da temperatura (tabela A1): Tabela A1 – Propriedades Água-Propilenoglicol (64,5/35,5) - ASHRAE Prop/Temp [°C] -15 -10 -5 0 5 10 15 20 25 ρ [kg/m³] 1045,8 1044,3 1042,6 1040,8 1038,9 1036,9 1034,7 1032,5 1030,0 C [kJ/kg.K] 3,672 3,687 3,702 3,717 3,732 3,747 3,763 3,778 3,793 k [W/m.K] 0,3858 0,3920 0,3978 0,4034 0,4087 0,4137 0,4185 0,4230 0,4272 µ [cP] 23,28 16,95 12,60 9,547 7,371 5,792 4,629 3,759 3,100 Fonte: SecCool, 2007 Propriedades termofísicas da solução água-glicerol, com 40% de concentração mássica de glicerol (PC real = -15,4 °C), em função da temperatura (tabela A2): Tabela A2 – Propriedades Água-Glicerol (60/40) - MELINDER Prop/Temp [°C] -15 -10 -5 0 5 10 15 20 25 ρ [kg/m³] 1111,5 1110,2 1108,8 1107,2 1105,4 1103,5 1101,4 1099,3 1097,0 C [kJ/kg.K] 3,281 3,305 3,329 3,353 3,377 3,400 3,424 3,448 3,471 k [W/m.K] 0,4200 0,4243 0,4287 0,4330 0,4374 0,4418 0,4462 0,4506 0,4549 µ [cP] 17,53 13,04 9,999 7,877 6,353 5,227 4,374 3,708 3,174 Fonte: SecCool, 2007 Propriedades termofísicas do glicerol puro, em função da temperatura (tabela A3): Tabela A3 – Propriedades Termofísicas do Glicerol Prop/Temp [°C] 0 5 10 15 20 25 ρ [kg/m³] 1276 1273 1270 1267 1264 1261 C [kJ/kg.K] 2,262 2,288 2,320 2,354 2,386 2,416 k [W/m.K] 0,2820 0,2835 0,2846 0,2856 0,2860 0,2860 Fonte: Çengel, 2007