A code is developed for the computation of three-dimensional aeroelastic problems such as wing flutter. The unsteady Navier-Stokes flow solver is based on a finite-volume approach with centered flux discretization and artificial diffusion. For the structural displacements a modal approach is applied. The temporal discretization is implicit for both the flow equations and the structural equations. An explicit dual-time method is used to integrate the coupled governing equations. A multigrid method is applied to advance the flow solution, and the computation is performed in parallel with a multiblock approach. A supercritical 2-D wing and the AGARD 445.6 wing serve as test cases for flutter investigations. Results for inviscid flow are compared with results obtained by solving the Navier-Stokes equations with the Baldwin-Lomax and k-ω turbulence models, respectively. Inclusion of viscous effects is critical for the 2-D wing. LCO of the 2-D wing is predicted, but with larger amplitude compared to experimental measurements. Predicted flutter boundary for the AGARD wing agrees well with experimental data in subsonic and transonic range but deviates significantly from experimental data in the supersonic range. Inclusion of viscous effects only slightly improves the result for this case.