ABSTRACT
Smoothness of the boundary conditions is established that will enable the iterates to asymptotically approximate the parabolic boundary layers of the solution. A domain decomposition technique based on the results of the asymptotic analysis is developed to numerically solve the preceding differential equation when the boundary functions are discontinuous. Convergence results of the numerical scheme are established, and a computational example is given. Divergence is computationally examined.
