Several complex processes such as electrochemical reaction, multiphase flow, mass, and heat transfer exist in the fuel cell. It is critical to explore the multi-coupling mechanism inside the fuel cell. Meanwhile, reasonably controlling the auxiliary system to make the fuel cell work under appropriate conditions is significant. Numerical models are powerful tools that can simulate, design, and optimize the fuel cell, thus reducing the cost and time of producing expensive prototypes. Generally, the distributed parameter model based on the partial differential conservation equation can explore the dynamic distribution of multiple physical fields and reveal the multiscale transmission phenomenon in the fuel cell. The lumped parameter model combined with the auxiliary subsystem model can reflect the impact on external dynamic changes of fuel cell performance. This makes it convenient to realize control strategy development and optimization. This chapter includes three parts, which take fuel cell modeling and simulation as the core. The state-of-the-art modeling and simulation concerning proton exchange membrane fuel cells are introduced in detail. Then, the process of a distributed parameter model and related simulation is presented, which aims to reveal the internal coupling law of polarization processes with different time constants and the dynamic response mechanism after a load change. Finally, a lumped parameter stack model is described in conjunction with the subsystem model. By this, different subsystem control strategies are investigated, which provide a guideline for controller design of the existing vehicular fuel cell system.