On 8/24/22 12:45, Lucas Myers wrote:
As a first guess I suspect I can just mark each of the subdomains with a
material_id or a user_index or something, and then only iterate over the
cells from a particular subdomain during the assembly and solve
normally. However, I am unsure how to mark the appropriate faces on the
outside of the subdomains as boundary faces.
Has anyone dealt with this kind of problem before? Does it seem like
this is the right direction to go? And how can I mark faces as belonging
to a boundary, even if the cells are internal to the domain?
The VectorTools::interpolate_boundary_values() function can't deal with
this, but it isn't very difficult to write your own replacement function
that simply loops over the faces of the cells that are part of a
subdomain, rather than the faces of the entire triangulation.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
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].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/adf75b27-0168-3b4e-5afe-8834fae298f7%40colostate.edu.
