skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Investigation of the NREL phase VI experiment with the incompressible CFD solver THETA

ORCiD logo [1];  [2];  [2]
  1. Deutsches Zentrum fur Luft- und Raumfahrt, Institute of Aerodynamics and Flow Technology, Braunschweig Germany
  2. Deutsches Zentrum fur Luft- und Raumfahrt, Institute of Aerodynamics and Flow Technology, Göttingen Germany
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Wind Energy
Additional Journal Information:
Journal Volume: 20; Journal Issue: 9; Related Information: CHORUS Timestamp: 2017-10-20 18:08:33; Journal ID: ISSN 1095-4244
Country of Publication:
United Kingdom

Citation Formats

Länger-Möller, A., Löwe, J., and Kessler, R. Investigation of the NREL phase VI experiment with the incompressible CFD solver THETA. United Kingdom: N. p., 2017. Web. doi:10.1002/we.2107.
Länger-Möller, A., Löwe, J., & Kessler, R. Investigation of the NREL phase VI experiment with the incompressible CFD solver THETA. United Kingdom. doi:10.1002/we.2107.
Länger-Möller, A., Löwe, J., and Kessler, R. 2017. "Investigation of the NREL phase VI experiment with the incompressible CFD solver THETA". United Kingdom. doi:10.1002/we.2107.
title = {Investigation of the NREL phase VI experiment with the incompressible CFD solver THETA},
author = {Länger-Möller, A. and Löwe, J. and Kessler, R.},
abstractNote = {},
doi = {10.1002/we.2107},
journal = {Wind Energy},
number = 9,
volume = 20,
place = {United Kingdom},
year = 2017,
month = 5

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on May 15, 2018
Publisher's Accepted Manuscript

Save / Share:
  • A multigrid solver for the steady incompressible Navier-Stokes equations on a curvilinear grid is constructed. The Cartesian velocity components are used in the discretization of the momentum equations. A staggered, geometrically symmetric distribution of velocity components is adopted which eliminates spurious pressure oscillations and facilitates the transformation between Cartesian and co-or contra-variant velocity components. The SCGS (symmetrical collective Gauss-Seidel) relaxation scheme proposed by Vanka on a Cartesian grid is extended to this case to serve as the smoothing procedure of the multigrid solver, in both [open quotes]box[close quotes] and [open quotes]box-line[close quotes] versions. Due to the symmetric distribution of velocitymore » components of this scheme, the convergence rate and numerical accuracy are not affected by grid orientation, in contrast to a scheme proposed in the literature in which difficulties arise when the grid lines turn 90 from the Cartesian coordinates. Some preliminary numerical experiences with this scheme are presented. 13 refs., 15 figs., 1 tab.« less
  • The development and applications of a parallel, time-dependent incompressible Navier-Stokes flow solver and a parallel multigrid elliptic kernel are described. The flow solver is based on a second-order projection method applied to a staggered finite-difference grid. The multigrid algorithms implemented in the parallel elliptic kernel, which is needed by the flow solver, and V-cycle and full V-cycle schemes. A grid-partition strategy is used in the parallel implementations of both the flow solver and the multigrid elliptic kernel on all fine and coarse grids. Numerical experiments and parallel performance tests show the parallel solver package is numerically stable, physically robust, andmore » computationally efficient. Both the multigrid elliptic kernel and the flow solver scale very well to a large number of processors on Intel Paragon and Cray T3D for computations with moderate granularity. The solver package has been carefully designed and coded so that it can be easily adapted to solving a variety of interesting two- and three-dimensional flow problems. The solver package is portable to parallel systems that support MPI, PVM, and NX for interprocessor communications. 13 refs., 23 figs., 3 tabs.« less
  • This paper presents a new forward-in-time advection method for nearly incompressible flow, MU, and its application to an adaptive multilevel flow solver for atmospheric flows. MU is a modification of Leonard et al.`s UTOPIA scheme. MU, like UTOPIA, is based on third-order accurate semi-Lagrangian multidimensional upwinding for constant velocity flows. for varying velocity fields, MU is a second-order conservative method. MU has greater stability and accuracy than UTOPIA and naturally decomposes into a monotone low-order method and a higher-order accurate correction for use with flux limiting. Its stability and accuracy make it a computationally efficient alternative to current finite-difference advectionmore » methods. We present a fully second-order accurate flow solver for the anelastic equations, a prototypical low Mach number flow. The flow solver is based on MU which is used for both momentum and scalar transport equations. This flow solver can also be implemented with any forward-in-time advection scheme. The multilevel flow solver conserves discrete global integrals of advected quantities and includes adaptive mesh refinements. Its second-order accuracy is verified using a nonlinear energy conservation integral for the anelastic equations. For a typical geophysical problem in which the flow is most rapidly varying in a small part of the domain, the multilevel flow solver achieves global accuracy comparable to uniform-resolution simulation for 10% of the computational cost. 36 refs., 10 figs.« less
  • This paper presents a parallel implementation of fractional solvers for the incompressible Navier-Stokes equations using an algebraic approach. Under this framework, predictor-corrector and incremental projection schemes are seen as sub-classes of the same class, making apparent its differences and similarities. An additional advantage of this approach is to set a common basis for a parallelization strategy, which can be extended to other split techniques or to compressible flows. The predictor-corrector scheme consists in solving the momentum equation and a modified 'continuity' equation (namely a simple iteration for the pressure Schur complement) consecutively in order to converge to the monolithic solution,more » thus avoiding fractional errors. On the other hand, the incremental projection scheme solves only one iteration of the predictor-corrector per time step and adds a correction equation to fulfill the mass conservation. As shown in the paper, these two schemes are very well suited for massively parallel implementation. In fact, when compared with monolithic schemes, simpler solvers and preconditioners can be used to solve the non-symmetric momentum equations (GMRES, Bi-CGSTAB) and to solve the symmetric continuity equation (CG, Deflated CG). This gives good speedup properties of the algorithm. The implementation of the mesh partitioning technique is presented, as well as the parallel performances and speedups for thousands of processors.« less
  • The purpose of this paper is to present a numerical study of flow fields for the NREL S805 and S809 airfoils using a spatially second-order symmetric total variational diminishing scheme. The steady two-dimensional flow is modeled as turbulent, viscous, and incompressible and is formulated in the pseudo-compressible form. The turbulent flow is closed by the Baldwind-Lomax algebraic turbulence model. Numerical solutions are obtained by the implicit approximate-factorization method. Numerical solutions are obtained by the implicit approximate-factorization method. The accuracy of the numerical results is compared with the Delft two-dimensional wind tunnel test data. For comparison, the Eppler code results aremore » also included. Numerical solutions of pressure and life coefficients show good agreement with the experimental data, but not the drag coefficients. To properly simulate the post-stall flow field, a better turbulence model should be used.« less