MINI-Element in 3D

mesh
1

Hi there,

so I’m implementing this finite element which is referred to as the MINI-Element. Using it to make a function space in 2D works fine. However in 3D it won’t work.

from dolfin import *

mesh     = UnitCubeMesh(10, 10, 10)
#mesh    = UnitSquareMesh(10, 10)

Qe      = FiniteElement("Lagrange", mesh.ufl_cell(), 1)
Be      = FiniteElement("Bubble",   mesh.ufl_cell(), 3)
Ve      = VectorElement(NodalEnrichedElement(Qe, Be))
element = MixedElement(Ve,Qe)
W       = FunctionSpace(mesh, element)

The error message is:

“RuntimeError: Bubble element of degree 3 and codimension 0 has no dofs.”

Any idea what to do here?

1 Like
2

See: https://fenicsproject.org/qa/7136/mini-element/

3 Likes