Accuracy aspects of conventional discretization methods for scalar transport with nonzero divergence velocity field arising from the energy balance equation

Research output: Contribution to journalArticleScientificpeer-review

1 Downloads (Pure)

Abstract

We are concerned with the numerical solution of a linear transport problem with nonzero divergence velocity field that originates from the spectral energy balance equation describing the evolution of wind waves and swells in coastal seas. The discretization error of the commonly used first-order upwind finite difference and first-order vertex-centered upwind finite volume schemes in one space dimension is analyzed. Smoothness of nondivergent velocity field plays a crucial role in this. No such analysis has been attempted to date for such problems. The two schemes studied differ in the manner in which they treat the scalar flux numerically. The finite difference variant is shock captured, whereas the vertex-centered finite volume approximation employs an arithmetic mean of the velocity and appears not to be flux conservative. The methods are subsequently extended to two dimensions on triangular meshes. Numerical experiments are provided to verify the convergence analysis. The main finding is that the finite difference scheme displays optimal rates of convergence and offers higher accuracy over the finite volume scheme, regardless the regularity of the velocity field. The latter scheme notably yields convergence rates of 0.5 and 0 in L2-norm and L-norm, respectively, when the velocity field is not smooth. A test case illustrating wave shoaling and refraction over submerged shoals is also presented and demonstrates the practical importance of flux conservation.

Original languageEnglish
Number of pages16
JournalInternational Journal for Numerical Methods in Fluids
DOIs
Publication statusPublished - 2020

Keywords

  • energy balance equation
  • finite difference method
  • finite volume method
  • flux conservation
  • nonzero divergence velocity
  • scalar transport equation

Fingerprint Dive into the research topics of 'Accuracy aspects of conventional discretization methods for scalar transport with nonzero divergence velocity field arising from the energy balance equation'. Together they form a unique fingerprint.

Cite this