Importing A Quadratic Mesh (GeometryFromMesh)

9 views (last 30 days)
ADSW121365
ADSW121365 on 1 Jun 2020
Edited: ADSW121365 on 1 Jun 2020
I'm using the PDE Toolbox and the geometryFromMesh function to import meshes made externally (gmsh). The PDE Toolbox supports both linear and quadratic meshes of tetrahedra. What I'm struggling with is successfully importing an external quadratic mesh.
So I have a geometry defined inside gmsh and can produce meshes of both orders on said geometry. If I create a linear mesh and import it via geometry from mesh as:
Single_coil_3v_2split
model = createpde(1)
nodes = msh.POS'; elements = msh.TETS; groupsID = elements(:,5);
elements = elements(:,1:4); elements =elements';
geometryFromMesh(model,nodes,elements,groupsID);
pdegplot(model,'CellLabels','on','FaceAlpha',0.5)
The geometry is detected correctly. However for a quadratic mesh on the same geometry, I end up with holes & distorted elements:
Single_coil_3v_splitquad
model = createpde(1)
nodes = msh.POS'; elements = msh.TETS; groupsID = elements(:,11);
elements = msh.TETS10; elements = elements(:,1:10); elements =elements';
geometryFromMesh(model,nodes,elements,groupsID);
pdegplot(model,'CellLabels','on','FaceAlpha',0.5)
This seems to be constistently a problem regardless of how I mesh (higher/lower resolutions still have the same fundamental issue). Ideally I'd import the geometry for the linear mesh, then import the quadratic mesh however this isn't supported by the toolbox from what I can tell.
Does anyone have any ideas, solutions or suggestions?

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by