This paper presents a second-order continuity
non-overlapping domain decomposition (DD) technique for
numerically solving second-order elliptic problems in
two-dimensional space. The proposed DD technique uses integrated Chebyshev polynomials to represent the solution in subdomains. The constants of integration are utilized to impose continuity of the second-order normal derivative of the solution at the interior points of subdomain interfaces. To also achieve a C2 (C squared) function
at the intersection of interfaces, two additional unknowns are introduced at each intersection point. Numerical results show that the present DD method yields a higher level of accuracy than conventional DD techniques based on differentiated Chebyshev polynomials
