We present a new method for solving the minimization problem in ferromagnetism. Our method is based on replacing the non-local non-convex total energy of magnetization by a new local non-convex energy of divergence-free fields. Such a general method works in all dimensions. However, for the two-dimensional case, since the divergence-free fields are equivalent to the rotated gradients, this new energy can be written as an integral functional of gradients and hence the minimization problem can be solved by some recent non-convex minimization procedures in the calculus of variations. We focus on the two-dimensional case in this paper and leave the three-dimensional situation to future work. Special emphasis is placed on the analysis of the existence/non-existence depending on the applied field and the physical domain.