DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: An assessment of atypical mesh topologies for low-Mach large-eddy simulation

Abstract

An implicit, low-dissipation, low-Mach, variable density control volume finite element formulation is used to explore foundational understanding of numerical accuracy for large-eddy simulation applications on hybrid meshes. Detailed simulation comparisons are made between low-order hexahedral, tetrahedral, pyramid, and wedge/prism topologies against a third-order, unstructured hexahedral topology. Using smooth analytical and manufactured low-Mach solutions, design-order convergence is established for the hexahedral, tetrahedral, pyramid, and wedge element topologies using a new open boundary condition based on energy-stable methodologies previously deployed within a finite-difference context. A wide range of simulations demonstrate that low-order hexahedral- and wedge-based element topologies behave nearly identically in both computed numerical errors and overall simulation timings. Moreover, low-order tetrahedral and pyramid element topologies also display nearly the same numerical characteristics. Although the superiority of the hexahedral-based topology is clearly demonstrated for trivial laminar, principally-aligned flows, e.g., a 1x2x10 channel flow with specified pressure drop, this advantage is reduced for non-aligned, turbulent flows including the Taylor–Green Vortex, turbulent plane channel flow (Reτ395), and buoyant flow past a heated cylinder. With the order of accuracy demonstrated for both homogenous and hybrid meshes, it is shown that solution verification for the selected complex flows can be established for all topology types. Althoughmore » the number of elements in a mesh of like spacing comprised of tetrahedral, wedge, or pyramid elements increases as compared to the hexahedral counterpart, for wall-resolved large-eddy simulation, the increased assembly and residual evaluation computational time for non-hexahedral is offset by more efficient linear solver times. Lastly, most simulation results indicate that modest polynomial promotion provides a significant increase in solution accuracy.« less

Authors:
 [1];  [1];  [2]
  1. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States). Computational Thermal and Fluid Mechanics Dept.
  2. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States). Aerosciences Dept.
Publication Date:
Research Org.:
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1487422
Alternate Identifier(s):
OSTI ID: 1636189
Report Number(s):
SAND-2018-9017J
Journal ID: ISSN 0045-7930; 667175
Grant/Contract Number:  
AC04-94AL85000; 17-SC-20-SC
Resource Type:
Accepted Manuscript
Journal Name:
Computers and Fluids
Additional Journal Information:
Journal Volume: 179; Journal Issue: C; Journal ID: ISSN 0045-7930
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; Low-Mach; Large-eddy simulation; Control volume finite element; Hybrid; Higher-order; Verification; Validation

Citation Formats

Domino, Stefan P., Sakievich, Philip, and Barone, Matthew F. An assessment of atypical mesh topologies for low-Mach large-eddy simulation. United States: N. p., 2018. Web. doi:10.1016/j.compfluid.2018.12.002.
Domino, Stefan P., Sakievich, Philip, & Barone, Matthew F. An assessment of atypical mesh topologies for low-Mach large-eddy simulation. United States. https://doi.org/10.1016/j.compfluid.2018.12.002
Domino, Stefan P., Sakievich, Philip, and Barone, Matthew F. Thu . "An assessment of atypical mesh topologies for low-Mach large-eddy simulation". United States. https://doi.org/10.1016/j.compfluid.2018.12.002. https://www.osti.gov/servlets/purl/1487422.
@article{osti_1487422,
title = {An assessment of atypical mesh topologies for low-Mach large-eddy simulation},
author = {Domino, Stefan P. and Sakievich, Philip and Barone, Matthew F.},
abstractNote = {An implicit, low-dissipation, low-Mach, variable density control volume finite element formulation is used to explore foundational understanding of numerical accuracy for large-eddy simulation applications on hybrid meshes. Detailed simulation comparisons are made between low-order hexahedral, tetrahedral, pyramid, and wedge/prism topologies against a third-order, unstructured hexahedral topology. Using smooth analytical and manufactured low-Mach solutions, design-order convergence is established for the hexahedral, tetrahedral, pyramid, and wedge element topologies using a new open boundary condition based on energy-stable methodologies previously deployed within a finite-difference context. A wide range of simulations demonstrate that low-order hexahedral- and wedge-based element topologies behave nearly identically in both computed numerical errors and overall simulation timings. Moreover, low-order tetrahedral and pyramid element topologies also display nearly the same numerical characteristics. Although the superiority of the hexahedral-based topology is clearly demonstrated for trivial laminar, principally-aligned flows, e.g., a 1x2x10 channel flow with specified pressure drop, this advantage is reduced for non-aligned, turbulent flows including the Taylor–Green Vortex, turbulent plane channel flow (Reτ395), and buoyant flow past a heated cylinder. With the order of accuracy demonstrated for both homogenous and hybrid meshes, it is shown that solution verification for the selected complex flows can be established for all topology types. Although the number of elements in a mesh of like spacing comprised of tetrahedral, wedge, or pyramid elements increases as compared to the hexahedral counterpart, for wall-resolved large-eddy simulation, the increased assembly and residual evaluation computational time for non-hexahedral is offset by more efficient linear solver times. Lastly, most simulation results indicate that modest polynomial promotion provides a significant increase in solution accuracy.},
doi = {10.1016/j.compfluid.2018.12.002},
journal = {Computers and Fluids},
number = C,
volume = 179,
place = {United States},
year = {Thu Dec 06 00:00:00 EST 2018},
month = {Thu Dec 06 00:00:00 EST 2018}
}

Journal Article:

Citation Metrics:
Cited by: 13 works
Citation information provided by
Web of Science

Save / Share: