Research Papers

Fluid–Structure Interaction Modeling of Vertical-Axis Wind Turbines

[+] Author and Article Information
Y. Bazilevs

Department of Structural Engineering,
University of California–San Diego,
La Jolla, CA 92093
e-mail: yuri@ucsd.edu

A. Korobenko, X. Deng, J. Yan

Department of Structural Engineering,
University of California–San Diego,
La Jolla, CA 92093

M. Kinzel, J. O. Dabiri

Department of Aerospace Engineering,
Division of Engineering and Applied Science,
California Institute of Technology,
Pasadena, CA 91125

Manuscript received March 25, 2014; final manuscript received April 14, 2014; accepted manuscript posted April 22, 2014; published online May 7, 2014. Assoc. Editor: Kenji Takizawa.

J. Appl. Mech 81(8), 081006 (May 07, 2014) (8 pages) Paper No: JAM-14-1135; doi: 10.1115/1.4027466 History: Received March 25, 2014; Revised April 14, 2014; Accepted April 22, 2014

Full-scale, 3D, time-dependent aerodynamics and fluid–structure interaction (FSI) simulations of a Darrieus-type vertical-axis wind turbine (VAWT) are presented. A structural model of the Windspire VAWT (Windspire energy, http://www.windspireenergy.com/) is developed, which makes use of the recently proposed rotation-free Kirchhoff–Love shell and beam/cable formulations. A moving-domain finite-element-based ALE-VMS (arbitrary Lagrangian–Eulerian-variational-multiscale) formulation is employed for the aerodynamics in combination with the sliding-interface formulation to handle the VAWT mechanical components in relative motion. The sliding-interface formulation is augmented to handle nonstationary cylindrical sliding interfaces, which are needed for the FSI modeling of VAWTs. The computational results presented show good agreement with the field-test data. Additionally, several scenarios are considered to investigate the transient VAWT response and the issues related to self-starting.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Fig. 1

Windspire VAWT structural model with dimensions included: (a) full model using isogeometric NURBS-based rotation-free shells and beams; (b) model cross section 1 showing attachment of the struts to the blades and tower shell; (c) model cross section 2 showing attachment of the struts and tower shell

Grahic Jump Location
Fig. 2

The VAWT aerodynamics computational domain in the reference configuration, including the inner cylindrical region, outer region, and sliding interface that is now allowed to move in space as a rigid object

Grahic Jump Location
Fig. 3

A 2D cross section of the computational mesh along the rotor axis. The view is from the top of the turbine, and the blades are numbered counterclockwise, which is the expected direction of rotation. The sliding interface may be seen along a circular curve where the mesh appears to be nonconforming.

Grahic Jump Location
Fig. 4

A 2D cross section of the blade boundary-layer mesh consisting of triangular prisms

Grahic Jump Location
Fig. 5

Time history of the aerodynamic torque for the pure aerodynamics simulations. (a) 8.0 m/s wind with experimental data from Ref. [25] and (b) 6.0 m/s wind with experimental data from Refs. [26,27].

Grahic Jump Location
Fig. 6

Time history of the rotor speed starting from 0 rad/s

Grahic Jump Location
Fig. 7

Time history of the rotor speed starting from 4 rad/s

Grahic Jump Location
Fig. 8

Time history of the rotor speed starting from 12 rad/s

Grahic Jump Location
Fig. 9

Vorticity isosurfaces at a time instant colored by velocity magnitude for the 4 rad/s case

Grahic Jump Location
Fig. 10

Vorticity isosurfaces of vorticity colored by velocity magnitude for the 4 rad/s case. Zoom on the rotor. From left to right: vorticity at 1.12 s, 1.24 s, 1.40 s, and 1.50 s.

Grahic Jump Location
Fig. 11

Turbine current configuration at two time instances for the 4 rad/s case. The tower centerline in the reference configuration is shown using the dashed line to illustrate the range of turbine motion during the cycle. The range of the tower tip displacement during the cycle is about 0.10–0.12 m.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In