Open AccessModeling advection in geophysical flows with particle level sets
Author(s) -
Samuel H.,
Evonuk M.
Publication year - 2010
Publication title -
geochemistry, geophysics, geosystems
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.928
H-Index - 136
ISSN - 1525-2027
DOI - 10.1029/2010gc003081
Subject(s) - classification of discontinuities , advection , eulerian path , level set method , geology , grid , geophysics , particle (ecology) , set (abstract data type) , scale (ratio) , mechanics , lagrangian , mathematics , computer science , statistical physics , physics , geodesy , mathematical analysis , oceanography , segmentation , artificial intelligence , image segmentation , thermodynamics , programming language , quantum mechanics
We have applied, for the first time in geodynamical flows, the Particle Level Set method for advecting compositional fields with sharp discontinuities. This robust and efficient Eulerian‐Lagrangian technique is based on the concept of Implicit Surfaces, which allows the use of high order accurate numerical schemes in the vicinity of discontinuities. We have tested the Particle Level Set method against the robust and popular Tracer‐in‐Cell method on well‐known 2D thermochemical benchmarks and typical 3D convective flows. The use of Lagrangian tracers in the Particle Level Set method yields accurate solutions of purely advective transport, where sub‐grid scale features can be resolved. In every case we ran we found that the Particle Level Set method accuracy equals or is better than the popular Tracer‐in‐Cell method, and can lead to significantly smaller computational cost, in particular in three‐dimensional flows, where the reduction of computational time for modeling advection processes is most needed.