Abstract
We consider the exact penalization of the incompressibility condition \(\text {div}(\mathbf {u})=0\) for the velocity field of a bi-viscous fluid in terms of the \(L^1\)–norm. This penalization procedure results in a nonsmooth optimization problem for which we propose an algorithm using generalized second-order information. Our method solves the resulting nonsmooth problem by considering the steepest descent direction and extra generalized second-order information associated to the nonsmooth term. This method has the advantage that the divergence-free property is enforced by the descent direction proposed by the method without the need of build-in divergence-free approximation schemes. The inexact penalization approach, given by the \(L^2\)-norm, is also considered in our discussion and comparison.
Similar content being viewed by others
References
Aposporidis, A., et al.: A mixed formulation of the Bingham fluid flow problem: analysis and numerical solution. Comput. Methods Appl. Mech. Eng. 200, 2434–2446 (2011)
Arrow, K.J., Azawa, H., Hurwicz, L., Uzawa, H.: Studies in Linear and Non-linear Programming. Stanford University Press, Stanford (1953)
Babušca, I.: Error-bounds for finite element method. Numerische Mathematik 16, 322–333 (1971)
Bercovier, M., Engelman, M.: A finite element method for incompressible non-Newtonian flows. J. Comput. Phys. 36, 313–326 (1980)
Burke, J.K.: Calmness and exact penalization. SIAM J. Control Optim. 29(2), 493–497 (1991)
Burke, J.K.: An exact penalization viewpoint of constrained optimization. SIAM J. Control Optim. 29(4), 968–998 (1991)
Ciarlet, P.G.: Linear and Nonlinear Functional Analysis with Applications. SIAM, New York (2013)
Clarke, F.: Functional Analysis, Calculus of Variations and Optimal Control. Springer, London (2013)
Bauschke, H.H., Combettes, P.L.: Convex Analysis and Monotone Operator Theory in Hilbert Spaces. Springer, New York (2011)
De los Reyes, J.C., González-Andrade, S.: Path following methods for steady laminar Bingham flow in cylindrical pipes. ESAIM Math. Model. Numer. Anal. 43, 81–117 (2009)
De los Reyes, J.C., González-Andrade, S.: Numerical simulation of two-dimensional Bingham fluid flow by semismooth Newton methods. J. Comput. Appl. Math. 235, 11–32 (2010)
De los Reyes, J.C., González-Andrade, S.: A combined BDF-semismooth Newton approach for time-dependent Bingham flow. Numer. Methods Partial Differ. Equ. 28, 834–860 (2012)
Dennis, J.E., Schnabel, R.B.: Numerical Methods for Unconstrained Optimization and Nonlinear Equations. Society for Industrial and Applied Mathematics, Philadelphia (1996)
Ekeland, I., Témam, R.: Convex Analysis and Variational Problems. SIAM, Philadelphia (1999)
Fusi, L., Farina, A., Rosso, F.: Retrieving the Bingham model from a bi-viscous model: some explanatory remarks. Appl. Math. Lett. 27, 11–14 (2014)
Gabay, D., Mercier, B.: A dual algorithm for the solution of nonlinear variational problems via finite element approximation. Comput. Math. Appl. 2, 17–40 (1976)
Girault, V., Raviart, P.A.: Finite Element Methods for Navier-Stokes Equations: Theory and Algorithms. Springer, Berlin (2012)
Giaquinta, M., Modica, G.: Mathematical Analysis: An Introduction to Functions of Several Variables. Birkhäuser, Boston (2010)
Glowinski, R.: On Alternating Direction Methods of Multipliers: a Historical Perspective. In Modeling Simulation and Optimization for Science and Technology, pp. 59–82. Springer, Berlin (2014)
González-Andrade, S.: Semismooth Newton and path-following methods for the numerical simulation of Bingham fluids, PhD thesis, EPN Quito (2008)
González-Andrade, S.: A preconditioned descent algorithm for variational inequalities of the second kind involving the p-Laplacian operator. Comput. Optim. Appl. 66, 123–162 (2017)
De los Reyes, J.C., González-Andrade, S.: Path following methods for steady laminar Bingham flow in cylindrical pipes. ESAIM Math. Modell. Numer. Anal. 43, 81–117 (2009)
González-Andrade, S., López-Ordóñez, S.: A multigrid optimization algorithm for the numerical solution of quasilinear variational inequalities involving the p-Laplacian. Comput. Math. Appl. 75, 1107–1127 (2018)
Huilgol, R.R., Nguyen, Q.D.: Variational principles and variational inequalities for the unsteady flows of a yield stress fluid. Int. J. Non-Linear Mech. 36, 49–67 (2001)
Huilgol, R.R., You, Z.: Application of the augmented Lagrangian method to steady pipe flows of Bingham. Casson and Herschel-Bulkley fluids. J. Non-Newtonian Fluid Mech. 128, 126–143 (2005)
Kelley, C.T.: Iterative Methods for Linear and Nonlinear Equations. SIAM, Philadelphia (1995)
Kikuchi, N., Oden, J.T.: Contact problems in elasticity: a study of variational inequalities and finite element methods. Studies in Applied Mathematics, SIAM, U.S.A. (1988)
Hintermüller, M., Ito, K., Kunisch, K.: The primal-dual active set strategy as a semismooth Newton method. SIAM J. Optim. 13, 865–888 (2002)
Laaber, P.: Numerical simulation of a three-dimensional Bingham fluid flow (2008)
Lions, J.L.: Optimal Control of Systems Governed by Partial Differential Equations. Springer, Germany (1971)
De los Reyes, J.C., Merino, P.: The second order method with enriched Hessian information for imaging composite sparse optimization problems. arXiv:2009.01878v3 (2021)
Peypouquet, J.: Convex Optimization in Normed Spaces: Theory, Methods and Examples. Springer, London (2015)
Stadler, G.: Elliptic optimal control problems with L 1-control cost and applications for the placement of control devices. Comput. Optim. Appl. 44, 159 (2009)
Tanner, R., Milthorpe, J.: Numerical simulation of the flow of fluids with yield stress. Num. Meth. Lam. Turb. Fl. 680–690 (1983)
O’Donovan, E.J., Tanner, R.I.: Numerical study of the Bingham squeeze film problem. J. Non-Newtonian Fluid Mech. 15, 75–83 (1984)
Temam, R.: Navier-Stokes Equations. Theory and Numerical Analysis. AMS Chelsea Publishing, New York (2001)
Treskatis, T., Moyers-González, M., Price, C.J.: An accelerated dual proximal gradient method for applications in viscoplasticity. J. Non-Newtonian Fluid Mech. 238, 115–130 (2016)
Treskatis, T.: Fast proximal algorithms for applications in viscoplasticity., PhD thesis, University of Canterbury (2016)
Ulbrich, M.: Semismooth Newton Methods for Variational Inequalities and Constrained Optimization Problems in Function Spaces. SIAM, Philadelphia (2011)
Wilbrandt, U.: Stokes–Darcy Equations Analytic and Numerical Analysis. Birkhäuser, Cham (2019)
Zowe, J., Kurcyusz, S.: Regularity and stability for the mathematical programming problem in banach spaces. Appl. Math. Optim. 5, 49–62 (1979)
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This research has been partially supported by Escuela Politécnica Nacional within the project PIGR-18-03 and Secretaría de Educación Superior, Ciencia, Tecnología e Innovación - Senescyt.
Appendix
Appendix
Lemma 8
Let \(\gamma \) and \(\sigma \) be two positive constants. The function \(\phi :\mathbb {R}\rightarrow \mathbb {R}\) defined by \(\phi (a):=\gamma \sigma \frac{ a}{\max (\sigma , \gamma |a|)}\) is Lipschitz continuous and semismooth.
Proof
Let us start by rewriting \(\phi (a)\) as \(\phi (a)=\gamma \sigma \frac{a}{\phi _m(a)}\), with \(\phi _m(a):=\max (\sigma ,\gamma |a|)\). Next, we notice that the max function is globally Lipschitz continuous with constant \(L_{max}\). This fact implies that
We conclude that \(\phi _m\) is Lipschitz continuous. Considering this result, we have that
Now, it is clear that \(0<\sigma \le \phi _m(a_2)\), thus \(\frac{1}{\phi _m(a_2)}\le \frac{1}{\sigma }\). By plugging this inequality in the above expression, we have that
Finally, since \(\left| \frac{a_1}{\phi _m(a_1)}\right| \le \frac{1}{\gamma }\), we conclude, thanks to (78), that
Regarding the semismoothness of \(\phi \), note that the absolute value \(|\cdot |: \mathbb {R} \rightarrow \mathbb {R}\) and the function \(\max (0, \cdot ) : \mathbb {R} \rightarrow \mathbb {R}\) are both semismooth (see [39, Sect. 2.5] and [28, Lemma 3.1] respectively). Then, since the composition of semismooth functions in \(\mathbb {R}^n\) is a semismooth function [39, Prop. 2.9], it follows that \(\phi (a)\) is semismooth.
Remark 4
The function \(\varphi _j:\mathbb {R}^m \rightarrow \mathbb {R}\) defined by \(\varphi _j(a):=\gamma \sigma \frac{ a_j}{\max (g, \beta |a|)}\) is also Lipschitz continuous and semismooth. The proof of this assertion is analogous to the one given in Lemma 8.
Lemma 9
Let \(\phi ({{\,\mathrm{div}\,}}\mathbf {u}(x))=\displaystyle \sigma \gamma \frac{ {{\,\mathrm{div}\,}}\mathbf {u}(x)}{\max (\sigma , \gamma |{{\,\mathrm{div}\,}}\mathbf {u}(x)|)}\) with \(\gamma \) and \(\sigma \) positive constants. A measurable selection \(M_{\phi }( \mathbf {u})\) of Clarke’s generalized Jacobian \( \partial \phi ({{\,\mathrm{div}\,}}\mathbf {u}) \) is :
a.e on \(\Omega \)
Proof
Let \(\phi _3=\phi _1 \circ \phi _2\), where \(\phi _1(z)= \max (0,z) + \sigma \) and \(\phi _2(y)= \gamma |y| - \sigma \). Then the following identity holds:
From [28, pp. 869] we have that \( M_{\phi _1}( \gamma |y| - \sigma ) \in \partial \phi _1(\gamma |y| - \sigma )\) given by
is a measurable selection of \(\partial \phi _1(\gamma |y| - \sigma )\). Next, since \(\phi _2\) involves the function \(| \cdot |\) evaluated at \(y\ne 0\). From [39, Exaple 2.5.1] we have that
Moreover, the chain rule for Clarke’s generalized Jacobian [39, Prop. 2.3] yields that:
Thus, since \(y \ne 0\),
Clearly, \(\phi (y)= \sigma \gamma \frac{y}{ \phi _3(y)}\). Then, from the composition of functions we obtain that
Then, from (80) the following cases can occur:
-
\(\gamma | y| > \sigma \). Then:
$$\begin{aligned} M_{\phi }(y) = \displaystyle \sigma \frac{1}{|y|} \, - \sigma \displaystyle \frac{y^2}{|y|^3} = 0. \end{aligned}$$ -
\(\gamma | y| \le \sigma \) gives:
$$\begin{aligned} M_{\phi }(y) = \gamma . \end{aligned}$$
Finally, by taking \(y={{\,\mathrm{div}\,}}\mathbf {u}(x)\) we have the desired result.
Remark 5
The measurable selection \(N_j( \mathbf { u}(x))\) of Clarke’s generalized Jacobian \( \partial \varphi _j(\mathcal {E} \mathbf {u}(x))\) is obtained by an analogous procedure to Lemma 9.
Rights and permissions
About this article
Cite this article
González-Andrade, S., López-Ordóñez, S. & Merino, P. Nonsmooth exact penalization second-order methods for incompressible bi-viscous fluids. Comput Optim Appl 80, 979–1025 (2021). https://doi.org/10.1007/s10589-021-00314-2
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10589-021-00314-2