Approaches to the numerical estimates of grid convergence of NSE in the presence of singularities

Chenguang Zhang, Krishnaswamy Nandakumar*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review


Evaluating the order of accuracy (order) is an integral part of the development and application of numerical algorithms. Apart from theoretical functional analysis to place bounds on error estimates, numerical experiments are often essential for nonlinear problems to validate the estimates in a reliable answer. The common workflow is to apply the algorithm using successively finer temporal/spatial grid resolutions δi, measure the error ∈i in each solution against the exact solution, the order is then obtained as the slope of the line that fits (log ∈i, log δi). We show that if the problem has singularities like divergence to infinity or discontinuous jump, this common workflow underestimates the order if solution at regions around the singularity is used. Several numerical examples with different levels of complexity are explored. A simple one-dimensional theoretical model shows it is impossible to numerically evaluate the order close to singularity on uniform grids.

Original languageEnglish
Pages (from-to)281-288
Number of pages8
JournalInternational Journal of Nonlinear Sciences and Numerical Simulation
Issue number3-4
StatePublished - 26 Jun 2018
Externally publishedYes


  • finite volume method
  • impact of singularity
  • numerical error estimates


Dive into the research topics of 'Approaches to the numerical estimates of grid convergence of NSE in the presence of singularities'. Together they form a unique fingerprint.

Cite this