Abstract:A low-dissipative solution framework integrating various types of high-order time scheme is proposed and implemented based on the open-source C++ library OpenFOAM. This framework aims to introduce different categories of low-dissipative time integration schemes into a unified solver convenient for comparison of scheme performance in finite volume computational fluid dynamics code, contributing to the development of low dissipation scheme appropriate for scale-resolving turbulence simulation. To demonstrate this general framework's ability of including a wide range of time integration method, in addition to typical Runge--Kutta family schemes of linear single-step method, two more complex linear multi-step method, Adams--Bashforth family and Adams--Bashforth--Moutton family schemes, are implemented with the projection algorithm, which increase the options of time discretization. The unified solver obtained by the solution framework select the specified time scheme from a variety of alternatives in runtime rather than maintaining multiple solvers with each compiled for a single scheme, while new scheme can be easily added according to the basic idea of this universal framework. Further research on the stability of the explicit scheme indicates that a multi-step method with an appropriate order may be an optimal choice when taking both the prediction accuracy and computational efficiency into account in unstable problems.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to develop a low - dissipation projection method framework that can integrate various high - order time integration schemes and be implemented based on the open - source C++ library OpenFOAM. Specifically, the objectives of the paper include: 1. **Provide a low - dissipation projection framework**: This framework can dynamically select different time integration schemes at runtime and is compatible with OpenFOAM's native architecture. Through this framework, researchers can easily compare the performance of different high - order time integration schemes in finite - volume computational fluid dynamics codes. 2. **Implement multiple high - order time integration schemes**: In addition to the common single - step methods of the Runge - Kutta family, linear multi - step methods of the Adams - Bashforth family and the Adams - Bashforth - Moulton family are also implemented. This increases the time discretization options, making the framework more general and flexible. 3. **Evaluate the advantages and disadvantages of the low - dissipation solver**: Through a series of benchmark test cases, comprehensively evaluate the advantages and disadvantages of the low - dissipation solver compared to traditional implicit solvers from the perspectives of computational accuracy and computational speed. 4. **Suggest how to better select low - dissipation schemes**: Through the study of the differences in accuracy and computational speed of different methods, propose how to better select low - dissipation schemes when performing scale - resolving turbulence simulations in industrial applications. ### Main contributions of the paper 1. **Unified low - dissipation projection solver**: Provides a unified low - dissipation projection solver, which provides more options for simulating practical engineering problems and helps to compare the real performance of different high - order methods in specific fluid problems. 2. **Convenience of the general framework**: Similar to OpenFOAM's native solver, it is convenient for researchers interested in numerical methods to implement multiple high - order low - dissipation schemes and carry out their own scheme designs to improve performance. 3. **Wide applicability**: Supports OpenFOAM's native source terms, turbulence models and other classes, and is applicable to a wider range of problems, which helps to evaluate various existing algorithms in specific fluid types. ### Implementation details The paper details the implementation details of the framework, including: - **Basic structure**: The basic structure of the framework is shown in Figure 2. By introducing multiple scheme processing modules in the native solver architecture, multi - scheme selection and invocation are realized. - **Composition and invocation of the projection method**: Describes in detail the implementation of the projection method, including the code blocks for the intermediate velocity calculation step and the Poisson equation solving step. - **Modular design**: Separates modules with different functions, facilitating the subsequent introduction of other low - dissipation high - order methods or the change of coupling algorithms. ### Conclusion Through this low - dissipation general solution framework, researchers can easily compare the performance of different time integration schemes, thereby finding the low - dissipation scheme most suitable for specific fluid problems. In addition, this framework also maintains support for OpenFOAM's native features, allowing other researchers to use it directly without a large number of compatibility modifications.

Development of low-dissipative projection method framework integrating various high-order time integration schemes using OpenFOAM

High-fidelity Simulation of Regular Waves Based on Multi-Moment Finite Volume Formulation and THINC Method

An Optimized Dispersion-Relation-Preserving Combined Compact Difference Scheme to Solve Advection Equations

A class of high-order finite difference schemes with minimized dispersion and adaptive dissipation for solving compressible flows

A family of TENOA-THINC-MOOD schemes based on diffuse-interface method for compressible multiphase flows

An Advection Scheme Using Paired Explicit Runge-Kutta Time Integration for Atmospheric Modeling

An Accurate and Efficient Multiphase Solver Based on THINC Scheme and Adaptive Mesh Refinement

A parallel solver framework for fully implicit monolithic fluid-structure interaction

A Sixth Order Hybrid Finite Difference Scheme Based on the Minimized Dispersion and Controllable Dissipation Technique

A subcycling/non-subcycling time advancement scheme-based DLM immersed boundary method framework for solving single and multiphase fluid–structure interaction problems on dynamically adaptive grids

A New High-Order Unstructured Numerical Scheme for Large Eddy Simulation

Implementation of Density-Based Solver for All Speeds in the Framework of OpenFOAM

A high-order diffuse-interface method with TENO-THINC scheme for compressible multiphase flows

A low dissipation finite difference nested multi-resolution WENO scheme for Euler/Navier-Stokes equations

Direct Modeling for Computational Fluid Dynamics and the Construction of High-order Compact Scheme for Compressible Flow Simulations

High order method for solving unsteady incompressible Navier-Stokes equations

A low-dissipation shock-capturing framework with flexible nonlinear dissipation control

An efficient multi-resolution SPH framework for multi-phase fluid-structure interactions

Development of an efficient large time step unsteady solver for incompressible flows using the IDEAL algorithm in OpenFOAM

A Scalable Fully Implicit Framework for Reservoir Simulation on Parallel Computers

Development of an accurate central finite-difference scheme with a compact stencil for the simulation of unsteady incompressible flows on staggered orthogonal grids