TY - JOUR

T1 - Convergence analysis of a numerical scheme for the porous medium equation by an energetic variational approach

AU - Duan, Chenghua

AU - Liu, Chun

AU - Wang, Cheng

AU - Yue, Xingye

N1 - Funding Information:
Acknowledgements The work of Yue is supported in part by NSF of China under the grants No. 11971342. Chun Liu and Cheng Wang are partially supported by NSF grants DMS-1216938, DMS-1418689, respectively.
Publisher Copyright:
© 2020 Global-Science Press.

PY - 2020/2

Y1 - 2020/2

N2 - The porous medium equation (PME) is a typical nonlinear degenerate parabolic equation. We have studied numerical methods for PME by an energetic variational approach in [C. Duan et al., J. Comput. Phys., 385 (2019), pp. 13-32], where the trajectory equation can be obtained and two numerical schemes have been developed based on different dissipative energy laws. It is also proved that the nonlinear scheme, based on f log f as the total energy form of the dissipative law, is uniquely solvable on an admissible convex set and preserves the corresponding discrete dissipation law. Moreover, under certain smoothness assumption, we have also obtained the second order convergence in space and the first order convergence in time for the scheme. In this paper, we provide a rigorous proof of the error estimate by a careful higher order asymptotic expansion and two step error estimates. The latter technique contains a rough estimate to control the highly nonlinear term in a discrete W1;1 norm and a refined estimate is applied to derive the optimal error order.

AB - The porous medium equation (PME) is a typical nonlinear degenerate parabolic equation. We have studied numerical methods for PME by an energetic variational approach in [C. Duan et al., J. Comput. Phys., 385 (2019), pp. 13-32], where the trajectory equation can be obtained and two numerical schemes have been developed based on different dissipative energy laws. It is also proved that the nonlinear scheme, based on f log f as the total energy form of the dissipative law, is uniquely solvable on an admissible convex set and preserves the corresponding discrete dissipation law. Moreover, under certain smoothness assumption, we have also obtained the second order convergence in space and the first order convergence in time for the scheme. In this paper, we provide a rigorous proof of the error estimate by a careful higher order asymptotic expansion and two step error estimates. The latter technique contains a rough estimate to control the highly nonlinear term in a discrete W1;1 norm and a refined estimate is applied to derive the optimal error order.

UR - http://www.scopus.com/inward/record.url?scp=85085964151&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85085964151&partnerID=8YFLogxK

U2 - 10.4208/NMTMA.OA-2019-0073

DO - 10.4208/NMTMA.OA-2019-0073

M3 - Article

AN - SCOPUS:85085964151

VL - 13

SP - 63

EP - 80

JO - Numerical Mathematics

JF - Numerical Mathematics

SN - 1004-8979

IS - 1

ER -