Spectral element applications in complex nuclear reactor geometries: Tet-to-hex meshing

Haomin Yuan, Mustafa A. Yildiz, Elia Merzari, Yiqi Yu, Aleksandr Obabko, Gerrit Botha, Giacomo Busco, Yassin A. Hassan, Duy Thien Nguyen

Research output: Contribution to journalArticle

Abstract

The spectral element code Nek5000 is an open-source, higher-order computational fluid dynamics code developed at Argonne National Laboratory. It is designed to solve incompressible Navier-Stokes equations, but it also has a low-Mach-number approximation feature available. Large eddy simulation is approached by explicit filtering of the velocity field (and other fields) to mimic the effect of dissipation due to the unresolved scale. The computational domain is decomposed into second-order hexahedral elements that conform to the boundaries. However, generating a high-quality pure-hexahedral mesh can be challenging for some problems. For simple geometries, traditional blocking methods can be used to decompose the domain into smaller blocks to generate a so-called structural mesh. A structural mesh can maintain good orthogonality but can have a highly skewed mesh to conform to the geometry, as well as unnecessary refinement in the far field. Moreover, for geometries with relative complexity, blocking the geometry becomes impossible. To address these issues, we adopted a tet-to-hex strategy to generate a pure hexahedral mesh for Nek5000. First, we generate a pure tetrahedral mesh for the geometry; then we divide one tetrahedral element into four hexahedral elements. A pure tetrahedral mesh could be easily generated for complex geometries by using many current meshing codes. In this paper, we use the commercial codes ANSYS meshing and ANSYS ICEM to generate the pure tetrahedral mesh and then convert it to a pure hexahedral mesh. Boundary layers are extruded in ANSYSICEM to maintain near-wall resolution.

Original languageEnglish (US)
Article number110422
JournalNuclear Engineering and Design
Volume357
DOIs
StatePublished - Feb 2020

Fingerprint

nuclear reactors
Nuclear reactors
mesh
geometry
Geometry
Navier-Stokes equations
large eddy simulation
Large eddy simulation
computational fluid dynamics
Navier Stokes equations
Mach number
nuclear reactor
dissipation
Computational fluid dynamics
Boundary layers
boundary layer
orthogonality
Navier-Stokes equation
far fields
code

All Science Journal Classification (ASJC) codes

  • Nuclear and High Energy Physics
  • Nuclear Energy and Engineering
  • Materials Science(all)
  • Safety, Risk, Reliability and Quality
  • Waste Management and Disposal
  • Mechanical Engineering

Cite this

Yuan, Haomin ; Yildiz, Mustafa A. ; Merzari, Elia ; Yu, Yiqi ; Obabko, Aleksandr ; Botha, Gerrit ; Busco, Giacomo ; Hassan, Yassin A. ; Nguyen, Duy Thien. / Spectral element applications in complex nuclear reactor geometries : Tet-to-hex meshing. In: Nuclear Engineering and Design. 2020 ; Vol. 357.
@article{72a9b14c86f24e4db127902795143f58,
title = "Spectral element applications in complex nuclear reactor geometries: Tet-to-hex meshing",
abstract = "The spectral element code Nek5000 is an open-source, higher-order computational fluid dynamics code developed at Argonne National Laboratory. It is designed to solve incompressible Navier-Stokes equations, but it also has a low-Mach-number approximation feature available. Large eddy simulation is approached by explicit filtering of the velocity field (and other fields) to mimic the effect of dissipation due to the unresolved scale. The computational domain is decomposed into second-order hexahedral elements that conform to the boundaries. However, generating a high-quality pure-hexahedral mesh can be challenging for some problems. For simple geometries, traditional blocking methods can be used to decompose the domain into smaller blocks to generate a so-called structural mesh. A structural mesh can maintain good orthogonality but can have a highly skewed mesh to conform to the geometry, as well as unnecessary refinement in the far field. Moreover, for geometries with relative complexity, blocking the geometry becomes impossible. To address these issues, we adopted a tet-to-hex strategy to generate a pure hexahedral mesh for Nek5000. First, we generate a pure tetrahedral mesh for the geometry; then we divide one tetrahedral element into four hexahedral elements. A pure tetrahedral mesh could be easily generated for complex geometries by using many current meshing codes. In this paper, we use the commercial codes ANSYS meshing and ANSYS ICEM to generate the pure tetrahedral mesh and then convert it to a pure hexahedral mesh. Boundary layers are extruded in ANSYSICEM to maintain near-wall resolution.",
author = "Haomin Yuan and Yildiz, {Mustafa A.} and Elia Merzari and Yiqi Yu and Aleksandr Obabko and Gerrit Botha and Giacomo Busco and Hassan, {Yassin A.} and Nguyen, {Duy Thien}",
year = "2020",
month = "2",
doi = "10.1016/j.nucengdes.2019.110422",
language = "English (US)",
volume = "357",
journal = "Nuclear Engineering and Design",
issn = "0029-5493",
publisher = "Elsevier BV",

}

