A rigorous unsteady state and spatially multidimensional model is presented and solved to describe the dynamic behavior of the primary and secondary drying stages of the lyophilization of a pharmaceutical product in vials for different operational policies. the results in this work strongly motivate the aggressive control of freeze drying and it is found that heat input control that runs the process close to the melting and scorch temperature constraints yields (i) faster drying times, and (ii) more uniform distributions of temperature and concentration of bound water at the end of the secondary drying stage. a rigorous unsteady state and spatially multidimensional model is presented and solved to describe the dynamic behavior of the primary and secondary drying stages of the lyophilization of a pharmaceutical product in vials for different operational policies. the results in this work strongly motivate the aggressive control of freeze drying and it is found that heat input control that runs the process close to the melting and scorch temperature constraints yields (i) faster drying times, and (ii) more uniform distributions of temperature and concentration of bound water at the end of the secondary drying stage.