The objective of this paper is to investigate the hydroelastic response of a flexible NACA66 hydrofoil in cavitating flows by combined experimental and numerical studies. Experimental results are presented for rigid/flexible NACA66 hydrofoils fixed at alpha(0) = 80 degrees for subcavitating (sigma = 8.0) and cavitating flows (sigma = 1.4). The high-speed video camera and Laser Doppler Vibrometer are applied to investigate the flow patterns and vibration characteristics. The multiphase flow is modeled with the incompressible and unsteady Reynolds Averaged Navier-Stokes (URANS) equations. The k - omega SST turbulence model with the turbulence viscosity correction and the Zwart cavitation model are introduced to the present simulations. The results showed that the cavitation has significant effect on the foil deformation and the unsteady characteristics of the hydroelastic response. The bending deformation is enhanced when the cavitation occurred. Meanwhile, the hydroelastic response has also affected the cavitation development and the vortex structure interactions. The cavity shedding frequency and vortex shedding and interacting frequency for the flexible hydrofoil are higher than that for the rigid hydrofoil. Compared to the periodic development of the hydrodynamic coefficients for the rigid hydrofoil, the hydrodynamic load coefficients of the flexible hydrofoil fluctuate more significantly, and the chaotic response of the flexible hydrofoil is mainly attributed to the disturbance caused by the flow-induced flutter and deformation of the foil. The evolution of the transient cavity shape and the corresponding hydrodynamic response can be divided into three stages: During the development of the attached cavity, the partial sheet cavity is formed and develops with the lift and drag coefficients increasing, while the maximum attached cavity formed on the suction side of the flexible hydrofoil is larger than that of the rigid hydrofoil, which is caused by the increase of the effective angle of attack due to the twist deformation. During the vortex structure interaction and cavity shedding process, the hydrodynamic loads for the flexible hydrofoil fluctuate because of the foil deformation, leading to a more complex cavitation pattern. During the residual cavity shedding and partial sheet cavity formation process, the cavities, together with the counter-rotational vortex structures, shed downstream totally and are followed by the formation of partial sheet cavities in next period, which is in advance for the flexible hydrofoil due to the larger effective angle of attack. (C) 2015 Elsevier Ltd. All rights reserved.