Physics-Compatible Numerical Methods for Flow through Porous Media

Abstract

In this thesis report, we investigate the difference in using physics-compatible elements and standard elements in the context of fluid flow through porous media. For this, we compare the H(div)-conforming Raviart-Thomas elements to Lagrange elements. Our model involves a channel flow over a porous media, resulting in a fluid and a porous subdomain. We distinguish between three approaches, as based on a research paper. The first is the one-domain approach called PE, which sets one equation on the entire domain. The other two approaches are two-domain approaches, which use different equations in the two subdomains, with interface conditions to couple them. The first of these is NSD, which couples Navier-Stokes with the Darcy equation. The second is NSF, which does the same but also incorporates the Forchheimer term, which models inertial effects.

In general, we were not able to make a completely successful model, as the PE case did not match the paper, and the results found when using Raviart-Thomas elements had some unexplained anomalies. We were also unable to carry out an accurate error analysis. We were able to see that Raviart-Thomas produced an exactly divergence-free solution, as expected from the theory, at the cost of more computing power and time needed.