The ability to understand and predict the flow behaviour of non-Newtonian fluids in swirling flow is industrially important for the efficient design and performance of processes which utilise fluids with complex rheological properties. In particular, fluids with elastic properties are not well described by non-Newtonian constitutive models, such that predictions using such models must be carefully validated. A benchmark problem is proposed here which provides a well defined geometry to study the swirling flow of non-Newtonian fluids as a test case for the validation of constitutive models. The confined swirling flow utilised is a torsionally driven cavity where the test fluid is confined in a cylinder with a rotating bottom lid, and stationary side walls and top lid. The flow field is three-dimensional and consists of both a primary motion, which is directed azimuthally, and a secondary motion, which is located in the radial and axial plane of the cylinder and driven by inertial and/or elastic forces.

The ability to understand and predict the flow behaviour of non-Newtonian fluids in swirling flow is industrially important for the efficient design and performance of processes which utilise fluids with complex rheological properties. In particular, fluids with elastic properties are not well described by non-Newtonian constitutive models, such that predictions using such models must be carefully validated. A benchmark problem is proposed here which provides a well defined geometry to study the swirling flow of non-Newtonian fluids as a test case for the validation of constitutive models. The confined swirling flow utilised is a torsionally driven cavity where the test fluid is confined in a cylinder with a rotating bottom lid, and stationary side walls and top lid. The flow field is three-dimensional and consists of both a primary motion, which is directed azimuthally, and a secondary motion, which is located in the radial and axial plane of the cylinder and driven by inertial and/or elastic forces.

In the recent years, researchers have discovered a multitude of applications using nanofibers in fields like composites, biotechnology, environmental engineering, defense, optics and electronics. This increase in nanofiber applications needs a higher rate of nanofiber production. Electrospinning has proven to be the best nanofiber manufacturing process because of simplicity and material compatibility. Study of effects of various electrospinning parameters is important to improve the rate of nanofiber processing. In addition, several applications demand well-oriented nanofibers. Researchers have experimentally tried to control the nanofibers using secondary external electric field. In the first study, the electrospinning process is modeled and the bending instability of a viscoelastic jet is simulated. For this, the existing discrete bead model is modified and the results are compared, qualitatively, with previous works in literature. In this study, an attempt is also made to simulate the effect of secondary electric field on electrospinning process and whipping instability. It is observed that the external secondary field unwinds the jet spirals, reduces the whipping instability and increases the tension in the fiber. Lattice Boltzmann method (LBM) has gained popularity in the past decade as the method is easy implement and can also be parallelized. In the second part of this thesis, a hybrid numerical scheme which couples lattice Boltzmann method with finite difference method for a Oldroyd-B viscoelastic solution is proposed. In this scheme, the polymer viscoelastic stress tensor is included in the equilibrium distribution function and the distribution function is updated using SRT-LBE model. Then, the local velocities from the distribution function are evaluated. These local velocities are used to evaluate local velocity gradients using a central difference method in space. Next, a forward difference scheme in time is used on the Maxwell Upper Convected model and the viscoelastic stress tensor is updated. Finally, using the proposed numerical method start-up Couette flow problem for Re = 0.5 and We = 1.1, is simulated. The velocity and stress results from these simulations agree very well with the analytical solutions.

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.

Viscoelastic fluid flow with immersed boundaries of complex geometry is widely found both in nature and engineering processes. Examples include haemocytes moving in human blood flow, self-propulsion of microscopic organisms in complex liquids, hydraulic fracturing with sand in oil flow, and suspension flow with viscoelastic medium. Computational modelling of such systems is important for understanding complex biological processes and assisting engineering designs. Conventional simulation methods use conformed meshes to resolve the boundaries of complex geometry. Dynamically updating the conformed mesh is computationally expensive and makes parallelization difficult. In comparison, Cartesian grid methods are more promising for large scale parallel simulation. Using Cartesian grid methods to simulate viscoelastic fluid flow with complex boundaries is a relatively unexplored area. In this thesis, a sharp interface Cartesian grid method (SICG) and a smoothed interface immersed boundary method (SIIB) are developed in order to simulate viscoelastic fluids in complex geometries. The SICG method shows a better prediction of the stress on stationary boundaries while the SIIB method shows reduced non-physical oscillations in the computation of drag and lift forces on moving boundaries. Parallel implementations of both solvers are developed. Convergence of the numerical schemes is shown and the implementations are validated with a few benchmark problems with both stationary and moving boundaries. This study also focuses on the simulation of flows past 2D cylindrical or 3D spherical particles. Lateral migration of particles induced by inertial and viscoelastic effects are investigated with different flow types. Equilibrium positions of inertia-induced migration are reported as a function of the particle Reynolds number and the blockage ratio. The migration in the viscoelastic fluid is simulated from zero elastic number to a finite elastic number. The inclusion of both inertial and viscoelastic effects on the lateral migration of a particle is the first of its kind. New findings are reported for the equilibrium positions of a spherical particle in square duct flow, which suggest the need for future experimental and computational work.

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.

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.

