Numerical simulation of non-isothermal viscoelastic fluid flows using a VMS stabilized finite element formulation

The effect of temperature in viscoelastic fluid flows is studied applying a stabilized finite element formulation based on both a standard and a log-conformation reformulation (LCR), and the Variational Multiscale (VMS) method as stabilization technique. The log-conformation reformulation turns out...

Descripción completa

Detalles Bibliográficos
Autores: Moreno Martínez, Laura, Codina, Ramon|||0000-0002-7412-778X, Baiges Aznar, Joan|||0000-0002-3940-5887
Tipo de recurso: artículo
Fecha de publicación:2021
País:España
Institución:Universitat Politècnica de Catalunya (UPC)
Repositorio:UPCommons. Portal del coneixement obert de la UPC
Idioma:inglés
OAI Identifier:oai:upcommons.upc.edu:2117/355902
Acceso en línea:https://hdl.handle.net/2117/355902
https://dx.doi.org/10.1016/j.jnnfm.2021.104640
Access Level:acceso abierto
Palabra clave:Fluid dynamics--Mathematical models
Non-isothermal fluid flow
Viscoelasticity
Log-conformation reformulation
Thermal coupling
Variational sub-grid scales
Dinàmica de fluids -- Càlcul numèric
Àrees temàtiques de la UPC::Física::Física de fluids::Flux de fluids
Descripción
Sumario:The effect of temperature in viscoelastic fluid flows is studied applying a stabilized finite element formulation based on both a standard and a log-conformation reformulation (LCR), and the Variational Multiscale (VMS) method as stabilization technique. The log-conformation reformulation turns out to be crucial to solve the cases with a high Weissenberg number. Regarding temperature coupling, a two-way coupling strategy is employed: on the one hand, the dependence of viscoelastic fluid parameters on temperature is established, together with the addition of a new term to the energy equation which takes into account the stress work. The formulations and the iterative algorithms are validated in the well-known flow past a cylinder benchmark. Besides, the extension 1:3 case is studied, in which several scenarios are explored varying the values of the main dimensionless numbers that characterize the problem to see how the flow pattern and temperature distribution change along the channel.