Skip to main content
Publications | Persons | Institutes | Projects
[ report an error in this record ]basket (0): add | show Print this page

Assessing Lagrangian schemes for simulating diffusion on non-flat isopycnal surfaces
Shah, S.H.A.M.; Heemink, A.W.; Deleersnijder, E. (2011). Assessing Lagrangian schemes for simulating diffusion on non-flat isopycnal surfaces. Ocean Modelling 39(3-4): 351-361. dx.doi.org/10.1016/j.ocemod.2011.05.008
In: Ocean Modelling. Elsevier: Oxford. ISSN 1463-5003; e-ISSN 1463-5011, more
Peer reviewed article  

Available in  Authors 

Keyword
    Marine/Coastal
Author keywords
    Lagrangian modelling; Non-flat isopycnal surfaces; Large-scale model;

Authors  Top 
  • Shah, S.H.A.M.
  • Heemink, A.W.
  • Deleersnijder, E., more

Abstract
    In large-scale ocean flows diffusion mostly occurs along the density surfaces and its representation resorts to the Redi isopycnal diffusivity tensor containing off-diagonal terms. This study focuses on the Lagrangian/particle framework for simulating such diffusive processes. A two-dimensional idealised test case for purely isopycnal diffusion on non-flat isopycnal surfaces is considered. Implementation of the higher order strong Euler, Milstein and order 1.5 Taylor schemes on our idealised test case shows that the higher order strong schemes produce the better pathwise approximations. The effective spurious diapycnal diffusivity is measured for each Lagrangian scheme under consideration. The propensity of the particles to move away from the isopycnal surface on which they were released is also measured. This shows that for non-flat isopycnals the order of convergence of the Euler scheme is not sufficient to achieve the desired accuracy. However, the Milstein scheme seems to be a good choice to achieve in an efficient way a fairly accurate result.

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