Search

We consider the mathematical analysis and numerical approximation of a system of nonlinear partial differential equations that arises in models that have relevance to steady isochoric flows of colloidal suspensions. The symmetric velocity gradient is assumed to be a monotone nonlinear function of the deviatoric part of the Cauchy stress tensor. We prove the existence of a unique weak solution to the problem, and under the additional assumption that the nonlinearity involved in the constitutive relation is Lipschitz continuous we also prove uniqueness of the weak solution. We then construct mixed finite element approximations of the system using both conforming and nonconforming finite element spaces. For both of these we prove the convergence of the method to the unique weak solution of the problem, and in the case of the conforming method we provide a bound on the error between the analytical solution and its finite element approximation in terms of the best approximation error from the finite element spaces. We propose first a Lions-Mercier type iterative method and next a classical fixed-point algorithm to solve the finite-dimensional problems resulting from the finite element discretisation of the system of nonlinear partial differential equations under consideration and present numerical experiments that illustrate the practical performance of the proposed numerical method., Comment: [v2] minor modifications; to appear in ESAIM: M2AN

A low-order finite element method is constructed and analysed for an incompressible non-Newtonian flow problem with power-law rheology. The method is based on a continuous piecewise linear approximation of the velocity field and piecewise constant approximation of the pressure. Stabilisation, in the form of pressure jumps, is added to the formulation to compensate for the failure of the inf-sup condition, and using an appropriate lifting of the pressure jumps a divergence-free approximation to the velocity field is built and included in the discretisation of the convection term. This construction allows us to prove the convergence of the resulting finite element method for the entire range $r>\frac{2 d}{d+2}$ of the power-law index $r$ for which weak solutions to the model are known to exist in $d$ space dimensions, $d \in \{2,3\}$., 23 pages, 1 figure