Este trabalho descreve uma solução híbrida analítico‐numérica aplicando a técnica da transformada integral generalizada (TTIG) na transferência de calor por convecção forçada de um escoamento laminar fluido dinamicamente desenvolvido e termicamente em desenvolvimento de um fluido Newtoniano em dutos de seção transversal com formato de setor de anel circular sob condições de contorno de Dirichlet, considerando um perfil de temperatura uniforme na entrada. Para facilitar o tratamento analítico e a aplicação das condições de contorno uma transformação conforme foi utilizada, visando alterar o domínio para um sistema de coordenadas mais apropriado. Feito isso, a TTIG foi aplicada na equação da energia para obtenção do campo de temperatura. Os resultados numéricos foram obtidos para parâmetros térmicos de interesse, tais como: temperatura média de mistura, números de Nusselt local e médio e comprimento de entrada térmica. Estes resultados foram comparados, quando possível, com os valores disponíveis na literatura e apresentaram uma ótima concordância.
This work describes a hybrid analytical‐numerical solution employing the Generalized Integral Transform Technique (GITT) to forced convection heat transfer in hydrodynamically fully developed and thermally non‐developed Newtonian laminar flow inside annular sector ducts under Dirichlet boundary conditions, considering uniform temperature entrance profile. In order to facilitate the analytical treatment and the application of the boundary conditions, a Conformal Transform was utilized to change the domain into a more suitable coordinate system. Thereafter, the GITT was applied on the energy equation to obtain the temperature field. Numerical results were obtained for quantities of practical interest, such as bulk mean temperature, local and average Nusselt number, and thermal entry length. These results were compared, as much as possible, with the parameter values available in the literature and they presented a good agreement.
Processos de transferência de calor para escoamento interno de fluidos em dutos representam uma classe de problemas difusivo‐convectivos de grande interesse na engenharia. Na concepção de um equipamento, dutos de seção transversal circular são largamente empregados em razão da sua simplicidade construtiva. Desta forma, inúmeros trabalhos e investigações relacionados com esta classe de problemas já foram realizados e uma vasta documentação é encontrada na literatura [1–3]. Por sua vez, o escoamento de fluidos em dutos que apresentam seção transversal de geometria não convencional encontra aplicações mais restritas, porém, do ponto de vista analítico‐numérico, estes problemas sempre despertam ampla atenção, pois, devido à sua maior complexidade, a busca de soluções fomenta o desenvolvimento de novas metodologias e a construção de novas ferramentas computacionais.
Neste contexto, o presente trabalho trata do cálculo de parâmetros térmicos da transferência de calor por convecção forçada para o problema de escoamento laminar de fluidos Newtonianos em regime permanente, termicamente em desenvolvimento, com perfil de temperatura de entrada uniforme e submetido a condições de contorno do primeiro tipo (condições de Dirichlet), em dutos de seção transversal com formato de setor de anel circular.
Geralmente, problemas com essa geometria são resolvidos utilizando‐se do sistema de coordenadas cilíndricas. Entretanto, para facilitar o tratamento analítico, utiliza‐se uma transformação conforme apropriada, com o objetivo de transformar o domínio original em um retângulo no novo sistema de coordenadas. Além de facilitar a aplicação das condições de contorno, esta transformação permite escrever a equação diferencial da energia de forma mais simples. Para a obtenção do campo de temperatura do escoamento aplica‐se a técnica da transformada integral generalizada (TTIG) [4], pois, para este caso, a equação diferencial parcial (EDP) da equação da energia é transformada em um sistema de equações diferenciais ordinárias (EDO) de primeira ordem, que é de fácil solução. Este procedimento vem sendo aplicado com sucesso para a obtenção da solução híbrida analítico‐numérica de problemas difusivos e difusivo‐convectivos com geometrias convencionais e não convencionais [5–25].
2Modelagem matemáticaPara a formulação do problema considera‐se o escoamento laminar de fluidos Newtonianos em regime permanente, termicamente não desenvolvido, com perfil de temperatura de entrada uniforme e temperatura prescrita na parede. Além disso, considera‐se ainda que as propriedades do fluido permanecem constantes em todo o domínio e que os efeitos da dissipação viscosa e da condução axial são desprezíveis. Desta forma, a equação da energia é expressa por:
sendo que, ρ é a massa específica, cp é o calor específico a pressão constante, k é condutividade térmica, w é a velocidade, T é a temperatura e Γ e Ω são o contorno e o domínio, respectivamente, da seção transversal do duto analisado, conforme ilustrado na figura 1.
As condições de entrada e de contorno são dadas por:
sendo que, T0 é a temperatura da entrada e TP é a temperatura da parede do duto.
O perfil de velocidade wx,y do escoamento no interior do duto de setor de anel circular é obtido através da TTIG, conforme apresentado detalhadamente em [26].
2.1AdimensionalizaçãoO potencial temperatura e os demais parâmetros físicos e geométricos foram adimensionalizados por:
sendo que, Dh é o diâmetro hidráulico, Pe é o número de Pèclet, AS é a área da seção transversal, Per é o perímetro do contorno analisado e wméd é a velocidade média do fluido.
Dessa forma, a equação da energia pode ser reescrita na forma adimensional, como:
E as condições de entrada e de contorno adimensionalizadas são:
2.2Transformação de coordenadas
Para facilitar o tratamento analítico do problema considera‐se a seguinte transformação conforme:
sendo que, Re=re/Dh, Z=X+iY e ω=u+iv.
Esta relação permite transformar o domínio do setor de anel circular no plano (X,Y) em um domínio retangular no plano (u,v) conforme ilustrado na figura 2.
As relações de transformação de coordenadas são dadas por:
A transformação conforme dada pela equação (14) satisfaz as condições de Cauchy‐Riemann. Portanto, os coeficientes métricos hu, hv e o Jacobiano J(u,v) desta transformação são dados por:
Como pode ser observado, o arco externo do setor de anel é dado pela reta v=0 no novo sistema de coordenadas e o arco interno pela reta v=v0, com v0=ln (re/ri). A coordenada u corresponde com a definição dada para o ângulo Θ do sistema de coordenadas cilíndricas.
Com isso, a equação da energia, equação (11), reescrita no novo sistema de coordenadas gerado pela transformação conforme é expressa por:
sendo que,
As condições de entrada e de contorno neste novo sistema de coordenadas são dadas por:
2.3Aplicação da técnica da transformada integral generalizada
Para a solução da equação da energia no novo sistema de coordenadas aplica‐se a TTIG que, na sua essência, procederá à remoção das derivadas parciais de segunda ordem. Para esta finalidade, o potencial θ(u,v,Z) é escrito em termos de uma expansão em autofunções normalizadas obtidas de problemas auxiliares de autovalor para cada coordenada espacial. Neste sentido, consideram‐se os seguintes problemas auxiliares de autovalor:
com:com:
Os autovalores e as autofunções associados a estes problemas são dados, respectivamente, por:
As autofunções acima definidas são ortogonais que permitem o desenvolvimento do seguinte par transformada‐inversa:
sendo que Ki(u) são as autofunções normalizadas dadas por:com:sendo que Zm(v) são as autofunções normalizadas dadas por:com:
Para a determinação do potencial transformado efetuando‐se o produto interno das autofunções normalizadas Ki(u) e Zm(v) com a equação da energia, obtém‐se o seguinte sistema de EDO:
Os parâmetros Bijmn são integráveis e, consequentemente, conhecidos. A solução deste sistema de EDO permite a obtenção do potencial transformado, quando submetido à condição de entrada transformada dada por:
Para fins computacionais, o potencial transformado pode ser determinado numericamente quando a expansão é truncada em uma dada ordem m=M e n=N. Portanto, utilizando‐se da fórmula de inversão, determina‐se o potencial temperatura adimensional:
Obviamente, quanto maior N e M maior será a precisão dos resultados, sendo que, do ponto de vista computacional, o duplo somatório poderia ser resolvido como proposto por [27], permitindo a redução deste custo através do reordenamento dos termos e substituição por um somatório simples. No entanto, pela simplicidade do problema em questão, isso não se torna necessário.
2.4Parâmetros térmicos de interesse2.4.1Temperatura média de misturaA temperatura média de mistura em sua forma adimensional é expressa por:
com:sendo que, AS∗=AS/Dh2 é a área da seção transversal adimensional. Portanto, no plano (u,v), θméd é dado por:2.4.2Número de Nusselt
Os números de Nusselt local e o médio podem ser definidos, respectivamente, por:
2.4.3Comprimento de entrada térmica
O comprimento de entrada térmica é definido como sendo a posição em que o número de Nusselt local é 5% maior que o número de Nusselt na região em que o fluido está termicamente desenvolvido [28]. Dessa forma, tem‐se:
3Resultados e discussão
Para a obtenção dos resultados numéricos, o método de quadratura de Gauss foi utilizado para o cálculo das integrais envolvida nos coeficientes Bijmn e nos parâmetros físicos de interesse. Desta forma, é necessária, também, a determinação das autofunções e do Jacobiano nestes pontos de quadratura. Feito isso, o sistema de equações diferenciais foi resolvido para a determinação dos coeficientes do potencial transformado utilizando‐se da rotina DIVPAG da biblioteca IMSL Fortran [29] e do truncamento da expansão dada pela equação (45) para diversas ordens M e N.
Foi verificado, através de uma análise de convergência, que na região de entrada do escoamento, Z<0,001, a convergência é relativamente lenta, sendo necessário considerar truncamentos da série com ordem M=N=25 para a obtenção de, pelo menos, 3 dígitos de precisão. Para a região onde o escoamento já se encontra termicamente desenvolvido é observado que a série que determina o potencial temperatura converge mais rapidamente, sendo necessário considerar apenas 20 termos para cada direção para uma precisão em torno de 4 dígitos.
Na tabela 1 são apresentados os resultados obtidos para a análise de convergência da temperatura adimensional média de mistura, em função do ângulo Θ e da ordem N=M de truncamento da série, para uma razão de aspecto r*=ri/re=0,20 e para uma coordenada axial adimensionalizada Z=0,001. A tabela 2 mostra a convergência da temperatura adimensional média de mistura em função da razão de aspecto e da ordem N=M de truncamento da série para um ângulo Θ=330° e para Z=0,001.
Convergência da temperatura adimensional média de mistura, θméd(Z), considerando uma razão de aspecto r*=0,20 em função da ordem de truncamento da série e do ângulo Θ para Z=10−3
| Θ | Ordem N=M de truncamento da série | |||||
|---|---|---|---|---|---|---|
| 10 | 15 | 20 | 25 | 30 | Valor convergido | |
| 5° | 0,9397 | 0,9403 | 0,9414 | 0,9416 | 0,9416 | 0,9416 |
| 10° | 0,9431 | 0,9436 | 0,9437 | 0,9433 | 0,9433 | 0,9433 |
| 15° | 0,9454 | 0,9455 | 0,9450 | 0,9440 | 0,9440 | 0,9440 |
| 45° | 0,9494 | 0,9478 | 0,9459 | 0,9444 | 0,9444 | 0,9444 |
| 60° | 0,9495 | 0,9475 | 0,9456 | 0,9444 | 0,9444 | 0,9444 |
| 90° | 0,9489 | 0,9464 | 0,9449 | 0,9441 | 0,9441 | 0,9441 |
| 150° | 0,9472 | 0,9446 | 0,9432 | 0,9426 | 0,9426 | 0,9426 |
| 180° | 0,9465 | 0,9438 | 0,9423 | 0,9417 | 0,9417 | 0,9417 |
| 210° | 0,9457 | 0,9431 | 0,9416 | 0,9409 | 0,9409 | 0,9409 |
| 270° | 0,9445 | 0,9420 | 0,9405 | 0,9396 | 0,9396 | 0,9396 |
| 300° | 0,9440 | 0,9416 | 0,9400 | 0,9390 | 0,9390 | 0,9390 |
| 350° | 0,9432 | 0,9409 | 0,9394 | 0,9382 | 0,9382 | 0,9382 |
Convergência da temperatura adimensional média de mistura, θméd(Z), considerando um ângulo Θ=330° em função da ordem de truncamento da série e da razão de aspecto r* para Z=10−3
| r* | Ordem N=M de truncamento da série | |||||
|---|---|---|---|---|---|---|
| 10 | 15 | 20 | 25 | 30 | Valor convergido | |
| 0,05 | 0,9507 | 0,9469 | 0,9433 | 0,9412 | 0,9412 | 0,9412 |
| 0,10 | 0,9474 | 0,9439 | 0,9414 | 0,9401 | 0,9401 | 0,9401 |
| 0,20 | 0,9435 | 0,9412 | 0,9396 | 0,9385 | 0,9385 | 0,9385 |
| 0,30 | 0,9411 | 0,9395 | 0,9382 | 0,9371 | 0,9371 | 0,9371 |
| 0,40 | 0,9377 | 0,9381 | 0,9370 | 0,9359 | 0,9359 | 0,9359 |
| 0,50 | 0,9393 | 0,9368 | 0,9358 | 0,9349 | 0,9349 | 0,9349 |
| 0,60 | 0,9361 | 0,9353 | 0,9346 | 0,9338 | 0,9338 | 0,9338 |
| 0,70 | 0,9345 | 0,9339 | 0,9333 | 0,9327 | 0,9327 | 0,9327 |
| 0,80 | 0,9329 | 0,9323 | 0,9319 | 0,9314 | 0,9314 | 0,9314 |
| 0,90 | 0,9313 | 0,9308 | 0,9304 | 0,9300 | 0,9300 | 0,9300 |
Nas figuras 3–5 são apresentados os comportamentos da temperatura adimensional média de mistura, dos números de Nusselt local e médio, respectivamente, em função da razão de aspecto r*=ri/re e para os diversos ângulos Θ (5°, 15°, 30°, 60°, 90°, 180°, 270° e 350°). Como esperado, os parâmetros térmicos de interesse decrescem com o aumento da coordenada axial adimensionalizada Z independentemente do ângulo Θ.
Nas tabelas 3 e 4, são apresentados os resultados obtidos para o número de Nusselt limite (Nu∞) e o comprimento de entrada térmica (Lth), respectivamente, considerando diversas configurações geométricas dos dutos de setor de anel circular.
Número de Nusselt limite, Nu∞, para diversas configurações dos dutos analisados
| r* | 5° | 15° | 30° | 60° | 90° | 180° | 270° | 350° |
|---|---|---|---|---|---|---|---|---|
| 0,05 | 1,567 | 2,061 | 2,489 | 2,911 | 3,094 | 3,284 | 3,457 | 3,639 |
| 0,10 | 1,715 | 2,219 | 2,626 | 2,979 | 3,102 | 3,308 | 3,604 | 3,875 |
| 0,15 | 1,878 | 2,384 | 2,757 | 3,024 | 3,095 | 3,381 | 3,782 | 4,112 |
| 0,20 | 2,055 | 2,553 | 2,879 | 3,043 | 3,090 | 3,492 | 3,978 | 4,348 |
| 0,25 | 2,246 | 2,724 | 2,981 | 3,042 | 3,054 | 3,632 | 4,186 | 4,582 |
| 0,30 | 2,453 | 2,891 | 3,055 | 3,032 | 3,133 | 3,798 | 4,403 | 4,812 |
| 0,35 | 2,671 | 3,029 | 3,095 | 3,030 | 3,194 | 3,985 | 4,626 | 5,037 |
| 0,40 | 2,901 | 3,180 | 3,098 | 3,045 | 3,287 | 4,191 | 4,852 | 5,257 |
| 0,45 | 3,136 | 3,279 | 3,072 | 3,087 | 3,413 | 4,414 | 5,080 | 5,478 |
| 0,50 | 3,372 | 3,330 | 3,033 | 3,162 | 3,573 | 4,652 | 5,312 | 5,685 |
| 0,55 | 3,596 | 3,322 | 3,000 | 3,277 | 3,769 | 4,903 | 5,538 | 5,898 |
| 0,60 | 3,790 | 3,258 | 2,991 | 3,437 | 4,005 | 5,165 | 5,779 | 6,106 |
| 0,65 | 3,928 | 3,158 | 3,024 | 3,651 | 4,284 | 5,439 | 6,003 | 6,304 |
| 0,70 | 3,978 | 3,056 | 3,118 | 3,937 | 4,607 | 5,727 | 6,223 | 6,493 |
| 0,75 | 3,891 | 2,988 | 3,296 | 4,279 | 4,976 | 5,974 | 6,453 | 6,682 |
| 0,80 | 3,649 | 3,004 | 3,596 | 4,718 | 5,401 | 6,316 | 6,683 | 6,857 |
| 0,85 | 3,303 | 3,180 | 4,074 | 5,251 | 5,865 | 6,608 | 6,899 | 7,028 |
| 0,90 | 3,009 | 3,673 | 4,817 | 5,924 | 6,378 | 6,914 | 7,105 | 7,192 |
| 0,95 | 3,225 | 4,848 | 5,944 | 6,666 | 6,932 | 7,207 | 7,293 | 7,326 |
| 1,00 | 7,619 | 7,619 | 7,619 | 7,619 | 7,619 | 7,619 | 7,619 | 7,619 |
Comprimento de entrada térmica, Lth, para diversas configurações dos dutos analisados
| r* | 5° | 15° | 30° | 60° | 90° | 180° | 270° | 350° |
|---|---|---|---|---|---|---|---|---|
| 0,05 | 0,3309 | 0,1494 | 0,0774 | 0,0497 | 0,0458 | 0,0520 | 0,0599 | 0,0650 |
| 0,10 | 0,3242 | 0,1372 | 0,0714 | 0,0481 | 0,0453 | 0,0530 | 0,0607 | 0,0646 |
| 0,15 | 0,3096 | 0,1260 | 0,0675 | 0,0469 | 0,0452 | 0,0542 | 0,0613 | 0,0632 |
| 0,20 | 0,2861 | 0,1163 | 0,0639 | 0,0458 | 0,0453 | 0,0554 | 0,0614 | 0,0607 |
| 0,25 | 0,2648 | 0,1085 | 0,0610 | 0,0448 | 0,0458 | 0,0569 | 0,0607 | 0,0564 |
| 0,30 | 0,2425 | 0,1015 | 0,0577 | 0,0444 | 0,0466 | 0,0581 | 0,0584 | 0,0511 |
| 0,35 | 0,2237 | 0,0954 | 0,0537 | 0,0442 | 0,0478 | 0,0589 | 0,0548 | 0,0454 |
| 0,40 | 0,2089 | 0,0901 | 0,0498 | 0,0446 | 0,0494 | 0,0584 | 0,0498 | 0,0394 |
| 0,45 | 0,1914 | 0,0828 | 0,0468 | 0,0453 | 0,0512 | 0,0567 | 0,0442 | 0,0324 |
| 0,50 | 0,1790 | 0,0727 | 0,0448 | 0,0467 | 0,0534 | 0,0532 | 0,0372 | 0,0277 |
| 0,55 | 0,1689 | 0,0620 | 0,0439 | 0,0485 | 0,0559 | 0,0482 | 0,0316 | 0,0227 |
| 0,60 | 0,1567 | 0,0546 | 0,0436 | 0,0511 | 0,0573 | 0,0416 | 0,0257 | 0,0193 |
| 0,65 | 0,1421 | 0,0478 | 0,0441 | 0,0536 | 0,0568 | 0,0337 | 0,0215 | 0,0166 |
| 0,70 | 0,1158 | 0,0448 | 0,0454 | 0,0554 | 0,0536 | 0,0256 | 0,0183 | 0,0143 |
| 0,75 | 0,0823 | 0,0429 | 0,0484 | 0,0567 | 0,0462 | 0,0201 | 0,0149 | 0,0120 |
| 0,80 | 0,0596 | 0,0423 | 0,0528 | 0,0514 | 0,0334 | 0,0165 | 0,0121 | 0,0105 |
| 0,85 | 0,0491 | 0,0464 | 0,0567 | 0,0394 | 0,0228 | 0,0130 | 0,0101 | 0,0093 |
| 0,90 | 0,0438 | 0,0538 | 0,0497 | 0,0205 | 0,0158 | 0,0100 | 0,0088 | 0,0084 |
| 0,95 | 0,0473 | 0,0534 | 0,0213 | 0,0122 | 0,0093 | 0,0084 | 0,0081 | 0,0081 |
| 1,00 | 0,0080 | 0,0080 | 0,0080 | 0,0080 | 0,0080 | 0,0080 | 0,0080 | 0,0080 |
Os comportamento do número de Nusselt limite e do comprimento de entrada térmica, parametrizados em Θ, são mostrados nas figuras 6 e 7, respectivamente. Pode ser notado que os resultados obtidos para os diversos dutos de setor de anel circular convergem para um mesmo valor quando r*→ 1, Nu∞≈7,6 e Lth=0,0080, valores estes, correspondentes ao problema térmico do escoamento laminar em canais de placas planas paralelas com temperatura uniforme em ambas as placas [30,31]. É interessante ainda observar que, para pequenos ângulos (Θ<30°), à medida que r* diminui a geometria do setor de anel circular vai se aproximando daquela correspondente à um duto de seção transversal quadrada, onde, por exemplo, o número de Nusselt limite passa por um mínimo local (Θ≈30°). Ressalta‐se que, o número de Nusselt limite para r*→ 0, corresponde à geometria do setor circular e os valores obtidos se aproximam daqueles referentes aos respectivos ângulos Θ [32].
Finalmente, destaca‐se que os resultados obtidos para o número de Nusselt limite e apresentados na tabela 3 e na figura 6, quando comparados com os mostrados em [32] apresentam uma ótima concordância, gerando uma diferença inferior a 1%.
4ConclusõesPara a solução do problema estudado, devido ao formato geométrico não convencional do duto, as dificuldades inerentes para a aplicação das condições de contorno foram removidas utilizando‐se uma transformação conforme adequada capaz de representar a geometria do duto de setor de anel circular em um domínio com contorno de formato retangular. A aplicação da TTIG sobre a equação da energia gerou um sistema acoplado e infinito de EDO para o potencial temperatura transformado que foi resolvido numericamente truncando‐se a expansão em um número finito de termos.
Observou‐se que a convergência do potencial temperatura é lenta na região de entrada do escoamento (Z<0,001), sendo necessária uma ordem de truncamento relativamente alta para o sistema de equações diferenciais do potencial transformado para a obtenção de, pelo menos, 3 dígitos de precisão.
Parâmetros térmicos de interesse, tais como, temperatura média de mistura, números de Nusselt local e médio e comprimento de entrada térmica, foram calculados e comparados, quando possível, com os resultados disponíveis na literatura, para dutos de diversas configurações geométricas, tendo sido observada uma ótima concordância.
Finamente, observa‐se que a TTIG é aplicada com sucesso para a obtenção de solução de problemas difusivos e difusivo‐convectivos multidimensionais, ressaltando‐se, aqui, o presente problema, o qual muitas vezes não admite solução pelas técnicas analíticas clássicas.