Miguel Moyers-Gonzalez

An accelerated dual proximal gradient method for applications in viscoplasticity

Abstract We present a very simple and fast algorithm for the numerical solution of viscoplastic flow problems without prior regularisation. Compared to the widespread alternating direction method of multipliers (ADMM / ALG2), the new method features three key advantages: firstly, it accelerates the worst-case convergence rate from O ( 1 / k ) to O(1/k), where k is the iteration counter. Secondly, even for nonlinear constitutive models like those of Casson or Herschel–Bulkley, no nonlinear systems of equations have to be solved in the subproblems of the algorithm. Thirdly, there is no need to augment the Lagrangian, which eliminates the difficulty of choosing a penalty parameter heuristically. In this paper, we transform the usual velocity-based formulation of viscoplastic flow problems to a dual formulation in terms of the stress. For the numerical solution of this dual problem we apply FISTA, an accelerated first-order optimisation algorithm from the class of so-called proximal gradient methods. Finally, we conduct a series of numerical experiments, focussing on stationary flow in two-dimensional square cavities. Our results confirm that Algorithm FISTA*, the new dual-based FISTA, outperforms state-of-the-art algorithms such as ADMM / ALG2 by several orders of magnitude. We demonstrate how this speedup can be exploited to identify the free boundary between yielded and unyielded regions with previously unknown accuracy. Since the accelerated algorithm relies solely on Stokes-type subproblems and nonlinear function evaluations, existing code based on augmented Lagrangians would require only few minor adaptations to obtain an implementation of FISTA*.

Linear stability analysis for plane-Poiseuille flow of an elastoviscoplastic fluid with internal microstructure for large Reynolds numbers

Abstract We study the linear stability of Plane Poiseuille flow of an elastoviscoplastic fluid using a revised version of the model proposed by Putz and Burghelea (A.M.V. Putz, T.I. Burghelea, Rheol. Acta 48 (2009) 673–689). The evolution of the microstructure upon a gradual increase of the external forcing is governed by a structural variable (the concentration of solid material elements) which decays smoothly from unity to zero as the stresses are gradually increased beyond the yield point. Stability results are in close conformity with the ones of a pseudo-plastic fluid. Destabilizing effects are related to the presence of an intermediate transition zone where elastic solid elements coexist with fluid elements. This region brings an elastic contribution which does modify the stability of the flow.