We emulate rough magnetic flux fronts and induction contours in a long hard type-II superconducting cylinder, using the elliptic critical-state model. Specifically, using a current density with spatial periodicity we perform numerical calculations to obtain magnetic induction distributions of the sample in screened and trapped flux regimes. We found that tuning the parameters of our model, it is possible to achieve different roughening levels at the flux front as well as at the flux profiles behind it, complex current circulation paths, the formation of Meissner state islands or cavities, and the so-called remagnetization fronts.