interface boundary condition of subdomains in parallel computing
 

I am working on non-overlapping parallel simulation of 2D black oil with matching grid using IMPES method (it is like solving Laplace equation). I have read many articles about domain decomposition. I know that pressure and mass flux must be equal on the interfaces of subdomains but I don’t know how to apply these boundary conditions between two subdomains. Can anyone help me?