SummaryThe ongoing trend in wind turbine development is towards larger rotors because of the resulting lower cost of energy. These large rotors lead to relatively flexible structures that are more susceptible to the unsteady aerodynamic loading occurring in normal operating conditions. An accurate prediction of these loadings is important for the design of an economically viable and technically reliable wind turbine. Simulation methods of wind turbine aerodynamics currently in use mainly fall into two categories: the first is the group of traditional low-fidelity engineering models and the second is the group of computationally expensive CFD methods based on the Navier-Stokes equations.For an engineering environment the search is for "medium fidelity" wind turbine simulation methods that bridge the gap between the computationally inexpensive low-fidelity methods and the computationally expensive CFD methods. The ultimate goal is a balanced mixture of higher accuracy of the representation of the physics and shorter simulation times for wind turbine aerodynamics simulation methods. This can be found in the combination of the theories for panel methods, integral boundary layer methods, strong viscous-inviscid coupling, and fluid structure interaction.The present study focuses on the development of the theory and the practical implementation of a fast multilevel integral transform in a computer program. We utilize this multilevel scheme in a low-order panel method. It is demonstrated that for the simulation of the wake flow of wind turbine rotors the computational burden is reduced from O(N 2 ) for a conventional panel method to O(N ) for the present method. This implies that the computational effort is reduced to grow linearly with problem size N , with N the number of panels.We consider the unsteady, incompressible flow around wind turbine rotors and assume the effects of viscosity to be confined to infinitesimal thin boundary layers and wake regions and assume irrotational flow elsewhere. These assumptions allows us to reduce the flow problem to the problem of solving the Laplace equation for a (scalar) velocity potential function. This makes it possible to reformulate the problem as an integral equation over the surface of the rotor and the wakes that emanate from the trailing edges of the rotating wind turbine blades.The mathematical model is discretized in the form of a low-order panel iii Multilevel Panel Method method. The implementation of the panel method is verified by considering the flow over a stationary ellipsoid in a uniform onset velocity field, the flow over a rotating ellipsoid in a fluid at rest, and by considering the flow over a high aspect ratio wing with elliptic planform with as cross-section a von Kármán-Trefftz airfoil. For the first two test cases the numerical results are compared with analytical solutions. The error in the velocity potential is shown to be O(h 2 ), with h a characteristic panel size. The third test case uses the analytical solution for the 2D von Kármán-Trefftz airfoil ...