Discrete analysis of domain decomposition approaches for mesh generation via the equidistribution principle

Moving mesh methods based on the equidistribution principle are powerful techniques for the space-time adaptive solution of evolution problems. Solving the resulting coupled system of equations, namely the original PDE and the mesh PDE, however, is challenging in parallel. Recently several Schwarz domain decomposition algorithms were proposed for this task and analyzed at the continuous level. However, after discretization, the resulting problems may not even be well posed, so the discrete algorithms requires a different analysis, which is the subject of this paper. We prove that when the number of grid points is large enough, the classical parallel and alternating Schwarz methods converge to the unique monodomain solution. Thus, such methods can be used in place of Newton’s method, which can suffer from convergence difficulties for challenging problems. The analysis for the nonlinear domain decomposition algorithms is based on M–function theory and is valid for an arbitrary number of subdomains. An asymptotic convergence rate is provided and numerical experiments illustrate the results.

[1]  Ronald D. Haynes,et al.  Alternating Schwarz methods for partial differential equation-based mesh generation , 2015, Int. J. Comput. Math..

[2]  Andrea Toselli,et al.  Domain decomposition methods : algorithms and theory , 2005 .

[3]  Weizhang Huang.  Practical aspects of formulation and solution of moving mesh partial differential equations , 2001 .

[4]  W. Rheinboldt,et al.  Error Estimates for Adaptive Finite Element Computations , 1978 .

[5]  Oliver Sander,et al.  Substructuring of a Signorini-type problem and Robin’s method for the Richards equation in heterogeneous soil , 2010, Comput. Vis. Sci..

[6]  M. Gander,et al.  Absorbing boundary conditions for the wave equation and parallel computing , 2004, Math. Comput..

[7]  Xiao-Chuan Cai,et al.  A class of parallel two-level nonlinear Schwarz preconditioned inexact Newton algorithms ☆ , 2007 .

[8]  Laurence Halpern,et al.  Optimized Schwarz Waveform Relaxation for Nonlinear Systems of Parabolic Type , 2014 .

[9]  Robert D. Russell,et al.  Convergence of de Boor's algorithm for the generation of equidistributing meshes , 2011 .

[10]  Weizhang Huang,et al.  Analysis Of Moving Mesh Partial Differential Equations With Spatial Smoothing , 1997 .

[11]  T. Politi,et al.  Some improvements for two-sided bounds on the inverse of diagonally dominant tridiagonal matrices , 2001 .

[12]  R. Kellogg.  A nonlinear alternating direction method , 1969 .

[13]  Peter Kuster,et al.  Domain Decomposition Methods In Science And Engineering Xix , 2016 .

[14]  Xiao-Chuan Cai,et al.  Additive Schwarz algorithms for parabolic convection-diffusion equations , 1991 .

[15]  L. Demkowicz.  One and two dimensional elliptic and Maxwell problems , 2006 .

[16]  Uri M. Ascher,et al.  DAEs That Should Not Be Solved , 2000 .

[17]  Lothar Collatz,et al.  Aufgaben monotoner Art , 1952 .

[18]  Ronald D. Haynes,et al.  Parallel stochastic methods for PDE based grid generation , 2013, Comput. Math. Appl..

[19]  Martin J. Gander,et al.  Alternating and Linearized Alternating Schwarz Methods for Equidistributing Grids , 2013, Domain Decomposition Methods in Science and Engineering XX.

[20]  Tarek P. Mathew,et al.  Domain Decomposition Methods for the Numerical Solution of Partial Differential Equations , 2008, Lecture Notes in Computational Science and Engineering.

[21]  Danping Yang,et al.  Additive Schwarz methods for parabolic problems , 2005, Appl. Math. Comput..

[22]  Frédéric Magoulès,et al.  A Schur Complement Method for Compressible Navier-Stokes Equations , 2013, Domain Decomposition Methods in Science and Engineering XX.

[23]  Martin J. Gander,et al.  Overlapping Schwarz Waveform Relaxation for Convection-Dominated Nonlinear Conservation Laws , 2005, SIAM J. Sci. Comput..

[24]  P. Bassanini,et al.  Elliptic Partial Differential Equations of Second Order , 1997 .

[25]  Martin J. Gander.  A Waveform Relaxation Algorithm with Overlapping Splitting for Reaction Diffusion Equations , 1999 .

[26]  Xiao-Chuan Cai,et al.  Domain Decomposition Methods for Monotone Nonlinear Elliptic Problems , 1994 .

[27]  Jinchao Xu,et al.  Domain Decomposition Methods in Scientific and Engineering Computing , 1994 .

[28]  Gérard Meurant,et al.  A Review on the Inverse of Symmetric Tridiagonal and Block Tridiagonal Matrices , 1992, SIAM J. Matrix Anal. Appl..

[29]  A. Ramage,et al.  Computational solution of two-dimensional unsteady PDEs using moving mesh methods , 2002 .

[30]  John A. Mackenzie,et al.  A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations , 2007, SIAM J. Sci. Comput..

[31]  Paul Andries Zegeling,et al.  Balanced monitoring of flow phenomena in moving mesh methods , 2009 .

[32]  Martin J. Gander,et al.  Optimized Schwarz Waveform Relaxation Methods for Advection Reaction Diffusion Problems , 2007, SIAM J. Numer. Anal..

[33]  Martin J. Gander,et al.  Optimized Schwarz Methods , 2006, SIAM J. Numer. Anal..

[34]  Arthur van Dam,et al.  A robust moving mesh finite volume method applied to 1D hyperbolic conservation laws from magnetohydrodynamics , 2006, J. Comput. Phys..

