Conservative domain decomposition procedure with unconditional stability and second-order accuracy (Q984378)

This article presents a conservative domain decomposition procedure for solving one-dimensional parabolic initial-boundary value problems. In this procedure, the fluxes at subdomain interfaces are calculated from the solution at the previous time level, then these fluxes serve as Neumann boundary data for an implicit, block-centered discretization in the subdomain. The unconditional stability and the second-order accuracy of the solution and the fluxes are proved. Numerical results examining the stability, accuracy and parallelism of the procedure are presented.