1. A method for reservoir simulation, which comprises: distributing new fluid pressures and compositions determined for the reservoir to one or more GPU's by a respectively linked CPU; performing a multi-phase flash calculation on each GPU that was distributed new fluid pressures and compositions with an unstable phase; and distributing a result for each multi-phase flash calculation to each respectively linked CPU for the reservoir simulation.