Three different models, relating to the study of microswimmers immersed in a low Reynolds number fluid, are presented. The underlying, mathematical concepts employed in each are developed using singularity methods of Stokes flow. The first topic concerns the motility of an artificial, three-sphere microswimmer with prescribed, non-reciprocal, internal forces. The swimmer progresses through a low Reynolds number, nonlinear, viscoelastic medium. The model developed illustrates that the presence of the viscoelastic rheology, when compared to a Newtonian environment, increases both the net displacement and swimming efficiency of the microswimmer. The second area concerns biological microswimming, modelling a sperm cell with a hyperactive waveform (vigorous, asymmetric beating), bound to the epithelial walls of the female, reproductive tract. Using resistive-force theory, the model concludes that, for certain regions in parameter space, hyperactivated sperm cells can induce mechanical forces that pull the cell away from the wall binding. This appears to occur via the regulation of the beat amplitude, wavenumber and beat asymmetry. The next topic presents a novel generalisation of slender-body theory that is capable of calculating the approximate flow field around a long, thin, slender body with circular cross sections that vary arbitrarily in radius along a curvilinear centre-line. New, permissible, slender-body shapes include a tapered flagellum and those with ribbed, wave-like structures. Finally, the detailed analytics of the generalised, slender-body theory are exploited to develop a numerical implementation capable of simulating a wider range of slender-body geometries compared to previous studies in the field.

Polímeros sintéticos ocupam uma posição de grande importância no estilo de vida moderno, servindo como matérias-primas para a construção de uma variedade de utensílios. Apesar do grande número de operações de processamento e produtos disponíveis, o planejamento de produtos e a otimização dos processos de produção raramente constituem-se de tarefas triviais. Isso deve-se ao fato da maioria das operações aplicadas na indústria de processamento de polímeros envolverem geometrias e padrões de escoamento complexos, além da dificuldade intrínseca relacionada ao comportamento reológico complexo de polímeros fundidos ou soluções poliméricas. Devido a estes fatores, o desenvolvimento de técnicas de dinâmica de fluido computacional (computational fluid dynamics – CFD) para a simulação de escoamentos de fluidos poliméricos e etapas de operações de processamento tem sido assunto de numerosos estudos durante as últimas décadas. Sob esta perspectiva, o solver viscoelasticFluidFoam, merece destaque. Ele é capaz de resolver simulações de escoamentos de fluidos viscoelásticos utilizando diferentes equações constitutivas. Contudo, apesar de resultados existentes na literatura apresentarem um bom potencial de aplicação, uma análise extensiva de seu desempenho numérico ainda não foi realizada. Neste contexto, a proposta do presente trabalho é a análise da influência de parâmetros de malha, numéricos e constitutivos no comportamento do solver. As bases para os testes compreendem uma geometria simples – escoamento laminar entre duas placas paralelas – o modelo constitutivo de Oldroyd-B e respectivas soluções analíticas para os campos de velocidade e tensão. Mesmo os testes demonstrando a inegável versatilidade do solver, eles revelam limitações em lidar com algumas configurações de malha e parâmetros constitutivos, principalmente com relação ao refinamento na direção perpendicular ao escoamento, diminuição do número de Reynolds e aumento do número de Weisenberg. Estas limitações podem ser parcialmente contornadas com escolha adequada de parâmetros de relaxação das variáveis e da razão de aspecto dos volumes de controle. Tais dificuldades não estão presentes em simulações de escoamentos de fluidos newtonianos em condições semelhantes, sugerindo que trabalhos futuros devem focar em implementações mais robustas do solver viscoelasticFluidFoam. / Synthetic polymers hold a position of great importance in modern lifestyle, serving as raw materials for the construction of a wide variety of appliances. Despite the large number of processing operations and products available, product planning and optimization of production processes rarely constitute a trivial task. This is due to the fact of operations applied in polymer processing industry involve complex geometries and flow patterns, plus the intrinsic difficulty related to the molten polymers or polymer solutions complex rheological behavior. Because of these factors, the development of techniques of computational fluid dynamics (CFD) for the simulation of flows of polymeric fluids and stages of processing operations has been the subject of numerous studies during the last decades. From this perspective, the viscoelasticFluidFoam solver deserves mention. The solver is capable of resolving simulations of viscoelastic fluid flows using different constitutive equations. However, despite the existing results in the literature present a great potential for application, an extensive analysis of their numerical performance has not been performed yet. The purpose of this paper is to examine the influence of mesh, numerical and constitutive parameters in the behavior of the solver. Bases for the tests comprise a simple geometry – laminar flow between two parallel plates – the constitutive model of Oldroyd-B and its analytical solutions for the velocity and stress fields. Although the tests show the undeniable versatility of the solver, they also reveal limitations in dealing with some mesh settings and constitutive parameters, particularly with respect to refinement in the direction perpendicular to the flow, decreasing in the Reynolds number and increasing in the Weisenberg number. This limitation can be partially circumvented with proper choice of variables relaxation parameters and aspect ratio of the control volumes. Such difficulties are not present in simulations of Newtonian fluids flows under similar conditions, suggesting that future works should focus on more robust implementations of the viscoelasticFluidFoam solver.

