In this work we use algebraic dual spaces with a domain decomposition method to solve the Darcy equations. We define the broken Sobolev spaces and their finite dimensional counterparts. A global trace space is defined that connects the solution between the broken spaces. Use of algebraic dual spaces results in a sparse, metric-free representation of the incompressibility constraint, the pressure gradient term, and on the continuity constraint between the sub domains. To demonstrate this, we solve two test cases: (i) a manufactured solution case, and (ii) an industrial benchmark reservoir modelling problem SPE10. The results demonstrate that the dual spaces can be used for domain decomposition formulation, and despite having more unknowns, requires less simulation time compared to the continuous Galerkin formulation, without compromising on the accuracy of the solution.

In ℝⁿ, let Λ^k(Ω) represent the space of smooth differential k-forms in Ω. The de Rham complex consists of a sequence of spaces, Λ^k(Ω), k = 0, 1…, n, connected by the exterior derivative, d: Λ^k(Ω) → Λ^k+1(Ω). Appropriately chosen B-spline spaces together with their associated dual B-spline spaces form a discrete double de Rham complex. In practical applications, this discrete double de Rham complex leads to very sparse systems. In this paper, this construction will be explained and illustrated by means of a non-trivial three-dimensional example.

Given a sequence of finite element spaces which form a de Rham sequence, we will construct dual representations of these spaces with associated differential operators which connect these spaces such that they also form a de Rham sequence. The dual representations also need to satisfy the de Rham sequence on the domain boundary. The matrix which converts primal representations to dual representations – the Hodge matrix – is the mass or Gram matrix. It will be shown that a bilinear form of a primal and a dual representation is equal to the vector inner product of the expansion coefficients (degrees of freedom) of both representations. This leads to very sparse system matrices, even for high order methods. The differential operators for dual representations will be defined. Vector operations, grad, curl and div, for primal and dual representations are both topological and do not depend on the metric, i.e. the size and shape of the mesh or the order of the numerical method. Derivatives are evaluated by applying sparse incidence and inclusion matrices to the expansion coefficients of the representations. As illustration of the use of dual representations, the method will be applied to (i) a mixed formulation for the Poisson problem in 3D, (ii) it will be shown that this approach allows one to preserve the equivalence between Dirichlet and Neumann problems in the finite dimensional setting and, (iii) the method will be applied to the approximation of grad–div eigenvalue problem on affine and non-affine meshes.

We present a hybrid mimetic spectral element formulation for Darcy flow. The discrete representations for (1) conservation of mass, and (2) inter-element continuity, are topological relations that lead to sparse matrix systems. These constraints are independent of the element size and shape, and thus invariant under mesh transformations. The resultant algebraic system is extremely sparse even for high degree polynomial basis. Furthermore, the system can be efficiently assembled and solved for each element separately.

In this paper, we will show that the equivalence of a div-grad Neumann problem and a grad-div Dirichlet problem can be preserved at the discrete level in 3-dimensional curvilinear domains if algebraic dual polynomial representations are employed. These representations will be introduced. Proof of the equivalence at the discrete level follows from the construction of the algebraic dual representations. A 3-dimensional test problem in curvilinear coordinates will illustrate this approach.

In this paper, we present a hybrid mimetic method which solves the mixed formulation of the Poisson problem on curvilinear quadrilateral meshes. The method is hybrid in the sense that the domain is decomposed into multiple disjoint elements and the interelement continuity is enforced using a Lagrange multiplier. The method is mimetic in the sense that the discrete divergence operator is exact. By using the mimetic basis functions and their dual representations, various metric-free discrete terms are obtained. The discrete system can be efficiently solved by first solving a reduced system for the Lagrange multiplier. Numerical experiments which validate the method are presented.

In this paper we will consider two curl-curl equations in two dimensions. One curl-curl problem for a scalar quantity F and one problem for a vector field E. For Dirichlet boundary conditions n× E= Ê_⊣ on E and Neumann boundary conditions n×curlF=Ê⊣, we expect the solutions to satisfy E = curl F. When we use algebraic dual polynomial representations, these identities continue to hold at the discrete level. Equivalence will be proved and illustrated with a computational example.

In this paper, we will use algebraic dual polynomials to set up a discrete Steklov-Poincaré operator for the mixed formulation of the Poisson problem. The method will be applied in curvilinear coordinates and to a test problem which contains a singularity. Exponential convergence of the trace variable in H 1 / 2 {H^{1/2}} -norm will be shown.

This chapter addresses the topological structure of steady, anisotropic, inhomogeneous diffusion problems. Differential operators are represented by sparse incidence matrices, while weighted mass matrices play the role of metric-dependent Hodge matrices. The resulting mixed formulation is point-wise divergence-free if the right hand side function f = 0. The method is inf-sup stable; no stabilization is required and the method displays optimal convergence on orthogonal and deformed grids.