An assessment of atypical mesh topologies for lowMach largeeddy simulation
Abstract
An implicit, lowdissipation, lowMach, variable density control volume finite element formulation is used to explore foundational understanding of numerical accuracy for largeeddy simulation applications on hybrid meshes. Detailed simulation comparisons are made between loworder hexahedral, tetrahedral, pyramid, and wedge/prism topologies against a thirdorder, unstructured hexahedral topology. Using smooth analytical and manufactured lowMach solutions, designorder convergence is established for the hexahedral, tetrahedral, pyramid, and wedge element topologies using a new open boundary condition based on energystable methodologies previously deployed within a finitedifference context. A wide range of simulations demonstrate that loworder hexahedral and wedgebased element topologies behave nearly identically in both computed numerical errors and overall simulation timings. Moreover, loworder tetrahedral and pyramid element topologies also display nearly the same numerical characteristics. Although the superiority of the hexahedralbased topology is clearly demonstrated for trivial laminar, principallyaligned flows, e.g., a 1x2x10 channel flow with specified pressure drop, this advantage is reduced for nonaligned, 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 topologymore »
 Authors:

 Sandia National Lab. (SNLNM), Albuquerque, NM (United States). Computational Thermal and Fluid Mechanics Dept.
 Sandia National Lab. (SNLNM), Albuquerque, NM (United States). Aerosciences Dept.
 Publication Date:
 Research Org.:
 Sandia National Lab. (SNLNM), Albuquerque, NM (United States)
 Sponsoring Org.:
 USDOE National Nuclear Security Administration (NNSA)
 OSTI Identifier:
 1487422
 Report Number(s):
 SAND20189017J
Journal ID: ISSN 00457930; 667175
 Grant/Contract Number:
 AC0494AL85000
 Resource Type:
 Journal Article: Accepted Manuscript
 Journal Name:
 Computers and Fluids
 Additional Journal Information:
 Journal Volume: 179; Journal Issue: C; Journal ID: ISSN 00457930
 Publisher:
 Elsevier
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; LowMach; Largeeddy simulation; Control volume finite element; Hybrid; Higherorder; Verification; Validation
Citation Formats
Domino, Stefan P., Sakievich, Philip, and Barone, Matthew F. An assessment of atypical mesh topologies for lowMach largeeddy 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 lowMach largeeddy simulation. United States. doi:10.1016/j.compfluid.2018.12.002.
Domino, Stefan P., Sakievich, Philip, and Barone, Matthew F. Thu .
"An assessment of atypical mesh topologies for lowMach largeeddy simulation". United States. doi: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 lowMach largeeddy simulation},
author = {Domino, Stefan P. and Sakievich, Philip and Barone, Matthew F.},
abstractNote = {An implicit, lowdissipation, lowMach, variable density control volume finite element formulation is used to explore foundational understanding of numerical accuracy for largeeddy simulation applications on hybrid meshes. Detailed simulation comparisons are made between loworder hexahedral, tetrahedral, pyramid, and wedge/prism topologies against a thirdorder, unstructured hexahedral topology. Using smooth analytical and manufactured lowMach solutions, designorder convergence is established for the hexahedral, tetrahedral, pyramid, and wedge element topologies using a new open boundary condition based on energystable methodologies previously deployed within a finitedifference context. A wide range of simulations demonstrate that loworder hexahedral and wedgebased element topologies behave nearly identically in both computed numerical errors and overall simulation timings. Moreover, loworder tetrahedral and pyramid element topologies also display nearly the same numerical characteristics. Although the superiority of the hexahedralbased topology is clearly demonstrated for trivial laminar, principallyaligned flows, e.g., a 1x2x10 channel flow with specified pressure drop, this advantage is reduced for nonaligned, 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 wallresolved largeeddy simulation, the increased assembly and residual evaluation computational time for nonhexahedral 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},
issn = {00457930},
number = C,
volume = 179,
place = {United States},
year = {2018},
month = {12}
}
Web of Science