Abstract
In this paper, a fast computational technique based on adaptive mesh generation for numerical solution of the obstacle problem is considered. The obstacle problem is an elliptic variational inequality problem, where its solution divides the domain into the contact and noncontact sets. The boundary between the contact and noncontact sets is a free boundary, which is priori unknown and the solution is not smooth on it. Due to lack of smoothness, numerical methods need a lot of mesh points in discretization for obtaining a numerical solution with a reasonable accuracy. In this paper, using an interpolating wavelet system and the fast wavelet transform, a multi-level algorithm for generating an appropriate-adapted mesh is presented. In each step of the algorithm, the semi-smooth Newton’s method or active set method is used for solving the discretized obstacle problem. We test the performance and accuracy of the proposed method by means of four numerical experiments. We show that the presented method significantly reduces CPU time in comparison with the full-grid algorithm and also can simultaneously capture the priori unknown free boundary.





Similar content being viewed by others
References
Adams DR, Hrynkiv V, Lenhart S (2010) Optimal control of a biharmonic obstacle problem. In: Around the research of vladimir Maz’ya III. Springer, pp 1–24
Akram M, Khan S (2013) Multilayered thresholding-based blood vessel segmentation for screening of diabetic retinopathy. Eng Comput 29(2):165–173
Alhadeff A, Leon S, Celes W, Paulino G (2016) Massively parallel adaptive mesh refinement and coarsening for dynamic fracture simulations. Eng Comput 32(3):533–552
Banz L, Schröder A (2015) Biorthogonal basis functions in hp-adaptive fem for elliptic obstacle problems. Comput Math Appl 70(8):1721–1742
Bertoluzza S (1996) Adaptive wavelet collocation method for the solution of burgers equation. Transp Theory Stat Phys 25(3–5):339–352
Braess D, Carstensen C, Hoppe RH (2007) Convergence analysis of a conforming adaptive finite element method for an obstacle problem. Numerische Mathematik 107(3):455–471
Burman E, Hansbo P, Larson MG, Stenberg R (2017) Galerkin least squares finite element method for the obstacle problem. Comput Methods Appl Mech Eng 313:362–374
Chalfant J, Patrikalakis N (2009) Three-dimensional object registration using wavelet features. Eng Comput 25(3):303–318
Chen X, Qi L, Sun D (1998) Global and superlinear convergence of the smoothing newton method and its application to general box constrained variational inequalities. Math Comput Am Math Soc 67(222):519–540
Comincioli V, Scapolla T, Naldi G, Venini P (2000) A wavelet-like galerkin method for numerical solution of variational inequalities arising in elastoplasticity. Commun Numer Methods Eng 16(2):133–144
De Marchi L, Franzè F, Baravelli E, Speciale N (2006) Wavelet-based adaptive mesh generation for device simulation. Solid State Electron 50(4):650–659
Dehghan M, Abbaszadeh M (2016) A finite element method for the numerical solution of Rayleigh stokes problem for a heated generalized second grade fluid with fractional derivatives. Eng Comput 33(3):587–605
Facchinei F, Pang JS (2007) Finite-dimensional variational inequalities and complementarity problems. Springer, Berlin
Friedman A (2010) Variational principles and free-boundary problems. Dover books on mathematics. Dover, USA
Grossmann C, Roos HG, Stynes M (2007) Numerical treatment of partial differential equations. Springer, Berlin
Hintermüller M, Ito K, Kunisch K (2002) The primal-dual active set strategy as a semismooth newton method. SIAM J Optim 13(3):865–888
Holmström M (1999) Solving hyperbolic pdes using interpolating wavelets. SIAM J Sci Comput 21(2):405–420
Kärkkäinen T, Kunisch K, Tarvainen P (2003) Augmented Lagrangian active set methods for obstacle problems. J Optim Theory Appl 119(3):499–533
Keshavarz E, Ordokhani Y, Razzaghi M (2014) Bernoulli wavelet operational matrix of fractional order integration and its applications in solving the fractional order differential equations. Appl Math Model 38(24):6038–6051
Kinderlehrer D, Stampacchia G (2000) An introduction to variational inequalities and their applications. Classics in applied mathematics. Society for Industrial and Applied Mathematics, Academic Press, New York
Li R, Liu W, Tang T (2000) Moving mesh finite element approximations for variational inequality I: static obstacle problem. Hong Kong Baptist University, Kowloon Tong
Paolucci S, Zikoski ZJ, Wirasaet D (2014) Wamr: An adaptive wavelet method for the simulation of compressible reacting flow. Part I. accuracy and efficiency of algorithm. J Comput Phys 272:814–841
Petrosyan A, Shahgholian H, Uraltseva NN, Shahgholian H (2012) Regularity of free boundaries in obstacle-type problems, vol 136. American Mathematical Society, Providence
Rastigejev YA, Paolucci S (2006) Wavelet-based adaptive multiresolution computation of viscous reactive flows. Int J Numer Methods Fluids 52(7):749–784
Reckinger SM, Vasilyev OV, Fox-Kemper B (2014) Adaptive wavelet collocation method on the shallow water model. J Comput Phys 271:342–359
Roberts D, Hegland M (2011) Solving variational inequalities using wavelet methods. ANZIAM J 52:949–966
Rometsch M, Urban K (2010) Adaptive wavelet methods for elliptic variational inequalitites I: analysis. University of Ulm, Ulm
Tremolieres R, Lions JL, Glowinski R (2011) Numerical analysis of variational inequalities. Elsevier, Amsterdam
Ulbrich M (2011) Semismooth Newton methods for variational inequalities and constrained optimization problems in function spaces, vol 11. SIAM, Seneca
Urban K (2009) Wavelet methods for elliptic partial differential equations. Oxford University Press, Oxford
Vasilyev OV, Kevlahan NKR (2005) An adaptive multilevel wavelet collocation method for elliptic problems. J Comput Phys 206(2):412–431
Vasilyev OV, Paolucci S (1997) A fast adaptive wavelet collocation algorithm for multidimensional PDEs. J Comput Phys 138(1):16–56
Wirasaet D (2007) Numerical solutions of multi-dimensional partial differential equations using an adaptive wavelet method. University of Notre Dame, Notre Dame
Zhang C (2005) Pricing American options by adaptive finite element method. Mathematics Department, University of Maryland
Zhao B (2015) The application of wavelet finite element method on multiple cracks identification of gear pump gear. Eng Comput 31(2):281–288
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Oshagh, M.Ke., Shamsi, M. & Dehghan, M. A wavelet-based adaptive mesh refinement method for the obstacle problem. Engineering with Computers 34, 577–589 (2018). https://doi.org/10.1007/s00366-017-0559-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00366-017-0559-1