We discuss the self-consistent time-dependent numerical boundary conditions on the basis of theory of characteristics for magnetohydrodynamics (MHD) simulations of solar plasma flows. The importance of using self-consistent boundary conditions is demonstrated by using an example of modeling coronal dynamic structures. This example demonstrates that the self-consistent boundary conditions assure the correctness of the numerical solutions. Otherwise, erroneous numerical solutions will appear.