Thermoelasticity solver

[aron2290]

The code as provided via the link in the yaml file is a thermoelasticity solver where several boundary conditions can be used. These conditions include Dirichlet and Neumann conditions for both the thermal and elastic solvers, and an option for convective cooling. These conditions are applied on one of the six sides of a rectangular prism. The output is given by 2 plots, and a vtk export for post-processing.

[gertjanvanzwieten]

Thanks! I have a few ideas for improving the script but my first question is about line 122, which is where you apply the Dirichlet boundary condition for the thermal problem, but instead of making it a hard constraint you add it as a term to the residual. Is there a reason for this? In the same script on line 162 you do use a hard constraint for the elastic problem.

[aron2290]

I did it like this to make it possible to add multiple constraints for one model. I assumed if you'd put multiple constraints, and add it as term to the residual, they would stack up giving the end result with all constraints combined.