A new approach for sensitivity analysis of chaotic dynamical systems is presented. Classical sensitivity analysis methods, used to understand how design and control variables of a given system affect an output of interest, bring to a cumulative error growth if they are applied to a chaotic system with trajectories which are long relative to the predictability time scale of the problem. This occurs because a finite amplitude perturbation in the initial condition leads the finite-time statistical quantities of the observable of interest to diverge exponentially from the reference trajectory. The purpose of this project is to validate a new sensitivity analysis algorithm to obtain a solution that remains bounded over the time interval considered.