A Posteriori Stabilized Sixth-Order Finite Volume Scheme with Adaptive Stencil Construction: Basics for the 1D Steady-State Hyperbolic Equations

doi:10.1007/s42967-021-00140-7

Communications on Applied Mathematics and Computation ›› 2023, Vol. 5 ›› Issue (2): 751-775.doi: 10.1007/s42967-021-00140-7

• ORIGINAL PAPERS • 上一篇下一篇

A Posteriori Stabilized Sixth-Order Finite Volume Scheme with Adaptive Stencil Construction: Basics for the 1D Steady-State Hyperbolic Equations

Gaspar J. Machado¹, Stéphane Clain¹, Raphaël Loubère²

1 Centre of Physics and Department of Mathematics, University of Minho, Campus of Azurém, 4800-058 Guimarães, Portugal;
2 CNRS and Institut de Mathématiques de Bordeaux(IMB), Université de Bordeaux, Talence, France

收稿日期:2020-09-21 修回日期:2021-04-13 出版日期:2023-06-20 发布日期:2023-05-26
通讯作者: Raphaël Loubère, raphael.loubere@u-bordeaux.fr;Gaspar J. Machado, gjm@math.uminho.pt;Stéphane Clain, clain@math.uminho.pt E-mail:raphael.loubere@u-bordeaux.fr;gjm@math.uminho.pt;clain@math.uminho.pt
基金资助:
S. Clain and G.J. Machado acknowledge the financial support by FEDER-Fundo Europeu de Desenvolvimento Regional, through COMPETE 2020-Programa Operational Fatores de Competitividade, and the National Funds through FCT-Fundação para a Ciência e a Tecnologia, project no. UID/FIS/04650/2019. S. Clain and G.J. Machado acknowledge the financial support by FEDER-Fundo Europeu de Desenvolvimento Regional, through COMPETE 2020-Programa Operacional Fatores de Competitividade, and the National Funds through FCT-Fundação para a Ciência e a Tecnologia, project no. POCI-01-0145-FEDER-028118. The material of this research has been partly built and discussed during the SHARK workshops taking place in Póvoa de Varzim, Portugal, http://www.SHARK-FV.eu/.

A Posteriori Stabilized Sixth-Order Finite Volume Scheme with Adaptive Stencil Construction: Basics for the 1D Steady-State Hyperbolic Equations

Gaspar J. Machado¹, Stéphane Clain¹, Raphaël Loubère²

1 Centre of Physics and Department of Mathematics, University of Minho, Campus of Azurém, 4800-058 Guimarães, Portugal;
2 CNRS and Institut de Mathématiques de Bordeaux(IMB), Université de Bordeaux, Talence, France

Received:2020-09-21 Revised:2021-04-13 Online:2023-06-20 Published:2023-05-26
Contact: Raphaël Loubère, raphael.loubere@u-bordeaux.fr;Gaspar J. Machado, gjm@math.uminho.pt;Stéphane Clain, clain@math.uminho.pt E-mail:raphael.loubere@u-bordeaux.fr;gjm@math.uminho.pt;clain@math.uminho.pt
Supported by:
S. Clain and G.J. Machado acknowledge the financial support by FEDER-Fundo Europeu de Desenvolvimento Regional, through COMPETE 2020-Programa Operational Fatores de Competitividade, and the National Funds through FCT-Fundação para a Ciência e a Tecnologia, project no. UID/FIS/04650/2019. S. Clain and G.J. Machado acknowledge the financial support by FEDER-Fundo Europeu de Desenvolvimento Regional, through COMPETE 2020-Programa Operacional Fatores de Competitividade, and the National Funds through FCT-Fundação para a Ciência e a Tecnologia, project no. POCI-01-0145-FEDER-028118. The material of this research has been partly built and discussed during the SHARK workshops taking place in Póvoa de Varzim, Portugal, http://www.SHARK-FV.eu/.

摘要/Abstract

摘要： We propose an adaptive stencil construction for high-order accurate finite volume schemes a posteriori stabilized devoted to solve one-dimensional steady-state hyperbolic equations. High accuracy (up to the sixth-order presently) is achieved, thanks to polynomial reconstructions while stability is provided with an a posteriori MOOD method which controls the cell polynomial degree for eliminating non-physical oscillations in the vicinity of discontinuities. We supplemented this scheme with a stencil construction allowing to reduce even further the numerical dissipation. The stencil is shifted away from troubles (shocks, discontinuities, etc.) leading to less oscillating polynomial reconstructions. Experimented on linear, Bürgers', and Euler equations, we demonstrate that the adaptive stencil technique manages to retrieve smooth solutions with optimal order of accuracy but also irregular ones without spurious oscillations. Moreover, we numerically show that the approach allows to reduce the dissipation still maintaining the essentially non-oscillatory behavior.

