Numerical Experiments for Darcy Flow on a Surface Using Mixed Exterior Calculus Methods

There are very few results on mixed finite element methods on surfaces. A theory for the study of such methods was given recently by Holst and Stern, using a variational crimes framework in the context of finite element exterior calculus. However, we are not aware of any numerical experiments where mixed finite elements derived from discretizations of exterior calculus are used for a surface domain. This short note shows results of our preliminary experiments using mixed methods for Darcy flow (hence scalar Poisson's equation in mixed form) on surfaces. We demonstrate two numerical methods. One is derived from the primal-dual Discrete Exterior Calculus and the other from lowest order finite element exterior calculus. The programming was done in the language Python, using the PyDEC package which makes the code very short and easy to read. The qualitative convergence studies seem to be promising.Comment: 14 pages, 11 figure

PyDEC: Software and Algorithms for Discretization of Exterior Calculus

This paper describes the algorithms, features and implementation of PyDEC, a Python library for computations related to the discretization of exterior calculus. PyDEC facilitates inquiry into both physical problems on manifolds as well as purely topological problems on abstract complexes. We describe efficient algorithms for constructing the operators and objects that arise in discrete exterior calculus, lowest order finite element exterior calculus and in related topological problems. Our algorithms are formulated in terms of high-level matrix operations which extend to arbitrary dimension. As a result, our implementations map well to the facilities of numerical libraries such as NumPy and SciPy. The availability of such libraries makes Python suitable for prototyping numerical methods. We demonstrate how PyDEC is used to solve physical and topological problems through several concise examples.Comment: Revised as per referee reports. Added information on scalability, removed redundant text, emphasized the role of matrix based algorithms, shortened length of pape

Well-Centered Triangulation

Meshes composed of well-centered simplices have nice orthogonal dual meshes (the dual Voronoi diagram). This is useful for certain numerical algorithms that prefer such primal-dual mesh pairs. We prove that well-centered meshes also have optimality properties and relationships to Delaunay and minmax angle triangulations. We present an iterative algorithm that seeks to transform a given triangulation in two or three dimensions into a well-centered one by minimizing a cost function and moving the interior vertices while keeping the mesh connectivity and boundary vertices fixed. The cost function is a direct result of a new characterization of well-centeredness in arbitrary dimensions that we present. Ours is the first optimization-based heuristic for well-centeredness, and the first one that applies in both two and three dimensions. We show the results of applying our algorithm to small and large two-dimensional meshes, some with a complex boundary, and obtain a well-centered tetrahedralization of the cube. We also show numerical evidence that our algorithm preserves gradation and that it improves the maximum and minimum angles of acute triangulations created by the best known previous method.Comment: Content has been added to experimental results section. Significant edits in introduction and in summary of current and previous results. Minor edits elsewher

Discretizations and Solvers for Coupling Stokes-Darcy Flows With Transport

This thesis studies a mathematical model, in which Stokes-Darcy flow system is coupled with a transport equation. The objective is to develop stable and convergent numerical schemes that could be used in environmental applications. Special attention is given to discretization methods that conserve mass locally. First, we present a global saddle point problem approach, which employs the discontinuous Galerkin method to discretize the Stokes equations and the mimetic finite difference method to discretize the Darcy equation. We show how the numerical scheme can be formulated on general polygonal (polyhedral in three dimensions) meshes if suitable operators mapping from degrees of freedom to functional spaces are constructed. The scheme is analyzed and error estimates are derived. A hybridization technique is used to solve the system effectively. We ran several numerical experiments to verify the theoretical convergence rates and depending on the mesh type we observed superconvergence of the computed solution in the Darcy region.Another approach that we use to deal with the flow equations is based on non-overlapping domain decomposition. Domain decomposition enables us to solve the coupled Stokes-Darcy flow problem in parallel by partitioning the computational domain into subdomains, upon which families of coupled local problems of lower complexity are formulated. The coupling of the subdomain problems is removed through an iterative procedure. We investigate the properties of this method and derive estimates for the condition number of the associated algebraic system. Results from computer tests supporting the convergence analysis of the method are provided. To discretize the transport equation we use the local discontinuous Galerkin (LDG) method, which can be thought as a discontinuous mixed finite element method, since it approximates both the concentration and the diffusive flux. We develop stability and convergence analysis for the concentration and the diffusive flux in the transport equation. The numerical error is a combination of the LDG discretization error and the error from the discretization of the Stokes-Darcy velocity. Several examples verifying the theory and illustrating the capabilities of the method are presented

Nonstandard Finite Element Methods

[no abstract available

Variational formulation of the finite calculus equations in solid mechanics and diffusion-reaction problems

We present a variational formulation of the finite calculus (FIC) equations for problems in mechanics governed by differential equations with symmetric operators. Applications considered include solid mechanics, diffusion-transport and diffusion-reaction problems. The key of the variational formulation is the identification of the FIC governing equations with the classical differential equations of mechanics written in terms of modified non-local variables. A total potential energy (TPE) functional is found in terms of the modified variables. The FIC equations in the domain and the boundary are recovered as the Euler-Lagrange equations and the natural boundary condition of the TPE functional, respectively. Symmetric finite element equations are obtained after discretization of the TPE functional, therefore preserving the symmetry of the governing infinitesimal equations. The variational FIC expression is reinterpreted as a Petrov Galerkin weighted residual form of the original FIC equations with non-local weighting functions. The analogy of the variational FIC-FEM formulation with a discontinuous Galerkin method is recognized. Extensions to multidimensional linear elastostatics and diffusion-reaction problems are presented