Skip to main content

IMIS

A new integrated search interface will become available in the next phase of marineinfo.org.
For the time being, please use IMIS to search available data

 

[ report an error in this record ]basket (0): add | show Print this page

Why the Euler scheme in particle tracking is not enough: the shallow-sea pycnocline test case
Gräwe, U.; Deleersnijder, E.; Shah, S.H.A.M.; Heemink, A.W. (2012). Why the Euler scheme in particle tracking is not enough: the shallow-sea pycnocline test case. Ocean Dynamics 62(4): 501-514. http://dx.doi.org/10.1007/s10236-012-0523-y
In: Ocean Dynamics. Springer-Verlag: Berlin; Heidelberg; New York. ISSN 1616-7341; e-ISSN 1616-7228, more
Peer reviewed article  

Available in  Authors 

Keyword
    Marine/Coastal
Author keywords
    Particle tracking; Milstein scheme; Pycnocline; Diffusion; Residence time

Authors  Top 
  • Gräwe, U.
  • Deleersnijder, E., more
  • Shah, S.H.A.M.
  • Heemink, A.W.

Abstract
    During the last decades, the Euler scheme was the common “workhorse” in particle tracking, although it is the lowest-order approximation of the underlying stochastic differential equation. To convince the modelling community of the need for better methods, we have constructed a new test case that will show the shortcomings of the Euler scheme. We use an idealised shallow-water diffusivity profile that mimics the presence of a sharp pycnocline and thus a quasi-impermeable barrier to vertical diffusion. In this context, we study the transport of passive particles with or without negative buoyancy. A semi-analytic solutions is used to assess the performance of various numerical particle-tracking schemes (first- and second-order accuracy), to treat the variations in the diffusivity profile properly. We show that the commonly used Euler scheme exhibits a poor performance and that widely used particle-tracking codes shall be updated to either the Milstein scheme or second-order schemes. It is further seen that the order of convergence is not the only relevant factor, the absolute value of the error also is.

All data in the Integrated Marine Information System (IMIS) is subject to the VLIZ privacy policy Top | Authors