Accreting black holes are observed in a large variety of systems in astronomy: active galactic nuclei, X-ray binaries, tidal disruption events, gamma-ray bursts. While analytical one-dimensional accretion models have been enormously useful for understanding the physics of accretion flows, many effects, e.g., the formation of jets and winds, are beyond the scope of these models. Numerical general relativistic MHD simulations are able to include more physics than analytical models and can study gas dynamics in two and three dimensions. The talk will review current progress in this field.