Weak solutions and simulations to a square phase‐field crystal model with Neumann boundary conditions. (5th January 2022)