Introdução: O tratamento térmico é um dos métodos mais empregados na indústria de alimentos e a utilização de modelagem e simulação matemática na transferência de calor vêm se mostrando uma poderosa ferramenta, promovendo redução de custo e tempo de resposta no desenvolvimento de novos produtos, equipamentos e processos. Objetivo: O objetivo deste trabalho foi avaliar a transferência de calor, simulando sistemas reais, através da linguagem de Programação Python e software OpenFoam. Metodologia: Foram estudadas duas formas de calcular a condução de calor pelo método dos volumes finitos, primeiramente utilizando a linguagem de programação Python, comparando a otimização dos resultados através da complexidade das malhas e variando a quantidade de volumes de controle. Posteriormente, foi realizada a simulação através do software OpenFOAM, aplicando malhas com volumes elementares de diferentes quantidades e geometrias. Resultados: Para a simulação utilizando o Python, foi utilizado uma coluna contendo duas faces isoladas termicamente, uma terceira sob ação de fluxo de calor e uma quarta sujeita a convecção. Para resolução literal, a geometria foi dividida em 9 volumes finitos, comprometendo a precisão dos resultados, obtendo valores pontuais com elevada amplitude térmica. Ampliando o número de volumes para 1.750, foram plotadas curvas de escoamento de temperatura, sendo possível verificar como a condutividade térmica influencia na transferência de calor, indicando os locais de maior acúmulo térmico. Desta forma, comparou-se o escoamento em materiais com condutividade diferente, utilizando o mesmo código de programação. Para análise com OpenFoam, estudou-se o escoamento incompressível em um cotovelo contendo duas entradas, com diferentes velocidades, e uma saída, analisando-se três sistemas, o primeiro com malha não ortogonal e geometrias triangulares e o segundo e terceiro com malhas ortogonais divididas em blocos, diferenciando-as em números de volume. Após simulação, a geometria ortogonal apresentou melhores resultados de visualização, onde as camadas do regime laminar ficaram melhor definidas. Conclusão: A aplicação de um número maior de volumes de controle possibilita um resultado mais próximo às situações térmicas reais e as curvas de escoamento de temperatura se tornam mais claras e definidas com o uso de malhas ortogonais, sendo mais nítidas quando utilizados um maior número de elementos.