Mathematical solutions to two way coupling multi-phase flows through numerical integration and PE/FD techniques