In this study, a two-dimensional computational domain featuring gas and solid phases is computationally studied for Geldart-B-type particles. In addition to the baseline case of a uniform gas-phase injection velocity, three different inlet velocity profiles were simulated, and their effects on the fluidized bed hydrodynamics and bubble dynamics have been studied. An in-house computer program was developed to track the bubbles and determine the temporal evolution of their size and position prior to their breakup. This program also provides information on the location of bubble coalescence and breakup. The gas-solid interactions were simulated using a Two-Fluid Model (TFM) with Gidaspow’s drag model. The results reveal that the bed hydrodynamics feature a counter-rotating vortex pair for the solid phase, and bubble dynamics, such as coalescence and breakup, can be correlated with the vortices’ outer periphery and the local gradients in the vorticity.