3D thermoelasticity program

The main objective of the paper is to present a program to solve the 3D BVP for the problem of linear thermoelasticity. Numerical algorithms for data structures, element-by-element finite volume approximations, and iterative solution are given.

The purpose of the present program is to solve a mixed boundary value problem (BVP) for a system of equations of the linear thermoelasticity theory in the computational domain composed of parallelepipeds, using finite-volume approximation, element-by-element approach to get global matrix, and fast incomplete factorization solvers.