High-order spectral volume scheme for multi-component flows using non-oscillatory kinetic flux

2017 ◽  
Vol 152 ◽  
pp. 120-133
Author(s):  
Na Liu ◽  
Xihua Xu ◽  
Yibing Chen
Keyword(s):  
High Order ◽  
Volume Scheme ◽  
Spectral Volume

Higher Resolution Hybrid-Upwind Spectral Finite-Volume Methods, for Flow in Porous and Fractured Media on Unstructured Grids

2021 ◽  
Author(s):  
Yawei Xie ◽  
Michael G. Edwards
Keyword(s):  
Finite Volume Method ◽  
Finite Volume ◽  
Unstructured Grids ◽  
Concentration Field ◽  
Convective Transport ◽  
High Order ◽  
Spectral Volume ◽  
Flux Approximation ◽  
Triangular Cell ◽  
Volume Method

Abstract A novel higher resolution spectral volume method coupled with a control-volume distributed multi-Point flux approximation (CVD-MPFA) is presented on unstructured triangular grids for subsurface reservoir simulation. The flow equations involve an essentially hyperbolic convection equation coupled with an elliptic pressure equation resulting from Darcy’s law together with mass conservation. The spectral volume (SV) method is a locally conservative, efficient high-order finite volume method for convective flow. In 2D geometry, the triangular cell is subdivided into sub-cells, and the average state variables in the sub-cells are used to reconstruct a high-order polynomial in the triangular cell. The focus here is on an efficient strategy for reconstruction of both a higher resolution approximation of the convective transport flux and Darcy-flux approximation on sub-cell interfaces, which is also coupled with a discrete fracture model. The strategy involves coupling of the SV method and reconstructed CVD-MPFA fluxes at the faces of the spectral volume, to obtain an efficient finer scale higher resolution finite-volume method which solves for both the saturation and pressure. A limiting procedure based on a Barth-Jespersen type limiter is used to prevent non-physical oscillations on unstructured grids. The fine scale saturation/concentration field is then updated via the reconstructed finite volume approximation over the sub-cell control-volumes. Performance comparisons are presented for two phase flow problems on 2D unstructured meshes including fractures. The results demonstrate that the spectral-volume method achieves further enhanced resolution of flow and fronts in addition to that of achieved by the standard higher resolution method over first order upwind, while improving upon efficiency.

A High Order Sharp-Interface Method with Local Time Stepping for Compressible Multiphase Flows

2011 ◽  
Vol 9 (1) ◽  
pp. 205-230 ◽  
Author(s):  
Angela Ferrari ◽  
Claus-Dieter Munz ◽  
Bernhard Weigand
Keyword(s):  
Finite Volume ◽  
High Order ◽  
Indicator Function ◽  
Finite Volume Scheme ◽  
Riemann Solver ◽  
Time Step ◽  
Material Interface ◽  
Volume Scheme ◽  
Time Stepping ◽  
Local Time Stepping

AbstractIn this paper, a new sharp-interface approach to simulate compressible multiphase flows is proposed. The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high order path-conservative discontinuous Galerkin finite element scheme to evolve an indicator function that tracks the material interface. At the interface our method applies ghost cells to compute the numerical flux, as the ghost fluid method. However, unlike the original ghost fluid scheme of Fedkiw et al. [15], the state of the ghost fluid is derived from an approximate-state Riemann solver, similar to the approach proposed in [25], but based on a much simpler formulation. Our formulation leads only to one single scalar nonlinear algebraic equation that has to be solved at the interface, instead of the system used in [25]. Away from the interface, we use the new general Osher-type flux recently proposed by Dumbser and Toro [13], which is a simple but complete Riemann solver, applicable to general hyperbolic conservation laws. The time integration is performed using a fully-discrete one-step scheme, based on the approaches recently proposed in [5,7]. This allows us to evolve the system also with time-accurate local time stepping. Due to the sub-cell resolution and the subsequent more restrictive time-step constraint of the DG scheme, a local evolution for the indicator function is applied, which is matched with the finite volume scheme for the solution of the Euler equations that runs with a larger time step. The use of a locally optimal time step avoids the introduction of excessive numerical diffusion in the finite volume scheme. Two different fluids have been used, namely an ideal gas and a weakly compressible fluid modeled by the Tait equation. Several tests have been computed to assess the accuracy and the performance of the new high order scheme. A verification of our algorithm has been carefully carried out using exact solutions as well as a comparison with other numerical reference solutions. The material interface is resolved sharply and accurately without spurious oscillations in the pressure field.

Sign in / Sign up

Export Citation Format

Share Document