The aerodynamic shape optimization of airfoils using unconstrained trust region methods