Resolução de Uma Questão Sobre Transferência de Calor A Partir de Uma Superfície Aletada
Resolução de Uma Questão Sobre Transferência de Calor A Partir de Uma Superfície Aletada
Resolução de Uma Questão Sobre Transferência de Calor A Partir de Uma Superfície Aletada
INTRODUÇÃO
Uma equação diferencial parcial é aquela a qual envolve uma ou mais variáveis
independentes e derivadas parciais de uma função, ou seja, uma variável dependente. Essas
equações são um modelo matemático que visam representar algum determinado fenômeno da
natureza, geralmente de forma idealizada e simplificada, cujas as soluções do modelo devem
apresentar um comportamento compatível com as propriedades do problema modelado.
(FRANCO, 2006)
Além disso, pela complexidade dos problemas muitas vezes se utiliza softwares para
auxiliarem na resolução, analise e também em reduzir o tempo de resolução desses
problemas. Dentre esses softwares se encontra o Matlab, que é bastante versátil em cálculos
matemáticos, visualização e gráficos, modelagens e simulações, além de outras funções.
Dessa forma, esse software se torna uma ferramenta adequada para implementar algoritmos
necessários para a resolução de EDP’s. Dessa forma, o presente trabalho visa mostrar a
resolução de uma EDP’s em relação a uma aleta, utilizando o método das linhas no software
de computação Matlab. (GILAT, 2009; YEO, 2017)
2. METODOLOGIA
2.1. Método das Linhas
O método das linhas que tem como objetivo resolver EDP’s através da discretização
das variáveis espaciais, as quais são aproximadas por uma diferença finita, enquanto uma
dessas variáveis permanece continua, geralmente o tempo, assim transformando a EDP em
um sistema de EDO’s. (DAVIS, 2013; SAUCEZ et al, 2001; YEO, 2017)
O calculo das diferenças finitas permite que seja possível tomar uma equação
diferencial e integrá-la numericamente calculando os valores da função em um número do
discreto de pontos. Também sendo possível fazer o inverso, caso um conjunto de valores
finitos estiver disponível, como dados experimentais, estes podem ser diferenciados, ou
integrados, usando o cálculo de diferenças finitas. (CONSTANTINIDES; MOSTOUFI, 1999)
d
2 2 f ( x i )−5 f ( xi +1 ) +4 f ( x i+2 )−f ( x i+3 )
2
f ( x i )= 2
dx h
d
2 f ( x i +1 )−2 f ( x i )−f ( x i−1 )
2
f ( x i )= 2
dx h
d
2 2 f ( x i )−5 f ( xi−1 ) + 4 f ( x i−2 )−f ( x i−3 )
2
f ( x i )= 2
dx h
dT
Qcond =−kA
dx
Desse modo então
ρcp A
dT d
=
dt dx (
−kA
dT
dx )
+hc p (T −T ∞)
dT k d T
2
hc p
= 2
− (T −T ∞ )
dt ρ c p d x ρ cp A
Considerando
k
α=
ρcp
hc p
β=
ρcp A
Logo
2
dT d T
=α 2
−β (T −T ∞)
dt dx
d T hc p
2
2
−
kA
( T −T ∞ )=0
dx
2.3. Implantação no Matlab
3. DISCUSSÃO E RESULTADOS
No gráfico é possível observar que quanto maior o comprimento da aleta, nesse caso
quanto mais perto chega do valor de 0,1 mm, menor é a temperatura do processo, quando a
aleta tem 0,1mm de comprimento se tem o menor valor de temperatura apresentado, igual a
293K. Já quando não há a presença de uma aleta, ou seja, no valor de 0 temos o maior valor
de temperatura. Assim mostra que o resultado está correto fisicamente já que a função de uma
aleta é aumentar a transferência de calor.
4. CONCLUSÃO
DAVIS, Mark E. Numerical methods and modeling for chemical engineers. Courier
Corporation, 2013.
YEO, Yeong Koo. Chemical engineering computation with MATLAB®. CRC Press,
2017.