You do not have permission to edit this page, for the following reason:

You are not allowed to execute the action you have requested.


You can view and copy the source of this page.

x
 
1
==Estudo numérico do processo de fusão de PCMs em cavidade cilíndrica==
2
3
===Numerical study of PCMs melting process in cylindrical cavity===
4
5
'''J. Juarez Estrazulas''', '''R.C. Oliveski'''
6
7
FTEC - Caxias do Sul, RS BRAZIL
8
9
10
11
==Resumo==
12
13
Inúmeras aplicações residenciais, comerciais e industriais tem seu rendimento aumentado quando um sistema de armazenamento de energia térmica é incorporado. Os PCMs (Phase Change Materials), devido ao seu alto calor latente de fusão, são materiais que representam uma alternativa viável à implementação de sistemas de armazenamento de energia térmica. Este trabalho apresenta um estudo numérico do processo de fusão de PCMs da família RT (Rubitherm Technologies GmbH) em cavidades cilíndricas. O estudo foi realizado através de simulação numérica por CFD (Computational Fluid Dynamics), com o software ANSYS Fluent. O modelo numérico adotado é bidimensional e foi validado com resultados numéricos e experimentais da literatura, obtendo-se boa aproximação. Para os PCMs aqui estudados, aumentando-se a temperatura de 10 <sup>o</sup>C para 20 <sup>o</sup>C e de 10 <sup>o</sup>C para 30 <sup>o</sup>C, acima da temperatura de mudança de fase, para frações líquidas entre 0,4 e 0,8, a redução média dos tempos de fusão é de, aproximadamente, 55,8% e 71,8%, respectivamente.
14
15
'''Palavras-chave''': Armazenamento térmico, phase change materials, PCM, CFD.
16
17
18
==Abstract==
19
20
Several residential, commercial and industrial applications have their efficiency increased when a thermal energy storage system is incorporated. The PCMs (Phase Change Materials), due to their high latent heat of fusion, are materials that represent a viable alternative to the implementation of thermal energy storage systems. This paper presents a numerical study of RT (Rubitherm Technologies GmbH) PCMs family fusion process in cylindrical cavities. The study was conducted through a CFD ''(''Computational Fluid Dynamics) numerical simulation, with ANSYS Fluent software. The numeric model adopted is two-dimensional and has been validated with numerical and experimental results of literature, achieving a good approximation. For PCMs here studied, increasing the temperature from 10 ''<sup>o</sup>C'' to 20 ''<sup>o</sup>C'' and 10 ''<sup>o</sup>C'' to 30 ''<sup>o</sup>C'', above the phase change temperature, for liquid fractions between 0.4 and 0.8, the average reduction of fusion time is approximately 55.8% and 71.8 %, respectively.
21
22
'''Keywords''': Thermal storage, phase change materials, PCM, CFD. 
23
24
25
==1 Introdução==
26
27
O contínuo crescimento da emissão de gases na atmosfera, decorrente do aumento da população mundial e os altos custos dos combustíveis fósseis e da energia elétrica são os motivos básicos para uma utilização mais eficiente da energia em todas as áreas da atividade humana. Edifícios comerciais e apartamentos residenciais são os consumidores de energia mais intensivos. Nos Estados Unidos, 41,1% de toda a energia primária consumida no ano de 2010, foi demandada por consumidores residenciais e comerciais, sendo o aquecimento, refrigeração e iluminação as maiores demandas [1]. Já na Europa, os edifícios residenciais consomem 40% de toda a energia [2]. No Brasil, estes setores consomem 13,7% de toda a energia primária [3]. A busca de novas tecnologias para diminuir a crescente preocupação com os problemas ambientais atuais, com a possível carência de energia, seu alto custo e os pesados investimentos necessários à criação de novas centrais de geração de energia tem constituído um dos principais problemas científicos nas últimas três décadas. O aspecto central deste problema é a necessidade de armazenamento de energia, térmica ou elétrica, com a finalidade de cobrir as lacunas existentes entre a geração de energia e o seu respectivo consumo. Estas lacunas podem ser de ordem temporal e de demanda [4].
28
29
Os métodos usados para armazenamento térmico de energia (TES-Thermal Energy Storage) são classificados em armazenamento térmico de energia por calor sensível (SHTES-Sensible Heat Thermal Energy Storage) e armazenamento térmico de energia por calor latente (LHTES - Latent Heat Thermal Energy Storage). Os SHTES dependem do aumento ou redução do calor sensível de um material. Já os LHTES dependem da capacidade do material em armazenar energia em uma condição praticamente isotérmica, com estreita faixa de variação. Estes tipos de sistemas utilizam PCMs (Phase Change Materials) e oferecem como vantagem, além da temperatura constante durante os ciclos de armazenamento e liberação da energia, maior capacidade de armazenamento e menor volume quando comparados aos SHTES [5]. Os LHTES constituem uma técnica eficaz, tanto para fins de aquecimento quanto refrigeração, porque eles são capazes de prover uma alta densidade de armazenamento de energia. Quando comparados com os SHTES, os LHTES requerem menor peso e volume de material para uma determinada quantidade de energia [4].
30
31
De acordo com Zhou et al. [6], a maioria dos PCMs sofrem de um problema comum que é a baixa condutividade térmica, com valores ao redor de 0,2 W/mK para os materiais parafínicos e 0,5 W/mK para os hidratos de sal e materiais eutéticos. Esta baixa condutividade aumenta o tempo dos processos de carga e descarga de um LHTES que utiliza PCM como meio de armazenamento térmico. Al-Abidi et al. [7] destacam que a maioria das técnicas propostas para a melhoria na condutividade térmica tem sido baseadas na introdução de aletas metálicas, nas mais variadas formas e configurações, nos tanques que armazenam os PCMs. Este fato é devido à simplicidade deste tipo de implementação, com facilidade de fabricação e baixo custo. Outras estratégias abordadas na pesquisa de Agyenim et al. [4] abrangem o uso de múltiplos tubos metálicos, a agitação através da inserção de bolhas, o uso de anéis metálicos, o uso de multitubos e escovas de carbono, o encapsulamento do PCM em material polimérico (microencapsulamento), o uso de uma matriz metálica ou de grafite, o uso de esferas metálicas ou de material poliolefínico, a adição de flocos de grafite e a compactação do PCM em painéis planos ou na forma de barras. Em alguns casos, como no acréscimo de anéis metálicos ao PCM, a condutividade térmica foi aumentada em 10 vezes.
32
33
Outro problema apresentado por alguns PCMs é o subresfriamento, ou seja, quando o PCM não se solidifica imediatamente após sua temperatura estar abaixo da temperatura de mudança de fase. Neste caso a solidificação do PCM só é conseguida levando-se a sua temperatura bem abaixo da temperatura de mudança de fase, para iniciar a cristalização do material e liberar o calor latente nele armazenado. A técnica mais comum para solucionar este problema é a incorporação de aditivos especiais, denominados ¨nucleadores¨, aos PCMs. A limitação atualmente existente é que estes ¨nucleadores¨ somente são eficazes numa faixa de temperatura entre 10 a 20  K acima ou abaixo da temperatura de mudança de fase [8, 9].
34
35
Como mencionado anteriormente, os LHTES são utilizados para o armazenamento de energia térmica e no controle e manutenção da temperatura. Dentre as aplicações dos PCMs nestes campos, destacam-se: atenuação do aquecimento de baterias de íons de lítio e de dispositivos móveis; tubos de calor e dissipadores de calor, para refrigeração de componentes eletrônicos; refrigeração comercial (caminhões) e residencial (freezers); embalagens para produtos farmacêuticos ou derivados do sangue; roupas para uso em regiões com temperaturas extremas; conforto térmico residencial (através do uso de PCM microencapsulado na argamassa de paredes, placas de PCM instaladas no piso, forro ou sótão); fogões solares e secadores de grãos; aquecimento solar de água.
36
37
Diversos pesquisadores focaram seus esforços no estudo de geometrias cilíndricas. Murray e Groulx [10], por exemplo, analisaram experimentalmente o processo de transferência de calor durante consecutivos ciclos carga/descarga de PCM visando aplicação em sistema de aquecimento solar de água. Eles demonstraram que o aumento da taxa de escoamento do líquido de transferência térmica (HTF), durante o processo de carga, resulta em fusão significativamente mais rápida. Por outro lado, durante o processo de descarga, um aumento do escoamento não afeta o tempo requerido para a solidificação do PCM.
38
39
Longeon et al. [11] estudaram, de forma numérica e experimental, a influência do ponto de entrada do fluido de transferência térmica (HTF) no tempo de fusão/solidificação do PCM RT 35 em um trocador de calor tipo casco-tubo. Os autores concluíram que o melhor rendimento do sistema era obtido executando a entrada do HTF pelo topo do trocador de calor durante o processo de carga e pela base durante o processo de descarga.
40
41
Zeng et al. [12] analisaram experimentalmente a fusão de PCM com propriedades térmicas incrementadas através da adição de nanopartículas (NePCM). Estes concluíram que o aumento de condutividade térmica obtido com a adição de nanopartículas reduz drasticamente os processos de convecção no PCM, não se traduzindo em redução no tempo de fusão do mesmo.
42
43
Al-Abidi et al. [13] investigaram numericamente os processos de solidificação do PCM RT 82 em um trocador de calor composto por três tubos verticais concêntricos, sendo o PCM armazenado no tubo intermediário, dotado de diferentes configurações de aletas e espessura destas aletas. Eles concluíram que o número de aletas e o comprimento delas influenciam diretamente na redução dos tempos de solidificação, ao passo que a espessura das aletas não tem efeito sobre o processo de solidificação.
44
45
Tay et al. [14]) executaram análise experimental do processo de fusão de PCM em tanque cilíndrico vertical dotado de tubos de cobre que executavam a circulação do fluido de transferência térmica e também da recirculação do PCM já liquefeito, processo ao qual estes autores deram o nome de “dynamic melting”. Com a adoção deste tipo de técnica, os autores demonstraram que a condutividade térmica pode ser duplicada, sendo este efeito comparável ao conseguido com a adição de aletas ao sistema.
46
47
Kalaiselvam et al. [15] avaliaram, tanto de forma experimental como numérica, os processos de fusão e solidificação de três PCMs parafínicos. Na etapa de solidificação, os processos numéricos resultaram em um tempo 11% inferior aos obtidos experimentalmente. Já durante o processo de fusão, os tempos obtidos numericamente foram 16,11% inferiores aos coletados nos testes.
48
49
Observa-se que pesquisadores tem feito esforços no sentido de apresentar estratégias de solução de problemas envolvendo PCMs. No entanto, devido à diversidade de PCMs disponíveis, vários destes materiais não foram ainda estudados detalhadamente, seja de forma numérica ou experimental.
50
51
Apesar da potencialidade da utilização dos PCMs na área de armazenamento térmico de energia, eles ainda apresentam deficiências e lacunas que vem sendo analisadas pelos pesquisadores, resultando em técnicas que visam o aprimoramento das características destes materiais e maior previsibilidade nos cálculos de capacidade térmica, tempos de fusão, solidificação e demais parâmetros necessários à correta aplicação dos PCMs. Dentre as deficiências dos PCMs pode-se citar a baixa condutividade térmica, a instabilidade de alguns PCMs ao longo dos ciclos de carga/descarga e a ausência de pesquisas detalhando inúmeros PCMs comercialmente disponíveis no mercado mundial.
52
53
Este trabalho visa contribuir para a diminuição desta lacuna, estudando detalhadamente a mudança de fase de quatro diferentes PCMs da família RT [16] em cavidades cilíndricas. Assim, os objetivos específicos são: (a) implementar numericamente o modelo matemático de mudança de fase; (b) validar a metodologia e o modelo matemático com resultados da literatura; (c) avaliar e (d) comparar o processo de fusão de diferentes PCMs quando submetidos a diferentes temperaturas acima da sua temperatura de mudança de fase (T<sub>pc</sub>).
54
55
==2 Metodologia==
56
57
===2.1 Abordagem Numérica===
58
59
A dinâmica de fluidos computacional (CFD) é o conjunto de ferramentas numéricas e computacionais utilizadas para resolver, visualizar e interpretar a solução das equações de balanço de momento, massa e energia. Tal ferramenta utiliza técnicas numéricas para solucionar problemas de engenharia, física e muitas outras áreas, já que esse tipo de simulação computacional é capaz de predizer fenômenos físicos e físico-químicos que ocorrem em um dado escoamento. Essa ferramenta vem sendo cada vez mais utilizada no desenvolvimento de produtos de consumo, aplicações industriais, aeroespaciais e de defesa. Em qualquer aplicação que envolva o escoamento de gases, escoamento de líquidos ou transferência de calor, a análise por CFD possibilita a inovação e novos patamares de eficiência para produtos e processos. Atualmente existem diversos softwares comerciais de CFD disponíveis no mercado. Dentre os mais conhecidos destacam-se o ANSYS CFX, o ANSYS Fluent, o Star-CD e o Flow-3D, entre outros.
60
61
O presente estudo foi desenvolvido utilizando-se os softwares CFD ANSYS Fluent (para as simulações numéricas), ANSYS CFX (para o pós-processamento dos resultados) e ANSYS ICEM (para a criação da malha computacional), todos disponíveis no Laboratório de Simulação Numérica do Programa de Pós-Graduação em Engenharia Mecânica da UNISINOS.
62
63
===2.2 Definição do problema===
64
65
====2.2.1  Definição dos parâmetros geométricos====
66
67
A geometria estudada consiste de uma cavidade cilíndrica vertical, conforme apresentado na Figura 1, a qual também indica que a aceleração da gravidade está em sentido oposto ao eixo Y. Nesta geometria o diâmetro (D) é de 30 mm e a altura H1 (=170 mm) representa a quantidade de PCM dentro do cilindro no tempo inicial (t=0). A altura H2 (=40 mm) é a folga necessária para compensar a dilatação térmica do PCM. Assim, durante o processo de aquecimento, ocorre um deslocamento ascendente da interface PCM-ar, sem que haja escape de PCM. A fase sólida é mostrada em cor azul e o ar, na parte superior (aberta ao ambiente), é mostrado em cor vermelha. Na base do cilindro está representada sua condição adiabática.
68
69
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
70
71
[[Image:draft_Samper_486932730-image1.png|center|318px]]
72
</div>
73
74
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
75
'''Figura 1''' – Geometria utilizada. Para o tempo t = 0 s, a fase sólida é representada pela cor azul e o ar pela cor vermelha (adaptado de [17]).</div>
76
77
====2.2.2  Seleção do material====
78
79
A empresa alemã Rubitherm dispõe de uma grande variedade de PCMs, denominados RT. Cada RT apresenta uma temperatura de mudança de fase que varia de -9 <sup>o</sup>C até 100 <sup>o</sup>C. O valor da temperatura de mudança de fase é indicado depois das letras RT, conforme mostrado na Tabela 1.
80
81
As análises propostas para este trabalho compreendem o estudo numérico do processo de fusão dos PCMs relacionados na Tabela 1, quando submetidos a temperaturas de 10, 20 e 30 <sup>o</sup>C acima de suas respectivas temperaturas de mudança de fase (T<sub>pc</sub>), ou seja com &#x0394;T=10, 20 e 30 ºC.
82
83
A escolha dos PCMs a serem estudados numericamente foi feita visando possíveis aplicações de mercado. Assim sendo, os RTs selecionados foram: o RT 4, pois pode ser utilizado para aplicações de conservação e transporte de alimentos ou produtos farmacêuticos; o RT 35 e RT 55 são de interesse na área de aquecimento solar de água e, por fim, o RT 82 também foi acrescentado, pois pode ser utilizado em algumas aplicações industriais.
84
85
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
86
'''Tabela 1''' – Temperatura de fusão, calor latente e condutividade térmica dos PCMs.</div>
87
88
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
89
Adaptado de [16]</div>
90
91
{| style="margin: 1em auto 1em auto;border-top: 1pt solid black;border-bottom: 1pt solid black;border-collapse: collapse;"
92
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
93
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">PCM</span>
94
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Temperatura</span>
95
96
<span style="text-align: center; font-size: 75%;">de Fusão [°C]</span>
97
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Calor Latente de Fusão.</span>
98
99
<span style="text-align: center; font-size: 75%;">L±7,5% [kJ kg<sup>-1</sup></span><span style="text-align: center; font-size: 75%;">]</span>
100
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Condutiv. Térmica</span>
101
102
<span style="text-align: center; font-size: 75%;"> k [W m<sup>-1</sup></span><span style="text-align: center; font-size: 75%;"> K<sup>-1</sup></span><span style="text-align: center; font-size: 75%;">]</span>
103
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
104
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">RT 4</span>
105
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2 – 4</span>
106
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">182</span>
107
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,2</span>
108
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
109
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">RT 35</span>
110
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">29 – 36</span>
111
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">170</span>
112
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,2</span>
113
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
114
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">RT 55</span>
115
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">51 – 57</span>
116
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">172</span>
117
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,2</span>
118
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
119
| style="border-top: none;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">RT 82</span>
120
| style="border-top: none;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">77 - 85</span>
121
| style="border-top: none;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">176</span>
122
| style="border-top: none;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,2</span>
123
|}
124
125
Na Tabela 2 são apresentados as simulações e detalhes do tipo de PCM, constante C (conhecida como ¨mushy zone¨, que corresponde à região do PCM que está em processo de mudança de fase) e o ∆T investigados. Para estas simulações serão analisadas a influência da temperatura sobre a velocidade do processo de fusão dos PCMs. Vale destacar que a constante C é detalhada no modelo matemático, Seção 3.3.
126
127
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
128
'''Tabela 2''' – Simulações, tipo de PCM, T<sub>pc</sub>, C e &#x0394;T.</div>
129
130
{| style="margin: 1em auto 1em auto;border: 1pt solid black;border-collapse: collapse;"
131
|- style="border: 1pt solid black;"
132
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<br/>Simulação
133
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|PCM
134
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|T<sub>pc</sub> [<sup>o</sup>C]
135
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|C [kg m<sup>-3</sup> s<sup>-1</sup>]
136
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|∆T [<sup>o</sup>C]
137
|- style="border: 1pt solid black;"
138
| style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;"|1- 3
139
| style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;"|RT 4
140
| style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;"|2-4
141
| style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;"|10<sup>7</sup>
142
| style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;"|10, 20 e 30
143
|- style="border: 1pt solid black;"
144
| style="border: none;text-align: center;"|4- 6
145
| style="border: none;text-align: center;"|RT 35
146
| style="border: none;text-align: center;"|29-36
147
| style="border: none;text-align: center;"|10<sup>7</sup>
148
| style="border: none;text-align: center;"|10, 20 e 30
149
|- style="border: 1pt solid black;"
150
| style="border: none;text-align: center;"|7- 9
151
| style="border: none;text-align: center;"|RT 55
152
| style="border: none;text-align: center;"|51-57
153
| style="border: none;text-align: center;"|10<sup>7</sup>
154
| style="border: none;text-align: center;"|10, 20 e 30
155
|- style="border: 1pt solid black;"
156
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|10- 12
157
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|RT 82
158
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|77-85
159
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|10<sup>7</sup>
160
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|10, 20 e 30
161
|}
162
163
As simulações foram feitas utilizando as especificações técnicas de cada PCM (temperaturas de fusão e solidificação, calor latente de fusão, massas específicas nos estados sólido e líquido, condutividade térmica e viscosidade dinâmica). Um detalhamento das propriedades termodinâmicas e de como a massa específica de cada PCM é interpolada consta no Apêndice B.
164
165
====2.2.3  Modelo matemático====
166
167
As equações de conservação utilizadas no modelamento de um sistema com PCM, através da técnica de modelagem entalpia-porosidade, são [17]:
168
169
{| style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" class="formulaSCP" 
170
|-
171
| 
172
{| style="text-align: center; margin:auto;width: 100%;"
173
|-
174
| <math display="inline">\frac{\partial \rho }{\partial t}+\nabla .\left(\rho {V}\right)=0,\,</math>
175
|}
176
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(1)</span>
177
|-
178
| 
179
{| style="text-align: center; margin:auto;width: 100%;"
180
|-
181
| <math display="inline">\, \, \rho \frac{D\rightarrow {V}}{Dt}=-\nabla p+\mu {\nabla }^{\, 2}\rightarrow {V}+\, \rho \rightarrow {g}+\rightarrow {S}</math>
182
|}
183
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(2)</span>
184
|-
185
| 
186
{| style="text-align: center; margin:auto;width: 100%;"
187
|-
188
| <math display="inline">\, \rho \frac{Dh}{Dt}=k{\nabla \, }^{2}T</math>
189
|}
190
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(3)</span>
191
|-
192
| 
193
{| style="text-align: center; margin:auto;width: 100%;"
194
|-
195
|  <math display="inline">\frac{\partial \left({\alpha }_{n}\right)}{\partial t}+\nabla .\left({\gamma }_{n}{V}\right)=0</math>
196
|}
197
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(4)</span>
198
|}
199
200
onde <math display="inline">{\alpha }_{n}</math> é a fração volumétrica do PCM na enésima célula computacional, podendo assumir três condições: se <math display="inline">{\alpha }_{n}=0,\,</math> a enésima célula não contém o PCM; se <math display="inline">{\alpha }_{n}=1,</math> a enésima célula está cheia de PCM; se <math display="inline">{0<\alpha }_{n}<1,</math> a enésima célula contém a interface entre o PCM e um ou mais fluidos existentes no domínio computacional; <math display="inline">\rho</math>  é a massa específica, k é a condutividade térmica, μ é a viscosidade dinâmica, <math display="inline">\rightarrow {S}</math> é o termo fonte da quantidade de movimento, <math display="inline">\rightarrow {V}</math> é o vetor velocidade, T é a temperatura, h é entalpia específica e p é a pressão.
201
202
A entalpia do material é definida como a soma da entalpia sensível <math display="inline">{h}_{s\, }(Equa\mbox{ç}\tilde{a}o\, \, 5</math>) e da variação de entalpia devida à mudança de fase (γL), sendo, ''h''<sub>ref</sub> a entalpia de referência na temperatura de referência (T<sub>ref</sub>), C<sub>p</sub> o calor específico, L a entalpia de fusão (calor latente do material) e <math display="inline">\gamma</math>  a fração líquida durante a mudança de fase que ocorre em uma faixa de temperatura T<sub>s</sub> < T < T<sub>l</sub>, definida pela seguinte relação [17]:
203
204
{| style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" class="formulaSCP" 
205
|-
206
| 
207
{| style="text-align: center; margin:auto;width: 100%;"
208
|-
209
| <math display="inline">{\, h}_{s}=\, {h}_{ref}+\, \int_{{T}_{ref}}^{T}{c}_{p}\, dT,</math>
210
|}
211
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(5)</span>
212
|-
213
| 
214
{| style="text-align: center; margin:auto;width: 100%;"
215
|-
216
| <math display="inline">\gamma =\left\{0 \ se \ T<T_s,\right\}</math>
217
|}
218
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(6)</span>
219
|}
220
221
A técnica entalpia-porosidade, analisada por Brent e Voller (1988), trata a região de mudança de fase (parcialmente solidificada), como se fosse um meio poroso. A porosidade em cada célula é definida como igual à fração líquida naquela célula. Em regiões totalmente solidificadas a porosidade é igual a zero, o que extingue as velocidades nestas regiões. O termo fonte ( <math display="inline">\rightarrow {S}</math>) da equação da quantidade de movimento (Equação 2) toma a seguinte forma:
222
223
{|
224
|-
225
| style="width: 55.0157mm;vertical-align: top;" | <big>
226
| style="" | </big> <math display="inline">\rightarrow {S}\, =\, -A\left(\gamma \right)\, \rightarrow {V}</math>
227
| style="" | 
228
| style="" | (7)
229
|}
230
231
onde A( <math display="inline">\gamma )</math> é a “função porosidade” definida por Brent e Voller (1988).
232
233
{|
234
|-
235
| style="width: 55.0157mm;vertical-align: top;" | <big>
236
| style="" | </big> <math display="inline">A(\gamma )\, =\frac{{C\, \left(1-\gamma \right)}^{\, 2}}{\left({\gamma }^{3}+\epsilon \right)},</math>
237
| style="" | (8)
238
|}
239
240
onde <math display="inline">\epsilon \,</math> é uma pequena constante para evitar a divisão por zero, recebendo valor de <math display="inline">0,001\,</math>  e C é a constante que reflete o processo cinético e a morfologia da região de mudança de fase. Entre outros, esta constante foi estudada por Brent e Voller [18] e Shmueli et al. [17]. Estes pesquisadores demonstraram que C controla, efetivamente, o grau de penetração de massa na região de mudança de fase e, desta maneira, conduz a um processo de fusão (ou solidificação) mais rápido ou mais lento, dependendo do valor adotado. Normalmente estes valores variam entre 10<sup>5</sup> à 10<sup>10</sup> e dependem da geometria utilizada. A constante C, também conhecida como “mushy zone”, foi assumida igual a 10<sup>7</sup>, mesmo valor utilizado por Shmueli et al. [17].
241
242
====2.2.4  Validação numérica====
243
244
Para a validação do modelo numérico utilizou-se os resultados de Shmueli et al. [17]. Estes autores investigaram a fusão do PCM RT 27 em cavidade cilíndrica vertical, mostrada na Figura 1, e compararam diversas simulações com os resultados experimentais de Katsman et al. [19].
245
246
O problema consiste em manter a temperatura da superfície externa 10 <sup>o</sup>C acima da temperatura média de fusão do PCM. A parte superior da geometria cilíndrica é exposta ao ar ambiente e a parte inferior foi considerada adiabática. As propriedades do RT 27, cuja temperatura de mudança de fase é 27 <sup>o</sup>C, são apresentadas na Tabela 3.
247
248
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
249
'''Tabela 3''' – Propriedades do RT 27 (adaptado de [17]).</div>
250
251
{| style="width: 100%;border-top: 1pt solid black;border-bottom: 1pt solid black;border-collapse: collapse;"
252
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
253
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">Estado:</span>
254
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">C<sub>p</sub></span>
255
256
<span style="text-align: center; font-size: 75%;">[J kg<sup>-1</sup></span><span style="text-align: center; font-size: 75%;">]</span>
257
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">k</span>
258
259
<span style="text-align: center; font-size: 75%;">[W m<sup>-1</sup></span><span style="text-align: center; font-size: 75%;"> K<sup>-1</sup></span><span style="text-align: center; font-size: 75%;">]</span>
260
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">&#x03c1;</span>
261
262
<span style="text-align: center; font-size: 75%;">[kg m<sup>-3</sup></span><span style="text-align: center; font-size: 75%;">]</span>
263
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">L</span>
264
265
<span style="text-align: center; font-size: 75%;">[kJ kg<sup>-1</sup></span><span style="text-align: center; font-size: 75%;">]</span>
266
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">μ</span>
267
268
<span style="text-align: center; font-size: 75%;">[Pa s]</span>
269
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
270
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">Sólido</span>
271
272
273
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">2500</span>
274
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,2</span>
275
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">870</span>
276
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">170</span>
277
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">-</span>
278
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
279
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">Líquido</span>
280
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">2500</span>
281
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,2</span>
282
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">870 em 299K</span>
283
284
<span style="text-align: center; font-size: 75%;">781,5 em 310K</span>
285
286
<span style="text-align: center; font-size: 75%;">750 em 343K</span>
287
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;"|<span style="text-align: center; font-size: 75%;">---</span>
288
| style="border-top: 1pt solid black;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">-1,137439x10<sup>-8</sup></span><span style="text-align: center; font-size: 75%;">T<sup>3</sup></span><span style="text-align: center; font-size: 75%;"> + 1,1788188x10<sup>-5</sup></span><span style="text-align: center; font-size: 75%;">T<sup>2</sup></span><span style="text-align: center; font-size: 75%;"> – 0,004111388T + 0,4857203</span>
289
|}
290
291
Na simulação executada, a temperatura inicial de todo o sistema foi definida como 22 <sup>o</sup>C, ou seja, o PCM está ligeiramente sub-resfriado, de forma a garantir que esteja em estado sólido no início da simulação. Foi assumido que ambas as fases, sólida e líquida, do PCM são homogêneas e isotrópicas e que o processo de fusão é simétrico em relação ao eixo central (Figura 1). Adicionalmente foi considerado que a fase líquida do PCM e o ar são fluidos Newtonianos, incompressíveis e o regime de escoamento laminar foi assumido para ambos.
292
293
====:2.2.5  Definição dos Parâmetros Computacionais====
294
295
Para realizar as simulações foi utilizado o software ANSYS Fluent. As condições de contorno adotadas foram as mesmas especificadas no trabalho de Shmueli et al. [17]. Para descrever a interface entre o PCM e o ar, a qual se desloca com o tempo, sem interpenetração dos dois fluidos, o modelo de “volume de fluido” ou VOF foi utilizado na simulação.
296
297
As condições de contorno para as velocidades são de não-deslizamento e não-penetração em todos os contornos sólidos. A parede externa da geometria cilíndrica foi mantida a uma temperatura constante de 38<sup>o</sup>C (10 <sup>o</sup>C acima da temperatura máxima de fusão do RT 27, que é de 28 <sup>o</sup>C) e na superfície inferior foi considerada a condição adiabática. No contorno superior, exposto à atmosfera ambiente, a condição de pressão de saída (Pressure Outlet) foi utilizada e a temperatura ambiente assumida foi de 38 <sup>o</sup>C.
298
299
A geração da malha computacional constitui uma etapa determinante no processo de análise numérica. Uma malha pobre, ou seja, com poucos elementos, faz com que o processo numérico não represente com a acurácia necessária os diferentes processos físicos que ocorrem durante a mudança de fase do PCM. Por outro lado, deve-se ter sempre em conta que, quanto maior o número de volumes de controle de uma malha, tanto maior será o custo computacional, em termos de tempo de processamento e de capacidade de equipamentos. A malha computacional mostrada na Figura 2 foi criada seguindo as recomendações do trabalho de Shmueli et al. [17], tendo cada um de seus elementos 0,5 mm de lado, com um total de 13.652 elementos.
300
301
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
302
303
[[Image:draft_Samper_486932730-image2.png|center|438px]]
304
</div>
305
306
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
307
'''Figura 2''' - Malha computacional utilizada na validação numérica: (i) aspecto geral da malha, (ii) detalhe da extremidade superior e (iii) detalhe da interface PCM-Ar.</div>
308
309
Foi adotado um refinamento em duas regiões distintas: junto à parede externa do cilindro e junto à interface entre o PCM e o ar. Estes refinamentos têm como objetivo, respectivamente, capturar de forma mais detalhada os efeitos de condução e convecção nos instantes iniciais da simulação e melhor acompanhar o deslocamento da interface, a qual ocorre em função do aumento de volume do PCM.
310
311
Na validação numérica foi utilizado passo de tempo (&#x0394;t) de 0,01 s e os critérios de convergência (CC) para continuidade e as componentes de velocidade foram de 10<sup>-3</sup>, sendo de 10<sup>-6</sup> para a energia, ambos calculados com dupla precisão. Vale destacar que estes critérios de convergência foram os mesmos utilizados por Shmueli et al. [17].
312
313
Visando confirmar que o ∆t e os critérios de convergência (CC) adotados não implicavam em erro para o processo numérico, foram executadas as simulações detalhadas na Tabela 4, variando estes parâmetros.
314
315
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
316
'''Tabela 4''' – Simulações com variação de ∆t e CC.</div>
317
318
{| style="border-collapse: collapse;"
319
|-
320
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|Simulação
321
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|Constante
322
323
C [kg m<sup>-3</sup> s<sup>-1</sup>]
324
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|Passo de
325
326
tempo ∆t [s]
327
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|CC
328
329
Contin. e Veloc.
330
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|CC
331
332
Energia
333
|-
334
| style="border-top: 1pt solid black;text-align: center;"|A
335
| style="border-top: 1pt solid black;text-align: center;"|10<sup>7</sup>
336
| style="border-top: 1pt solid black;text-align: center;"|0,01
337
| style="border-top: 1pt solid black;text-align: center;"|10<sup>-3</sup>
338
| style="border-top: 1pt solid black;text-align: center;"|10<sup>-6</sup>
339
|-
340
| style="text-align: center;"|B
341
| style="text-align: center;"|10<sup>7</sup>
342
| style="text-align: center;"|0,01
343
| style="text-align: center;"|10<sup>-6</sup>
344
| style="text-align: center;"|10<sup>-8</sup>
345
|-
346
| style="text-align: center;"|C
347
| style="text-align: center;"|10<sup>7</sup>
348
| style="text-align: center;"|0,002
349
| style="text-align: center;"|10<sup>-3</sup>
350
| style="text-align: center;"|10<sup>-6</sup>
351
|-
352
| style="border-bottom: 1pt solid black;text-align: center;"|D
353
| style="border-bottom: 1pt solid black;text-align: center;"|10<sup>7</sup>
354
| style="border-bottom: 1pt solid black;text-align: center;"|0,002
355
| style="border-bottom: 1pt solid black;text-align: center;"|10<sup>-6</sup>
356
| style="border-bottom: 1pt solid black;text-align: center;"|10<sup>-8</sup>
357
|}
358
359
Na Tabela 5 são apresentados os resultados das simulações de fração líquida para as 4 simulações da Tabela 4. Os resultados indicam que, mesmo com o uso de passo de tempo e critérios de convergência mais refinados como os usados nas simulações B, C e D, os resultados são praticamente os mesmos obtidos com os critérios adotados na simulação A. Os desvios médios obtidos para os resultados das frações líquidas das simulações B, C e D, em relação aos obtidos na simulação A foram, respectivamente, 0,35%, 0,63% e 0,26%. Se forem tomados apenas os resultados obtidos após 2 min estes mesmos desvios ficam em 0,06%, 0,05% e 0,02%. Tais constatações permitem afirmar que os parâmetros utilizados para o ∆t e CC nesta validação numérica (os mesmos da simulação A) são válidos e resultam em velocidades de fusão que reproduzem, com bom grau de precisão, os dados experimentais de Katsman [19].
360
361
'''Tabela 5''' – Resultados das simulações com variação de ∆t e CC.
362
363
{| style="width: 100%;border-collapse: collapse;"
364
|-
365
| style="border-top: 1pt solid black;text-align: center;"|
366
| style="border-top: 1pt solid black;"|<span style="text-align: center; font-size: 75%;">Simulação A</span>
367
| colspan='2'  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Simulação   B</span>
368
| colspan='2'  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Simulação   C</span>
369
| colspan='2'  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Simulação   D</span>
370
|-
371
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Tempo [min]</span>
372
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Fração</span>
373
374
<span style="text-align: center; font-size: 75%;">Líquida</span>
375
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Fração</span>
376
377
<span style="text-align: center; font-size: 75%;">Líquida</span>
378
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Desvio</span>
379
380
<span style="text-align: center; font-size: 75%;">[%]</span>
381
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Fração</span>
382
383
<span style="text-align: center; font-size: 75%;">Líquida</span>
384
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Desvio</span>
385
386
<span style="text-align: center; font-size: 75%;">[%]</span>
387
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Fração</span>
388
389
<span style="text-align: center; font-size: 75%;">Líquida</span>
390
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Desvio</span>
391
392
<span style="text-align: center; font-size: 75%;">[%]</span>
393
|-
394
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,167</span>
395
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,044239</span>
396
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,045346</span>
397
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2,50</span>
398
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,046819</span>
399
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">5,51</span>
400
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,045670</span>
401
| style="text-align: center;vertical-align: top;"|3,13
402
|-
403
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2</span>
404
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,164532</span>
405
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,166829</span>
406
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1,40</span>
407
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,168182</span>
408
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">2,17</span>
409
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,164612</span>
410
| style="text-align: center;vertical-align: top;"|0,05
411
|-
412
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4</span>
413
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,236794</span>
414
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,236769</span>
415
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0,01</span>
416
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,236921</span>
417
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,05</span>
418
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,236669</span>
419
| style="text-align: center;vertical-align: top;"|-0,05
420
|-
421
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6</span>
422
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,297303</span>
423
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,297392</span>
424
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,03</span>
425
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,297170</span>
426
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">-0,04</span>
427
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,297167</span>
428
| style="text-align: center;vertical-align: top;"|-0,05
429
|-
430
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">8</span>
431
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,336182</span>
432
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,336224</span>
433
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,01</span>
434
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,336340</span>
435
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,05</span>
436
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,336330</span>
437
| style="text-align: center;vertical-align: top;"|0,04
438
|-
439
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">10</span>
440
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,375469</span>
441
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,374736</span>
442
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0,20</span>
443
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,374837</span>
444
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">-0,17</span>
445
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,375104</span>
446
| style="text-align: center;vertical-align: top;"|-0,10
447
|-
448
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12</span>
449
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,409504</span>
450
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,411161</span>
451
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,40</span>
452
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,411253</span>
453
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,43</span>
454
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,410346</span>
455
| style="text-align: center;vertical-align: top;"|0,21
456
|-
457
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">14</span>
458
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,445273</span>
459
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,444849</span>
460
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0,10</span>
461
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,445006</span>
462
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">-0,06</span>
463
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,445186</span>
464
| style="text-align: center;vertical-align: top;"|-0,02
465
|-
466
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">16</span>
467
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,478774</span>
468
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,480461</span>
469
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,35</span>
470
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,479022</span>
471
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,05</span>
472
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,478816</span>
473
| style="text-align: center;vertical-align: top;"|0,01
474
|-
475
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">18</span>
476
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,511701</span>
477
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,511500</span>
478
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0,04</span>
479
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,511609</span>
480
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">-0,02</span>
481
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,511792</span>
482
| style="text-align: center;vertical-align: top;"|0,02
483
|-
484
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">20</span>
485
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,541824</span>
486
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,542962</span>
487
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,21</span>
488
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,543028</span>
489
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,22</span>
490
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,542620</span>
491
| style="text-align: center;vertical-align: top;"|0,15
492
|-
493
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">22</span>
494
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,570371</span>
495
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,570318</span>
496
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0,01</span>
497
| style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,570325</span>
498
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">-0,01</span>
499
| style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,570367</span>
500
| style="text-align: center;vertical-align: top;"|0,00
501
|-
502
| style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">24</span>
503
| style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,6004645</span>
504
| style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,600446</span>
505
| style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,00</span>
506
| style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,600478</span>
507
| style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,00</span>
508
| style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0,600549</span>
509
| style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0,01
510
|}
511
512
A presente validação consistiu em uma análise do tipo transiente em um espaço planar bidimensional. Para solução do problema foram adotados os seguintes modelos: volume de fluido (VOF) como modelo multifásico, escoamento laminar como modelo de viscosidade e modelo de solidificação e derretimento (solidification & melting). Os resultados parciais da validação foram salvos a cada 600 passos de tempo e o número máximo de iterações para atingir a convergência, em cada passo de tempo, foi definido como 1000.
513
514
Além da malha padrão, foram executadas simulações com malhas adicionais, contendo um número de elementos 30% e 60% acima da malha padrão. As curvas mostradas na Figura 3 representam a fração líquida do PCM para as três malhas citadas. Os resultados obtidos com a malha padrão apresentaram um desvio máximo de 0,5% em relação aos obtidos com a malha mais refinada, significando que ela está adequada ao uso nas simulações deste trabalho.
515
516
Destaca-se que os resultados apresentados na Figura 3 referem-se a ∆T = 30 <sup>o</sup>C. Esta condição implica em processos térmicos e fluidodinâmicos mais intensos do que aqueles obtidos com ∆T = 10 <sup>o</sup>C, exigindo assim malhas computacionais mais refinadas. Portanto, uma vez feita a validação da malha computacional com o processo mais intenso, os menos intensos também estão assegurados.
517
518
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
519
520
[[Image:draft_Samper_486932730-image3.png|center|456px]]
521
</div>
522
523
'''Figura 3''' – Validação da malha computacional com ∆T = 30 <sup>o</sup>C e C = 10<sup>7</sup>.
524
525
==3 Resultados e Discussão==
526
527
Nesta seção serão apresentados os resultados da validação numérica e do estudo dos PCMs selecionados.
528
529
===3.1 Validação numérica===
530
531
A Figura 4 mostra vetores velocidade em diferentes instantes de tempo e permite confirmar as observações feitas por Sparrow e Broadbent [20], ratificadas no trabalho de Shmueli et al. [17]. No início do processo a transferência de calor ocorre por condução da parede externa da geometria cilíndrica para a fase sólida do PCM, através de uma finíssima camada de PCM líquido. Na medida em que o processo de fusão avança, a convecção natural no PCM líquido se torna dominante, o que pode ser observado através dos vetores velocidade mostrados, os quais indicam intensidade da velocidade local. A evolução destes processos vai alterando o formato do PCM sólido, inicialmente cilíndrico, para cônico, e, paulatinamente, vai reduzindo o volume da fase sólida.
532
533
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
534
535
[[Image:draft_Samper_486932730-image4.png|center|504px]]
536
</div>
537
538
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
539
'''Figura 4''' – Vetores velocidade da validação numérica.</div>
540
541
Na Figura 5 são apresentados os resultados de fração líquida quando o RT 27 é submetido a temperaturas de 10<sup>o</sup>C e 30<sup>o</sup>C acima de sua T<sub>pc</sub>. Observa-se boa convergência entre os resultados numéricos de Shmueli et al. [17] e os obtidos na validação numérica do presente trabalho, tanto para ∆T = 10 <sup>o</sup>C quanto para ∆T = 30 <sup>o</sup>C. Comparando-se os resultados numéricos, do presente trabalho e de Shmueli et al. [17]) com os resultados experimentais de Katsman [19], pode se observar que: (a) para ∆T = 30 <sup>o</sup>C, os três resultados são muito semelhantes, qualitativa e quantitativamente, até a fração líquida de, aproximadamente, 0,8. Acima deste valor, a diferença média é de, aproximadamente, 7 % e (b) para ∆T = 10 <sup>o</sup>C, os resultados numéricos e experimentais apresentam valores muito próximos até fração líquida de 0,7. Acima deste valor, os resultados experimentais e numéricos, tanto de Shmueli quanto do presente trabalho, apresentaram, aproximadamente, 12% de desvio máximo entre eles. No entanto, apesar das diferenças máximas de 7 e 12 %, deve-se lembrar que este é um processo transiente e que estas diferenças ocorrem apenas em curto intervalo de tempo. Se considerado o tempo total do processo, a diferença média é de 2,5  e 6,1% para ∆T de 10 e 30 <sup>o</sup>C, respectivamente, o que pode ser considerado bastante razoável para processos de mudança de fase, principalmente para pequenos ∆T (10 ºC), quando o processo de investigação, tanto numérico quanto experimental, se torna cada vez mais exigente . Shmueli et al. [17] testou vários valores para a constante C. O uso da constante C=10<sup>7</sup> foi o que apresentou resultados mais próximos dos obtidos experimentalmente.
542
543
Ainda em relação à Figura 5, tomando-se como parâmetro de comparação o tempo necessário para atingir uma fração líquida de 0,6 com ∆T = 10<sup>o</sup>C (24 min), constata-se que o aumento para ∆T = 20<sup>o</sup>C reduz este tempo em 53,9% (11,06 min) e, usando-se ∆T = 30<sup>o</sup>C, este mesmo tempo é reduzido em 70,8% (7,02 min). Outra observação importante é que o tempo de fusão total, observado no processo numérico, foi reduzido em 66,7%, quando se eleva o &#x0394;T = 10 °C para 30 °C. Este valor ficou muito próximo dos 63,6%, observados nos experimentos de Katsman [19].
544
545
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
546
547
[[Image:draft_Samper_486932730-image5.png|center|456px]]
548
</div>
549
550
'''Figura 5''' – Fração líquida vs. tempo, para o RT 27, variando &#x0394;T.
551
552
Na Figura 6 são apresentados os perfis de fluxo de calor do RT 27 para &#x0394;T de 10 <sup>o</sup>C, 20 <sup>o</sup>C e 30 <sup>o</sup>C, com C=10<sup>7</sup>. Observa-se que o fluxo de calor é intenso nos instantes iniciais do processo de fusão, quando a troca térmica é dominada pelo processo de condução, e decresce exponencialmente na medida em que há um incremente do processo de convecção no PCM já liquefeito. Além disso, no tempo de 8 min, por exemplo, também pode-se observar que, quanto maior é o &#x0394;T maior é o fluxo de calor, correspondendo assim às expectativas para este tipo de processo.
553
554
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
555
<big>
556
[[Image:draft_Samper_486932730-image6.png|center|456px]]
557
</big></div>
558
559
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
560
'''Figura 6''' – Fluxo de calor vs. tempo, para o RT 27, variando &#x0394;T.</div>
561
562
===3.2 Resultados obtidos com a análise dos PCMs===
563
564
Os resultados apresentados nesta seção referem-se aos 12 tipos de PCMs citados apresentados na Tabela 2. Foi avaliada a influência da temperatura em relação aos tempos de fusão e também sobre fluxo de calor, considerando &#x0394;T de 10 <sup>o</sup>C, 20 <sup>o</sup>C e 30 <sup>o</sup>C acima da T<sub>pc</sub> de cada PCM. Para a comparação das simulações são investigados os tempos necessários para atingir as frações líquidas de 0,4, 0,6 , 0,8 e 1,0.
565
566
A Figura 7 apresenta a variação da fração líquida em função do tempo para o PCM RT 4, quando submetido a &#x0394;T de 10 <sup>o</sup>C, 20 <sup>o</sup>C e 30 <sup>o</sup>C. Tomandose como referência inicial o tempo necessário para que o RT 4 atinja a fração líquida de 0,6 (~18 min), com &#x0394;T = 10 <sup>o</sup>C, observa-se que o aumento do &#x0394;T para 20 <sup>o</sup>C e 30 <sup>o</sup>C faz com que este tempo seja reduzido em 57% e 72,7%, totalizando, aproximadamente, 7,74 min e 4,93 min, respectivamente.
567
568
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
569
570
[[Image:draft_Samper_486932730-image7.png|center|456px]]
571
</div>
572
573
'''Figura 7''' – Fração líquida vs. tempo, para o RT 4, variando &#x0394;T.
574
575
A variação da fração líquida com o tempo, para o PCM RT 35, é mostrada na Figura 8 quando este é submetido a &#x0394;T de 10 <sup>o</sup>C, 20 <sup>o</sup>C e 30 <sup>o</sup>C. O tempo necessário para atingir a fração líquida de 0,6 é de 19,4 min. 7,8 min e 4,9 min, apresentando redução de 59,8% e de 74,6%, respectivamente, em relação ao tempo obtido com &#x0394;T de 10 <sup>o</sup>C.
576
577
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
578
579
[[Image:draft_Samper_486932730-image8.png|center|456px]]
580
</div>
581
582
'''Figura 8''' – Fração líquida vs. tempo, para o RT 35, variando &#x0394;T.
583
584
Na Tabela 6 é apresentado um sumário dos tempos de fusão para a fração líquida de 0,6, assim como outras características dos PCM estudados. Observa-se que o incremento do &#x0394;T de 10 <sup>o</sup>C para 20 <sup>o</sup>C reduz o tempo de fusão, em média, em 56,6%. Já o incremento no &#x0394;T de 10 <sup>o</sup>C para 30 <sup>o</sup>C, reduz em 72,1% o tempo de fusão médio. <br/>
585
586
'''Tabela 6''' – Sumário dos tempos de fusão até a fração líquida de 0,6.
587
588
{| style="width: 100%;border: 1pt solid black;border-collapse: collapse;"
589
|-
590
| rowspan='3' style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;"|PCM
591
| rowspan='3' style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;"|L
592
593
[kJ kg<sup>-1</sup>]
594
| colspan='3'  style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;vertical-align: top;"|Tempo para atingir α=0,6 [min]
595
| colspan='2'  style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;vertical-align: top;"|Redução de Tempo [%]
596
|-
597
| rowspan='2' style="border: none;text-align: center;"|DT=10<sup>o</sup>C
598
| rowspan='2' style="border: none;text-align: center;"|DT=20<sup>o</sup>C
599
| rowspan='2' style="border: none;text-align: center;"|DT=30<sup>o</sup>C
600
| style="border: none;text-align: center;vertical-align: top;"|de DT=10<sup>o</sup>C
601
| style="border: none;text-align: center;vertical-align: top;"|de DT=10<sup>o</sup>C
602
|-
603
| style="border: none;text-align: center;vertical-align: top;"|para DT=20<sup>o</sup>C
604
| style="border: none;text-align: center;vertical-align: top;"|para DT=30<sup>o</sup>C
605
|-
606
| style="border: none;text-align: center;vertical-align: top;"|RT  4
607
| style="border: none;text-align: center;vertical-align: top;"|182
608
| style="border: none;text-align: center;vertical-align: top;"|18,05
609
| style="border: none;text-align: center;vertical-align: top;"|7,74
610
| style="border: none;text-align: center;vertical-align: top;"|4,93
611
| style="border: none;text-align: center;vertical-align: top;"|57,1%
612
| style="border: none;text-align: center;vertical-align: top;"|72,7%
613
|-
614
| style="border: none;text-align: center;vertical-align: top;"|RT 27
615
| style="border: none;text-align: center;vertical-align: top;"|179
616
| style="border: none;text-align: center;vertical-align: top;"|24,00
617
| style="border: none;text-align: center;vertical-align: top;"|11,06
618
| style="border: none;text-align: center;vertical-align: top;"|7,02
619
| style="border: none;text-align: center;vertical-align: top;"|53,9%
620
| style="border: none;text-align: center;vertical-align: top;"|70,8%
621
|-
622
| style="border: none;text-align: center;vertical-align: top;"|RT 35
623
| style="border: none;text-align: center;vertical-align: top;"|170
624
| style="border: none;text-align: center;vertical-align: top;"|19,40
625
| style="border: none;text-align: center;vertical-align: top;"|7,80
626
| style="border: none;text-align: center;vertical-align: top;"|4,90
627
| style="border: none;text-align: center;vertical-align: top;"|59,8%
628
| style="border: none;text-align: center;vertical-align: top;"|74,7%
629
|-
630
| style="border: none;text-align: center;vertical-align: top;"|RT 55
631
| style="border: none;text-align: center;vertical-align: top;"|172
632
| style="border: none;text-align: center;vertical-align: top;"|23,70
633
| style="border: none;text-align: center;vertical-align: top;"|10,75
634
| style="border: none;text-align: center;vertical-align: top;"|6,60
635
| style="border: none;text-align: center;vertical-align: top;"|54,6%
636
| style="border: none;text-align: center;vertical-align: top;"|72,2%
637
|-
638
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|RT 82
639
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|176
640
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|17,50
641
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|7,40
642
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|5,20
643
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|57,7%
644
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|70,3%
645
|}
646
647
Na Figura 9 são apresentadas as reduções médias obtidas nos tempos de fusão dos diversos PCMs para as frações líquidas de 0,4, 0,6, 0,8 e 1,0 (fusão completa). Estes resultados mostram que, para frações líquidas entre 0,4 e 0,8, as reduções percentuais obtidas nos tempos de fusão, situam-se em patamares muito próximos, sendo que as reduções médias são de 55,8% e 71,8%, quando incrementado o &#x0394;T= 10 <sup>o</sup>C para 20 <sup>o</sup>C e de 20 <sup>o</sup>C para 30 <sup>o</sup>C, respectivamente. Para α=1,0 e &#x0394;T=30 <sup>o</sup>C a redução média de 63,6% aproximou-se bastante do observado experimentalmente para o RT 27 (66,7%), conforme Figura 5.
648
649
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
650
651
[[Image:draft_Samper_486932730-image9.png|center|438px]]
652
</div>
653
654
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
655
'''Figura 9''' – Reduções médias dos tempos de fusão, variando α e &#x0394;T.</div>
656
657
No que diz respeito ao fluxo de calor (q¨), o comportamento verificado para os PCMs estudados foi similar ao mostrado na Figura 6, ou seja, o fluxo de calor é intenso nos instantes iniciais, quando a troca térmica é dominada pelo processo de condução e decai exponencialmente na medida em que o efeito convectivo avança. A Figura 10 apresenta o comportamento do fluxo de calor do RT 35 em função do tempo, para diferentes &#x0394;T.
658
659
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
660
661
[[Image:draft_Samper_486932730-image10.png|center|456px]]
662
</div>
663
664
'''Figura 10''' – Fluxo de calor vs. tempo, para o RT 35, variando &#x0394;T.
665
666
Na Figura 11 é apresentado o resultado do comportamento do fluxo de calor em função da fração líquida variando &#x0394;T, para o RT 35. Vale destacar que também aqui constatou-se que há similaridade do comportamento do fluxo de calor para todos os PCMs estudados, ou seja o fluxo de calor decresce exponencialmente na medida em que o PCM se liquefaz.
667
668
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
669
670
[[Image:draft_Samper_486932730-image11.png|center|456px]]
671
</div>
672
673
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
674
'''Figura 11''' – Fluxo de calor vs. fração líquida, para o RT 35, variando &#x0394;T.</div>
675
676
Na Tabela 7 são apresentados os valores do fluxo de calor dos PCMs estudados (em t=10 s) e os incrementos obtidos quando o &#x0394;T é aumentado de 10 °C para 20 °C e também para 30 °C. Os incrementos no fluxo de calor variaram entre 51,7% e 74,1% (aumentando o &#x0394;T de 10 °C para 20 °C) e entre 98,1% a 124,3% (aumentando o &#x0394;T de 10 °C para 30 °C).
677
678
<br/>
679
680
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
681
'''Tabela 7''' – Sumário dos incrementos obtidos em q¨ (10 s), variando &#x0394;T.</div>
682
683
{| style="width: 100%;border: 1pt solid black;border-collapse: collapse;"
684
|-
685
| rowspan='3' style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;"|PCM
686
| rowspan='3' style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;"|L
687
688
[kJ kg<sup>-1</sup>]
689
| colspan='3'  style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;vertical-align: top;"|q”<sub>máx</sub> [kW m<sup>-2</sup>]
690
| colspan='2'  style="border-top: 1pt solid black;border-left: none;border-bottom: none;border-right: none;text-align: center;vertical-align: top;"|Redução de Tempo [%]
691
|-
692
| rowspan='2' style="border: none;text-align: center;"|DT=10<sup>o</sup>C
693
| rowspan='2' style="border: none;text-align: center;"|DT=20<sup>o</sup>C
694
| rowspan='2' style="border: none;text-align: center;"|DT=30<sup>o</sup>C
695
| style="border: none;text-align: center;vertical-align: top;"|de DT=10<sup>o</sup>C
696
| style="border: none;text-align: center;vertical-align: top;"|de DT=10<sup>o</sup>C
697
|-
698
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|para DT=20<sup>o</sup>C
699
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|para DT=30<sup>o</sup>C
700
|-
701
| style="border-left: none;border-bottom: none;border-right: none;text-align: center;vertical-align: top;"|RT  4
702
| style="border-left: none;border-bottom: none;border-right: none;text-align: center;vertical-align: top;"|182
703
| style="border-left: none;border-bottom: none;border-right: none;text-align: center;vertical-align: top;"|4,65
704
| style="border-left: none;border-bottom: none;border-right: none;text-align: center;vertical-align: top;"|7,21
705
| style="border-left: none;border-bottom: none;border-right: none;text-align: center;vertical-align: top;"|9,21
706
| style="border-left: none;border-bottom: none;border-right: none;text-align: center;vertical-align: top;"|55,1%
707
| style="border-left: none;border-bottom: none;border-right: none;text-align: center;vertical-align: top;"|98,1%
708
|-
709
| style="border: none;text-align: center;vertical-align: top;"|RT 35
710
| style="border: none;text-align: center;vertical-align: top;"|170
711
| style="border: none;text-align: center;vertical-align: top;"|3,70
712
| style="border: none;text-align: center;vertical-align: top;"|6,44
713
| style="border: none;text-align: center;vertical-align: top;"|8,30
714
| style="border: none;text-align: center;vertical-align: top;"|74,1%
715
| style="border: none;text-align: center;vertical-align: top;"|124,3%
716
|-
717
| style="border: none;text-align: center;vertical-align: top;"|RT 55
718
| style="border: none;text-align: center;vertical-align: top;"|172
719
| style="border: none;text-align: center;vertical-align: top;"|4,18
720
| style="border: none;text-align: center;vertical-align: top;"|6,34
721
| style="border: none;text-align: center;vertical-align: top;"|8,43
722
| style="border: none;text-align: center;vertical-align: top;"|51,7%
723
| style="border: none;text-align: center;vertical-align: top;"|101,7%
724
|-
725
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|RT 82
726
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|176
727
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|3,83
728
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|6,66
729
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|8,47
730
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|73,9%
731
| style="border-top: none;border-left: none;border-bottom: 1pt solid black;border-right: none;text-align: center;vertical-align: top;"|121,2%
732
|}
733
734
<span id='_GoBack'></span>
735
Em termos de características de escoamento, todos os processos de fusão analisados são muito similares. Estes são representados na Figura 12(a-b), para o RT 4, com &#x0394;T=10 ºC no t=10 min. Com o aquecimento da parede vertical, o PCM, inicialmente sólido, muda de fase e passa a aumentar sua temperatura para valores acima da temperatura de fusão. Em decorrência disso, forças de empuxo se tornam presentes, decorrentes da variação da massa específica.
736
737
738
[[Image:draft_Samper_486932730-image12.png|center|588px]]
739
740
'''Figura 12''' – RT 4, &#x0394;T=10 ºC e t=10 min: (a) Campos de velocidade e (b) de linhas de corrente.
741
742
Nas Figuras 12(a-b) são apresentados os campos de velocidade e de linhas de corrente. Com o aumento da temperatura junto à parede vertical, e consequente redução da massa específica, ocorre um escoamento ascendente do PCM líquido. Esta condição pode ser verificada na Figura 12(a), sendo melhor visualizada nos Detalhes 1 e 2. Nestes detalhes também é possível identificar o escoamento descente junto à parcela sólida, em consequência da conservação da massa. Também é possível verificar a região cisalhante entre os escoamentos ascendente e descendente, onde a velocidade é zero. A possibilidade de visualização deste detalhe indica que o grau de refinamento da malha computacional está adequado para representar o processo. Já na Figura 12(b), onde estão representadas linhas de corrente, além da região referente à parcela sólida, pode-se observar duas regiões distintas com recirculação: a do PCM líquido (junto à parede lateral) e a do ar (no topo do domínio computacional). Neta última observa-se duas recirculações bem definidas, as quais são decorrentes do arrasto provocado pela recirculação do PCM líquido na interface entre os dois fluidos, indicando que ali também a malha computacional está adequada para o processo analisado.
743
744
==4 Conclusões==
745
746
Este trabalho teve como objetivo investigar numericamente o processo de fusão de diferentes tipos de PCMs em cavidades cilíndricas, inicialmente validando a metodologia e o modelo matemático com resultados da literatura e, posteriormente, avaliando e comparando o processo de fusão de diferentes PCMs quando submetidos a diferentes temperaturas acima da sua temperatura de mudança de fase (T<sub>pc</sub>).
747
748
A implementação do modelo numérico foi validada com resultados para PCM RT 27. Estes foram confrontados com resultados numéricos e experimentais da literatura, apresentando boa aproximação. Assim sendo, a metodologia, modelo matemático, malha computacional e demais parâmetros utilizados foram considerados aptos para os estudos das simulações.
749
750
Foram simulados os processos de fusão de quatro diferentes PCMs considerando-se três &#x0394;T (10, 20 e 30 ºC) acima da temperatura de fusão destes PCMs, totalizando 12 diferentes simulações. Para estas simulações as conclusões obtidas foram que o incremento do &#x0394;T, de 10 <sup>o</sup>C para 20 <sup>o</sup>C, ocasiona uma redução média geral nos tempos de fusão de 55,8% para atingir frações líquidas entre 0,4 e 0,8. Já o aumento do &#x0394;T de 10 <sup>o</sup>C para 30 <sup>o</sup>C resulta em uma redução média geral do tempo de fusão de 71,8%, para a mesma faixa de frações líquidas. As reduções de tempo médias obtidas para atingir a fusão completa dos PCMS (fração líquida 1,0) foram de 46,9% (de &#x0394;T de 10 <sup>o</sup>C para 20 <sup>o</sup>C) e de 63,6% (de &#x0394;T de 10 <sup>o</sup>C para 30 <sup>o</sup>C). Este último resultado mostrou-se similar ao obtido experimentalmente por Katsman [19], no qual a redução de tempo para fusão completa do RT 27 foi de 66,6% (de &#x0394;T de 10 <sup>o</sup>C para 30 <sup>o</sup>C).
751
752
No que diz respeito ao fluxo de calor (q¨), verificou-se que a utilização de &#x0394;T mais elevado propicia um forte incremento no fluxo de calor inicial, contribuindo no processo de troca térmica e reduzindo o tempo de fusão do PCM.
753
754
<br/>
755
756
==Apêndice A==
757
758
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
759
'''Tabela A1''' – Variáveis utilizadas e suas respectivas unidades.</div>
760
761
{| style="width: 100%;border: 1pt solid black;border-collapse: collapse;"
762
|-
763
| style="text-align: center;vertical-align: top;"|<br/>Símbolo
764
| style="text-align: center;vertical-align: top;"|Significado
765
| style="text-align: center;vertical-align: top;"|Unidade
766
|-
767
| style="text-align: center;vertical-align: top;"|C
768
| style="vertical-align: top;"|Constante Mushy Zone
769
| style="text-align: center;vertical-align: top;"|kg m<sup>-3</sup> s<sup>-1</sup>
770
|-
771
| style="text-align: center;vertical-align: top;"|CC
772
| style="vertical-align: top;"|Critério de convergência
773
| style="text-align: center;vertical-align: top;"|-
774
|-
775
| style="text-align: center;vertical-align: top;"|<big>''c<sub>p</sub>''</big>
776
| style="vertical-align: top;"|Calor específico do fluido a pressão constante
777
| style="text-align: center;vertical-align: top;"|J kg<sup>-1</sup> K<sup>-1</sup>
778
|-
779
| style="text-align: center;vertical-align: top;"|<big>''D''</big>
780
| style="vertical-align: top;"|Diâmetro do cilindro de PCM
781
| style="text-align: center;vertical-align: top;"|mm
782
|-
783
| style="text-align: center;vertical-align: top;"|D<sub>m</sub>
784
| style="vertical-align: top;"|Diâmetro médio
785
| style="text-align: center;vertical-align: top;"|mm
786
|-
787
| style="text-align: center;vertical-align: top;"|g
788
| style="vertical-align: top;"|Aceleração da gravidade
789
| style="text-align: center;vertical-align: top;"|m s<sup>-2</sup>
790
|-
791
| style="text-align: center;vertical-align: top;"|''h''
792
| style="vertical-align: top;"|Entalpia específica
793
| style="text-align: center;vertical-align: top;"|J kg<sup>-1</sup>
794
|-
795
| style="text-align: center;vertical-align: top;"|''h<sub>s</sub>''
796
| style="vertical-align: top;"|Entalpia sensível
797
| style="text-align: center;vertical-align: top;"|J kg<sup>-1</sup>
798
|-
799
| style="text-align: center;vertical-align: top;"|''h<sub>ref</sub>''
800
| style="vertical-align: top;"|Entalpia de referência
801
| style="text-align: center;vertical-align: top;"|J kg<sup>-1</sup>
802
|-
803
| style="text-align: center;vertical-align: top;"|''H<sub>1</sub>''
804
| style="vertical-align: top;"|Altura da cavidade estudada contendo PCM
805
| style="text-align: center;vertical-align: top;"|mm
806
|-
807
| style="text-align: center;vertical-align: top;"|''H<sub>2</sub>''
808
| style="vertical-align: top;"|Altura da cavidade estudada contendo ar
809
| style="text-align: center;vertical-align: top;"|mm
810
|-
811
| style="text-align: center;vertical-align: top;"|''k''
812
| style="vertical-align: top;"|Condutividade térmica
813
| style="text-align: center;vertical-align: top;"|W m<sup>-1</sup> K<sup>-1</sup>
814
|-
815
| style="text-align: center;vertical-align: top;"|''L''
816
| style="vertical-align: top;"|Calor latente de fusão
817
| style="text-align: center;vertical-align: top;"|J kg<sup>-1</sup>
818
|-
819
| style="text-align: center;vertical-align: top;"|p
820
| style="vertical-align: top;"|Pressão
821
| style="text-align: center;vertical-align: top;"|Pa
822
|-
823
| style="text-align: center;vertical-align: top;"|''q''
824
| style="vertical-align: top;"|Taxa de transferência de calor
825
| style="text-align: center;vertical-align: top;"|W
826
|-
827
| style="text-align: center;vertical-align: top;"|''q"''
828
| style="vertical-align: top;"|Fluxo de calor
829
| style="text-align: center;vertical-align: top;"|W m<sup>-</sup>²]
830
|-
831
| style="text-align: center;vertical-align: top;"| <math>\rightarrow {S}</math>
832
| style="vertical-align: top;"|Termo fonte de momento
833
| style="text-align: center;vertical-align: top;"|Pa m<sup>-1</sup>
834
|-
835
| style="text-align: center;vertical-align: top;"|''t''
836
| style="vertical-align: top;"|Tempo
837
| style="text-align: center;vertical-align: top;"|s
838
|-
839
| style="text-align: center;vertical-align: top;"|''T''
840
| style="vertical-align: top;"|Temperatura
841
| style="text-align: center;vertical-align: top;"|<sup>o</sup>C, K
842
|-
843
| style="text-align: center;vertical-align: top;"|''T<sub>s</sub>''
844
| style="vertical-align: top;"|Temperatura do PCM na fase sólida
845
| style="text-align: center;vertical-align: top;"|<sup>o</sup>C, K
846
|-
847
| style="text-align: center;vertical-align: top;"|''T<sub>l</sub>''
848
| style="vertical-align: top;"|Temperatura do PCM na fase líquida
849
| style="text-align: center;vertical-align: top;"|<sup>o</sup>C, K
850
|-
851
| style="text-align: center;vertical-align: top;"|''T<sub>pc</sub>''
852
| style="vertical-align: top;"|Temperatura de mudança de fase
853
| style="text-align: center;vertical-align: top;"|<sup>o</sup>C, K
854
|-
855
| style="text-align: center;vertical-align: top;"|''V''
856
| style="vertical-align: top;"|Volume
857
| style="text-align: center;vertical-align: top;"|m<sup>3</sup>
858
|-
859
| style="text-align: center;vertical-align: top;"| <math>\rightarrow {V}</math>
860
| style="vertical-align: top;"|Vetor velocidade
861
| style="text-align: center;vertical-align: top;"|m/s
862
|-
863
| style="text-align: center;vertical-align: top;"|''&#x03b1;<sub>n</sub>''
864
| style="vertical-align: top;"|Fração líquida na enésima célula computacional
865
| style="text-align: center;vertical-align: top;"|-
866
|-
867
| style="text-align: center;vertical-align: top;"|''&#x03b1;(t)''
868
| style="vertical-align: top;"|Fração líquida do PCM em um instante ¨t¨
869
| style="text-align: center;vertical-align: top;"|-
870
|-
871
| style="text-align: center;vertical-align: top;"|
872
| style="vertical-align: top;"|Variação
873
| style="text-align: center;vertical-align: top;"|-
874
|-
875
| style="text-align: center;vertical-align: top;"| <math>\epsilon</math> 
876
| style="vertical-align: top;"|Constante numérica (0,001)
877
| style="text-align: center;vertical-align: top;"|-
878
|-
879
| style="text-align: center;vertical-align: top;"|ρ
880
| style="vertical-align: top;"|Massa específica
881
| style="text-align: center;vertical-align: top;"|kg m<sup>-3</sup>
882
|-
883
| style="text-align: center;vertical-align: top;"|''&#x03bc;''
884
| style="vertical-align: top;"|Viscosidade dinâmica
885
| style="text-align: center;vertical-align: top;"|Pa s
886
|-
887
| style="text-align: center;vertical-align: top;"| <math>\gamma</math> 
888
| style="vertical-align: top;"|Fração líquida
889
| style="text-align: center;vertical-align: top;"|-
890
|}
891
892
<br/>
893
894
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
895
'''Apêndice B'''</div>
896
897
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
898
''' PROPRIEDADES DOS PCMs ESTUDADOS:'''</div>
899
900
Para a execução deste trabalho foram consideradas as características técnicas disponíveis no website da empresa Rubitherm Technologies GmbH ([http://www.rubitherm.com www.rubitherm.com]), abaixo transcritas.
901
902
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
903
'''Tabela B1''' – Propriedades termodinâmicas dos PCMs estudados</div>
904
905
{| style="width: 100%;border-top: 1pt solid black;border-bottom: 1pt solid black;border-collapse: collapse;"
906
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
907
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Propriedade'''
908
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''RT 4'''
909
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''RT 35'''
910
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''RT 55'''
911
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''RT 82'''
912
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
913
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|Temperatura de Fusão [<sup>o</sup>C]
914
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2-4
915
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|29-36
916
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|51-57
917
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|77-85
918
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
919
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|L [kJ kg<sup>-1</sup>]
920
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|182 (-5 a 10 <sup>o</sup>C)
921
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|170 (26 a 41 <sup>o</sup>C)
922
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|172 (48 a 63 <sup>o</sup>C)
923
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|176 (70 a 85 <sup>o</sup>C)
924
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
925
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|c<sub>p</sub> [kJ kg<sup>-1</sup> K<sup>-1</sup>]
926
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2
927
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2
928
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2
929
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2
930
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
931
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|K [W m<sup>-1</sup> K<sup>-1</sup>]
932
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0,2
933
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0,2
934
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0,2
935
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0,2
936
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
937
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|Expansão volumétrica [%]
938
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|12,5
939
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|12,5
940
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|14
941
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|12,5
942
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
943
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|Flash point [<sup>o</sup>C]
944
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|110
945
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|167
946
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|>200
947
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|>200
948
|}
949
950
Através de contatos feitos via e-mail, com a direção da empresa Rubitherm, foram coletados, adicionalmente, os seguintes dados:
951
952
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
953
'''Tabela B2''' – Massas específicas dos PCMs estudados</div>
954
955
{| style="width: 100%;border-top: 1pt solid black;border-bottom: 1pt solid black;border-collapse: collapse;"
956
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
957
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;vertical-align: top;"|'''Propriedade'''
958
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;vertical-align: top;"|'''RT 4'''
959
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;vertical-align: top;"|'''RT 35'''
960
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;vertical-align: top;"|'''RT 55'''
961
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;vertical-align: top;"|'''RT 82'''
962
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
963
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|&#x0072; (sólido) [kg m<sup>-3</sup>]
964
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|880,0 (-15 <sup>o</sup>C)
965
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|860,0 (15<sup> o</sup>C)
966
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|880,0 (15 <sup>o</sup>C)
967
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|880,0 (15 <sup>o</sup>C)
968
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
969
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|&#x0072; (início da fusão) [kg m<sup>-3</sup>]
970
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|870,0 (2 <sup>o</sup>C)
971
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|850,0 (29 <sup>o</sup>C)
972
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|870,0 (51 <sup>o</sup>C)
973
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|870,0 (77 <sup>o</sup>C)
974
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
975
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|&#x0072; (término da fusão) [kg m<sup>-3</sup>]
976
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|870,0 (4 <sup>o</sup>C)
977
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|781,5 (36 <sup>o</sup>C)
978
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|870,0 (57 <sup>o</sup>C)
979
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|870,0 (85 <sup>o</sup>C)
980
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
981
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|&#x0072; (líquido) [kg m<sup>-3</sup>]
982
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|770,0 (15 <sup>o</sup>C)
983
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|770,0 (45 <sup>o</sup>C)
984
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|770,0 (80 <sup>o</sup>C) 
985
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|770,0 (90 <sup>o</sup>C)
986
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
987
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|&#x0072; (dado adicional) [kg m<sup>-3</sup>]
988
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|700,0 (40 <sup>o</sup>C)
989
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|700,0 (70 <sup>o</sup>C)
990
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|700,0 (100 <sup>o</sup>C)
991
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|700,0 (100 <sup>o</sup>C)
992
|}
993
994
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
995
'''Tabela B3''' – Viscosidade dos PCMs estudados</div>
996
997
{| style="width: 100%;margin: 1em auto 1em auto;border-top: 1pt solid black;border-bottom: 1pt solid black;border-collapse: collapse;"
998
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
999
| rowspan='2' style="border-top: 1pt solid black;border-bottom: none;text-align: center;"|'''PCM'''
1000
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;vertical-align: top;"|Temp. de Referência
1001
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;vertical-align: top;"|Viscosidade Cinética
1002
| style="border-top: 1pt solid black;border-bottom: none;text-align: center;vertical-align: top;"|Viscosidade Dinâmica
1003
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
1004
| style="border-top: none;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|[<sup>o</sup>C]
1005
| style="border-top: none;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|[mm<sup>2</sup> s<sup>-1</sup>]
1006
| style="border-top: none;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|[mPa s]
1007
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
1008
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''RT 4'''
1009
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|40
1010
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2,6
1011
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1,82
1012
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
1013
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''RT 35'''
1014
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|70
1015
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3,3
1016
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2,31
1017
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
1018
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''RT 55'''
1019
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|100
1020
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3,8
1021
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2,66
1022
|- style="border-top: 1pt solid black;border-bottom: 1pt solid black;"
1023
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''RT 82'''
1024
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|100
1025
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|4,5
1026
| style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3,15
1027
|}
1028
1029
==Referências==
1030
1031
{| style="width: 100%;"
1032
|-
1033
| style="text-align: right;vertical-align: top;"|[1]
1034
| style="vertical-align: top;"|U.S. DEPARTMENT OF ENERGY.Building Energy Data Book. Disponível em [http://buildingsdatabook.eren.doe.gov/TableView.aspx?table=1.1.3 http://buildingsdatabook.eren.doe.gov/TableView.aspx?table=1.1.3]. Acesso em 22 de março de 2015.
1035
|-
1036
| style="text-align: right;vertical-align: top;"|[2]
1037
| style="vertical-align: top;"|EUROPEAN COMISSION. Statistical Pocketbook 2014. Disponível em: [http://ec.europa.eu/transport/facts-fundings/statistics/pocketbook-2014_en.htm http://ec.europa.eu/transport/facts-fundings/statistics/pocketbook-2014_en.htm]. Acesso em 22 de março de 2015.
1038
|-
1039
| style="text-align: right;vertical-align: top;"|[3]
1040
| style="vertical-align: top;"|MINISTÉRIO DE MINAS E ENERGIA. Balanço Energético Nacional – Ano Base 2013. Disponível em: [https://ben.epe.gov.br/BENRelatorioFinal2013.aspx https://ben.epe.gov.br/BENRelatorioFinal2013.aspx]. Acesso em 22 de março de 2015.
1041
|-
1042
| style="text-align: right;vertical-align: top;"|[4]
1043
| style="vertical-align: top;"|AGYENIM, F.; HEWITT, N.; EAMES, P.; SMYTH, M. A review of materials, heat transfer and phase change problem formulation for latent heat thermal energy storage systems (LHTESS). Renewable and Sustainable Energy Reviews, v. 14, p. 615-628, 2010.
1044
|-
1045
| style="text-align: right;vertical-align: top;"|[5]
1046
| style="vertical-align: top;"|PIELICHOWSKA, K.; PIELICHOWSKI, K.; Phase change materials for termal energy storaga. Progress in Materials Science, v. 65, p. 67-123, 2014.
1047
|-
1048
| style="text-align: right;vertical-align: top;"|[6]
1049
| style="vertical-align: top;"|ZHOU, D.; ZHAO, C.Y.; TIAN, Y. Review on termal energy storage with phase change materials in building applications. Applied Energy, v. 92, p. 593-605, 2012.
1050
|-
1051
| style="text-align: right;vertical-align: top;"|[7]
1052
| style="vertical-align: top;"|AL-ABIDI, A.A.; MAT, S.B.; SOPIAN, K.; SULAIMAN, M.Y.; LIM, C.H.; MOHAMMAD, A.T. A. Review of thermal energy storage for air conditioning systems. Renewable and Sustainable Energy Reviews, v. 16, p. 5802-5819, 2012.
1053
|-
1054
| style="text-align: right;vertical-align: top;"|[8]
1055
| style="vertical-align: top;"|ORÓ, E.; GRACIA, A.; CASTELL, A. FARID, M.M.; CABEZA, L.F. Review on phase change materials (PCMs) for cold termal energy storage applications. Applied Energy, v. 99, p. 513-533, 2012.
1056
|-
1057
| style="text-align: right;vertical-align: top;"|[9]
1058
| style="vertical-align: top;"|MEHLING, H.; CABEZA, L. F. Heat and cold storage with PCM: An up to date introduction into basics and applications. 1. ed. Berlin: Springer, 2008.
1059
|-
1060
| style="text-align: right;vertical-align: top;"|[10]
1061
| style="vertical-align: top;"|MURRAY, R. E.; GROULX, D. Experimental study of the phase change and energy characteristics inside a cylindrical latent heat energy storage system: Part 1 consecutive charging and discharging. Renewable Energy, v. 62, p. 571-581, 2014.
1062
|-
1063
| style="text-align: right;vertical-align: top;"|[11]
1064
| style="vertical-align: top;"|LONGEON, M.; SOUPART A.; FOURMIGUÉ, J.F.; BRUCH, A.; MARTY, P. Experimental and numerical study of annular PCM storage in the presence of natural convection. Applied Energy, v. 112, p. 175-184, 2013.
1065
|-
1066
| style="text-align: right;vertical-align: top;"|[12]
1067
| style="vertical-align: top;"|ZENG, Y.; FAN, L.W.; XIAO, Y.Q.; YU, Z.T.; CEN, K.F. An experimental investigation of melting of nanoparticle-enhanced phase change materials (NePCMs) in a bottom-heated vertical cylindrical cavity. International Journal of Heat and Mass Transfer, v. 66, p. 111-117, 2013.
1068
|-
1069
| style="text-align: right;vertical-align: top;"|[13]
1070
| style="vertical-align: top;"|AL-ABIDI, A.A.; MAT, S.B.; SOPIAN, K.; SULAIMAN, M.Y.; MOHAMMAD, A.T.. Numerical study of PCM solidification in a triplex tube heat exchanger with internal and external fins. International Journal of Heat and Mass Transfer, v. 61, p. 684-695, 2013.
1071
|-
1072
| style="text-align: right;vertical-align: top;"|[14]
1073
| style="vertical-align: top;"|TAY, N.H.; BRUNO F.; BELUSKO, M. Experimental investigation of dynamic melting in a tube-in-tank PCM system. Applied Energy, v. 104, p. 137-148, 2013.
1074
|-
1075
| style="text-align: right;vertical-align: top;"|[15]
1076
| style="vertical-align: top;"|KALAISELVAM, S.; VEERAPPAN, M.; AARON, A.A.; INIYAN, S. Experimental and analytical investigation of solidification and melting of PCMs inside cylindrical encapsulation. International Journal of Thermal Sciences, v. 47, p. 858-874, 2008.
1077
|-
1078
| style="text-align: right;vertical-align: top;"|[16]
1079
| style="vertical-align: top;"|RUBITHERM. Disponível em [http://www.rubitherm.de/english/index.htm http://www.rubitherm.de/english/index.htm]. Acesso em: 07 de maio de 2014.
1080
|-
1081
| style="text-align: right;vertical-align: top;"|[17]
1082
| style="vertical-align: top;"|SHMUELI, H.; ZISKIND, G.; LETAN, R. Melting in a vertical cylindrical tube: Numerical investigation and comparison with experiments.  International Journal of Heat and Mass Transfer, v. 53, p. 4082-4091, 2010.
1083
|-
1084
| style="text-align: right;vertical-align: top;"|[18]
1085
| style="vertical-align: top;"|BRENT, A.D.; VOLLER, V.R. Enthalpy-porosity technique for modeling convection-diffusion phase change – application to the melting of a pure metal.  Numerical Heat Transfer, v. 13, p. 297-318, 1988.
1086
|-
1087
| style="text-align: right;vertical-align: top;"|[19]
1088
| style="vertical-align: top;"|KATSMAN, L. Investigation of phase change in cylindrical geometry with internal fins. M.Sc. Thesis, Heat Transfer Laboratory, Department of Mechanical Engineering, Ben-Gurion University of the Negev, Beer-Sheva, Israel, 2006.
1089
|-
1090
| style="text-align: right;vertical-align: top;"|[20]
1091
| style="vertical-align: top;"|SPARROW, E.M.; BROADBENT, J.A. Inward melting in a vertical tube which allows free expansion of the phase-change medium, ASME Journal of Heat Transfer, v. 104, p. 309-315, 1982. 
1092
|}
1093

Return to Estrazulas Oliveski 2016a.

Back to Top

Document information

Published on 03/01/18
Accepted on 19/09/16
Submitted on 11/01/16

Volume 34, Issue 1, 2018
DOI: 10.23967/j.rimni.2017.5.001
Licence: CC BY-NC-SA license

Document Score

0

Times cited: 1
Views 299
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?