关键词: Finite volume, MOOD, Adaptive stencil, Steady-state solution, Euler equations, High order

Abstract: We propose an adaptive stencil construction for high-order accurate finite volume schemes a posteriori stabilized devoted to solve one-dimensional steady-state hyperbolic equations. High accuracy (up to the sixth-order presently) is achieved, thanks to polynomial reconstructions while stability is provided with an a posteriori MOOD method which controls the cell polynomial degree for eliminating non-physical oscillations in the vicinity of discontinuities. We supplemented this scheme with a stencil construction allowing to reduce even further the numerical dissipation. The stencil is shifted away from troubles (shocks, discontinuities, etc.) leading to less oscillating polynomial reconstructions. Experimented on linear, Bürgers', and Euler equations, we demonstrate that the adaptive stencil technique manages to retrieve smooth solutions with optimal order of accuracy but also irregular ones without spurious oscillations. Moreover, we numerically show that the approach allows to reduce the dissipation still maintaining the essentially non-oscillatory behavior.

Key words: Finite volume, MOOD, Adaptive stencil, Steady-state solution, Euler equations, High order

中图分类号:

Gaspar J. Machado, Stéphane Clain, Raphaël Loubère. A Posteriori Stabilized Sixth-Order Finite Volume Scheme with Adaptive Stencil Construction: Basics for the 1D Steady-State Hyperbolic Equations[J]. Communications on Applied Mathematics and Computation, 2023, 5(2): 751-775.

参考文献

1. Beam, R.M., Warming, R.F.: An implicit finite-difference algorithm for hyperbolic systems in conservation law form. J. Comput. Phys. 22(1), 87–110 (1976)
2. Beam, R.M., Warming, R.F.: Upwind second-order difference schemes and applications in aerodynamic flows. AIAA J. 14(9), 1241–1249 (1976)
3. Clain, S., Diot, S., Loubère, R.: A high-order finite volume method for systems of conservation laws— multi-dimensional optimal order detection (MOOD). J. Comput. Phys. 230(10), 4028–4050 (2011)
4. Clain, S., Loubère, R., Machado, G.J.: A posteriori stabilized sixth-order finite volume scheme for one-dimensional steady-state hyperbolic equations. Adv. Comput. Math. 44(2), 571–607 (2018)
5. Diot, S., Clain, S., Loubère, R.: Improved detection criteria for the multi-dimensional optimal order detection (MOOD) on unstructured meshes with very high-order polynomials. Comput. Fluids 64, 43–63 (2012)
6. Diot, S., Loubère, R., Clain, S.: The MOOD method in the three-dimensional case: very-high-order finite volume method for hyperbolic systems. Int. J. Numer. Methods Fluids 73, 362–392 (2013)
7. Einfeldt, B., Munz, C.D., Roe, P.L., Sjögreen, B.: On Godunov-type methods near low densities. J. Comput. Phys. 92, 273–295 (1991)
8. Harten, A., Engquist, B., Osher, S., Chakravarthy, S.R.: Uniformly high order accurate essentially nonoscillatory schemes III. J. Comput. Phys. 71, 231–303 (1987)
9. Harten, A., Lax, P.D., van Leer, B.: On upstream differencing and Godunov-type schemes for hyperbolic conservation laws. SIAM Rev. 25(1), 35–61 (1983)
10. Jiang, G.-S., Shu, C.-W.: Efficient implementation of weighted ENO schemes. J. Comput. Phys. 126, 202–228 (1996)
11. Liu, X.D., Osher, S., Chan, T.: Weighted essentially non-oscillatory schemes. J. Comput. Phys. 115, 200–212 (1994)
12. Pulliam, T.H.: Early development of implicit methods for computational fluid dynamics at NASA Ames. Comput. Fluids 38, 491–495 (2009)
13. Pulliam, T.H., Lomax, H.: Simulation of three-dimensional compressible viscous flow on the Illiac IV computer. AIAA paper 79-0206 (1979)
14. Pulliam, T.H., Steger, J.L.: Implicit finite-difference simulations of three dimensional compressible flow. AIAA J. 18(2), 159 (1980)
15. Salas, M.D., Abarbanel, S., Gottlieb, D.: Multiple steady states for characteristic initial value problems. Appl. Numer. Math. 2, 193–210 (1986)
16. Shu, C.-W.: Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws. In: Lecture Notes in Mathematics, pp. 325–432. Springer, Berlin (1998)
17. Steger, J.L.: Implicit finite difference simulation of flow about arbitrary geometries with application to airfoils. AIAA J. 16(7), 679–686 (1978)
18. Venkatakrishnan, V., Mavriplis, D.: Implicit solvers for unstructured meshes. J. Comput. Phys. 105(1), 83–91 (1993)