[35]  Herbert B. Keller,et al.  Space-time domain decomposition for parabolic problems , 2002, Numerische Mathematik.

[36]  W. Hackbusch,et al.  On the nonlinear domain decomposition method , 1997 .

[37]  Xiao-Chuan Cai,et al.  Multiplicative Schwarz Methods for Parabolic Problems , 1994, SIAM J. Sci. Comput..

[38]  Zi-Cai Li,et al.  Schwarz Alternating Method , 1998 .

[39]  L. Petzold.  A description of dassl: a differential/algebraic system solver , 1982 .

[40]  H. H. Rachford,et al.  The Numerical Solution of Parabolic and Elliptic Differential Equations , 1955 .

[41]  Xuecheng Tai,et al.  Rate of Convergence of Some Space Decomposition Methods for Linear and Nonlinear Problems , 1998 .

[42]  Robert D. Russell,et al.  Adaptivity with moving grids , 2009, Acta Numerica.

[43]  HuangWeizhang,et al.  Moving mesh partial differential equations (MMPDES) based on the equidistribution principle , 1994 .

[44]  S. H. Lui,et al.  On linear monotone iteration and Schwarz methods for nonlinear elliptic PDEs , 2002, Numerische Mathematik.

[45]  Robert D. Russell,et al.  A Schwarz Waveform Moving Mesh Method , 2007, SIAM J. Sci. Comput..

[46]  Maciej Paszyński,et al.  Computing with hp-ADAPTIVE FINITE ELEMENTS: Volume II Frontiers: Three Dimensional Elliptic and Maxwell Problems with Applications , 2007 .

[47]  David E. Keyes,et al.  Nonlinearly Preconditioned Inexact Newton Algorithms , 2002, SIAM J. Sci. Comput..

[48]  Olaf Steinbach,et al.  FETI Methods for the Simulation of Biological Tissues , 2013, Domain Decomposition Methods in Science and Engineering XX.

[49]  John Mackenzie,et al.  An analysis of stability and convergence of a finite-difference discretization of a model parabolic PDE in 1D using a moving mesh , 2006 .

[50]  Martin J. Gander,et al.  Domain Decomposition Approaches for Mesh Generation via the Equidistribution Principle , 2012, SIAM J. Numer. Anal..

[51]  Baodong Zheng,et al.  Spectral radius and infinity norm of matrices , 2008 .

[52]  Xiao-Chuan Cai,et al.  One‐level Newton–Krylov–Schwarz algorithm for unsteady non‐linear radiation diffusion problem , 2004, Numer. Linear Algebra Appl..

[53]  Weizhang Huang,et al.  A Simple Adaptive Grid Method in Two Dimensions , 1994, SIAM J. Sci. Comput..

[54]  Martin J. Gander,et al.  A homographic best approximation problem with application to optimized Schwarz waveform relaxation , 2009, Math. Comput..

[55]  Xiao-Chuan Cai,et al.  A parallel nonlinear additive Schwarz preconditioned inexact Newton algorithm for incompressible Navier-Stokes equations , 2005 .

[56]  Bevan K. Youse,et al.  Introduction to real analysis , 1972 .

[57]  Robert D. Russell,et al.  Adaptive Moving Mesh Methods , 2010 .

[58]  D. Keyes,et al.  Non‐linear additive Schwarz preconditioners and application in computational fluid dynamics , 2002 .

[59]  W. Rheinboldt.  On M-functions and their application to nonlinear Gauss-Seidel iterations and to network flows☆ , 1970 .

[60]  Tao Tang,et al.  Moving Mesh Methods for Computational Fluid Dynamics , 2022 .

[61]  Weizhang Huang,et al.  Moving mesh partial differential equations (MMPDES) based on the equidistribution principle , 1994 .

[62]  C. D. Boor,et al.  Good approximation by splines with variable knots. II , 1974 .

[63]  Ruo Li,et al.  Moving Mesh Finite Element Methods for the Incompressible Navier-Stokes Equations , 2005, SIAM J. Sci. Comput..

[64]  Robert D. Russell,et al.  A Moving Mesh Method for Time—dependent Problems Based on Schwarz Waveform Relaxation , 2008 .

[65]  Martin J. Gander,et al.  A waveform relaxation algorithm with overlapping splitting for reaction diffusion equations , 1999, Numer. Linear Algebra Appl..

[66]  James M. Ortega,et al.  Iterative solution of nonlinear equations in several variables , 2014, Computer science and applied mathematics.

[67]  Martin J. Gander,et al.  Space-Time Continuous Analysis of Waveform Relaxation for the Heat Equation , 1998, SIAM J. Sci. Comput..

[68]  Oliver Sander.  Coupling Geometrically Exact Cosserat Rods and Linear Elastic Continua , 2013, Domain Decomposition Methods in Science and Engineering XX.

[69]  Martin J. Gander,et al.  Optimal Schwarz Waveform Relaxation for the One Dimensional Wave Equation , 2003, SIAM J. Numer. Anal..

[70]  Weizhang Huang,et al.  A moving collocation method for solving time dependent partial differential equations , 1996 .

[71]  S. Lui,et al.  On Monotone and Schwarz Alternating Methods for Nonlinear Elliptic PDEs , 2001 .

[72]  Anthony Michel,et al.  Newton-Schwarz Optimised Waveform Relaxation Krylov Accelerators for Nonlinear Reactive Transport , 2013, Domain Decomposition Methods in Science and Engineering XX.

[73]  H. G. Burchard,et al.  Splines (with optimal knots) are better , 1974 .