Javier, On Thursday, July 28, 2016 at 9:34:15 AM UTC-4, Jose Javier Muñoz Criollo wrote: > > > I have a problem that involves the consideration of a mesh that varies in > size over time. In particular it requires the creation and removal of a > layer of elements on top of a cubic domain. This is, at time t0 the extra > layer doesn't exist and only the cubic domain is considered. At time t1, > the layer is created and remains over the cubic domain until time t2 when > is removed. The information contained in the layer's degrees of freedom can > be safely destroyed during removal. > > From step 49 I can see that it is possible to merge two meshes, that might > solve the problem of creation of the layer. But, how could I erase the > layer after this? > if your code is serial, by far the simplest method is to use hp::DoFHandler with FE_Nothing (see http://dealii.org/developer/doxygen/deal.II/step_46.html) and to activate the elements only between t1 and t2. So you can work with only one mesh and one dof handler and there is no projection. If you need MPI, you can remove cells from a mesh using create_triangulation_with_removed_cells (http://dealii.org/developer/doxygen/deal.II/namespaceGridGenerator.html#ada140ece81bf38a23e738e0e57f89e97) or just save the triangulation before merging.
> > Also, I would like to proyect the information contained in the old mesh > correspoding to the cubic domain to the new mesh while initializing the > extra layer with a user defined value. Is this possible? > > That's the hard part. You will need to map the cells from the old mesh to the new mesh and then do the projection yourself (you will probably need set_dof_values and get_dof_values to do it http://dealii.org/developer/doxygen/deal.II/classDoFCellAccessor.html#a72c2a8569468c2b249c65d0b87e591a4). Best, Bruno -- 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
