The aim of this contribution is to present a new Newton-type solver for yield stress fluids, for instance for viscoplastic Bingham fluids. In contrast to standard globally defined (‘outer’) damping strategies, we apply weighting strategies for the different parts inside of the resulting Jacobian matrices (after discretizing with FEM), taking into account the special properties of the partial operators which arise due to the differentiation of the corresponding nonlinear viscosity function. Moreover, we shortly discuss the corresponding extension to fluids with a pressure-dependent yield stress which are quite common for modelling granular material. From a numerical point of view, the presented method can be seen as a generalized Newton approach for non-smooth problems.
[1]
David G. Schaeffer,et al.
Instability in the evolution equations describing incompressible granular flow
,
1987
.
[2]
B. Cazacliu,et al.
Analysis of a Regularized Bingham Model with Pressure-Dependent Yield Stress
,
2015
.
[3]
Olivier Pouliquen,et al.
A constitutive law for dense granular flows
,
2006,
Nature.
[4]
Stefan Turek,et al.
An investigation of frictional and collisional powder flows using a unified constitutive equation
,
2010
.
[5]
Stefan Turek,et al.
Efficient Solvers for Incompressible Flow Problems - An Algorithmic and Computational Approach
,
1999,
Lecture Notes in Computational Science and Engineering.