Comparison of OpenFOAM and EllipSys3D for neutral atmospheric flow over complex terrain

Research output: Contribution to journalJournal articleResearchpeer-review

1139 Downloads (Pure)

Abstract

The flow solvers OpenFOAM and EllipSys3D are compared in the case of neutral atmospheric flow over terrain using the test cases of Askervein and Bolund hills. Both solvers are run using the steady-state Reynolds-averaged Navier–Stokes k– turbulence model. One of the main modeling differences between the two solvers is the wall-function approach. The Open-FOAM v.1.7.1 uses a Nikuradse’s sand roughness model, while EllipSys3D uses a model based on the atmospheric roughness length. It is found that Nikuradse’s model introduces an error dependent on the near-wall cell height. To mitigate this error the near-wall cells should be at least 10 times larger than the surface roughness. It is nonetheless possible to obtain very similar results between EllipSys3D and OpenFOAM v.1.7.1. The more recent OpenFOAM v.2.2.1, which includes the atmospheric roughness length wall-function approach, has also been tested and compared to the results of OpenFOAM v.1.7.1 and EllipSys3D. The numerical results obtained using the same wall-modeling approach in both EllipSys3D and Open-FOAM v.2.1.1 proved to be almost identical. Two meshing strategies are investigated using HypGrid and SnappyHexMesh. The performance of Open-FOAM on SnappyHexMesh-based low-aspect-ratio unstructured meshes is found to be almost an order of magnitude faster than on HypGrid-based structured and high-aspect-ratio meshes. However, proper control of boundary layer resolution is found to be very difficult when the SnappyHexMesh tool is utilized for grid generation purposes. The OpenFOAM is generally found to be 2–6 times slower than EllipSys3D in achieving numerical results of the same order of accuracy on similar or identical computational meshes, when utilization of EllipSys3D default grid sequencing procedures is included.
Original languageEnglish
JournalWind Energy Science
Volume1
Pages (from-to)55–70
ISSN2366-7443
DOIs
Publication statusPublished - 2016

Bibliographical note

© Author(s) 2016. CC Attribution 3.0 License.

Cite this

@article{fa740d5a500b4f92ba64d1eb98ccf360,
title = "Comparison of OpenFOAM and EllipSys3D for neutral atmospheric flow over complex terrain",
abstract = "The flow solvers OpenFOAM and EllipSys3D are compared in the case of neutral atmospheric flow over terrain using the test cases of Askervein and Bolund hills. Both solvers are run using the steady-state Reynolds-averaged Navier–Stokes k– turbulence model. One of the main modeling differences between the two solvers is the wall-function approach. The Open-FOAM v.1.7.1 uses a Nikuradse’s sand roughness model, while EllipSys3D uses a model based on the atmospheric roughness length. It is found that Nikuradse’s model introduces an error dependent on the near-wall cell height. To mitigate this error the near-wall cells should be at least 10 times larger than the surface roughness. It is nonetheless possible to obtain very similar results between EllipSys3D and OpenFOAM v.1.7.1. The more recent OpenFOAM v.2.2.1, which includes the atmospheric roughness length wall-function approach, has also been tested and compared to the results of OpenFOAM v.1.7.1 and EllipSys3D. The numerical results obtained using the same wall-modeling approach in both EllipSys3D and Open-FOAM v.2.1.1 proved to be almost identical. Two meshing strategies are investigated using HypGrid and SnappyHexMesh. The performance of Open-FOAM on SnappyHexMesh-based low-aspect-ratio unstructured meshes is found to be almost an order of magnitude faster than on HypGrid-based structured and high-aspect-ratio meshes. However, proper control of boundary layer resolution is found to be very difficult when the SnappyHexMesh tool is utilized for grid generation purposes. The OpenFOAM is generally found to be 2–6 times slower than EllipSys3D in achieving numerical results of the same order of accuracy on similar or identical computational meshes, when utilization of EllipSys3D default grid sequencing procedures is included.",
author = "Dalibor Cavar and Pierre-Elouan R{\'e}thor{\'e} and Andreas Bechmann and S{\o}rensen, {Niels N.} and Benjamin Martinez and Frederik Zahle and Jacob Berg and Kelly, {Mark C.}",
note = "{\circledC} Author(s) 2016. CC Attribution 3.0 License.",
year = "2016",
doi = "10.5194/wes-2016-3",
language = "English",
volume = "1",
pages = "55–70",
journal = "Wind Energy Science",
issn = "2366-7443",
publisher = "Copernicus GmbH",

}

