Path following methods for steady laminar Bingham flow in cylindrical pipes

Juan Carlos De Los Reyes; González, Sergio

ESAIM: Modélisation mathématique et analyse numérique, Tome 43 (2009) no. 1, pp. 81-117.

This paper is devoted to the numerical solution of stationary laminar Bingham fluids by path-following methods. By using duality theory, a system that characterizes the solution of the original problem is derived. Since this system is ill-posed, a family of regularized problems is obtained and the convergence of the regularized solutions to the original one is proved. For the update of the regularization parameter, a path-following method is investigated. Based on the differentiability properties of the path, a model of the value functional and a correspondent algorithm are constructed. For the solution of the systems obtained in each path-following iteration a semismooth Newton method is proposed. Numerical experiments are performed in order to investigate the behavior and efficiency of the method, and a comparison with a penalty-Newton-Uzawa-conjugate gradient method, proposed in [Dean et al., J. Non-newtonian Fluid Mech. 142 (2007) 36-62], is carried out.

[1] J. Alberty, C. Carstensen and S. Funken, Remarks around 50 lines of $\mathrm{M}$atlab: short finite element implementation. Numer. Algorithms 20 (1999) 117-137. | MR | Zbl

[2] H.W. Alt, Lineare Funktionalanalysis. Springer-Verlag (1999). | Zbl

[3] A. Bensoussan and J. Frehse, Regularity Results for Nonlinear Elliptic Systems and Applications, Applied Mathematical Sciences 151. Springer-Verlag (2002). | MR | Zbl

[4] H. Brézis, Monotonicity methods in $\mathrm{H}$ilbert spaces and some applications to nonlinear partial differential equations, in Contributions to Non-linear Functional Analysis, E. Zarantonello Ed., Acad. Press (1971) 101-156. | MR | Zbl

[5] J.C. De Los Reyes and K. Kunisch, A semi-smooth $\mathrm{N}$ewton method for control constrained boundary optimal control of the $\mathrm{N}$avier-$\mathrm{S}$tokes equations. Nonlinear Anal. 62 (2005) 1289-1316. | MR | Zbl

[6] E.J. Dean, R. Glowinski and G. Guidoboni, On the numerical simulation of $\mathrm{B}$ingham visco-plastic flow: Old and new results. J. Non-Newtonian Fluid Mech. 142 (2007) 36-62. | Zbl

[7] G. Duvaut and J.L. Lions, Inequalities in Mechanics and Physics. Springer-Verlag, Berlin (1976). | MR | Zbl

[8] I. Ekeland and R. Temam, Convex Analysis and Variational Problems. North-Holland Publishing Company, The Netherlands (1976). | MR | Zbl

[9] M. Fuchs and G. Seregin, Some remarks on non-Newtonian fluids including nonconvex perturbations of the $\mathrm{B}$ingham and $\mathrm{P}$owell-$\mathrm{E}$yring model for viscoplastic fluids. Math. Models Methods Appl. Sci. 7 (1997) 405-433. | MR | Zbl

[10] M. Fuchs and G. Seregin, Regularity results for the quasi-static $\mathrm{B}$ingham variational inequality in dimensions two and three. Math. Z. 227 (1998) 525-541. | MR | Zbl

[11] M. Fuchs, J.F. Grotowski and J. Reuling, On variational models for quasi-static Bingham fluids. Math. Methods Appl. Sci. 19 (1996) 991-1015. | MR | Zbl

[12] R. Glowinski, Numerical Methods for Nonlinear Variational Problems, Springer Series in Computational Physics. Springer-Verlag (1984). | MR | Zbl

[13] R. Glowinski, J.L. Lions and R. Tremolieres, Analyse numérique des inéquations variationnelles1976). | Zbl

[14] M. Hintermüller and K. Kunisch, Path-following methods for a class of constrained minimization problems in function spaces. SIAM J. Optim. 17 (2006) 159-187. | MR | Zbl

[15] M. Hintermüller and K. Kunisch, Feasible and non-interior path-following in constrained minimization with low multiplier regularity. SIAM J. Contr. Opt. 45 (2006) 1198-1221. | MR | Zbl

[16] M. Hintermüller and G. Stadler, An infeasible primal-dual algorithm for TV-based inf-convolution-type image restoration. SIAM J. Sci. Comput. 28 (2006) 1-23. | MR | Zbl

[17] M. Hintermüller, K. Ito and K. Kunisch, The primal-dual active set strategy as a semi-smooth $\mathrm{N}$ewton method. SIAM J. Optim. 13 (2003) 865-888. | MR | Zbl

[18] R.R. Huilgol and Z. You, Application of the augmented $\mathrm{L}$agrangian method to steady pipe flows of $\mathrm{B}$ingham, $\mathrm{C}$asson and $\mathrm{H}$erschel-$\mathrm{B}$ulkley fluids. J. Non-Newtonian Fluid Mech. 128 (2005) 126-143.

[19] K. Ito and K. Kunisch, Augmented $\mathrm{L}$agrangian methods for nonsmooth, convex optimization in $\mathrm{H}$ilbert spaces. Nonlinear Anal. 41 (2000) 591-616. | MR | Zbl

[20] K. Ito and K. Kunisch, Semi-smooth $\mathrm{N}$ewton methods for variational inequalities of the first kind. ESAIM: M2AN 37 (2003) 41-62. | Numdam | MR | Zbl

[21] J.-L. Lions, Optimal Control of Systems Governed by Partial Differential Equations. Springer-Verlag (1971). | MR | Zbl

[22] P.P. Mosolov and V.P. Miasnikov, Variational methods in the theory of the fluidity of a viscous-plastic medium. J. Appl. Math. Mech. (P.M.M.) 29 (1965) 468-492. | Zbl

[23] T. Papanastasiou, Flows of materials with yield. J. Rheology 31 (1987) 385-404. | Zbl

[24] G. Stadler, Infinite-dimensional Semi-smooth $\mathrm{N}$ewton and Augmented $\mathrm{L}$agrangian Methods for Friction and Contact Problems in Elasticity. Ph.D. thesis, Karl-Franzens University of Graz, Graz, Austria (2004).

[25] G. Stadler, Path-following and augmented Lagrangian methods for contact problems in linear elasticity. J. Comp. Appl. Math. 203 (2007) 533-547. | MR | Zbl

[26] D. Sun and J. Han, Newton and quasi-$\mathrm{N}$ewton methods for a class of nonsmooth equations and related problems. SIAM J. Optim. 7 (1997) 463-480. | MR | Zbl

[27] M. Ulbrich, Nonsmooth $\mathrm{N}$ewton-like methods for variational inequalities and constrained optimization problems in function spaces. Habilitation thesis, Technische Universität München, Germany (2001-2002).