A hybrid incremental projection method for thermal-hydraulics applications
| Autor: | Mark A. Christon, Alan K. Stagg, Yidong Xia, Hong Luo, Jozsef Bakosi, Markus Berndt, Balasubramanya T. Nadiga, Marianne M. Francois |
|---|---|
| Rok vydání: | 2016 |
| Předmět: |
Numerical Analysis
Mathematical optimization Physics and Astronomy (miscellaneous) Discretization Applied Mathematics Multiphysics Solver 01 natural sciences Finite element method Projection (linear algebra) 010305 fluids & plasmas Computer Science Applications Computational science 010101 applied mathematics Computational Mathematics Modeling and Simulation 0103 physical sciences Projection method 0101 mathematics Dykstra's projection algorithm Interpolation Mathematics |
| Zdroj: | Journal of Computational Physics. 317:382-404 |
| ISSN: | 0021-9991 |
| Popis: | 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-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 ? C F L ? 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. A new second-order hybrid finite-element/finite-volume projection algorithm for transient viscous flow has been introduced.The hybrid discretization prevents pressure modes without using Rhie-Chow interpolation or a Petrov-Galerkin formulation.A monotonicity-preserving advection method shown to deliver second-order accuracy on mixed element topology meshes.Verification studies demonstrate hybrid projection solver accuracy and temporal stability for super-CFL conditions. |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|
Full Text from ScienceDirect