A method to compute sensitivities for use in gradient-based shape optimization applied to unsteady aeroelasticity problems is presented. The method takes into account the coupling between all three fundamental aspects of computational aeroelasticity, namely, unsteady flow equations, time varying structural response to aerodynamic loads, and dynamic meshes that accommodate geometric deformations...