Comparison of OpenFOAM and EllipSys3D for neutral atmospheric flow over complex terrain. / Cavar, Dalibor; Réthoré, Pierre-Elouan; Bechmann, Andreas; Sørensen, Niels N.; Martinez, Benjamin; Zahle, Frederik; Berg, Jacob; Kelly, Mark C.

In: Wind Energy Science, Vol. 1, 2016, p. 55–70.

Research output: Contribution to journalJournal articleResearchpeer-review

TY - JOUR

T1 - Comparison of OpenFOAM and EllipSys3D for neutral atmospheric flow over complex terrain

AU - Cavar, Dalibor

AU - Réthoré, Pierre-Elouan

AU - Bechmann, Andreas

AU - Sørensen, Niels N.

AU - Martinez, Benjamin

AU - Zahle, Frederik

AU - Berg, Jacob

AU - Kelly, Mark C.

N1 - © Author(s) 2016. CC Attribution 3.0 License.

PY - 2016

Y1 - 2016

N2 - The flow solvers OpenFOAM and EllipSys3D are compared in the case of neutral atmospheric flow over terrain using the test cases of Askervein and Bolund hills. Both solvers are run using the steady-state Reynolds-averaged Navier–Stokes k– turbulence model. One of the main modeling differences between the two solvers is the wall-function approach. The Open-FOAM v.1.7.1 uses a Nikuradse’s sand roughness model, while EllipSys3D uses a model based on the atmospheric roughness length. It is found that Nikuradse’s model introduces an error dependent on the near-wall cell height. To mitigate this error the near-wall cells should be at least 10 times larger than the surface roughness. It is nonetheless possible to obtain very similar results between EllipSys3D and OpenFOAM v.1.7.1. The more recent OpenFOAM v.2.2.1, which includes the atmospheric roughness length wall-function approach, has also been tested and compared to the results of OpenFOAM v.1.7.1 and EllipSys3D. The numerical results obtained using the same wall-modeling approach in both EllipSys3D and Open-FOAM v.2.1.1 proved to be almost identical. Two meshing strategies are investigated using HypGrid and SnappyHexMesh. The performance of Open-FOAM on SnappyHexMesh-based low-aspect-ratio unstructured meshes is found to be almost an order of magnitude faster than on HypGrid-based structured and high-aspect-ratio meshes. However, proper control of boundary layer resolution is found to be very difficult when the SnappyHexMesh tool is utilized for grid generation purposes. The OpenFOAM is generally found to be 2–6 times slower than EllipSys3D in achieving numerical results of the same order of accuracy on similar or identical computational meshes, when utilization of EllipSys3D default grid sequencing procedures is included.

AB - The flow solvers OpenFOAM and EllipSys3D are compared in the case of neutral atmospheric flow over terrain using the test cases of Askervein and Bolund hills. Both solvers are run using the steady-state Reynolds-averaged Navier–Stokes k– turbulence model. One of the main modeling differences between the two solvers is the wall-function approach. The Open-FOAM v.1.7.1 uses a Nikuradse’s sand roughness model, while EllipSys3D uses a model based on the atmospheric roughness length. It is found that Nikuradse’s model introduces an error dependent on the near-wall cell height. To mitigate this error the near-wall cells should be at least 10 times larger than the surface roughness. It is nonetheless possible to obtain very similar results between EllipSys3D and OpenFOAM v.1.7.1. The more recent OpenFOAM v.2.2.1, which includes the atmospheric roughness length wall-function approach, has also been tested and compared to the results of OpenFOAM v.1.7.1 and EllipSys3D. The numerical results obtained using the same wall-modeling approach in both EllipSys3D and Open-FOAM v.2.1.1 proved to be almost identical. Two meshing strategies are investigated using HypGrid and SnappyHexMesh. The performance of Open-FOAM on SnappyHexMesh-based low-aspect-ratio unstructured meshes is found to be almost an order of magnitude faster than on HypGrid-based structured and high-aspect-ratio meshes. However, proper control of boundary layer resolution is found to be very difficult when the SnappyHexMesh tool is utilized for grid generation purposes. The OpenFOAM is generally found to be 2–6 times slower than EllipSys3D in achieving numerical results of the same order of accuracy on similar or identical computational meshes, when utilization of EllipSys3D default grid sequencing procedures is included.

U2 - 10.5194/wes-2016-3

DO - 10.5194/wes-2016-3

M3 - Journal article

VL - 1

SP - 55

EP - 70

JO - Wind Energy Science

JF - Wind Energy Science

SN - 2366-7443

ER -