A hybrid incremental projection method for thermal-hydraulics applications

Christon, Mark A.; Bakosi, Jozsef; Nadiga, Balasubramanya T.; Berndt, Markus; Francois, Marianne M.; Stagg, Alan K.; Xia, Yidong; Luo, Hong

doi:10.1016/j.jcp.2016.04.061

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

Full Record
Other Related Research

Abstract

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 »« less

Authors:

Christon, Mark A. ^[1]; Bakosi, Jozsef ^[2]; Nadiga, Balasubramanya T. ^[2];

^[2]; Francois, Marianne M. ^[2];

^[3];

^[4]; Luo, Hong ^[5]

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)

Publication Date:: Fri Jul 01 00:00:00 EDT 2016

Research Org.:: Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)

Sponsoring Org.:: USDOE

OSTI Identifier:: 1256295

Alternate Identifier(s):: OSTI ID: 1261272; OSTI ID: 1347629

Report Number(s):: LA-UR-16-22436
Journal ID: ISSN 0021-9991

Grant/Contract Number:: AC05-00OR22725; AC52-06NA25396

Resource Type:: Accepted Manuscript

Journal Name:: Journal of Computational Physics

Additional Journal Information:: Journal Volume: 317; Journal Issue: C; Journal ID: ISSN 0021-9991

Publisher:: Elsevier

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; FVM; FEM; Incompressible flow; Monotonicity-preserving advection; Projection method; Mixed-topology meshes; Thermal-hydraulics

Citation Formats


                    Christon, Mark A., Bakosi, Jozsef, Nadiga, Balasubramanya T., Berndt, Markus, Francois, Marianne M., Stagg, Alan K., Xia, Yidong, and Luo, Hong. A hybrid incremental projection method for thermal-hydraulics applications.  United States: N. p., 2016. 
Web.  doi:10.1016/j.jcp.2016.04.061.

Copy to clipboard


                    Christon, Mark A., Bakosi, Jozsef, Nadiga, Balasubramanya T., Berndt, Markus, Francois, Marianne M., Stagg, Alan K., Xia, Yidong, & Luo, Hong. A hybrid incremental projection method for thermal-hydraulics applications.  United States.  https://doi.org/10.1016/j.jcp.2016.04.061

Copy to clipboard


                    Christon, Mark A., Bakosi, Jozsef, Nadiga, Balasubramanya T., Berndt, Markus, Francois, Marianne M., Stagg, Alan K., Xia, Yidong, and Luo, Hong. Fri .  
"A hybrid incremental projection method for thermal-hydraulics applications".  United States.  https://doi.org/10.1016/j.jcp.2016.04.061.  https://www.osti.gov/servlets/purl/1256295.

