We integrate the MHD ideal equations of a slender flux tube to simulate the internal plasma dynamics of coronal post-flare loops. We study the onset and evolution of the internal plasma instability to compare with observations and to gain insight into physical processes and characteristic parameters associated with flaring events. The numerical approach uses a finite-volume Harten-Yee TVD scheme to integrate the 1D1/2 MHD equations specially designed to capture supersonic flow discontinuities. We could reproduce the observational sliding down and upwardly propagating of brightening features along magnetic threads of an event occurred on October 1st, 2001. We show that high--speed downflow perturbations, usually interpreted as slow magnetoacoustic waves, could be better interpreted as slow magnetoacoustic shock waves. This result was obtained considering adiabaticity in the energy balance equation. However, a time--dependent forcing from the basis is needed to reproduce the reiteration of the event which resembles observational patterns -commonly known as quasi--periodic pulsations (QPPs)- which are related with large scale characteristic longitudes of coherence. This result reinforces the interpretation that the QPPs are a response to the pulsational flaring activity.