[1]	Johannes Markert, Gregor Gassner, Stefanie Walch. A Sub-element Adaptive Shock Capturing Approach for Discontinuous Galerkin Methods[J]. Communications on Applied Mathematics and Computation, 2023, 5(2): 679-721.
[2]	J. C. González-Aguirre, J. A. González-Vázquez, J. Alavez-Ramírez, R. Silva, M. E. Vázquez-Cendón. Numerical Simulation of Bed Load and Suspended Load Sediment Transport Using Well-Balanced Numerical Schemes[J]. Communications on Applied Mathematics and Computation, 2023, 5(2): 885-922.
[3]	Paula Luna, Arturo Hidalgo. Numerical Approach of a Coupled Pressure-Saturation Model Describing Oil-Water Flow in Porous Media[J]. Communications on Applied Mathematics and Computation, 2023, 5(2): 946-964.
[4]	M. Semplice, E. Travaglia, G. Puppo. One- and Multi-dimensional CWENOZ Reconstructions for Implementing Boundary Conditions Without Ghost Cells[J]. Communications on Applied Mathematics and Computation, 2023, 5(1): 143-169.
[5]	Bao-Shan Wang, Wai Sun Don, Alexander Kurganov, Yongle Liu. Fifth-Order A-WENO Schemes Based on the Adaptive Diffusion Central-Upwind Rankine-Hugoniot Fluxes[J]. Communications on Applied Mathematics and Computation, 2023, 5(1): 295-314.
[6]	Xucheng Meng, Yaguang Gu, Guanghui Hu. A Fourth-Order Unstructured NURBS-Enhanced Finite Volume WENO Scheme for Steady Euler Equations in Curved Geometries[J]. Communications on Applied Mathematics and Computation, 2023, 5(1): 315-342.
[7]	G. Puppo, M. Semplice, G. Visconti. Quinpi: Integrating Conservation Laws with CWENO Implicit Methods[J]. Communications on Applied Mathematics and Computation, 2023, 5(1): 343-369.
[8]	R. Abgrall. A Combination of Residual Distribution and the Active Flux Formulations or a New Class of Schemes That Can Combine Several Writings of the Same Hyperbolic Problem: Application to the 1D Euler Equations[J]. Communications on Applied Mathematics and Computation, 2023, 5(1): 370-402.
[9]	Zepeng Liu, Yan Jiang, Mengping Zhang, Qingyuan Liu. High Order Finite Difference WENO Methods for Shallow Water Equations on Curvilinear Meshes[J]. Communications on Applied Mathematics and Computation, 2023, 5(1): 485-528.
[10]	Hendrik Ranocha, Lisandro Dalcin, Matteo Parsani, David I. Ketcheson. Optimized Runge-Kutta Methods with Automatic Step Size Control for Compressible Computational Fluid Dynamics[J]. Communications on Applied Mathematics and Computation, 2022, 4(4): 1191-1228.
[11]	Hendrik Ranocha, Gregor J. Gassner. Preventing Pressure Oscillations Does Not Fix Local Linear Stability Issues of Entropy-Based Split-Form High-Order Schemes[J]. Communications on Applied Mathematics and Computation, 2022, 4(3): 880-903.
[12]	Andreas Dedner, Robert Klöfkorn. Extendible and Efcient Python Framework for Solving Evolution Equations with Stabilized Discontinuous Galerkin Methods[J]. Communications on Applied Mathematics and Computation, 2022, 4(2): 657-696.
[13]	Yu Pan, Zhen-Guo Yan, Joaquim Peiró, Spencer J. Sherwin. Development of a Balanced Adaptive Time-Stepping Strategy Based on an Implicit JFNK-DG Compressible Flow Solver[J]. Communications on Applied Mathematics and Computation, 2022, 4(2): 728-757.
[14]	Vít Dolejší, Filip Roskovec. Goal-Oriented Anisotropic hp-Adaptive Discontinuous Galerkin Method for the Euler Equations[J]. Communications on Applied Mathematics and Computation, 2022, 4(1): 143-179.
[15]	Giacomo Albi, Lorenzo Pareschi. High Order Semi-implicit Multistep Methods for Time-Dependent Partial Diferential Equations[J]. Communications on Applied Mathematics and Computation, 2021, 3(4): 701-718.

A Posteriori Stabilized Sixth-Order Finite Volume Scheme with Adaptive Stencil Construction: Basics for the 1D Steady-State Hyperbolic Equations

A Posteriori Stabilized Sixth-Order Finite Volume Scheme with Adaptive Stencil Construction: Basics for the 1D Steady-State Hyperbolic Equations

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价