<p>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 [<CitationRef CitationID="CR12">12</CitationRef>], 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 <InlineEquation ID="IEq1"> <EquationSource Format="TEX">\(\mathbb {R}^d\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mrow> <mi mathvariant="double-struck">R</mi> </mrow> <mi>d</mi> </msup> </math></EquationSource> </InlineEquation>, <InlineEquation ID="IEq2"> <EquationSource Format="TEX">\(d\ge 2\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>d</mi> <mo>≥</mo> <mn>2</mn> </mrow> </math></EquationSource> </InlineEquation>. 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.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

A Structure-Preserving Numerical Method for a Radially Symmetric Solution of Porous Medium Equation

  • Qingqing Zhang,
  • Chenghua Duan,
  • Chun Liu,
  • Cheng Wang,
  • Xingye Yue

摘要

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 \(\mathbb {R}^d\) R d , \(d\ge 2\) 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.