We employed a quantum simulation approach to investigate the magnetic properties of monolayer square nanodisks with Dzyaloshinsky-Moriya (DM) interaction. The computational program converged very quickly, and generated chiral spin structures on the disk planes with good symmetry. When the DM interaction is sufficiently strong, multi-domain structures appears, their sizes or average distance between each pair of domains can be approximately described by a modified grid theory. We further found that the external magnetic field and uniaxial magnetic anisotropy both normal to the disk plane lead to reductions of the total free energy and total energy of the nanosystems, thus are able to stabilize and/or induce the vortical structures, however, the chirality of the vortex is still determined by the sign of the DM interaction parameter. Moreover, the geometric shape of the nanodisk affects the spin configuration on the disk plane as well.