Finite Volume Methods

Finite volume methods discretize conservation laws by dividing the domain into control volumes and updating the average value in each from the fluxes across its boundaries, conserving mass, momentum, and energy by construction.

Definition

A finite volume method is a discretization of a conservation law that stores the average of the solution over each control volume and evolves these averages by balancing the numerical fluxes through the volume boundaries, so that the discrete scheme is locally and globally conservative.

Scope

This topic covers the integral (conservation) form of PDEs, cell-averaged unknowns and numerical flux functions, the Godunov method and approximate Riemann solvers, high-resolution schemes with slope limiters that suppress spurious oscillations near shocks, and the role of finite volume methods in computational fluid dynamics.

Core questions

• Why does working from the integral conservation form make the method inherently conservative?
• How are numerical fluxes at cell interfaces defined, and what makes a flux consistent and stable?
• How do Godunov-type schemes and Riemann solvers capture discontinuities such as shocks?
• How do high-resolution methods avoid the oscillations of high-order schemes near discontinuities?

Key theories

Conservation and the numerical flux

By updating cell averages using a single numerical flux shared between adjacent cells, the method conserves the underlying quantity exactly; consistency of the flux with the true flux and a suitable stability condition yield convergence to weak solutions of the conservation law.

Godunov methods and Riemann solvers

Godunov's approach treats each cell interface as a local Riemann problem whose (exact or approximate) solution defines the flux, allowing the scheme to capture shocks and contact discontinuities sharply and correctly.

High-resolution schemes and limiters

To overcome the first-order accuracy of basic Godunov schemes without introducing spurious oscillations, high-resolution methods reconstruct higher-order interface states and apply slope or flux limiters that enforce a total-variation-diminishing property near discontinuities.

Mechanisms

Integrating the conservation law over a control volume converts spatial derivatives into surface fluxes, so the rate of change of a cell average equals the net flux through its faces. Because adjacent cells share each face flux, whatever leaves one cell enters its neighbour and the total quantity is conserved exactly. The interface flux is computed by solving, exactly or approximately, the Riemann problem posed by the differing states on either side; this captures wave structure and discontinuities. High-resolution variants first reconstruct a higher-order profile within each cell and limit it to prevent new extrema, achieving second-order accuracy in smooth regions while remaining oscillation-free at shocks.

Clinical relevance

Finite volume methods are the standard discretization in computational fluid dynamics and are central to simulating compressible and incompressible flows, aerodynamics, shock and detonation phenomena, shallow-water and atmospheric and oceanic flows, and porous-media and reservoir simulation, precisely because they respect physical conservation laws and handle discontinuities robustly.

History

The conservative, Riemann-problem-based approach originated with Godunov's 1959 scheme; the development of high-resolution, total-variation-diminishing methods and limiters by van Leer, Harten, and others in the 1970s and 1980s made finite volume methods the dominant framework for compressible flow and other hyperbolic conservation laws.

Key figures

• Sergei Godunov
• Peter Lax
• Bram van Leer
• Randall J. LeVeque
• Eleuterio Toro

Seminal works

• leveque2002
• toro2009

Frequently asked questions

Why are finite volume methods conservative?

They update cell averages using fluxes shared between neighbouring cells, so any quantity leaving one cell exactly enters the adjacent one. Summing over the whole domain, the total quantity changes only through the domain boundary, which mirrors the physical conservation law.

Why are limiters needed near shocks?

Plain high-order schemes produce spurious oscillations (overshoots and undershoots) near discontinuities. Slope and flux limiters detect steep gradients and locally reduce the reconstruction order, keeping the solution monotone and oscillation-free while preserving accuracy in smooth regions.