Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Issues with meshing domain with many lines

Please login with a confirmed email address before reporting spam

I'm trying to solve a problem where the domain is a unit cube with lines generated via a periodic voronoi tesselation. I require these lines for my physics. I am very often running into issues when I try to mesh this domain. Some tesselations succeed and some do not. Often the error is either "failed to insert point", "intersecting edges" or similar. Sometimes when the mesh does succeed the solution will fail due to inf and NaN in the stiffness matrix. Since this solution error only appears some times I'm assuming it is also caused by the mesh.

How should I buid my mesh procedure to effectively mesh my domain? I've tried meshing along the edges and boundaries first but this does not help any. I've tried collapsing the short edges and removing details but these do not improve the success rate either.

I have attached an example of what the tesselation may look like. Due to the sometimes very tight intersections between lines I do not expect every tesselation to be possible, but currently I have around a 10% success rate which I too low.


1 Reply Last Post Feb 18, 2020, 2:00 a.m. EST
Bjorn Bretz COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Feb 18, 2020, 2:00 a.m. EST
Updated: 4 years ago Feb 18, 2020, 2:41 a.m. EST

Hi Daniel,

In general, it is recommended that you mesh the entire geometry with one operation if possible. For this model this means that you should only use one Free Tetrahedral node, perhaps with Size nodes defining appropriate sizing on the lines inside the domain, instead of first using Edge and/or Free Triangular nodes, creating mesh on edges and faces, followed by a Free Tetrahedral node to mesh the interior (domain). The motivation for only using one operation can be found in the following blog.

https://www.comsol.com/blogs/best-practices-for-meshing-domains-with-different-size-settings/

If you need further assistance I suggest you contact COMSOL support.

https://www.comsol.com/support

Regards,
Bjorn Bretz
COMSOL Development

Hi Daniel, In general, it is recommended that you mesh the entire geometry with one operation if possible. For this model this means that you should only use one Free Tetrahedral node, perhaps with Size nodes defining appropriate sizing on the lines inside the domain, instead of first using Edge and/or Free Triangular nodes, creating mesh on edges and faces, followed by a Free Tetrahedral node to mesh the interior (domain). The motivation for only using one operation can be found in the following blog. https://www.comsol.com/blogs/best-practices-for-meshing-domains-with-different-size-settings/ If you need further assistance I suggest you contact COMSOL support. https://www.comsol.com/support Regards, Bjorn Bretz COMSOL Development

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.