Numerical scheme for simulation of transient flows of non-Newtonian fluids characterised by a non-monotone relation between the symmetric part of the velocity gradient and the Cauchy stress tensor

Abstract

We propose a numerical scheme for simulation of transient flows of incompressible non-Newtonian fluids characterised by a non-monotone relation between the symmetric part of the velocity gradient (shear rate) and the Cauchy stress tensor (shear stress). The main difficulty in dealing with the governing equations for flows of such fluids is that the nonmonotone constitutive relation allows several values of the stress to be associated with the same value of the symmetric part of the velocity gradient. This issue is handled via a reformulation of the governing equations. The equations are reformulated as a system for the triple pressure-velocity-apparent viscosity, where the apparent viscosity is given by a scalar implicit equation. We prove that the proposed numerical scheme has-on the discrete level-a solution, and using the proposed scheme we numerically solve several flow problems.1 The norm of a tensorial quantity A is defined as the standard Frobenius norm, A = def (Tr (AA ⊺ ))) 1 2 .

Thermodynamical Extension of a Symplectic Numerical Scheme with Half Space and Time Shifts Demonstrated on Rheological Waves in Solids

“…which we have already met in (9), as the thermodynamical requirement (6) at the continuum level, and which also induces, via (57),…”

“…As expected, these conditions prove to be equivalent to the ones obtained via Jury's criteria -we omit the details to avoid redundant repetition. 9 Following Jury [25], a matrix is positive innerwise if the determinant of the matrix and its all inners are positive. Here, inner ∆ m−2 of an m × m matrix is formed by deleting its first and mth rows and columns, inner ∆ m−4 is the inner of ∆ m−2 , and the procedure is continued until ∆ 1 or ∆ 2 is reached.…”

Dynamic and Renormalization-Group Extensions of the Landau Theory of Critical Phenomena

“…This potential was found in [ 24 ] and is based on a work by Le Roux and Rajagopal [ 25 ]. The rate thermodynamic analysis of the qualitative implied by the potential was then confirmed by numerical simulations in [ 26 ].…”

Proposal of an open-source computational toolbox for solving PDEs in the context of chemical reaction engineering using FEniCS and complementary components

“…The FEniCS Project has been used within a wide range of research fields, such as tidal energy, geoscience, fluid mechanics, theoretical biology, strain gradient elasticity, biophysics and metamodelling ( Abali, 2019 ; Epanchintsev et al., 2016 ; Funke et al., 2019 ; Goodwin et al., 2019 ; Haagenson et al., 2020 ; Janečka et al., 2019 ; Lejeune, 2020 ; Murray and Young, 2020 ; Phunpeng and Baiz, 2015 ; Zhu and Yan, 2019 ). Abali (2017) , for instance, demonstrated several modelling examples with different engineering applications through a continuum mechanics approach.…”