We have quantitatively investigated the radiation belt's dynamic variations of 1.5-6.0 MeV electrons during 54 CME (coronal mass ejection)-driven storms from 1993 to 2003 and 26 CIR (corotating interaction region)-driven recurrent storms in 1995 by utilizing case and statistical studies based on the data from the SAMPEX satellite. It is found that the boundaries determined by fitting an exponential to the flux as a function of L shell obtained in this study agree with the observed outer and inner boundaries of the outer radiation belt. Furthermore, we have constructed the Radiation Belt Content (RBC) index by integrating the number density of electrons between those inner and outer boundaries. According to the ratio of the maximum RBC index during the recovery phase to the pre-storm average RBC index, we conclude that CME-driven storms produce more relativistic electrons than CIR-driven storms in the entire outer radiation belt, although the relativistic electron fluxes during CIR-related storms are much higher than those during CME-related storms at geosynchronous orbit. The physical radiation belt model STEERB is based on the three-dimensional Fokker-Planck equation and includes the physical processes of local wave-particle interactions, radial diffusion, and adiabatic transport. Due to the limitation of numerical schemes, formal radiation belt models do not include the cross diffusion term of local wave-particle interactions. The numerical experiments of STEERB have shown that the energetic electron fluxes can be overestimated by a factor of 5 or even several orders (depending on the pitch angle) if the cross diffusion term is ignored. This implies that the cross diffusion term is indispensable for the evaluation of radiation belt electron fluxes. Formal radiation belt models often adopt dipole magnetic field; the time varying Hilmer-Voigt geomagnetic field was adopted by the STEERB model, which self-consistently included the adiabatic transport process. The test simulations clearly indicate that the adiabatic