Hamed,
Am Dienstag, 15. November 2016 19:09:18 UTC+1 schrieb Hamed Babaei:
>
> Hi all,
>
> Following step-45, I am going to implement periodic boundary condition on
> my code for a Thermo-elastic problem, in which the thermal and elastic
> parts are solved separately (I have two different dof_handelers for
> temprature and displacement fields without any block structure). I am
> wondering how to apply periodic constrains, made in setup_system, in the
> thermal field assemble_system for which I don't use
> "constraints.distribute_local_to_global ()" function. Thermal
> Assemble_system is the following:
>
[...]
>
> for (unsigned int i=0; i<dofs_per_cell_eta; ++i)
> {
> for (unsigned int j=0; j<dofs_per_cell_eta; ++j)
> {
> mass_matrix. add (local_dof_indices[i],
> local_dof_indices[j],
> cell_mass_matrix(i,j));
>
> laplace_matrix.add (local_dof_indices[i],
> local_dof_indices[j],
> cell_laplace_matrix(i,j));
>
> nl_matrix. add (local_dof_indices[i],
> local_dof_indices[j],
> cell_nl_matrix(i,j));
>
> }
> nl_term (local_dof_indices[i]) += cell_nl_term(i);
> }
> I guessed perhaps I can replace the red part with the following,
>
> constraints.distribute_local_to_global (cell_mass_matrix,
>
> local_dof_indices,
> mass_matrix);
> constraints.distribute_local_to_global (cell_laplace_matrix,
>
> local_dof_indices,
> laplace_matrix);
> constraints.distribute_local_to_global (cell_nl_matrix,
>
> local_dof_indices,
> nl_matrix);
> constraints.distribute_local_to_global (cell_nl_term,
>
> local_dof_indices,
> nl_term);
>
Can you confirm that these approaches are the same if no constraints are
specified? This should definitely hold.
>
> but the newton_rophson solver for thermal part diverged when I did so!!
>
It is difficult to help you with just this information. How do the
constraints you are using actually look like?
Are they correct?
Best,
Daniel
--
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.
