Análise da qualidade de tensões obtidas na simulação de escoamentos de fluidos viscoelásticos usando a formulação log-conformação

Abstract

Uma das mais recentes abordagens propostas na literatura para tratar o problema do alto número de Weissenberg (We) é a Formulação Log-Conformação (FLC). Nesta formulação, a equação constitutiva viscoelástica utilizada é reescrita em termos de uma variável Ψ, que é o logaritmo do tensor conformação. Apesar do potencial de aplicação da FLC, pouca atenção tem sido dirigida para análise da acurácia da solução obtida para o campo de tensões quando se utiliza esta formulação. Assim, o objetivo do presente trabalho foi estudar a acurácia da solução obtida pela FLC na análise de escoamentos de fluidos viscoelásticos usando duas geometrias padrão de estudo: placas paralelas e cavidade quadrada com tampa móvel. Primeiramente, a FLC foi implementada no pacote de CFD OpenFOAM. Em seguida foram verificados os limites do número de Weissenberg na formulação numérica padrão (Welim,P), onde para a geometria de placas paralelas foi encontrado Welim,P = 0,3 e para a geometria da cavidade quadrada com tampa móvel foi encontrado Welim,P = 0,8. Depois o código implementado foi aplicado em ambas as geometrias, comparando-se a solução obtida pela FLC com aquela da formulação padrão na faixa de We < Welim,P. Os resultados obtidos na geometria de placas paralelas apresentaram boa concordância com a solução padrão e solução analítica. Para a geometria da cavidade quadrada com tampa móvel, que não possui solução analítica, boa concordância dos resultados também foi observada em comparação com a solução padrão. Posteriormente foram comparados os resultados obtidos pela FLC na faixa de We > Welim,P. Na geometria de placas paralelas, além da boa concordância com a solução analítica, obteve-se convergência em todos os casos estudados neste trabalho, com o maior número de Weissenberg utilizado sendo igual a 8 Os resultados da geometria da cavidade quadrada com tampa móvel também apresentaram boa concordância em comparação com dados da literatura, porém a convergência foi obtida até para We = 2. Com respeito à comparação das formulações numéricas com a solução analítica, feita apenas na geometria de placas paralelas, foi observado um erro máximo de 7,57% na solução padrão e de 12,33% na FLC. Em relação à análise da qualidade das tensões usando os resíduos da equação constitutiva viscoelástica como critério de acurácia, foi verificado nas duas geometrias que os valores de tensão obtidos usando a FLC são menos acurados que aqueles obtidos pela formulação explícita no tensor das tensões nos casos em que esta última converge. Também foi observado que a acurácia diminui com o aumento do We. Esse efeito pôde ser melhor notado na geometria de placas paralelas. Uma razão para a perda de acurácia da tensão provavelmente ocorre devido à natureza matemática da transformação algébrica inversa de Ψxx para τxx. O novo solver implementado neste trabalho apresentou convergência e soluções corretas para as duas geometrias, logo foi implementado corretamente. Ele também potencializa o solver de partida viscoelastiFluidFoam ao estender simulações para uma faixa maior do número de Weissenberg.

A recent approach proposed in the literature to deal with the High Weissenberg Number Problem is the Log-Conformation formulation (LCF). In this formulation the viscoelastic constitutive equation is rewritten in terms of the logarithm of the conformation tensor Ψ. Despite the great potential application of the LCF, little attention has been given in the literature to the accuracy of the obtained stress fields. The purpose of this work was to study the solution obtained by LCF in the analysis of viscoelastic flows using two benchmark geometries: parallel plates and lid driven cavity. Firstly, the LCF was implemented in the OpenFOAM CFD package. Then, the limits of Weissenberg number for the standard numerical formulation (Welim,P) were verified, obtaining Welim,P = 0.3 for the parallel plates and Welim,P = 0.8 for the lid driven cavity. When comparing the solution obtained by the LCF with that of the standard formulation in a range of We < Welim,P, the results obtained for the parallel plates geometry showed good agreement with the standard solution and the analytical solution. For the lid driven cavity geometry, for which there is not analytical solution, good agreement with the standard solution was also observed. For We > Welim,P in the parallel plates geometry, in addition to the good agreement with the analytical solution, it was possible to obtain convergence in all the cases studied in this work, with the largest number of Weissenberg used being equal to 8 The results of the lid driven cavity geometry also presented good agreement in comparison with literature data, but convergence was obtained up to We = 2. With respect to the comparison of the numerical formulations with the analytical solution for the parallel plates geometry, a maximum error of 7.57% was observed in the standard solution and of 12.33% in the LCF. When using the residues of the viscoelastic constitutive equation as a criterion of accuracy, it was verified that for the two geometries the stress values obtained using the LCF were less accurate than those obtained by the explicit formulation in the stress tensor. It has also been observed that accuracy decreases with increasing of We. One reason for the loss of stress accuracy probably occurs because of the mathematical nature of the inverse algebraic transformation from Ψxx to τxx. The new solver implemented in this work presented convergence and correct solutions for the two geometries, so it was implemented correctly. It also potentiates the viscoelastiFluidFoam starting solver by extending simulations to a larger range of Weissenberg number.