skip to main content

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

In this paper, a new second-order accurate, hybrid, incremental projection method for time-dependent incompressible viscous flow is introduced in this paper. The hybrid finite-element/finite-volume discretization circumvents the well-known Ladyzhenskaya–Babuška–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. The Hydra toolkit is written in C++ and provides a rich suite of extensible and fully-parallel components that permitmore » rapid application development, supports multiple discretization techniques, provides I/O 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
 [1] ;  [2] ;  [2] ; ORCiD logo [2] ;  [2] ; ORCiD logo [3] ; ORCiD logo [4] ;  [5]
  1. Computational Sciences International, Los Alamos, NM (United States); Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
  2. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
  3. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
  4. Idaho National Lab. (INL), Idaho Falls, ID (United States)
  5. North Carolina State Univ., Raleigh, NC (United States)
Publication Date:
OSTI Identifier:
Report Number(s):
Journal ID: ISSN 0021-9991
Grant/Contract Number:
AC05-00OR22725; AC52-06NA25396
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 317; Journal Issue: C; Journal ID: ISSN 0021-9991
Research Org:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org:
Country of Publication:
United States
97 MATHEMATICS AND COMPUTING; FVM, FEM, incompressible flow, monotonicity-preserving advection, projection method, mixed-topology meshes, thermal-hydraulics; FVM; FEM; Incompressible flow; Monotonicity-preserving advection; Projection method; Mixed-topology meshes; Thermal-hydraulics