Polímeros sintéticos ocupam uma posição de grande importância no estilo de vida moderno, servindo como matérias-primas para a construção de uma variedade de utensílios. Apesar do grande número de operações de processamento e produtos disponíveis, o planejamento de produtos e a otimização dos processos de produção raramente constituem-se de tarefas triviais. Isso deve-se ao fato da maioria das operações aplicadas na indústria de processamento de polímeros envolverem geometrias e padrões de escoamento complexos, além da dificuldade intrínseca relacionada ao comportamento reológico complexo de polímeros fundidos ou soluções poliméricas. Devido a estes fatores, o desenvolvimento de técnicas de dinâmica de fluido computacional (computational fluid dynamics – CFD) para a simulação de escoamentos de fluidos poliméricos e etapas de operações de processamento tem sido assunto de numerosos estudos durante as últimas décadas. Sob esta perspectiva, o solver viscoelasticFluidFoam, merece destaque. Ele é capaz de resolver simulações de escoamentos de fluidos viscoelásticos utilizando diferentes equações constitutivas. Contudo, apesar de resultados existentes na literatura apresentarem um bom potencial de aplicação, uma análise extensiva de seu desempenho numérico ainda não foi realizada. Neste contexto, a proposta do presente trabalho é a análise da influência de parâmetros de malha, numéricos e constitutivos no comportamento do solver. As bases para os testes compreendem uma geometria simples – escoamento laminar entre duas placas paralelas – o modelo constitutivo de Oldroyd-B e respectivas soluções analíticas para os campos de velocidade e tensão. Mesmo os testes demonstrando a inegável versatilidade do solver, eles revelam limitações em lidar com algumas configurações de malha e parâmetros constitutivos, principalmente com relação ao refinamento na direção perpendicular ao escoamento, diminuição do número de Reynolds e aumento do número de Weisenberg. Estas limitações podem ser parcialmente contornadas com escolha adequada de parâmetros de relaxação das variáveis e da razão de aspecto dos volumes de controle. Tais dificuldades não estão presentes em simulações de escoamentos de fluidos newtonianos em condições semelhantes, sugerindo que trabalhos futuros devem focar em implementações mais robustas do solver viscoelasticFluidFoam. / Synthetic polymers hold a position of great importance in modern lifestyle, serving as raw materials for the construction of a wide variety of appliances. Despite the large number of processing operations and products available, product planning and optimization of production processes rarely constitute a trivial task. This is due to the fact of operations applied in polymer processing industry involve complex geometries and flow patterns, plus the intrinsic difficulty related to the molten polymers or polymer solutions complex rheological behavior. Because of these factors, the development of techniques of computational fluid dynamics (CFD) for the simulation of flows of polymeric fluids and stages of processing operations has been the subject of numerous studies during the last decades. From this perspective, the viscoelasticFluidFoam solver deserves mention. The solver is capable of resolving simulations of viscoelastic fluid flows using different constitutive equations. However, despite the existing results in the literature present a great potential for application, an extensive analysis of their numerical performance has not been performed yet. The purpose of this paper is to examine the influence of mesh, numerical and constitutive parameters in the behavior of the solver. Bases for the tests comprise a simple geometry – laminar flow between two parallel plates – the constitutive model of Oldroyd-B and its analytical solutions for the velocity and stress fields. Although the tests show the undeniable versatility of the solver, they also reveal limitations in dealing with some mesh settings and constitutive parameters, particularly with respect to refinement in the direction perpendicular to the flow, decreasing in the Reynolds number and increasing in the Weisenberg number. This limitation can be partially circumvented with proper choice of variables relaxation parameters and aspect ratio of the control volumes. Such difficulties are not present in simulations of Newtonian fluids flows under similar conditions, suggesting that future works should focus on more robust implementations of the viscoelasticFluidFoam solver.