skip to main content

DOE PAGESDOE PAGES

Title: A hybrid incremental projection method for thermal-hydraulics applications

A new second-order accurate, hybrid, incremental projection method for time-dependent incompressible viscous flow is introduced. The hybrid finite-element/finite-volume discretization circumvents the well-known Ladyzhenskaya–Babuska–Brezzi conditions for stability, and does not require special treatment to filter pressure modes by either Rhie–Chow interpolation or by using a Petrov–Galerkin finite element formulation. The use of a co-velocity with a high-resolution advection method and a linearly consistent edge-based treatment of viscous/diffusive terms yields a robust algorithm for a broad spectrum of incompressible flows. The high-resolution advection method is shown to deliver second-order spatial convergence on mixed element topology meshes, and the implicit advective treatment significantly increases the stable time-step size. The algorithm is robust and extensible, permitting the incorporation of features such as porous media flow, RANS and LES turbulence models, and semi-/fully-implicit time stepping. A series of verification and validation problems are used to illustrate the convergence properties of the algorithm. The temporal stability properties are demonstrated on a range of problems with 2≤CFL≤100. The new flow solver is built using the Hydra multiphysics toolkit. Finally, the Hydra toolkit is written in C++ and provides a rich suite of extensible and fully-parallel components that permit rapid application development, supports multiple discretization techniques, provides I/Omore » interfaces, dynamic run-time load balancing and data migration, and interfaces to scalable popular linear solvers, e.g., in open-source packages such as HYPRE, PETSc, and Trilinos.« less
Authors:
 [1] ;  [2] ;  [2] ; ORCiD logo [2] ;  [3] ; ORCiD logo [4] ; ORCiD logo [5] ;  [6]
  1. Computational Sciences International, Los Alamos, NM (United States); Los Alamos National Lab. (LANL), Los Alamos, NM (United States). Computational Physics and Methods Group (CCS-2). Guest Scientist
  2. Los Alamos National Lab. (LANL), Los Alamos, NM (United States). Computational Physics and Methods Group (CCS-2)
  3. Los Alamos National Lab. (LANL), Los Alamos, NM (United States). Methods and Algorithms Group (XCP-4)
  4. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Modeling and Simulation Group
  5. Idaho National Lab. (INL), Idaho Falls, ID (United States). Dept. of Energy Resource Recovery and Sustainability
  6. North Carolina State Univ., Raleigh, NC (United States). Dept. of Mechanical and Aerospace Engineering
Publication Date:
Report Number(s):
LA-UR--16-22436
Journal ID: ISSN 0021-9991; NT0304000; NEAF343
Grant/Contract Number:
AC05-00OR22725; AC52-06NA25396
Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 317; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Research Org:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States); Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org:
USDOE
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING FVM; FEM; Incompressible flow; Monotonicity-preserving advection; Projection method; Mixed-topology meshes; Thermal-hydraulics
OSTI Identifier:
1261272
Alternate Identifier(s):
OSTI ID: 1261272