The energy of chemical reaction is visualized in real space using the electronic energy density n(E)((r) over right arrow) associated with the electron density n(<(<(r)over right arrow>)over right arrow>). The electronic energy density n(E)((r) over right arrow) is decomposed into the kinetic energy density n(T)((r) over right arrow), the external potential energy density n(V)((r) over right arrow), and the interelectron potential energy density n(W)((r) over right arrow). Using the electronic energy density n(E)((r) over right arrow) we can pick up any point in a chemical reaction system and find how the electronic energy E is assigned to the selected point. We can then integrate the electronic energy density n(E)((r) over right arrow) in any region R surrounding the point and find out the regional electronic energy E-R to the global E. The kinetic energy density n(T)((r) over right arrow) is used to identify the intrinsic shape of the reactants, the electronic transition state, and the reaction products along the course of the chemical reaction coordinate. The intrinsic shape is identified with the electronic interface S that discriminates the region R-D of the electronic drop from the region R-A of the electronic atmosphere in the density distribution of the electron gas. If the R spans the whole space, then the integral gives the total E. The regional electronic energy E-R in thermodynamic ensemble is realized in electrochemistry as the intrinsic Volta electric potential phi (R) and the intrinsic Herring-Nichols work function Phi (R). We have picked up first a hydrogen-like atom for which we have analytical exact expressions of the relativistic kinetic energy density n(TM)((r) over right arrow) and its nonrelativistic version n(T)((r) over right arrow). These expressions are valid for any excited bound states as well as the ground state. Second, we have selected the following five reaction systems and show the figures of the n(T)((r) over right arrow) as well as the other energy densities along the intrinsic reaction coordinates: a protonation reaction to He, addition reactions of HF to C2H4 and C2H2, hydrogen abstraction reactions of NH3+ from HF and NH3. Valence electrons possess their unique delocalized drop region remote from those heavily localized drop regions adhered to core electrons. The kinetic energy density n(T)((r) over right arrow) and the tension density <()over right arrow>(S)((r) over right arrow) can vividly demonstrate the formation of the chemical bond. Various basic chemical concepts in these chemical reaction systems have been clearly visualized in real three-dimensional space.