@article{fdi:010078177,
  title = {3-{D} numerical modelling of crustal polydiapirs with volume-of-fluid methods},
  author = {{L}ouis-{N}apoleon, {A}. and {G}erbault, {M}uriel and {B}onometti, {T}. and {T}hieulot, {C}. and {M}artin, {R}. and {V}anderhaeghe, {O}.},
  editor = {},
  language = {{ENG}},
  abstract = {{G}ravitational instabilities exert a crucial role on the {E}arth dynamics and in particular on its differentiation. {T}he {E}arth's crust can be considered as a multilayered fluid with different densities and viscosities, which may become unstable in particular with variations in temperature. {W}ith the specific aim to quantify crustal scale polydiapiric instabilities, we test here two codes, {JADIM} and {O}pen{FOAM}, which use a volume-of-fluid ({VOF}) method without interface reconstruction, and compare them with the geodynamics community code {ASPECT}, which uses a tracking algorithm based on compositional fields. {T}he {VOF} method is well-known to preserve strongly deforming interfaces. {B}oth {JADIM} and {O}pen{FOAM} are first tested against documented two and three-layer {R}ayleigh-{T}aylor instability configurations in 2-{D} and 3-{D}. 2-{D} and 3-{D} results show diapiric growth rates that fit the analytical theory and are found to be slightly more accurate than those obtained with {ASPECT}. {W}e subsequently compare the results from {VOF} simulations with previously published {R}ayleigh-{B}enard analogue and numerical experiments. {W}e show that the {VOF} method is a robust method adapted to the study of diapirism and convection in the {E}arth's crust, although it is not computationally as fast as {ASPECT}. {O}pen{FOAM} is found to run faster than, and conserve mass as well as {JADIM}. {F}inally, we provide a preliminary application to the polydiapiric dynamics of the orogenic crust of {N}axos {I}sland ({G}reece) at about 16 {M}yr, and propose a two-stages scenario of convection and diapirism. {T}he timing and dimensions of the modelled gravitational instabilities not only corroborate previous estimates of timing and dimensions associated to the dynamics of this hot crustal domain, but also bring preliminary insight on its rheological and tectonic contexts.},
  keywords = {{G}eomechanics ; {N}umerical modelling ; {C}rustal structure ; {D}iapirism ; {D}ynamics: gravity and tectonics ; {R}heology: crust and lithosphere},
  booktitle = {},
  journal = {{G}eophysical {J}ournal {I}nternational},
  volume = {222},
  numero = {1},
  pages = {474--506},
  ISSN = {0956-540{X}},
  year = {2020},
  DOI = {10.1093/gji/ggaa141},
  URL = {https://www.documentation.ird.fr/hor/fdi:010078177},
}