Quantitative evaluation of numerical integration schemes for Lagrangian particle dispersion models

3Citations
Citations of this article
6Readers
Mendeley users who have this article in their library.

Abstract

A rigorous methodology for the evaluation of integration schemes for Lagrangian particle dispersion models (LPDMs) is presented. A series of one-dimensional test problems are introduced, for which the Fokker-Planck equation is solved numerically using a finite-difference discretisation in physical space and a Hermite function expansion in velocity space. Numerical convergence errors in the Fokker-Planck equation solutions are shown to be much less than the statistical error associated with a practical-sized ensemble ( Combining double low line 106) of LPDM solutions; hence, the former can be used to validate the latter. The test problems are then used to evaluate commonly used LPDM integration schemes. The results allow for optimal time-step selection for each scheme, given a required level of accuracy. The following recommendations are made for use in operational models. First, if computational constraints require the use of moderate to long time steps, it is more accurate to solve the random displacement model approximation to the LPDM rather than use existing schemes designed for long time steps. Second, useful gains in numerical accuracy can be obtained, at moderate additional computational cost, by using the relatively simple "small-noise" scheme of Honeycutt.

Cite

CITATION STYLE

APA

Mohd Ramli, H., & Gavin Esler, J. (2016). Quantitative evaluation of numerical integration schemes for Lagrangian particle dispersion models. Geoscientific Model Development, 9(7), 2411–2457. https://doi.org/10.5194/gmd-9-2441-2016

Register to see more suggestions

Mendeley helps you to discover research relevant for your work.

Already have an account?

Save time finding and organizing research with Mendeley

Sign up for free