An Essentially Nonoscillatory Spectral Deferred Correction Method for Hyperbolic Problems

2016 ◽  
Vol 13 (03) ◽  
pp. 1650017 ◽  
Author(s):  
Samet Y. Kadioglu

We present a computational method based on the spectral deferred corrections (SDC) time integration technique and the essentially nonoscillatory (ENO) finite volume method for hyperbolic problems. The SDC technique is used to advance the solutions in time with high-order of accuracy. The ENO method is used to define high-order cell edge quantities that are then used to evaluate numerical fluxes. The coupling of the SDC method with a high-order finite volume method (piece-wise parabolic method (PPM)) is first carried out by Layton et al. [J. Comput. Phys. 194(2) (2004) 697]. Issues about this approach have been addressed and some improvements have been added to it in Kadioglu et al. [J. Comput. Math. 1(4) (2012) 303]. Here, we investigate the implications when the PPM method is replaced with the well-known ENO method. We note that the SDC-PPM method is fourth-order accurate in time and space. Therefore, we kept the order of accuracy of the ENO procedure as fourth-order in order to be able to make a consistent comparison between the two approaches (SDC-ENO versus SDC-PPM). We have tested the new SDC-ENO technique by solving smooth and nonsmooth hyperbolic problems. Our numerical results indicate that the fourth-order of accuracy in both space and time has been achieved for smooth problems. On the other hand, the new method performs very well when it is applied to nonlinear problems that involve discontinuous solutions. In other words, we have obtained highly resolved discontinuous solutions with essentially no-oscillations at or around the discontinuities.

2021 ◽  
Author(s):  
Yawei Xie ◽  
Michael G. Edwards

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.


2015 ◽  
Vol 123 ◽  
pp. 202-217 ◽  
Author(s):  
Stephen M. Guzik ◽  
Xinfeng Gao ◽  
Landon D. Owen ◽  
Peter McCorquodale ◽  
Phillip Colella

Sign in / Sign up

Export Citation Format

Share Document