The objective of this work is to determine the optimal shape, gains and duration of an electric pulse applied to a Peltier cell, together with the length of the thermoelectric to maximize cooling while minimizing electric consumption. For this purpose, a fully coupled, multiphysics, dynamic finite-element model, which solves for the thermal, electric and mechanical fields is used. Because of the demanding computing requirements of the optimization process, a special mesh is designed and a convergence analysis is carried out before using the multiphysics model. The highly nonlinear optimization is done by simulated annealing, a heuristic algorithm in the Markov chain Monte-Carlo family. A preliminary parametric investigation is presented, analyzing the impact of some of the parameters. The results of this preliminary analysis help to understand the effect of the different shapes in the evolution of the cold face temperature. Some of these results are expected and have already been discussed elsewhere, but others can only be explained after further analysis and a full system modeling. Pulse optimization is multiobjective and multiparametric, i.e., it can consider several targets such as maximizing the cooling temperature, the cooling duration or others. The trade-offs between the different targets are studied. In all cases, stresses inside the thermoelement are examined at all points, and the pulses must meet the restriction that an equivalent stress is not above the allowable value.