To investigate the interplay between particle charging and hydrodynamics in fluidized beds, models for triboelectric charging and electrostatic forces were built into a computational fluid dynamics‐discrete element method model. Charge transfer was governed by the difference in effective work function between contacting materials as well as the electric field at the point of contact. Monodisperse particles were fluidized with an effective work function difference between the particles and the conducting walls. For smaller work function differences, hydrodynamics were not changed significantly as compared with an uncharged case. In these simulations, the average charge saturated at a value much lower than the value anticipated based on the work function difference, and a unimodal distribution of charges was observed. For larger work function differences, particles stuck to walls and bed height oscillations due to slugging were less pronounced. For these cases, a bimodal distribution of charges emerged due to effects from strong electric fields. © 2016 American Institute of Chemical Engineers AIChE J, 63: 1872–1891, 2017