Numerical method for Darcy flow derived using Discrete Exterior Calculus

Mathematics – Numerical Analysis

Scientific paper

Rate now

  [ 0.00 ] – not rated yet Voters 0   Comments 0

Details

Added numerical experiment for flow on a surface. Other small changes in meshing related comments

Scientific paper

We derive a numerical method for Darcy flow, hence also for Poisson's equation in mixed (first order) form, based on discrete exterior calculus (DEC). Exterior calculus is a generalization of vector calculus to smooth manifolds and DEC is one of its discretizations on simplicial complexes such as triangle and tetrahedral meshes. DEC is a coordinate invariant discretization, in that it does not depend on the embedding of the simplices or the whole mesh. We start by rewriting the governing equations of Darcy flow using the language of exterior calculus. This yields a formulation in terms of flux differential form and pressure. The numerical method is then derived by using the framework provided by DEC for discretizing differential forms and operators that act on forms. We also develop a discretization for spatially dependent Hodge star that varies with the permeability of the medium. This also allows us to address discontinuous permeability. The matrix representation for our discrete non-homogeneous Hodge star is diagonal, with positive diagonal entries. The resulting linear system of equations for flux and pressure are saddle type, with a diagonal matrix as the top left block. The performance of the proposed numerical method is illustrated on many standard test problems. These include patch tests in two and three dimensions, comparison with analytically known solution in two dimensions, layered medium with alternating permeability values, and a test with a change in permeability along the flow direction. We also show numerical evidence of convergence of the flux and the pressure. A convergence experiment is also included for Darcy flow on a surface. A short introduction to the relevant parts of smooth and discrete exterior calculus is included in this paper. We also include a discussion of the boundary condition in terms of exterior calculus.

No associations

LandOfFree

Say what you really think

Search LandOfFree.com for scientists and scientific papers. Rate them and share your experience with other people.

Rating

Numerical method for Darcy flow derived using Discrete Exterior Calculus does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Numerical method for Darcy flow derived using Discrete Exterior Calculus, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Numerical method for Darcy flow derived using Discrete Exterior Calculus will most certainly appreciate the feedback.

Rate now

     

Profile ID: LFWR-SCP-O-265066

  Search
All data on this website is collected from public sources. Our data reflects the most accurate information available at the time of publication.