A positivity-preserving unigrid method for elliptic PDEs

While constraints arise naturally in many physical models, their treatment in mathematical and numerical models varies widely, depending on the nature of the constraint and the availability of simulation tools to enforce it. In this paper, we consider the solution of discretized PDE models that have a natural constraint on the positivity (or non-negativity) of the solution. While discretizations of such models often offer analogous positivity properties on their exact solutions, the use of approximate solution algorithms (and the unavoidable effects of floating -- point arithmetic) often destroy any guarantees that the computed approximate solution will satisfy the (discretized form of the) physical constraints, unless the discrete model is solved to much higher precision than discretization error would dictate. Here, we introduce a class of iterative solution algorithms, based on the unigrid variant of multigrid methods, where such positivity constraints can be preserved throughout the approximate solution process. Numerical results for one- and two-dimensional model problems show both the effectiveness of the approach and the trade-off required to ensure positivity of approximate solutions throughout the solution process.

[1]  Paul J. Schweitzer,et al.  A Survey of Aggregation-Disaggregation in Large Markov Chains , 2021 .

[2]  Albert Y. Zomaya,et al.  Partial Differential Equations , 2005, Explorations in Numerical Analysis.

[3]  Michael Griebel,et al.  Stochastic subspace correction methods and fault tolerance , 2018, Math. Comput..

[4]  P. Oswald,et al.  Greedy and Randomized Versions of the Multiplicative Schwarz Method , 2012 .

[5]  Andrea L. Bertozzi,et al.  A numerical scheme for particle-laden thin film flow in two dimensions , 2011, J. Comput. Phys..

[6]  Weizhang Huang,et al.  Adaptive Moving Mesh Methods , 2010 .

[7]  Irad Yavneh,et al.  Square and stretch multigrid for stochastic matrix eigenproblems , 2010, Numer. Linear Algebra Appl..

[8]  Thomas A. Manteuffel,et al.  Algebraic Multigrid for Markov Chains , 2010, SIAM J. Sci. Comput..

[9]  Thomas A. Manteuffel,et al.  Smoothed Aggregation Multigrid for Markov Chains , 2010, SIAM J. Sci. Comput..

[10]  B. Jovanovic,et al.  Increasing of the Accuracy in the Computation of the Concentrations in Diffusion Models with Localized Chemical Reactions , 2008 .

[11]  Thomas A. Manteuffel,et al.  Projection Multilevel Methods for Quasilinear Elliptic Partial Differential Equations: Numerical Results , 2006, SIAM J. Numer. Anal..

[12]  Andrea L. Bertozzi,et al.  Positivity-Preserving Numerical Schemes for Lubrication-Type Equations , 1999, SIAM J. Numer. Anal..

[13]  Udo R. Krieger,et al.  On a two-level multigrid solution method for finite Markov chains , 1995 .

[14]  Wei Wu,et al.  Numerical Experiments with Iteration and Aggregation for Markov Chains , 1992, INFORMS J. Comput..

[15]  Kyle W. Kindle,et al.  An iterative aggregation-disaggregation algorithm for solving linear equations , 1986 .

[16]  S. F. McCormick,et al.  Unigrid for Multigrid Simulation , 1983 .

[17]  W. R. Holland,et al.  Unigrid methods for boundary value problems with nonrectangular domains , 1982 .

[18]  Robert J. Plemmons,et al.  Nonnegative Matrices in the Mathematical Sciences , 1979, Classics in Applied Mathematics.

[19]  J. A. White,et al.  On Selection of Equidistributing Meshes for Two-Point Boundary-Value Problems , 1979 .

[20]  A. Bertozzi THE MATHEMATICS OF MOVING CONTACT LINES IN THIN LIQUID FILMS , 1998 .

[21]  U. Krieger Numerical Solution of Large Finite Markov Chains by Algebraic Multigrid Techniques , 1995 .

[22]  William L. Briggs,et al.  A multigrid tutorial , 1987 .