3D multi-phase simulation of metal bipolar plate proton exchange membrane fuel cell stack with cooling flow field