Hi Behrooz, To complement Simon’s suggestion, I can think of two more ways to achieve what you’re wanting but will affect the mesh itself. Whether or the second suggestion is a reasonable approach for you depends on the complexity of the geometry that you’re creating. 1. Use an external mesh generator. Generating disconnected regions in the domain is normally a near trivial thing to do with these tools. 2. You can use the deal.II mesh tools to create a mesh that has a obvious physical disconnect between the surfaces of interest. Maybe you’d do this by creating the coarse triangulation by hand, or two triangulations that you "merge together <https://www.dealii.org/current/doxygen/deal.II/namespaceGridGenerator.html#a7cd88e7eacd46697dee80ad2b8438d54>" (watch for the merge tolerance, which defines how close vertices must be to one another to be considered coincident). You can then use the mesh movement functions (see step-49 <https://www.dealii.org/current/doxygen/deal.II/step_49.html#ModifyingaMesh>) to move the grid such that you make the surfaces that you want coincident. When you move the mesh, the cell connectivity remains as is and so you’ll end up with two boundaries at this interface with no connection between the two.
Best, Jean-Paul > On 02 Oct 2020, at 14:30, simon...@gmail.com <simonsti...@gmail.com> wrote: > > Hi, > if you have a single triangulation, you can put the solution to the "left" > and "right" of the interface into different components of a vector valued > element. For example: > > FESystem<dim> left_element(FE_Q<dim>(1), FE_Nothing<dim>()); > FESystem<dim> right_element(FE_Nothing<dim>(), FE_Q<dim>(1)); > > hp::FECollection<dim> fe_collection(left_element, right_element); > > Then you can select which cells should use which element as in step-46: > www.dealii.org/current/doxygen/deal.II/step_46.html > <http://www.dealii.org/current/doxygen/deal.II/step_46.html> > > Best, > Simon > > > On Friday, October 2, 2020 at 11:51:50 AM UTC+2 Behrooz Karami wrote: > Dear All, > > Is it possible to duplicate nodes/vertices in deal.ii? > I want to duplicate nodes at the interface of two domains so that the could > get fully separated. > Any ideas are appreciated. > > Best, > Behrooz > > > > > > -- > The deal.II project is located at http://www.dealii.org/ > <http://www.dealii.org/> > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > <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 > <mailto:dealii+unsubscr...@googlegroups.com>. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/7affe3d6-fd0e-4986-96a6-ee5d67420072n%40googlegroups.com > > <https://groups.google.com/d/msgid/dealii/7affe3d6-fd0e-4986-96a6-ee5d67420072n%40googlegroups.com?utm_medium=email&utm_source=footer>. -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/C48AECA2-B1E8-4FE8-A131-48431CBC575D%40gmail.com.
