[Libmesh-users] Question about hybrid mesh

[Viviana P. B. <vpa...@uc...>]

Hi,

I am trying to assemble the system's matrix for a mesh composed of TET10
(bulk) bounded by TRI6 (surface) elements but I get an error related to an
unsupported type of elements. The problem solves for a single variable.

I create my meshes using CUBIT and export it in in exodus format: 1. Create
the volume and mesh. 2. Create 2 blocks, one for the volume and one for the
surface. 3. Give an element type to each block.

In my program, libmesh recognizes both subdomains with the ids I gave in
CUBIT. When I try to assemble the matrix, I get the error:
*"ERROR: Unsupported type: 4."* I don't understand if I should create
separate quadrature rules because there are two types of elements. But then
the elements are attached to the same variable, so how would I be able to
distinguish between them?

Thanks in advance. Best,

Viviana.