We report on an investigation of a method that applies simultaneously two different mathematical models in order to optimize the design of a Zeeman Slower towards the implementation of ultra cold atoms in solid state physics. We introduce the implementation of a finite element simulation that allows us to predict with great accuracy the magnetic field intensity profile generated by the proposed design. Through the prediction of the behavior of the Zeeman Slower a greater control is acquired, which allows the optimization of the different experimental variables. We applied the method in the design of a multilayer solenoidal Spin-Flip Zeeman Slower for strontium atoms. The magnetic intensity profile generated by the Zeeman Slower is in agreement with the magnetic field strength profile necessary for the atom cooling and tends to zero in both end sides. The latter terms are essential in order to optimize the amount of trapped and cooled atoms.