Raviart-Thomas for quadrilaterals

Gigio · May 12, 2021, 8:01pm

Is there any availability estimate for using Raviart-Thomas elements for quadrilaterals? Is there some available alternative way?

dokken · May 12, 2021, 8:52pm

I guess you are referring to either:
https://defelement.com/elements/qcurl.html
or
https://defelement.com/elements/qdiv.html

Both of these elements are implemented and tested in dolfin-x:

github.com

FEniCS/dolfinx/blob/bc894b5ab586617caac779a4eba59193a939b7f1/python/test/unit/fem/test_fem_pipeline.py#L358-L363


@pytest.mark.parametrize("family", ["NCE", "NCF"])
@pytest.mark.parametrize("degree", [1, 2, 3])
def test_NC_hex(family, degree, cell_type, datadir):
    mesh = get_mesh(cell_type, datadir)
    V = FunctionSpace(mesh, (family, degree))
    run_vector_test(mesh, V, degree - 1)

github.com

FEniCS/dolfinx/blob/bc894b5ab586617caac779a4eba59193a939b7f1/python/test/unit/function/test_interpolation.py#L208-L212


def test_NCE_interpolation(cell_type, order):
    random.seed(8)
    mesh = one_cell_mesh(cell_type)
    V = FunctionSpace(mesh, ("NCE", order))
    run_vector_test(V, order - 1)

Topic		Replies	Views
Raviart Thomas glyph visualisation in paraview I/O	6	134	January 17, 2024
Quad Mesh Fails during evaluation on coordinates	6	1107	May 1, 2020
Custom quadrature rule for Quadrature element	1	672	September 17, 2021
Interpolating Test Function between Mismatched Meshes variational formulation	1	655	July 21, 2020
Approximation by piecwise linear functions, Accuracy Table. Triangle General	2	187	July 10, 2023

Raviart-Thomas for quadrilaterals

Related Topics