A semi-exact primal-dual method for steady viscoplastic flows
The SMAI Journal of computational mathematics, Volume 12 (2026), pp. 135-170

Constitutive laws for viscoplastic materials involve a multivalued nonlinearity. Duality methods for such equations are known to converge, but the convergence is slow, requiring a large number of iterations. We introduce here a new iterative method of implicit primal-dual type for a class of variational problems, with a particular asymmetrical form in terms of the primal and dual unknowns, with an exact resolution of one relation and an approximate resolution of the other (semi-exact method). Two fast algorithms are proposed. The first includes an adaptive choice of a variable parameter, and the second includes a Newton-like linearized correction procedure. Applying the semi-exact method to a steady viscoplastic problem leads at each iteration to the resolution of a Laplace type equation, thus equivalent in terms of cost per iteration to the well-known Augmented Lagrangian or dual FISTA methods. Numerical tests show that the semi-exact method achieves a faster convergence in terms of number of iterations, compared to the augmented Lagrangian method or to the FISTA* method with acceleration, for Bingham or Herschel–Bulkley laws. This is particularly true when a large zero order term is present, as is the case when solving a time dependent problem.

Published online:
DOI: 10.5802/jcm.145
Classification: 35J20, 76A05, 65K15, 74S05
Keywords: Viscoplastic materials, semi-exact method, implicit primal-dual algorithm, adaptive parameter, Newton linearized correction, steady flows

François Bouchut  1 ; David Maltese  1

