Hello everyone,
I am modifying the Step-44 (Quasi-Static Finite-Strain Compressible
Elasticity) tutorial in deal.II to incorporate thermo-mechanical equations.
However, I am facing an issue with applying temperature Dirichlet boundary
conditions.
Problem Description:
I am applying the temperature boundary condition using:
VectorTools::interpolate_boundary_values(
dof_handler_ref,
boundary_id,
Functions::ConstantFunction<dim>(current_applied_temperature_0, dim +
1),
constraints,
fe.component_mask(Temperature)
);
After that, before solving the system I apply:
constraints.condense(tangent_matrix, system_rhs);
Despite this, my solution converges, but the computed nodal temperatures
are negative and completely incorrect.
However, when I apply boundary conditions using:
MatrixTools::apply_boundary_values(boundary_values,tangent_matrix,solution_n,system_rhs);
I obtain the correct nodal temperatures.
Could someone clarify what am I doing wrong in my first approach? Thanks
for your support
--
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.
To view this discussion visit
https://groups.google.com/d/msgid/dealii/fbc849e3-1f25-4f8b-8e4a-772a870d7eabn%40googlegroups.com.
