On 11/27/18 6:56 PM, Oleg Kmechak wrote:
> 
>     But you can't set boundary indicators on faces
>     that are not actually at the boundary.
> 
> 
> You mean that I can't set Dirichlet Conditions on such structure as the tree 
> above?
> Cos this feature is critical for my work

What I'm saying is that you can't set boundary indicators on faces between two 
cells (i.e., on *interior* faces). That's because these are not actually on 
the *boundary*.

If you want to set conditions on degrees of freedom in the interior of the 
domain, you need a different approach. The easiest way is to make sure that 
the faces where you want to set boundary conditions are really at the boundary 
-- i.e., these faces have exactly one adjacent cell; but you can have two 
faces each bounding one cell at the physically same location: these are then 
boundary faces.

Other approaches are to copy the logic of interpolate_boundary_values() and 
make it work on interior faces that you mark through some other approach.

Best
  W.

-- 
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 bange...@colostate.edu
                            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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply via email to