Spectral element applications in complex nuclear reactor geometries : Tet-to-hex meshing. / Yuan, Haomin; Yildiz, Mustafa A.; Merzari, Elia; Yu, Yiqi; Obabko, Aleksandr; Botha, Gerrit; Busco, Giacomo; Hassan, Yassin A.; Nguyen, Duy Thien.

In: Nuclear Engineering and Design, Vol. 357, 110422, 02.2020.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Spectral element applications in complex nuclear reactor geometries

T2 - Tet-to-hex meshing

AU - Yuan, Haomin

AU - Yildiz, Mustafa A.

AU - Merzari, Elia

AU - Yu, Yiqi

AU - Obabko, Aleksandr

AU - Botha, Gerrit

AU - Busco, Giacomo

AU - Hassan, Yassin A.

AU - Nguyen, Duy Thien

PY - 2020/2

Y1 - 2020/2

N2 - The spectral element code Nek5000 is an open-source, higher-order computational fluid dynamics code developed at Argonne National Laboratory. It is designed to solve incompressible Navier-Stokes equations, but it also has a low-Mach-number approximation feature available. Large eddy simulation is approached by explicit filtering of the velocity field (and other fields) to mimic the effect of dissipation due to the unresolved scale. The computational domain is decomposed into second-order hexahedral elements that conform to the boundaries. However, generating a high-quality pure-hexahedral mesh can be challenging for some problems. For simple geometries, traditional blocking methods can be used to decompose the domain into smaller blocks to generate a so-called structural mesh. A structural mesh can maintain good orthogonality but can have a highly skewed mesh to conform to the geometry, as well as unnecessary refinement in the far field. Moreover, for geometries with relative complexity, blocking the geometry becomes impossible. To address these issues, we adopted a tet-to-hex strategy to generate a pure hexahedral mesh for Nek5000. First, we generate a pure tetrahedral mesh for the geometry; then we divide one tetrahedral element into four hexahedral elements. A pure tetrahedral mesh could be easily generated for complex geometries by using many current meshing codes. In this paper, we use the commercial codes ANSYS meshing and ANSYS ICEM to generate the pure tetrahedral mesh and then convert it to a pure hexahedral mesh. Boundary layers are extruded in ANSYSICEM to maintain near-wall resolution.

AB - The spectral element code Nek5000 is an open-source, higher-order computational fluid dynamics code developed at Argonne National Laboratory. It is designed to solve incompressible Navier-Stokes equations, but it also has a low-Mach-number approximation feature available. Large eddy simulation is approached by explicit filtering of the velocity field (and other fields) to mimic the effect of dissipation due to the unresolved scale. The computational domain is decomposed into second-order hexahedral elements that conform to the boundaries. However, generating a high-quality pure-hexahedral mesh can be challenging for some problems. For simple geometries, traditional blocking methods can be used to decompose the domain into smaller blocks to generate a so-called structural mesh. A structural mesh can maintain good orthogonality but can have a highly skewed mesh to conform to the geometry, as well as unnecessary refinement in the far field. Moreover, for geometries with relative complexity, blocking the geometry becomes impossible. To address these issues, we adopted a tet-to-hex strategy to generate a pure hexahedral mesh for Nek5000. First, we generate a pure tetrahedral mesh for the geometry; then we divide one tetrahedral element into four hexahedral elements. A pure tetrahedral mesh could be easily generated for complex geometries by using many current meshing codes. In this paper, we use the commercial codes ANSYS meshing and ANSYS ICEM to generate the pure tetrahedral mesh and then convert it to a pure hexahedral mesh. Boundary layers are extruded in ANSYSICEM to maintain near-wall resolution.

UR - http://www.scopus.com/inward/record.url?scp=85074641090&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85074641090&partnerID=8YFLogxK

U2 - 10.1016/j.nucengdes.2019.110422

DO - 10.1016/j.nucengdes.2019.110422

M3 - Article

AN - SCOPUS:85074641090

VL - 357

JO - Nuclear Engineering and Design

JF - Nuclear Engineering and Design

SN - 0029-5493

M1 - 110422

ER -