The interaction of solar wind discontinuities with reflected solar wind particles upstream of Earth's bow shock leads to the formation of large scale transient phenomena such as Foreshock Bubbles (FBs) and hot flow anomalies. These transient phenomena play an important role in accelerating and energizing plasma and could have global impacts on the near-Earth environment. Direct derivations from the Vlasov-Maxwell equation provide the equations that describe the temporal evolution of the kinetic and thermal energy. In this ongoing study, we investigate the behavior of the fluid energy terms that directly affect the evolution of the kinetic and thermal energy associated with these transients, with a particular focus on FBs. Specifically, we analyze the behavior of these energy terms in different sub-regions of the FB, including its core, sheath, and the shock created by its expansion at different stages in the lifetime of an FB. We employ a 2D global hybrid-Vlasov simulation performed with the Vlasiator model and compare the numerical results with FBs observed by the Magnetospheric MultiScale (MMS) mission. We discuss the role of FBs in accelerating, heating the plasma and producing or annihilating magnetic energy.