Kyusik, [...] > Additional Information: > The input data for creating a triangulation contained information about a > line with indices 272 and 100that is supposed to have boundary indicator . > However, this is an internal line not located on the boundary. You cannot > assign a boundary indicator to it. > > If this happened at a place where you call > Triangulation::create_triangulation() yourself, you need to check the > SubCellData object you pass to this function. > > If this happened in a place where you are reading a mesh from a file, then > you need to investigate why such a line ended up in the input file. A > typical case is a geometry that consisted of multiple parts and for which > the mesh generator program assumes that the interface between two parts is > a boundary when that isn't supposed to be the case, or where the mesh > generator simply assigns 'geometry indicators' to lines at the perimeter of > a part that are not supposed to be interpreted as 'boundary indicators'. > -------------------------------------------------------- > > Aborting! > ---------------------------------------------------- > > But, there is no problem when I only add the physical line on the outer > one.(circle) > > Could you let me know what the problem is? > In my opinion the error message is quite clear, isn't it? Assigning boundary indicators is only allowed for faces on the actual boundary while adding a physical line on the inner boundary in gmsh results in the attempt to set boundary indicators on faces in the interior. If you really want to treat the interface between inner and outer mesh as a boundary, have a look at the implementation of GridGenerator::merge_triangulations [1]. Essentially, you would do the same but neglect the call to GridTools::delete_duplicated_vertices. This way you end up with a triangulation in which inner and outer part don't share any vertex. Is this what you want?
Best, Daniel [1] https://www.dealii.org/8.4.0/doxygen/deal.II/grid__generator_8cc_source.html#l03559 -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.