Over het archief
Het OWA, het open archief van het Waterbouwkundig Laboratorium heeft tot doel alle vrij toegankelijke onderzoeksresultaten van dit instituut in digitale vorm aan te bieden. Op die manier wil het de zichtbaarheid, verspreiding en gebruik van deze onderzoeksresultaten, alsook de wetenschappelijke communicatie maximaal bevorderen.
Dit archief wordt uitgebouwd en beheerd volgens de principes van de Open Access Movement, en het daaruit ontstane Open Archives Initiative.
Basisinformatie over ‘Open Access to scholarly information'.
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, meer
|  |
| Trefwoord |
|
| Author keywords |
Particle tracking; Milstein scheme; Pycnocline; Diffusion; Residence time |
| Auteurs | | Top |
- Gräwe, U.
- Deleersnijder, E., meer
- 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. |
IMIS is ontwikkeld en wordt gehost door het VLIZ.
