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

DOI10.1016/J.JCP.2018.02.010zbMATH Open1391.76734arXiv1709.02581OpenAlexW2752498785MaRDI QIDQ1637564FDOQ1637564

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)

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.

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

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

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)

Cites Work

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

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)