Journal
GEOCHEMISTRY GEOPHYSICS GEOSYSTEMS
Volume 16, Issue 6, Pages 2015-2023Publisher
AMER GEOPHYSICAL UNION
DOI: 10.1002/2015GC005824
Keywords
conservative interpolation; divergence free; PIC method; viscosity contrast
Categories
Funding
- EU [264517]
- NERC [NE/J008028/1]
- European Research Council (ERC) [StG 279828]
- NERC [NE/J007854/1, NE/J008028/1] Funding Source: UKRI
- Natural Environment Research Council [NE/J007854/1, NE/J008028/1] Funding Source: researchfish
Ask authors/readers for more resources
The particle-in-cell method is generally considered a flexible and robust method to model the geodynamic problems with chemical heterogeneity. However, velocity interpolation from grid points to particle locations is often performed without considering the divergence of the velocity field, which can lead to significant particle dispersion or clustering if those particles move through regions of strong velocity gradients. This may ultimately result in cells void of particles, which, if left untreated, may, in turn, lead to numerical inaccuracies. Here we apply a two-dimensional conservative velocity interpolation (CVI) scheme to steady state and time-dependent flow fields with strong velocity gradients (e.g., due to large local viscosity variation) and derive and apply the three-dimensional equivalent. We show that the introduction of CVI significantly reduces the dispersion and clustering of particles in both steady state and time-dependent flow problems and maintains a locally steady number of particles, without the need for ad hoc remedies such as very high initial particle densities or reseeding during the calculation. We illustrate that this method provides a significant improvement to particle distributions in common geodynamic modeling problems such as subduction zones or lithosphere-asthenosphere boundary dynamics.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available