Analisa-se o desempenho de armazenadores térmicos por calor latente com geometria cilíndrica de eixo vertical. O problema é formulado em coordenadas cilíndricas com as equações de conservação de movimento acopladas a equação da energia para se considerar o efeito da convecção natural nos processos de fusão e solidificação. Essas equações são resolvidas numericamente utilizando-se a técnica dos volumes de controle finitos desenvolvida por Patankar e o modelo proposto por Bonacina para simular o efeito do calor latente de fusão. Desenvolve-se um programa computacional para resolver os sistemas de equações lineares resultantes da discretização. Os resultados para o calor armazenado, para o fluxo local de transferência de calor e para a taxa global de transferência de calor são apresentados numa forma adimensional para diversos regimes de operação do armazenador e comparados com resultados experimentais. São analisadas ainda as influências de parâmetros geométricos e de transferência de calor no seu comportamento.