TY - JOUR

T1 - A high-order finite volume method for solving one-dimensional convection and diffusion equations

AU - Pan, Dartzi

N1 - Publisher Copyright:
© 2017 Taylor & Francis.
Copyright:
Copyright 2017 Elsevier B.V., All rights reserved.

PY - 2017/6/3

Y1 - 2017/6/3

N2 - A spatially high-order finite volume method for solving convection and diffusion equations is developed and tested in this work. The method performs a high-order piecewise polynomial reconstruction of the local flow field based on the relationship between Taylor’s series expansion and the volume-averaged flow quantities. A 5 × 5 matrix inversion for each cell is done to compute the cell-center variables and derivatives up to fourth order. While a fixed symmetric grid stencil is maintained in smooth flow regions, a detector for large change in linear data slopes is developed to trigger the use of ENO stencil around flow discontinuities. Regular time integration scheme such as the four-stage Runge–Kutta method or the Euler implicit method is used for time integration. The present finite volume method is shown to be spatially fifth-order accurate for the linear convection equation, fourth-order accurate for the linear diffusion equation, and fourth-order accurate for the linear convection–diffusion equation. The shocks captured in solving the inviscid Burger’s equation are sharp and oscillation free. For the system of Euler equations, a characteristic limiter is further developed to limit the growth of total variation of the solution. Test examples solving shock-tube problems and the interactions of two blast waves show that various flow discontinuities are captured sharply without spurious oscillations.

AB - A spatially high-order finite volume method for solving convection and diffusion equations is developed and tested in this work. The method performs a high-order piecewise polynomial reconstruction of the local flow field based on the relationship between Taylor’s series expansion and the volume-averaged flow quantities. A 5 × 5 matrix inversion for each cell is done to compute the cell-center variables and derivatives up to fourth order. While a fixed symmetric grid stencil is maintained in smooth flow regions, a detector for large change in linear data slopes is developed to trigger the use of ENO stencil around flow discontinuities. Regular time integration scheme such as the four-stage Runge–Kutta method or the Euler implicit method is used for time integration. The present finite volume method is shown to be spatially fifth-order accurate for the linear convection equation, fourth-order accurate for the linear diffusion equation, and fourth-order accurate for the linear convection–diffusion equation. The shocks captured in solving the inviscid Burger’s equation are sharp and oscillation free. For the system of Euler equations, a characteristic limiter is further developed to limit the growth of total variation of the solution. Test examples solving shock-tube problems and the interactions of two blast waves show that various flow discontinuities are captured sharply without spurious oscillations.

UR - http://www.scopus.com/inward/record.url?scp=85020457346&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85020457346&partnerID=8YFLogxK

U2 - 10.1080/10407790.2017.1326769

DO - 10.1080/10407790.2017.1326769

M3 - Article

AN - SCOPUS:85020457346

VL - 71

SP - 533

EP - 548

JO - Numerical Heat Transfer, Part B: Fundamentals

JF - Numerical Heat Transfer, Part B: Fundamentals

SN - 1040-7790

IS - 6

ER -