Abstract
AbstractIn this work, we present the construction of two distinct finite element approaches to solve the porous medium equation (PME). In the first approach, we transform the PME to a log-density variable formulation and construct a continuous Galerkin method. In the second approach, we introduce additional potential and velocity variables to rewrite the PME into a system of equations, for which we construct a mixed finite element method. Both approaches are first-order accurate, mass conserving, and proved to be unconditionally energy stable for their respective energies. The mixed approach is shown to preserve positivity under a CFL condition, while a much stronger property of unconditional bound preservation is proved for the log-density approach. A novel feature of our schemes is that they can handle compactly supported initial data without the need for any perturbation techniques. Furthermore, the log-density method can handle unstructured grids in any number of dimensions, while the mixed method can handle unstructured grids in two dimensions. We present results from several numerical experiments to demonstrate these properties.
Funder
Division of Mathematical Sciences
Publisher
Springer Science and Business Media LLC
Reference42 articles.
1. Aronson, D.G., Caffarelli, L.A., Kamin, S.: How an initially stationary interface begins to move in porous medium flow. SIAM J. Math. Anal. 14, 639–658 (1983)
2. Baines, M., Hubbard, M., Jimack, P.: A moving mesh finite element algorithm for the adaptive solution of time-dependent partial differential equations with moving boundaries. Appl. Numer. Math. 54, 450–469 (2005). (Selected papers from the 16th Chemnitz Finite Element Symposium (2003))
3. Baranger, J., Maitre, J.-F., Oudin, F.: Connection between finite volume and mixed finite element methods. ESAIM Math. Model. Numer. Anal. 30, 445–465 (1996)
4. Barenblatt, I.S.: On some unsteady motions of a liquid and gas in a porous medium. Akad. Nauk SSSR. Prikl. Mat. Mekh. 16, 67–78 (1952)
5. Barth, T.J.: Aspects of unstructured grids and finite-volume solvers for the Euler and Navier–Stokes equations, AGARD, special course on unstructured grid methods for advection dominated flows (1992)