Multiphysics Contact
Multiphysics contact problems often exhibit a high degree of nonlinearity, which leads to convergence problems for the nonlinear solver. As ab example, consider heat transfer through the contact area, where initially only a small spot is in contact. The solution for the temperature field is then extremely sensitive to the size of the contact area. If at the same time, the solid deforms due to thermal expansion, there may be large changes in the contact area between each iteration,
It is important to resolve the size of the contact area accurately, that is, to use a very fine mesh in the contact area when modeling fully coupled multiphysics problems.
If the contact area is larger, a fine mesh is not required because then the temperature solution is not that sensitive to the size of the contact area. If possible, start with an initial configuration where the contact area is not very small.
You can use the contact variables (gap and contact pressure) in expressions for quantities in other physics interfaces. As an example, a thermal resistance in the contact region can depend on the contact pressure.
In many cases, the penalty method is preferred in multiphysics contact problems because of its better stability and less restrictive requirements on solver selections. If the contact conditions depend strongly on the contact pressure, use the augmented Lagrangian method because if its higher accuracy.