Four-Field Finite Element Solver and Sensitivities for Quasi-Newtonian Flows

Abstract

A computationally efficient finite element algorithm for power-law fluid is elaborated in view of extensive direct and inverse simulations. We adopt a splitting technique to simplify the nonlinear structure of the fluids equations and derive a four-field saddle-point formulation for which we prove the existence of a solution. The resolution of the corresponding variational inequalities is based on an augmented Lagrangian method and a mixed finite element discretization. The resulting iterative solver proves to be fast and robust with low memory consumption. The time-saving provided by the algorithm compared to the standard algorithms of fixed point and Newton increases with the number of degrees of freedom and the nonlinearity of the problem. It is therefore well suited for the solution of large problems with a great number of elements and for corresponding adjoint-based computations. Bidimensional numerical experiments are performed on two realistic situations of gravity flows: an experimental viscoplast...