Importing A Quadratic Mesh (GeometryFromMesh)
7 views (last 30 days)
Show older comments
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)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/307782/image.png)
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?
0 Comments
Answers (0)
See Also
Categories
Find more on Geometry and Mesh in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!