Re: [deal.II] Impose values inside a material/ not on the boundary

2022-02-17 Thread Wolfgang Bangerth
On 2/15/22 07:32, 'Markus Mehnert' via deal.II User Group wrote: Thank you for your response. I identified the correct dofs by trial and error which worked (9 dofs for a scalar quantity with quadratic polyorder). However, when I use /"fe_face.system_to_base_index().first.first" /to identify,

Re: [deal.II] Impose values inside a material/ not on the boundary

2022-02-15 Thread 'Markus Mehnert' via deal.II User Group
Dear Wolfgang, Thank you for your response. I identified the correct dofs by trial and error which worked (9 dofs for a scalar quantity with quadratic polyorder). However, when I use *"fe_face.system_to_base_index().first.first" *to identify, which base element these dofs belong to, it does no

Re: [deal.II] Impose values inside a material/ not on the boundary

2022-02-14 Thread Wolfgang Bangerth
On 2/11/22 12:12, 'Markus Mehnert' via deal.II User Group wrote: *** Caution: EXTERNAL Sender *** Dear Wolfgang, Thank you for your quick response. The dofs per face are defined as / / / / /const unsigned int dofs_per_face = fe_cell.dofs_per_face;/ / / where fe_cell is the/FESystem/ that consis

Re: [deal.II] Impose values inside a material/ not on the boundary

2022-02-11 Thread Wolfgang Bangerth
Markus, how do you define dofs_per_face? Do you adjust it for the polynomial degree? Best W. On 2/11/22 10:51, 'Markus Mehnert' via deal.II User Group wrote: *** Caution: EXTERNAL Sender *** Hi everyone, Instead of prescribing classical Dirichlet boundary conditions (i.e. the values for t

[deal.II] Impose values inside a material/ not on the boundary

2022-02-11 Thread 'Markus Mehnert' via deal.II User Group
Hi everyone, Instead of prescribing classical Dirichlet boundary conditions (i.e. the values for the solution field on the boundary) I am trying to impose the values for a specific degree of freedom inside my material. My idea was to iterate over all cells and their faces. Then I use the locati