Numerical artifacts in the generalized porous medium equation: why harmonic averaging itself is not to blame

DOI10.1016/J.JCP.2018.02.010MaRDI QIDQ1637564zbMATH OpenOpenAlexFDO

Authors Danielle C. Maddix, Luiz E. B. Sampaio, Margot Gerritsen

Publication date 8 June 2018

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1709.02581

zbMATH Keywords

degenerate nonlinear parabolic equationstemporal oscillationsgeneralized porous medium equationharmonic and arithmetic averagingmodified equation analysisself-sharpening

Mathematics Subject Classification ID

Degenerate parabolic equations (35K65) Flows in porous media; filtration; seepage (76S05)

Abstract: The degenerate parabolic Generalized Porous Medium Equation (GPME) poses numerical challenges due to self-sharpening and its sharp corner solutions. For these problems, we show results for two subclasses of the GPME with differentiable $k(p)$ with respect to $p$, namely the Porous Medium Equation (PME) and the superslow diffusion equation. Spurious temporal oscillations, and nonphysical locking and lagging have been reported in the literature. These issues have been attributed to harmonic averaging of the coefficient $k(p)$ for small $p$, and arithmetic averaging has been suggested as an alternative. We show that harmonic averaging is not solely responsible and that an improved discretization can mitigate these issues. Here, we investigate the causes of these numerical artifacts using modified equation analysis. The modified equation framework can be used for any type of discretization. We show results for the second order finite volume method. The observed problems with harmonic averaging can be traced to two leading error terms in its modified equation. This is also illustrated numerically through a Modified Harmonic Method (MHM) that can locally modify the critical terms to remove the aforementioned numerical artifacts.

Recommendations

• Numerical artifacts in the discontinuous generalized porous medium equation: how to avoid spurious temporal oscillations
• A numerical procedure for the porous media equation
• Numerical methods with interface estimates for the porous medium equation
• On numerical analyses in the presence of unstable saturated porous materials
• Numerical approximations to interface curves for a porous media equation
• scientific article; zbMATH DE number 1023448
• scientific article; zbMATH DE number 1281313
• Numerical methods for porous medium equation by an energetic variational approach
• Numerical Homogenization of Well Singularities in the Flow Transport through Heterogeneous Porous Media
• scientific article; zbMATH DE number 617812

Cites work

• scientific article; zbMATH DE number 3124687 (Why is no real title available?)
• scientific article; zbMATH DE number 3151882 (Why is no real title available?)
• scientific article; zbMATH DE number 2196537 (Why is no real title available?)
• scientific article; zbMATH DE number 3294428 (Why is no real title available?)
• A stability technique for evolution partial differential equations. A dynamical systems approach
• A study on moving mesh finite element solution of the porous medium equation
• Exact solutions and invariant subspaces of nonlinear partial differential equations in mechanics and physics.
• Interfaces in multidimensional diffusion equations with absorption terms
• METHODS OF CONSTRUCTING APPROXIMATE SELF-SIMILAR SOLUTIONS OF NONLINEAR HEAT EQUATIONS. I
• Numerical artifacts in the discontinuous generalized porous medium equation: how to avoid spurious temporal oscillations
• Numerical artifacts in the generalized porous medium equation: why harmonic averaging itself is not to blame
• Numerical simulation for porous medium equation by local discontinuous Galerkin finite element method
• Regularity Propeties of Flows Through Porous Media
• Systems of conservation laws
• The mimetic finite difference method for elliptic and parabolic problems with a staggered discretization of diffusion coefficient
• The porous medium equation. Mathematical theory
• Total variation diminishing Runge-Kutta schemes

Cited in

(5)

• A decoupled and positivity-preserving DDFVS scheme for diffusion problems on polyhedral meshes
• Regularization errors introduced by the one-fluid formulation in the solution of two-phase elliptic problems
• Numerical artifacts in the discontinuous generalized porous medium equation: how to avoid spurious temporal oscillations
• Numerical artifacts in the generalized porous medium equation: why harmonic averaging itself is not to blame
• Learning physical models that can respect conservation laws

This page was built for publication: Numerical artifacts in the generalized porous medium equation: why harmonic averaging itself is not to blame

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1637564)