skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Scalability of OpenFOAM Density-Based Solver with Runge–Kutta Temporal Discretization Scheme

Journal Article · · Scientific Programming
DOI:https://doi.org/10.1155/2020/9083620· OSTI ID:1604122
ORCiD logo [1]; ORCiD logo [2]; ORCiD logo [3]
  1. Department of Mechanical and Industrial Engineering, University of Illinois at Chicago, Chicago, IL, USA
  2. Department of Mechanical and Industrial Engineering, University of Illinois at Chicago, Chicago, IL, USA, Computational Science Division and Leadership Computing Facility, Argonne National Laboratory, Lemont, IL, USA
  3. Intel Corporation, Schaumburg, IL, USA

Compressible density-based solvers are widely used in OpenFOAM, and the parallel scalability of these solvers is crucial for large-scale simulations. In this paper, we report our experiences with the scalability of OpenFOAM’s native rhoCentralFoam solver, and by making a small number of modifications to it, we show the degree to which the scalability of the solver can be improved. The main modification made is to replace the first-order accurate Euler scheme in rhoCentralFoam with a third-order accurate, four-stage Runge-Kutta or RK4 scheme for the time integration. The scaling test we used is the transonic flow over the ONERA M6 wing. This is a common validation test for compressible flows solvers in aerospace and other engineering applications. Numerical experiments show that our modified solver, referred to as rhoCentralRK4Foam, for the same spatial discretization, achieves as much as a 123.2% improvement in scalability over the rhoCentralFoam solver. As expected, the better time resolution of the Runge–Kutta scheme makes it more suitable for unsteady problems such as the Taylor–Green vortex decay where the new solver showed a 50% decrease in the overall time-to-solution compared to rhoCentralFoam to get to the final solution with the same numerical accuracy. Finally, the improved scalability can be traced to the improvement of the computation to communication ratio obtained by substituting the RK4 scheme in place of the Euler scheme. All numerical tests were conducted on a Cray XC40 parallel system, Theta, at Argonne National Laboratory.

Sponsoring Organization:
USDOE
Grant/Contract Number:
#ANL 4J-30361-0030A
OSTI ID:
1604122
Journal Information:
Scientific Programming, Journal Name: Scientific Programming Vol. 2020; ISSN 1058-9244
Publisher:
Hindawi Publishing CorporationCopyright Statement
Country of Publication:
Egypt
Language:
English
Citation Metrics:
Cited by: 3 works
Citation information provided by
Web of Science

References (18)

Application of a fractional-step method to incompressible Navier-Stokes equations journal June 1985
The numerical simulation of two-dimensional fluid flow with strong shocks journal April 1984
Scalability of OpenFOAM for bio-medical flow simulations journal November 2014
The Tau Parallel Performance System journal May 2006
A New Flux Splitting Scheme journal July 1993
On the implementation of low-dissipative Runge–Kutta projection methods for time dependent flows using OpenFOAM® journal April 2014
Numerical solution of unsteady Navier–Stokes equations on curvilinear meshes journal June 2012
Factorization methods for the numerical approximation of Navier–Stokes equations journal July 2000
Implementation of semi-discrete, non-staggered central schemes in a colocated, polyhedral, finite volume framework, for high-speed viscous flows journal January 2009
Communication Optimization for Multiphase Flow Solver in the Library of OpenFOAM journal October 2018
Implementation of density-based solver for all speeds in the framework of OpenFOAM journal October 2014
New High-Resolution Central Schemes for Nonlinear Conservation Laws and Convection–Diffusion Equations journal May 2000
Modeling of Ice Accretion over Aircraft Wings Using a Compressible OpenFOAM Solver journal June 2019
Theta: Rapid installation and acceptance of an XC40 KNL system: Theta: Rapid installation and acceptance of an XC40 KNL system journal September 2017
A low-dissipative solver for turbulent compressible flows on unstructured meshes, with OpenFOAM implementation journal July 2017
A tensorial approach to computational continuum mechanics using object-oriented techniques journal January 1998
reactingFoam-SCI: An open source CFD platform for reacting flow simulation journal August 2019
PT-Scotch: A tool for efficient parallel graph ordering journal July 2008

Similar Records

Modeling of Ice Accretion over Aircraft Wings Using a Compressible OpenFOAM Solver
Journal Article · Mon Jun 03 00:00:00 EDT 2019 · International Journal of Aerospace Engineering · OSTI ID:1604122

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs
Journal Article · Mon Mar 14 00:00:00 EDT 2022 · SIAM Journal on Scientific Computing · OSTI ID:1604122

Stage-parallel fully implicit Runge–Kutta solvers for discontinuous Galerkin fluid simulations
Journal Article · Fri Jan 27 00:00:00 EST 2017 · Journal of Computational Physics · OSTI ID:1604122

Related Subjects