A wavelet-based adaptive mesh refinement method for the obstacle problem

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.

[1]  Martin Stynes,et al.  Numerical Treatment of Partial Differential Equations , 2007 .

[2]  Carsten Carstensen,et al.  Convergence analysis of a conforming adaptive finite element method for an obstacle problem , 2007, Numerische Mathematik.

[3]  Bin Zhao.  The application of wavelet finite element method on multiple cracks identification of gear pump gear , 2013, Engineering with Computers.

[4]  S. Lenhart,et al.  Optimal Control of a Biharmonic Obstacle Problem , 2010 .

[5]  A. Petrosyan,et al.  Regularity of Free Boundaries in Obstacle-type Problems , 2012 .

[6]  Silvia Bertoluzza.  Adaptive wavelet collocation method for the solution of Burgers equation , 1996 .

[7]  WENBIN LIU,et al.  MOVING MESH FINITE ELEMENT APPROXIMATIONS FOR VARIATIONAL INEQUALITY I: STATIC OBSTACLE PROBLEM , 2000 .

[8]  Nicholas M. Patrikalakis,et al.  Three-dimensional object registration using wavelet features , 2009, Engineering with Computers.

[9]  O. Vasilyev,et al.  A Fast Adaptive Wavelet Collocation Algorithm for Multidimensional PDEs , 1997 .

[10]  Kazufumi Ito,et al.  The Primal-Dual Active Set Strategy as a Semismooth Newton Method , 2002, SIAM J. Optim..

[11]  R. Glowinski,et al.  Numerical Analysis of Variational Inequalities , 1981 .

[12]  Oleg V. Vasilyev,et al.  Adaptive wavelet collocation method on the shallow water model , 2014, J. Comput. Phys..

[13]  Glaucio H. Paulino,et al.  Massively parallel adaptive mesh refinement and coarsening for dynamic fracture simulations , 2015, Engineering with Computers.

[14]  A wavelet‐like Galerkin method for numerical solution of variational inequalities arising in elastoplasticity , 2000 .

[15]  Mats Holmström,et al.  Solving Hyperbolic PDEs Using Interpolating Wavelets , 1999, SIAM J. Sci. Comput..

[16]  Nicholas K.-R. Kevlahan,et al.  An adaptive multilevel wavelet collocation method for elliptic problems , 2005 .

[17]  Solving variational inequalities using wavelet methods , 2011 .

[18]  Andreas Schröder,et al.  Biorthogonal basis functions in hp-adaptive FEM for elliptic obstacle problems , 2015, Comput. Math. Appl..

[19]  D. Kinderlehrer,et al.  An introduction to variational inequalities and their applications , 1980 .

[20]  Damrongsak Wirasaet,et al.  Numerical solutions of multi-dimensional partial differential equations using an adaptive wavelet method , 2007 .

[21]  Peter Hansbo,et al.  Galerkin least squares finite element method for the obstacle problem , 2016, 1609.03431.

[22]  Nicolò Speciale,et al.  Wavelet-based adaptive mesh generation for device simulation , 2005 .

[23]  M. Usman Akram,et al.  Multilayered thresholding-based blood vessel segmentation for screening of diabetic retinopathy , 2011, Engineering with Computers.

[24]  Karsten Urban,et al.  Wavelet Methods for Elliptic Partial Differential Equations , 2008 .

[25]  Samuel Paolucci,et al.  WAMR: An adaptive wavelet method for the simulation of compressible reacting flow. Part I. Accuracy and efficiency of algorithm , 2014, J. Comput. Phys..

[26]  Mehdi Dehghan,et al.  A finite element method for the numerical solution of Rayleigh–Stokes problem for a heated generalized second grade fluid with fractional derivatives , 2017, Engineering with Computers.

[27]  Samuel Paolucci,et al.  Wavelet‐based adaptive multiresolution computation of viscous reactive flows , 2006 .

[28]  Karl Kunisch,et al.  Augmented Lagrangian Active Set Methods for Obstacle Problems , 2003 .

[29]  Yadollah Ordokhani,et al.  Bernoulli wavelet operational matrix of fractional order integration and its applications in solving the fractional order differential equations , 2014 .

[30]  Xiaojun Chen,et al.  Global and superlinear convergence of the smoothing Newton method and its application to general box constrained variational inequalities , 1998, Math. Comput..

[31]  Michael Ulbrich,et al.  Semismooth Newton Methods for Variational Inequalities and Constrained Optimization Problems in Function Spaces , 2011, MOS-SIAM Series on Optimization.

[32]  F. Facchinei,et al.  Finite-Dimensional Variational Inequalities and Complementarity Problems , 2003 .

[33]  A. Friedman.  Variational principles and free-boundary problems , 1982 .

[34]  M. Rometsch.  ADAPTIVE WAVELET METHODS FOR ELLIPTIC VARIATIONAL INEQUALITIES I: ANALYSIS , 2010 .