1 Laboratoire d’Analyse et de Mathématiques Appliquées (UMR 8050), CNRS, Univ. Gustave Eiffel, UPEC, F-77454, Marne-la-Vallée, France
@article{SMAI-JCM_2026__12__135_0,
     author = {Fran\c{c}ois Bouchut and David Maltese},
     title = {A semi-exact primal-dual method for steady viscoplastic flows},
     journal = {The SMAI Journal of computational mathematics},
     pages = {135--170},
     year = {2026},
     publisher = {Soci\'et\'e de Math\'ematiques Appliqu\'ees et Industrielles},
     volume = {12},
     doi = {10.5802/jcm.145},
     language = {en},
     url = {https://smai-jcm.centre-mersenne.org/articles/10.5802/jcm.145/}
}
TY  - JOUR
AU  - François Bouchut
AU  - David Maltese
TI  - A semi-exact primal-dual method for steady viscoplastic flows
JO  - The SMAI Journal of computational mathematics
PY  - 2026
SP  - 135
EP  - 170
VL  - 12
PB  - Société de Mathématiques Appliquées et Industrielles
UR  - https://smai-jcm.centre-mersenne.org/articles/10.5802/jcm.145/
DO  - 10.5802/jcm.145
LA  - en
ID  - SMAI-JCM_2026__12__135_0
ER  - 
%0 Journal Article
%A François Bouchut
%A David Maltese
%T A semi-exact primal-dual method for steady viscoplastic flows
%J The SMAI Journal of computational mathematics
%D 2026
%P 135-170
%V 12
%I Société de Mathématiques Appliquées et Industrielles
%U https://smai-jcm.centre-mersenne.org/articles/10.5802/jcm.145/
%R 10.5802/jcm.145
%G en
%F SMAI-JCM_2026__12__135_0
François Bouchut; David Maltese. A semi-exact primal-dual method for steady viscoplastic flows. The SMAI Journal of computational mathematics, Volume 12 (2026), pp. 135-170. doi: 10.5802/jcm.145

[1] Vincent Acary; Paul Armand; Hoang Minh Nguyen; Maksym Shpakovych Second-order cone programming for frictional contact mechanics using interior point algorithm, Optim. Methods Softw., Volume 39 (2024) no. 3, pp. 634-663 | Zbl | DOI | MR

[2] Kenneth J. Arrow; Leonid Hurwicz; Hirofumi Uzawa Studies in linear and non-linear programming, Stanford Mathematical Studies in the Social Sciences, 2, Stanford University Press, 1958, 229 pages (with contributions by H. B. Chenery, S. M. Johnson, S. Karlin, T. Marschak, R. M. Solow) | Zbl

[3] Amir Beck; Marc Teboulle A fast iterative shrinkage-thresholding algorithm for linear inverse problems, SIAM J. Imaging Sci., Volume 2 (2009) no. 1, pp. 183-202 | Zbl | DOI | MR

[4] Amir Beck; Marc Teboulle A fast dual proximal gradient algorithm for convex minimization and applications, Oper. Res. Lett., Volume 42 (2014) no. 1, pp. 1-6 | DOI | MR | Zbl

[5] Clément Berger Fluides à seuil : interactions modèles et données, Ph. D. Thesis, Ecole Normale Supérieure de Lyon, France (2024)

[6] Alfredo Bermudez; Carlos Moreno Duality methods for solving variational inequalities, Comput. Math. Appl., Volume 7 (1981), pp. 43-58 | Zbl | DOI | MR

[7] François Bouchut; Carsten Carstensen; Alexandre Ern H 1 regularity of the minimizers for the inviscid total variation and Bingham fluid problems for H 1 data, Nonlinear Anal., Theory Methods Appl., Volume 258 (2025), 113809, 21 pages | Zbl | MR

[8] François Bouchut; Robert Eymard; Alain Prignet Convergence of conforming approximations for inviscid incompressible Bingham fluid flows and related problems, J. Evol. Equ., Volume 14 (2014) no. 3, pp. 635-669 | Zbl | DOI | MR

[9] Miroslav Bulíček; Erika Maringová; Josef Málek On nonlinear problems of parabolic type with implicit constitutive equations involving flux, Math. Models Methods Appl. Sci., Volume 31 (2021) no. 10, pp. 2039-2090 | DOI | Zbl | MR

[10] Rénald Chalayer; Laurent Chupin; Thierry Dubois A bi-projection method for incompressible Bingham flows with variable density, viscosity, and yield stress, SIAM J. Numer. Anal., Volume 56 (2018) no. 4, pp. 2461-2483 | Zbl | DOI | MR

[11] Antonin Chambolle; Thomas Pock A first-order primal-dual algorithm for convex problems with applications to imaging, J. Math. Imaging Vis., Volume 40 (2011) no. 1, pp. 120-145 | Zbl | DOI | MR

[12] Yunmei Chen; Guanghui Lan; Yuyuan Ouyang Optimal primal-dual methods for a class of saddle point problems, SIAM J. Optim., Volume 24 (2014) no. 4, pp. 1779-1814 | Zbl | DOI | MR

[13] Wei Deng; Wotao Yin On the global and linear convergence of the generalized alternating direction method of multipliers, J. Sci. Comput., Volume 66 (2016) no. 3, pp. 889-916 | Zbl | DOI | MR

[14] Ivar Ekeland; Roger Témam Convex analysis and variational problems, Classics in Applied Mathematics, 28, Society for Industrial and Applied Mathematics, 1999, xiv+402 pages (unabridged, corrected republication of the 1976 English original) | Zbl | DOI

[15] Maurice Fabien; Johnny Guzmán; Michael Neilan; Ahmed Zytoon Low-order divergence-free approximations for the Stokes problem on Worsey–Farin and Powell–Sabin splits, Comput. Methods Appl. Mech. Eng., Volume 390 (2022), 114444, 21 pages | Zbl | MR

[16] Enrique D. Fernández-Nieto; José M. Gallardo; Paul Vigneaux Efficient numerical schemes for viscoplastic avalanches. Part 2: The 2D case, J. Comput. Phys., Volume 353 (2018), pp. 460-490 | Zbl | DOI

[17] François Févotte; Ari Rappaport; Martin Vohralík Adaptive regularization, discretization, and linearization for nonsmooth problems based on primal-dual gap estimators, Comput. Methods Appl. Mech. Eng., Volume 418B (2024), 116558, 33 pages | Zbl | MR

[18] Martin Fuchs; Gregory Seregin Variational methods for problems from plasticity theory and for generalized Newtonian fluids, Lecture Notes in Mathematics, 1749, Springer, 2000, vi+269 pages | Zbl | MR

[19] Roland Glowinski Numerical methods for nonlinear variational problems, Scientific Computation, Springer, 2008, xviii+493 pages | Zbl | MR

[20] Roland Glowinski On alternating direction methods of multipliers: a historical perspective, Modeling, Simulation and Optimization for Science and Technology (Computational Methods in Applied Sciences), Volume 34, Springer, 2014, pp. 59-82 | DOI | Zbl

[21] Roland Glowinski; Jacques-Louis Lions; Raymond Trémolières Numerical analysis of variational inequalities, Studies in Mathematics and its Applications, 8, North-Holland, 1981, xxix+776 pages | Zbl | MR

[22] Friedrich Hecht New development in FreeFem++, J. Numer. Math., Volume 20 (2012) no. 3-4, pp. 251-265 | Zbl | MR

[23] Duc Nguyen Hoai Analysis and approximation of compressible viscoplastic models with general nonlinearity for granular flows, Ph. D. Thesis, Univ. Gustave Eiffel (2020) (https://theses.hal.science/tel-03078670)

[24] Nicolas Roquet; Pierre Saramito An adaptive finite element method for Bingham fluid flows around a cylinder, Comput. Methods Appl. Mech. Eng., Volume 192 (2003) no. 31-32, pp. 3317-3341 | Zbl | DOI | MR

[25] Pierre Saramito A damped Newton algorithm for computing viscoplastic fluid flows, J. Non-Newton. Fluid Mech., Volume 238 (2016), pp. 6-15 | DOI | MR

[26] Pierre Saramito; Anthony Wachs Progress in numerical simulation of yield stress fluid flows, Rheol. Acta, Volume 56 (2017), pp. 211-230 | DOI

[27] Timm Treskatis; Miguel A. Moyers-González; Chris J. Price An accelerated dual proximal gradient method for applications in viscoplasticity, J. Non-Newton. Fluid Mech., Volume 238 (2016), pp. 115-130 | DOI | MR

[28] Timm Treskatis; Ali Roustaei; Ian Frigaard; Anthony Wachs Practical guidelines for fast, efficient and robust simulations of yield-stress flows without regularisation: A study of accelerated proximal gradient and augmented Lagrangian methods, J. Non-Newton. Fluid Mech., Volume 262 (2018), pp. 149-164 | DOI | MR

Cited by Sources: