How to Use Isoparametric Elements with P2-P1 in FEniCSx for DFG 2D-2 Benchmark?

It does, as i have already mentioned the mesh is higher order, made possible in the line:

Generating the mesh#

We are now ready to generate the mesh. > However, we have to decide if our mesh should consist of triangles or quadrilaterals. In this demo, to match the DFG 2D-3 benchmark, we use second order quadrilateral elements.

    gmsh.model.mesh.generate(gdim)
    gmsh.model.mesh.setOrder(2)
    gmsh.model.mesh.optimize("Netgen")
1 Like