With the gradual upscaling of offshore wind turbines, the blade size becomes large and the blade tip velocity reaches almost a Mach number of 0.3, the limit of the incompressible flow assumption.To simulate wind turbine flows, an incompressible solver was often a good choice as it gave good predictions and the Mach number of the wind turbine flows were small. When the Mach number is increased to be close to 0.3, the validity of incompressible flow may be questionable and thus needs to be checked. In this paper, an accurate wind turbine flow solution method is proposed to study the compressibility effects. First, the National Renewable Energy Laboratory (NREL) Phase VI wind turbine is used as an example to verify the method at a small Mach number. Through a comparative analysis between the compressible, the incompressible numerical results, and experimental data, the proposed method is found with a good simulation accuracy. Next, the flow fields of the large wind turbine IEA-15 MW, under different wind conditions are simulated. The local Mach number, density and pressure coefficient at different blade cross-sections of the compressible results are analyzed and compared to its incompressible ones, and the influences of compressibility on the thrust, torque and power of the wind turbine are quantified. The results show that at the rated wind speed, the thrust of the compressible flow simulation is 1.40 % higher than that of the incompressible one, and the torque of the compressible flow simulation is nearly 11.07 % higher than that of the incompressible one. As a result, the compressibility effects are recommended to be added in future studies of large wind turbines.