×

An inexact Newton method for solving complementarity problems in hydrodynamic lubrication. (English) Zbl 1384.76045

Summary: We present an iterative procedure based on a damped inexact Newton iteration for solving linear complementarity problems. We introduce the method in the framework of a popular problem arising in mechanical engineering: the analysis of cavitation in lubricated contacts. In this context, we show how the perturbation and the damping parameter are chosen in our method and we prove the global convergence of the entire procedure. A Fortran implementation of the method is finally analyzed. First, we validate the procedure and analyze all its components, performing also a comparison with a recently proposed technique based on the Fischer-Burmeister-Newton iteration. Then, we solve a 2D problem and provide some insights on an efficient implementation of the method exploiting routines of the Lapack and of the PETSc packages for the solution of inner linear systems.

MSC:

76M25 Other numerical methods (fluid mechanics) (MSC2010)
76D08 Lubrication theory
65H10 Numerical computation of solutions to systems of equations
65K05 Numerical mathematical programming methods
90C51 Interior-point methods
76B10 Jets and cavities, cavitation, free-streamline theory, water-entry problems, airfoil and hydrofoil theory, sloshing

Software:

LAPACK; SuperLU; PETSc; MUMPS
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Almqvist, A., Wall, P.: Modelling cavitation in (elasto)hydrodynamic lubrication. In: Darji, P.H. (ed.) Advances in Tribology (2016) · Zbl 0551.76007
[2] Amestoy, P.R., Duff, I.S., Koster, J., L’Excellent, J.Y.: A fully asynchronous multifrontal solver using distributed dynamic scheduling. SIAM J. Matrix. Anal. A 23(1), 15-41 (2001) · Zbl 0992.65018 · doi:10.1137/S0895479899358194
[3] Amestoy, P.R., Guermouche, A., L’Excellent, J.Y., Pralet, S.: Hybrid scheduling for the parallel solution of linear systems. Parallel Comput. 32(2), 136-156 (2006) · doi:10.1016/j.parco.2005.07.004
[4] Anderson, E., Bai, Z., Bischof, C., Blackford, S., Demmel, J., Dongarra, J., Du Croz, J., Greenbaum, A., Hammarling, S., McKenney, A., Sorensen, D.: LAPACK Users’ Guide, 3rd edn. Society for Industrial and Applied Mathematics, Philadelphia, PA (1999) · Zbl 0934.65030
[5] Argaez, M., Tapia, R., Velazquez, L.: Numerical comparisons of path-following strategies for a primal-dual interior-point method for nonlinear programming. J. Optim. Theory App. 114, 255-272 (2002) · Zbl 1026.90095 · doi:10.1023/A:1016047200413
[6] Armijo, L.: Minimization of functions having Lipschitz-continuous first partial derivatives. Pac. J. Math. 16, 1-3 (1966) · Zbl 0202.46105 · doi:10.2140/pjm.1966.16.1
[7] Balay, S.; Gropp, WD; McInnes, LC; Smith, BF; Bruaset, AM (ed.); Langtangen, HP (ed.); Arge, E. (ed.), Efficient management of parallelism in object oriented numerical software libraries, 163-202 (1997), Cambridge · doi:10.1007/978-1-4612-1986-6_8
[8] Balay, S., Abhyankar, S., Adams, M.F., Brown, J., Brune, P., Buschelman, K., Dalcin, L., Eijkhout, V., Gropp, W.D., Kaushik, D., Knepley, M.G., McInnes, L.C., Rupp, K., Smith, B.F., Zampini, S., Zhang, H., Zhang, H.: PETSc users manual. Tech. Rep. ANL-95/11 - Revision 3.8, Argonne National Laboratory, (2017). http://www.mcs.anl.gov/petsc · Zbl 0478.65030
[9] Benzi, M.; Wathen, AJ; Schilders, WHA (ed.); Vorst, HA (ed.); Rommes, J. (ed.), Some preconditioning techniques for saddle point problems, No. 13, 195-211 (2008), Berlin · doi:10.1007/978-3-540-78841-6_10
[10] Bernstein, D.S.: Matrix Mathematics: Theory, Facts and Formulas, 2nd edn. Princeton University Press, Princeton (2009) · Zbl 1183.15001
[11] Bertocchi, L., Dini, D., Giacopini, M., Fowell, M., Baldini, A.: Fluid film lubrication in the presence of cavitation: a mass-conserving two-dimensional formulation for compressible, piezoviscous and non-Newtonian fluids. Tribol. Int. 67, 61-71 (2013) · doi:10.1016/j.triboint.2013.05.018
[12] Bonettini, S., Galligani, E., Ruggiero, V.: Inner solvers for interior point methods for large scale nonlinear programming. Comput. Optim. Appl. 37, 1-34 (2007) · Zbl 1146.90068 · doi:10.1007/s10589-007-9012-5
[13] Bonneau, D., Hajjam, M.: Modélisation de la rupture et de la reformation des films lubrifiants dans les contacts élastohydrodynamiques. Eur. J. Comput. Mech. 10, 679-704 (2001) · Zbl 1107.76346
[14] Capriz, G., Cimatti, G.: Free boundary problems in the theory of hydrodynamic lubrication: A survey. Tech. Rep. Nota Informatica C81-7, Istituto di Matematica, University of Pisa (1981) · Zbl 0557.76038
[15] Chapiro, G., Gutierrez, A.E., Herskovits, J., Mazorche, S.R., Pereira, W.S.: Numerical solution of a class of moving boundary problems with a nonlinear complementarity approach. J. Optim. Theory Appl. 168(2), 534-550 (2016) · Zbl 1337.65115 · doi:10.1007/s10957-015-0816-7
[16] Cimatti, G., Menchi, O.: On the numerical solution of a variational inequality connected with the hydrodynamica lubrication of a complete journal bearing. Calcolo 15, 249-258 (1978) · Zbl 0432.76036 · doi:10.1007/BF02575916
[17] Cottle, R.W., Dantzig, G.B.: Complementarity pivot theory of mathematical programming. Linear Algebra Appl. 1, 103-125 (1968) · Zbl 0155.28403 · doi:10.1016/0024-3795(68)90052-9
[18] Cottle, R.W., Pang, J.S., Stone, R.E.: The Linear Complementarity Problem. Classics in Applied Mathematics, SIAM (2009) · Zbl 1192.90001
[19] Cryer, C.: The numerical solution of a degenerate variational inequality. In: Parter, S.V. (ed.) Numerical Methods for Partial Differential Equations (1979) · Zbl 0442.73082
[20] Cryer, C., Dempster, A.: Equivalence of linear complementarity problems and linear programs in vector lattice Hilbert spaces. SIAM J. Control Optim. 18(1), 76-90 (1980) · Zbl 0428.90077 · doi:10.1137/0318005
[21] Dembo, R., Eisenstat, S., Steihaug, T.: Inexact Newton methods. SIAM J. Numer. Anal. 19, 400-408 (1982) · Zbl 0478.65030 · doi:10.1137/0719025
[22] Demmel, J.W., Eisenstat, S.C., Gilbert, J.R., Li, X.S., Liu, J.W.H.: A supernodal approach to sparse partial pivoting. SIAM J. Matrix Anal. A 20(3), 720-755 (1999) · Zbl 0931.65022 · doi:10.1137/S0895479895291765
[23] Durazzi, C.: On the Newton interior-point method for nonlinear programming problems. J. Optim. Theory Appl. 104, 73-90 (2000) · Zbl 0960.90093 · doi:10.1023/A:1004624721836
[24] Eisenstat, S., Walker, H.: Globally convergent inexact Newton methods. SIAM J. Optim. 4, 393-422 (1994) · Zbl 0814.65049 · doi:10.1137/0804022
[25] El-Bakri, A., Tapia, R., Tsuchiya, T., Zhang, Y.: On the formulation and theory of the Newton interior point method for nonlinear programming. J. Optim. Theory Appl. 89, 507-541 (1996) · Zbl 0851.90115 · doi:10.1007/BF02275347
[26] Fischer, A.: A special Newton-type optimization method. Optimization 24, 269-284 (1992) · Zbl 0814.65063 · doi:10.1080/02331939208843795
[27] Fokkema, D., Sleijpen, G., Van der Vost, H.: Accelerated inexact Newton schemes for large systems of nonlinear equations. SIAM J. Sci. Comput. 19, 657-674 (1998) · Zbl 0916.65050 · doi:10.1137/S1064827595296148
[28] Galligani, E.: Analysis of the convergence of an inexact Newton method for solving Karush-Kuhn-Tucker systems. In: Atti del Seminario Matematico e Fisico dell’Università di Modena e Reggio Emilia LII, pp. 331-368 (2004) · Zbl 1330.90137
[29] Giacopini, M., Fowell, M., Dini, D., Strozzi, A.: A mass-conserving complementarity formulation to study lubricant films in the presence of cavitation. J. Tribol. 132(4), 041702 (2010) · Zbl 0931.65022
[30] Horn, R., Johnson, C.: Matrix Analysis. Cambridge University Press, Cambridge (1985) · Zbl 0576.15001 · doi:10.1017/CBO9780511810817
[31] Kelley, C.: Iterative Methods for Linear and Nonlinear Equations. Frontiers in Applied Mathematics. SIAM, Philadelphia (1995) · Zbl 0832.65046 · doi:10.1137/1.9781611970944
[32] Khonsari, M., Booser, E.: Applied Tribology: Bearing, Design and Lubrication, 2nd edn. Wiley, Chichester (2008) · doi:10.1002/9780470059456
[33] Kostreva, M.: Elasto-hydrodynamic lubrication: a nonlinear complementarity problem. Int. J. Numer. Meth. Fl 4, 377-397 (1984) · Zbl 0551.76007 · doi:10.1002/fld.1650040407
[34] Lakshmikantham, V., Trigiante, D.: Theory of Difference Equations: Numerical Methods and Applications, 2nd edn. Marcel Dekker Inc., New York (2002) · Zbl 1014.39001 · doi:10.1201/9780203910290
[35] Laratta, A., Menchi, O.: Approssimazione della soluzione di una disequazione variazionale. Applicazione ad un problema di frontiera libera. Calcolo 11, 243-267 (1974) · Zbl 0308.65063 · doi:10.1007/BF02575828
[36] Li, X., Demmel, J., Gilbert, J., Grigori, L., Shao, M., Yamazaki, I.: SuperLU Users’ Guide. Tech. Rep. LBNL-44289, Lawrence Berkeley National Laboratory, http://crd.lbl.gov/ xiaoye/SuperLU/ (1999)
[37] McAllister, G., Rohde, S.: An optimization problem in hydrodynamic lubrication theory. Appl. Math. Optim. 2, 223-235 (1976) · doi:10.1007/BF01464268
[38] Morini, B., Simoncini, V., Tani, M.: A comparison of reduced and unreduced kkt systems arising from interior point methods. Comput. Optim. Appl. 68(1), 1-27 (2017) · Zbl 1406.90088 · doi:10.1007/s10589-017-9907-8
[39] Nocedal, J., Wright, S.: Numerical Optimization. Springer, New York (1999) · Zbl 0930.65067 · doi:10.1007/b98874
[40] Oh, K., Goenka, P.: The elastohydrodynamic solution of journal bearings under dynamic loading. J. Tribol. 107(3), 389-395 (1985) · doi:10.1115/1.3261088
[41] Rheinboldt, W.: Methods for Solving Systems of Nonlinear Equations, 2nd edn. SIAM, Philadelphia (1998) · Zbl 0906.65051 · doi:10.1137/1.9781611970012
[42] Siefert, C., De Sturler, E.: Preconditioners for generalized saddle-point problems. SIAM J. Numer. Anal. 44(3), 1275-1296 (2006) · Zbl 1120.65061 · doi:10.1137/040610908
[43] Süli, E., Mayers, D.: An Introduction to Numerical Analysis. Cambridge University Press, Cambridge (2003) · Zbl 1033.65001 · doi:10.1017/CBO9780511801181
[44] Varga, R.: Matrix Iterative Analysis, 2nd edn. Springer, Berlin (2000) · Zbl 0998.65505 · doi:10.1007/978-3-642-05156-2
[45] Woloszynski, T., Podsiadlo, P., Stachowiak, G.W.: Efficient solution to the cavitation problem in hydrodynamic lubrication. Tribol. Lett. 58, 1-11 (2015) · doi:10.1007/s11249-015-0487-4
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.