EDIT:
I fixed the problem that involved inhomogeneous constraints (i had not
initialized the bc matrix). However I am still struggling since boundary
conditions are only being interpolated on the external face of my boundary
cells.
Le mercredi 26 juin 2024 à 14:00:46 UTC+3, Tomas Fernandez Bouvier a écrit :
> In the context of a parallel implementation of step-26, where I inspired
> myself in step-32, I am trying to set a different constant in the initial
> and boundary conditions that I apply to my system. Below is an image of the
> result when I set the initial system to 1 instead of 0. I tried to
> investigate the implementation of a rhs. Somehow the conditions appear to
> be inhomogeneous and therefore the rhs is modified. The only improvement
> that I got was when I muted this lines in my code:
>
> if (constraints.is_inhomogeneously_constrained(
> data.local_dof_indices[i]))
> {
> for (unsigned int j = 0; j < dofs_per_cell; ++j)
> data.matrix_for_bc(j, i) +=0 ;
> /*
> (scratch.phi_T[i] * scratch.phi_T[j] *
> (use_bdf2_scheme ? ((2 * time_step + old_time_step)
> /
> (time_step + old_time_step)) :
> 1.) ) *
> scratch.fe_values.JxW(q);
> */
> }
>
> but I am still not able to set the entire system initial temperature to 1
> or any other initial value.
>
> Can you help me?
>
> Thank you!
>
>
>
> [image: Screenshot from 2024-06-26 13-55-07.png]
--
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/efb4046c-8d37-4243-8dba-bc959af2b0dan%40googlegroups.com.
