The cutoff method for the numerical computation of nonnegative solutions of parabolic PDEs with application to anisotropic diffusion and lubrication-type equations

Abstract

The cutoff method, which cuts off the values of a function less than a given number, is studied for the numerical computation of nonnegative solutions of parabolic partial differential equations. A convergence analysis is given for a broad class of finite difference methods combined with cutoff for linear parabolic equations. Two applications are investigated, linear anisotropic diffusion problems satisfying the setting of the convergence analysis and nonlinear lubrication-type equations for which it is unclear if the convergence analysis applies. The numerical results are shown to be consistent with the theory and in good agreement with existing results in the literature. The convergence analysis and applications demonstrate that the cutoff method is an effective tool for use in the computation of nonnegative solutions. Cutoff can also be used with other discretization methods such as collocation, finite volume, finite element, and spectral methods and for the computation of positive solutions.

arXiv comment19 pages, 41 figures

References in corpus (5)

• Preserving Monotonicity in Anisotropic Diffusion
• An anisotropic mesh adaptation method for the finite element solution of heterogeneous anisotropic diffusion problems
• Maximum principle in linear finite element approximations of anisotropic diffusion-convection-reaction problems
• Maximum principle for the finite element solution of time dependent anisotropic diffusion problems
• Discrete maximum principle and a Delaunay-type mesh condition for linear finite element approximations of two-dimensional anisotropic diffusion problems

Cited by in corpus (10)

• A new Lagrange multiplier approach for constructing structure preserving schemes, I. positivity preserving
• A Note on why Enforcing Discrete Maximum Principles by a simple a Posteriori Cutoff is a Good Idea
• A moving-boundary model of reactive settling in wastewater treatment. Part 2: Numerical scheme
• Anisotropic Mesh Adaptation for Finite Element Solution of Anisotropic Porous Medium Equation
• Arbitrarily high-order exponential cut-off methods for preserving maximum principle of parabolic equations
• Moving mesh finite difference solution of non-equilibrium radiation diffusion equations
• A Relaxed Lagrange Multiplier Approach for Phase Field Models
• A positivity-preserving finite volume element method for anisotropic diffusion problems on quadrilateral meshes
• $Q$-tensor gradient flow with quasi-entropy and discretizations preserving physical constraints
• On mesh restrictions to satisfy comparison principles, maximum principles, and the non-negative constraint: Recent developments and new results