Numerical method for solving joint thermo-diffusive problems in an infinite combined domain with thin resistant interphase.

This work deals with a class of Boundary Value Problems describing joint thermo-diffussive fields in an infinite combined domain, which consists of two subdomains, matched by a thin intermediate layer. The main problem is reduced to an equivalent one given in the bounded subdomain, with non-local boundary condition on the transmission surface. Such a condition incorporates all the information about the infinite subdomain and the intermediate layer. The equivalent problem is solved by means of Finite Element Method in frames of Matlab package. As it is not possible to introduce the non-local boundary conditions along a part of the boundary directly into FEM code, a dedicated iterative subroutine is constructed. Effectiveness of the method has been checked on selected benchmarks. Accuracy and convergence of the procedure have been addressed in the analysis.