The importance of wind energy has progressed rapidly in the last years. Although Horizontal Axis Wind Turbines (HAWT) are most well-spread, there is an increasing interest in Vertical Axis Wind Turbines (VAWT), especially in the H-Darrieus concept, as these rotors are omni-directional and affordable. However, the physics of these rotors is more complex; they can only be analyzed using transient CFD simulations. Due to the finite aspect ratio of the rotors, a wingtip vortex is created, which generates losses. Optimizing the wingtip geometry could be advantageous for increasing the efficiency of the rotors: this can only be achieved with three-dimensional turbulent transient simulations. For the optimization of winglets, the whole process (mesh generation, CFD computation, post-processing) has to be automated. This is achieved using the OPtimization Algorithm Library++ (OPAL++), a custom C++ code for the description of blended and canted winglets, coupled with a CD-Adapco StarCCM+ JAVA script for the automatization of the mesh generation and CFD computations. To check the viability of the present concept, two parameters have been varied in the simulations. As shown in what follows, an efficient automatic optimization of wind turbine wingtips can be implemented in this manner.