A hybrid incremental projection method for thermal-hydraulics applications

Journal Article · Fri Jul 01 00:00:00 EDT 2016 · Journal of Computational Physics

DOI:https://doi.org/10.1016/j.jcp.2016.04.061· OSTI ID:1256295

Christon, Mark A. ^[1]; Bakosi, Jozsef ^[2]; Nadiga, Balasubramanya T. ^[2]; Berndt, Markus ^[2]; Francois, Marianne M. ^[2]; Stagg, Alan K. ^[3]; Xia, Yidong ^[4]; Luo, Hong ^[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)

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 permit 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.

• Research Organization: Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Los Alamos National Laboratory (LANL); Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
• Sponsoring Organization: USDOE
• Grant/Contract Number: AC05-00OR22725; AC52-06NA25396
• OSTI ID: 1256295
• Alternate ID(s): OSTI ID: 1347629; OSTI ID: 1261272
• Report Number(s): LA-UR--16-22436; LA-UR-16-22436
• Journal Information: Journal of Computational Physics, Journal Name: Journal of Computational Physics Journal Issue: C Vol. 317; ISSN 0021-9991
• Publisher: Elsevier
• Country of Publication: United States
• Language: English

Related Subjects

97 MATHEMATICS AND COMPUTING
FEM
FVM
FVM
FEM
incompressible flow
monotonicity-preserving advection
projection method
mixed-topology meshes
thermal-hydraulics
Incompressible flow
Mixed-topology meshes
Monotonicity-preserving advection
Projection method
Thermal-hydraulics