Premium
Accuracy aspects of conventional discretization methods for scalar transport with nonzero divergence velocity field arising from the energy balance equation
Author(s) -
Zijlema Marcel
Publication year - 2021
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.4935
Subject(s) - finite volume method , mathematics , mathematical analysis , discretization , finite difference , upwind scheme , vector field , finite difference method , rate of convergence , geometry , physics , mechanics , channel (broadcasting) , electrical engineering , engineering
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 L 2 ‐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.