Use of algebraic dual spaces in domain decomposition methods for Darcy flow in 3D domains

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...

A mass-, kinetic energy- and helicity-conserving mimetic dual-field discretization for three-dimensional incompressible Navier-Stokes equations, part I: Periodic domains

We introduce a mimetic dual-field discretization which conserves mass, kinetic energy and helicity for three-dimensional incompressible Navier-Stokes equations. The discretization makes use of a conservative dual-field mixed weak formulation where two evolution equations of velocity are employed and dual representations of the solution are...

Construction and application of algebraic dual polynomial representations for finite element methods on quadrilateral and hexahedral meshes

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...

A high order hybrid mimetic discretization on curvilinear quadrilateral meshes for complex geometries

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...

Discrete equivalence of adjoint neumann–dirichlet div-grad and grad-div equations in curvilinear 3d domains

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...

A conservative hybrid method for darcy flow

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...

Spectral mimetic least-squares method for div-curl systems

In this paper the spectral mimetic least-squares method is applied to a two-dimensional div-curl system. A test problem is solved on orthogonal and curvilinear meshes and both h- and p-convergence results are presented. The resulting solutions will be pointwise divergence-free for these test problems. For N> 1 optimal convergence rates on...