A hybrid incremental projection method for thermal-hydraulics applications
- Computational Sciences International, Los Alamos, NM (United States); Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
- Idaho National Lab. (INL), Idaho Falls, ID (United States)
- North Carolina State Univ., Raleigh, NC (United States)
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/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); 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: 1261272; OSTI ID: 1347629
- Report Number(s):
- LA-UR-16-22436
- Journal Information:
- Journal of Computational Physics, Vol. 317, Issue C; ISSN 0021-9991
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Similar Records
Improvement of Numerical Methods in Pronghorn
Interfacing MFIX with PETSC and HYPRE Linear Solver Libraries (Final Report)