Copy to clipboard


                    
@article{osti_1256295,

  title        = {A hybrid incremental projection method for thermal-hydraulics applications},

  author       = {Christon, Mark A. and Bakosi, Jozsef and Nadiga, Balasubramanya T. and Berndt, Markus and Francois, Marianne M. and Stagg, Alan K. and Xia, Yidong and Luo, Hong},

  abstractNote = {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.},

  doi          = {10.1016/j.jcp.2016.04.061},

  journal      = {Journal of Computational Physics},

  number       = C,

  volume       = 317,

  place        = {United States},

  year         = {Fri Jul 01 00:00:00 EDT 2016},

  month        = {Fri Jul 01 00:00:00 EDT 2016}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1016/j.jcp.2016.04.061

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 4 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Similar Records in DOE PAGES and OSTI.GOV collections:

High-order fully implicit solver for all-speed fluid dynamics: AUSM ride from nearly incompressible variable-density flows to shock dynamics

Journal Article Nourgaliev, R. ; Greene, P. ; Weston, B. ; ... - Shock Waves

We introduce a novel Newton–Krylov (NK)-based fully implicit algorithm for solving fluid flows in a wide range of flow conditions—from variable density nearly incompressible to supersonic shock dynamics. The key enabling feature of our all-speed solver is the ability to efficiently solve conservation laws by choosing a set of independent variables that produce a well-conditioned Jacobian matrix for the linear iterations of the global nonlinear iterative solver. In particular, instead of choosing to discretize the conservative variables (density, momentum, total energy), which is traditionally used in Eulerian high-speed compressible fluid dynamics, we demonstrate superior performance by discretizing the primitive variables—pressure–velocity–temperaturemore »« less
Cited by 2
https://doi.org/10.1007/s00193-018-0871-8

Full Text Available
Improvement of Numerical Methods in Pronghorn

Technical Report Lindsay, Alexander D. ; Tano Retamales, Mauricio Eduardo ; Giudicelli, Guillaume Louis ; ...

During the fiscal year of 2021, the finite-volume method (FVM) was deployed in Pronghorn to provide improved efficiency, stability, and accuracy for coarse-mesh, thermal-hydraulics problems. While the main goals of fiscal year 2021 were met, several issues emerged from the early deployment of the finite-volume method in Pronghorn. These issues were: The compressible and incompressible formulations that were implemented are inadequate for many nuclear reactor flow problems; Omission of terms accounting for the porosity and Darcy-Forchheimer body force discontinuities in the Rhie-Chow interpolation lead to oscillations in pressure and velocity at these discontinuities; The FVM lacks a correction for non-orthogonalmore »« less
https://doi.org/10.2172/2008350

Full Text Available
Interfacing MFIX with PETSC and HYPRE Linear Solver Libraries (Final Report)

Technical Report Krishnamoorthy, Gautham ; Thornock, Jeremy

A computational bottleneck during the solution to multiphase formulations of the incompressible Navier-Stokes equations is often during the implicit solution of the pressure correction equation that results from operator-splitting methods. Since density is a coefficient in the pressure-correction equation, large variations or discontinuities among the phase densities greatly increase the condition number of the pressure-correction matrix and retard the convergence of iterative methods employed in its solution. To alleviate this shortcoming, the open-source multiphase code MFiX is interfaced with the linear solver libraries PETSc and HYPRE. Through an appropriate mapping of matrix and vector data structures between the two softwaremore »« less
https://doi.org/10.2172/1592576

Full Text Available
A Comparison of Classical and Aggregation-Based Algebraic Multigrid Preconditioners for High-Fidelity Simulation of Wind Turbine Incompressible Flows

Journal Article Thomas, Stephen ; Ananthan, Shreyas ; Yellapantula, Shashank ; ... - SIAM Journal on Scientific Computing

This paper presents a comparison of parallel strong scaling performance of classical and aggregation algebraic multigrid (AMG) preconditioners in the context of wind turbine simulations. Fluid motion is governed by the incompressible Navier--Stokes equations, discretized in space with control-volume finite elements and in time with an inexact projection scheme using an implicit integrator. A discontinuous-Galerkin sliding-mesh algorithm captures rotor motion. The momentum equations are solved with iterative Krylov methods, preconditioned by symmetric Gauss--Seidel (SGS) in Trilinos andmore »« less
https://doi.org/10.1137/18M1179018
A Comparison of Classical and Aggregation-Based Algebraic Multigrid Preconditioners for High-Fidelity Simulation of Wind Turbine Incompressible Flows

Journal Article Thomas, Stephen J. ; Ananthan, Shreyas ; Yellapantula, Shashank ; ... - SIAM Journal on Scientific Computing

This paper presents a comparison of parallel strong scaling performance of classical and aggregation algebraic multigrid (AMG) preconditioners in the context of wind turbine simulations. Fluid motion is governed by the incompressible Navier--Stokes equations, discretized in space with control-volume finite elements and in time with an inexact projection scheme using an implicit integrator. A discontinuous-Galerkin sliding-mesh algorithm captures rotor motion. The momentum equations are solved with iterative Krylov methods, preconditioned by symmetric Gauss--Seidel (SGS) in Trilinos andmore »« less
Cited by 4
https://doi.org/10.1137/18M1179018

Full Text Available

Similar Records