With Monte Carlo simulations, we investigate the relaxation dynamics with a domain wall for magnetic systems at the critical temperature. The dynamic scaling behavior is carefully analyzed, and a dynamic roughening process is observed. For comparison, similar analysis is applied to the relaxation dynamics with a free or disordered surface.In the past years, much effort of physicists has been devoted to the understanding of nonequilibrium dynamic processes. Phase ordering dynamics, spin glass dynamics, structural glass dynamics and interface growth etc are important examples. Since the pioneer work by Janssen et al [1], the universal dynamic scaling form in critical dynamics has been explored up to the macroscopic short-time regime [1-10], when the system is still far from equilibrium.Although the spatial correlation length is still short in the beginning of the time evolution, the short-time dynamic scaling form is induced by the divergent correlating time around a continuous phase transition. Based on the short-time dynamic scaling form, new methods for the determination of both dynamic and static critical exponents as well as the critical temperature have been developed [4,5,[7][8][9][10][11]. Since the measurements are carried out in the short-time regime, one does not suffer from critical slowing down.In understanding the dynamic scaling form far from equilibrium, we should keep in mind that it holds after a time scale t mic , which is sufficiently long in the microscopic sense, but still short in the macroscopic sense. More importantly, the macroscopic initial condition should be taken into account in the dynamic scaling form [1,9,12]. For the dynamic relaxation starting from an ordered state, i.e., a state with an initial magnetization m 0 = 1, for example, the magnetization decays by a power law [8,9,12]. If m 0 is smaller but close to 1, there emerge corrections to scaling. For the dynamic relaxation starting from a random state, i.e., a state without spatial correlations and with a small m 0 , however, the magnetization does not decay, and rather shows an initial increase in the macroscopic short-time regime. An independent critical exponent x 0 must be introduced to describe the scaling dimension of the initial magnetization [1,6,9,10]. If m 0 = 0, the magnetization naturally remains zero during the dynamic evolution, but x 0 is still needed to describe the auto-correlation function etc. This critical exponent also explains the power-law decay of the remanent magnetization in spin glasses [2,13,14]. On the other hand, the short-time dynamic scaling form is universal, in the sense that it does not depend on the microscopic details of the dynamic system, such as the lattice types, interactions, and updating schemes etc. Up to now, the dynamic relaxation with the ordered and random initial states has been systematically investigated.Recent progress in the non-equilibrium critical dynamics and its applications includes, for