The porous medium equation (PME) is a nonlinear degenerate equation with the finite speed and possible waiting time. Based on an energetic variational approach, a structure-preserving numerical scheme for the one-dimensional PME has been constructed in [
12
], which preserves the positivity of the solution and the mass conservation, the energy dissipation law, as well as an efficient calculation of the finite speed and the waiting time. In this paper, we use a similar approach to numerically solve a radially symmetric solution of PME in
R
d
,
d
≥
2
. It is proved that the numerical scheme is uniquely solvable on an admissible convex set and satisfies the original discrete energy dissipation law. Moreover, the optimal rate convergence analysis and error estimate is theoretically established. The two and three dimensional simulation results indicate that the numerical method is effective in solving the axisymmetric solutions of the PME equation, and compute the finite speed and the waiting time effectively. As another advantage, we give the convergence order of the radially symmetric solution with a support and the waiting time numerically in two dimension.
- A Structure-Preserving Numerical Method for a Radially Symmetric Solution of Porous Medium Equation
- Qingqing Zhang - Shanghai UniversityChenghua Duan - Shanghai UniversityChun Liu - Illinois Institute of TechnologyCheng Wang - University of Massachusetts DartmouthXingye Yue - Soochow University
- Journal of scientific computing, Vol.108(1)
- Springer US; NEW YORK
- 32
- NSFC 11901109 / National Natural Science Foundation of China (http://dx.doi.org/10.13039/501100001809) DMS-1950868?DMS-2118181; DMS-2012269 and DMS-2309548; NSFC 12371401 / National Science Foundation (http://dx.doi.org/10.13039/100000001)
- Department of Mathematics
- English
- Journal article
- https://doi.org/10.1007/s10915-026-03340-6
- 9914534306501301