Force method conception using transfer matrix to apply to multiphase flow by one-by-one corresponding Particle-Cartesian cell (P/CC) model

. This work presents a novel force method for multiphase flow in which displacements and forces progress while strains terminate on surfaces between two different phases, such as the surface between air and liquid. This report is part of a research project to apply Helmholtz decomposition (H-d) to the finite element method. State vectors continues to exist on the surface of multiphase. The conventional scheme uses displacement method and applies a model to set Dirichlet boundary conditions for air on the air-liquid interface. The proposed scheme allows simulating it without such modeling.