High-order reconstruction schemes for the solution of hyperbolic conservation laws in orthogonal curvilinear coordinates are revised in the finite volume approach. The formulation employs a piecewise polynomial approximation to the zone-average values to reconstruct left and right interface states from within a computational zone to arbitrary order of accuracy by inverting a Vandermonde-like linear system of equations with spatially varying coefficients. The approach is general and can be used on uniform and non-uniform meshes although explicit expressions are derived for polynomials from second to fifth degree in cylindrical and spherical geometries with uniform grid spacing. It is shown that, in regions of large curvature, the resulting expressions differ considerably from their Cartesian counterparts and that the lack of such corrections can severely degrade the accuracy of the solution close to the coordinate origin. Limiting techniques and monotonicity constraints are revised for conventional reconstruction schemes, namely, the piecewise linear method (PLM), third-order weighted essentially non-oscillatory (WENO) scheme and the piecewise parabolic method (PPM). The performance of the improved reconstruction schemes is investigated in a number of selected numerical benchmarks involving the solution of both scalar and systems of nonlinear equations (such as the equations of gas dynamics and magnetohydrodynamics) in cylindrical and spherical geometries in one and two dimensions. Results confirm that the proposed approach yields considerably smaller errors, higher convergence rates and it avoid spurious numerical effects at a symmetry axis.

High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates

Mignone, A.
2014-01-01

Abstract

High-order reconstruction schemes for the solution of hyperbolic conservation laws in orthogonal curvilinear coordinates are revised in the finite volume approach. The formulation employs a piecewise polynomial approximation to the zone-average values to reconstruct left and right interface states from within a computational zone to arbitrary order of accuracy by inverting a Vandermonde-like linear system of equations with spatially varying coefficients. The approach is general and can be used on uniform and non-uniform meshes although explicit expressions are derived for polynomials from second to fifth degree in cylindrical and spherical geometries with uniform grid spacing. It is shown that, in regions of large curvature, the resulting expressions differ considerably from their Cartesian counterparts and that the lack of such corrections can severely degrade the accuracy of the solution close to the coordinate origin. Limiting techniques and monotonicity constraints are revised for conventional reconstruction schemes, namely, the piecewise linear method (PLM), third-order weighted essentially non-oscillatory (WENO) scheme and the piecewise parabolic method (PPM). The performance of the improved reconstruction schemes is investigated in a number of selected numerical benchmarks involving the solution of both scalar and systems of nonlinear equations (such as the equations of gas dynamics and magnetohydrodynamics) in cylindrical and spherical geometries in one and two dimensions. Results confirm that the proposed approach yields considerably smaller errors, higher convergence rates and it avoid spurious numerical effects at a symmetry axis.
2014
270
784
814
https://www.sciencedirect.com/science/article/pii/S0021999114002538?via=ihub
Curvilinear geometry; Finite volume; Hydrodynamics; Magnetohydrodynamics (MHD); Methods: numerical; Reconstruction methods; Physics and Astronomy (miscellaneous); Computer Science Applications1707 Computer Vision and Pattern Recognition
Mignone, A.*
File in questo prodotto:
File Dimensione Formato  
2014.Mignone.JCP.pdf

Accesso riservato

Tipo di file: PDF EDITORIALE
Dimensione 1.15 MB
Formato Adobe PDF
1.15 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
1404.0537.pdf

Accesso aperto

Tipo di file: POSTPRINT (VERSIONE FINALE DELL’AUTORE)
Dimensione 842.33 kB
Formato Adobe PDF
842.33 kB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/2318/1693387
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 97
  • ???jsp.display-item.citation.isi??? 88
social impact