Viscoelastic fluid flow simulation using the contravariant deformation formulation
| Autor: | MA Martien Hulsen, Patrick D. Anderson, Markus Hütter, Mick A. Carrozza |
|---|---|
| Přispěvatelé: | Processing and Performance |
| Jazyk: | angličtina |
| Rok vydání: | 2019 |
| Předmět: |
Finite element method
General Chemical Engineering Multiplicative decomposition 01 natural sciences 010305 fluids & plasmas Physics::Fluid Dynamics Viscosity 0103 physical sciences Weissenberg number General Materials Science Tensor Couette flow Physics Finite volume method 010304 chemical physics Applied Mathematics Mechanical Engineering Mathematical analysis Numerical stability Condensed Matter Physics Flow (mathematics) Planar Couette flow Flow around a confined cylinder Shear flow Conformation-tensor based models |
| Zdroj: | Journal of Non-Newtonian Fluid Mechanics, 270, 23-35. Elsevier |
| ISSN: | 0377-0257 |
| Popis: | The new formulation for conformation-tensor based viscoelastic fluid models, written in terms of the contravariant deformation tensor as introduced by Hutter et al. (2018), has been approached numerically. It has been implemented in a finite element method framework with the DEVSS-G/SUPG method for stabilisation. A stress-implicit as well as a stress-explicit formulation is presented that allows to integrate the non-linear flow equations, with or without a solvent contribution, forward in time with second-order accuracy. The time-dependent stability in shear flow is maintained using the contravariant deformation, which is tested by solving a perturbed planar Couette flow of an upper-convected Maxwell (UCM) fluid. The stability and accuracy of our implementation has furthermore been tested by solving the flow around a cylinder with confining walls. The new formulation turns out to be more stable as compared to the formulation in terms of the conformation tensor. It enables to simulate the flow of a Giesekus fluid with non-linear parameter α = 0.01 and viscosity ratio ν = η s / η 0 = 0.59 way beyond the Weissenberg number at which the High Weissenberg Number Problem manifests itself using the standard formulation. The stability appears to be similar to the log-conformation representation. The computed stress profiles up to a Weissenberg number of Wi = 0.6 for an Oldroyd-B fluid with viscosity ratio ν = 0.59 , compare well with benchmark results from other studies, including other finite element, finite volume and spectral element methods. |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|
Full Text from ScienceDirect