The bursting and quiescent activity observed in many GRB remains to be fully understood, thus is significant to study the effects that a pulsed central engine has on the prompt GRB emission variability with numerical simulations. In this study, for the first time ever, numerous models of a variable jet drilling through the stellar envelope and then through the ISM is followed using 3D numerical simulations. Since there is a correlation between the variability and the observed peak isotropic luminosity, we calculate the light curve for each model and discuss how the variability in the central engine affects the photospheric emission. Finally, we also deliberate how our results compare with observations from Swift and other instruments.