It takes three to tango: 1. Simulating buoyancy-driven flow in the presence of large viscosity contrasts

Journal of Geophysical Research: Solid Earth

Volumn 115, Issue 7, 2010, Pages

  Suckale, Jenny ^a   Nave, Jean Christophe ^a   Hager, Bradford H ^a  

^a MASSACHUSETTS INSTITUTE OF TECHNOLOGY   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

ANALYTICAL METHOD; BUOYANCY; FLUID DYNAMICS; GEODYNAMICS; MULTIPHASE FLOW; NAVIER-STOKES EQUATIONS; NUMERICAL MODEL; VISCOELASTIC FLUID; GEOPHYSICS; STOKES FORMULA; VISCOSITY;

EID: 77955213944     PISSN: 21699313     EISSN: 21699356     Source Type: Journal    
DOI: 10.1029/2009JB006916     Document Type: Article

References (66)

1. Adalsteinsson, D., and J. A. Sethian (1999), The fast construction of extension velocities in level set methods, J. Comput. Phys., 148(1), 2-22, doi:10.1006/jcph.1998.6090.
2. Alley, K. M., and E. M. Parmentier (1998), Numerical experiments on thermal convection in a chemically stratified viscous fluid heated from below: Implications for a model of lunar evolution, Phys. Earth Planet. Inter., 108(1), 15-32, doi:10.1016/S0031-9201(98)00096-X. (Pubitemid 28441603)
3. Bergantz, G. W., and J. Ni (1999), A numerical study of sedimentation by dripping instabilities in viscous fluids, Int. J. Multiphase Flow, 25(2), 307-320, doi:10.1016/S0301-9322(98)00050-0. (Pubitemid 29347938)
4. Bourgouin, L., H. B. Mülhaus, A. J. Hale, and A. Arsac (2007), Studying the influence of a solid shell on lava dome growth and evolution using the level set method, Geophys. J. Int., 170(3), 1431-1438, doi:10.1111/ j.1365-246X.2007.03471.x.
5. Chern, I. L., and Y. C. Shu (2007), A coupling interface method for elliptic interface problems, J. Comput. Phys., 225(2), 2138-2174, doi:10.1016/j. jcp.2007.03.012.
6. Chopp, D. L. (2001), Some improvements of the fast marching method, SIAM J. Sci. Stat. Comput., 23(1), 230-244, doi:10.1137/S106482750037617X.
7. Chopp, D. L. (2009), Another look at velocity extensions in the level set method, SIAM J. Sci. Comput., 31(5), 3255-3273, doi:10.1137/ 070686329.
8. Christensen, U. R. (1992), An Eulerian technique for thermomechanical modeling of lithospheric extension, J. Geophys. Res., 97(B2), 2015-2036, doi:10.1029/91JB02642.
9. Conrad, C. P., and P. Molnar (1997), The growth of Rayleigh-Taylor-type instabilities in the lithosphere for various rheological and density structures, Geophys. J. Int., 129(1), 95-112, doi:10.1111/j.1365-246X.1997. tb00939.x.
10. Deubelbeiss, Y., and B. J. P. Kaus (2008), Comparison of Eulerian and Lagrangian numerical techniques for the Stokes equations in the presence of strongly varying viscosity, Phys. Earth Planet. Inter., 171(1-4), 92-111, doi:10.1016/j.pepi.2008.06.023.
11. Elkins-Tanton, L. T. (2007), Continental magmatism, volatile recycling, and a heterogeneous mantle caused by lithospheric gravitational instabilities, J. Geophys. Res., 112, B03405, doi:10.1029/2005JB004072. (Pubitemid 47386355)
12. Enright, D., R. Fedkiw, J. Ferziger, and I. Mitchell (2002), A hybrid particle level set method for improved interface capturing, J. Comput. Phys., 183(1), 83-116, doi:10.1006/jcph.2002.7166. (Pubitemid 36100551)
13. Enright, D., F. Losasso, and R. Fedkiw (2005), A fast and accurate semi-Lagrangian particle level set method, Comput. Struct., 83(6-7), 479-490, doi:10.1016/j.compstruc.2004.04.024. (Pubitemid 40219668)
14. Fedkiw, R. P., T. Aslam, B. Merriman, and S. Osher (1999), A nonoscillatory Eulerian approach to interfaces in multimaterial flows (the ghost fluid method), J. Comput. Phys., 152(2), 457-492, doi:10.1006/jcph. 1999.6236.
15. Fleitout, L., and C. Froidevaux (1982), Tectonics and topography for a lithosphere containing density heterogeneities, Tectonics, 1(1), 21-56, doi:10.1029/TC001i001p00021.
16. Gerya, T. V., and D. A. Yuen (2003), Characteristics-based marker-in-cell method with conservative finite-differences schemes for modeling geological flows with strongly variable transport properties, Phys. Earth Planet. Inter., 140(4), 293-318, doi:10.1016/j.pepi.2003.09.006. (Pubitemid 37502872)
17. Gerya, T. V., and D. A. Yuen (2007), Robust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems, Phys. Earth Planet. Inter., 163(1-4), 83-105, doi:10.1016/j.pepi. 2007.04.015.
18. Glimm, J., D. Marchesin, and O. McBryan (1981), A numerical method for two phase flow with an unstable interface, J. Comput. Phys., 39(1), 179-200, doi:10.1016/0021-9991(81)90144-3.
19. Hale, A. J. (2008), Lava dome growth and evolution with an independently deformable talus, Geophys. J. Int., 174(1), 391-417, doi:10.1111/j.1365-246X. 2008.03806.x.
20. Hale, A. J., L. Bourgouin, and H. B. Mühlhaus (2007), Using the level-set method to model endogenous lava dome growth, J. Geophys. Res., 112, B03213, doi:10.1029/2006JB004445. (Pubitemid 47386342)
21. Hoogenboom, T., and G. A. Houseman (2006), Rayleigh-Taylor instability as a mechanism for corona formation on Venus, Icarus, 180(2), 292-307, doi:10.1016/j.icarus.2005.11.001. (Pubitemid 43092578)
22. Houseman, G. A., D. P. McKenzie, and P. Molnar (1981), Convective instability of a thickened boundary layer and its relevance for the thermal evolution of continental convergent belts, J. Geophys. Res., 86(B7), 6115-6132, doi:10.1029/JB086iB07p06115.
23. Ismail-Zadeh, A., I. Tsepelev, C. Talbot, and A. Korotkii (2004), Threedimensional forward and backward modelling of diapirism: Numerical approach and its applicability to the evolution of salt structures in the Pricaspian basin, Tectonophysics, 387(1-4), 81-103, doi:10.1016/j.tecto. 2004.06.006.
24. Jiang, G.-S., and D. Peng (2000), Weighted ENO schemes for Hamilton Jacobi equat ions, SIAM J. Sci. Comput., 21(6), 2126-2143, doi:10.1137/ S106482759732455X. (Pubitemid 32076060)
25. Jiang, G.-S., and C.-W. Shu (1996), Efficient implementation of weighted ENO schemes, J. Comput. Phys., 126(1), 202-228, doi:10.1006/ jcph.1996.0130.
26. Kang, M., R. P. Fedkiw, and X. D. Liu (2000), A boundary condition capturing method for multiphase incompressible flow, J. Sci. Comput., 15(3), 323-360, doi:10.1023/A:1011178417620. (Pubitemid 32508315)
27. Keck, R. (2007), Reinitialization for level set methods, Master's thesis, Univ. of Kaiserslautern, Kaiserslautern, Germany.
28. Lenardic, A., and W. M. Kaula (1993), A numerical treatment of geodynamic viscous flow problems involving the advection of material interfaces, J. Geophys. Res., 98(B5), 8243-8260, doi:10.1029/92JB02858.
29. LeVeque, R. J., and Z. Li (1994), The immersed interface method for elliptic equations with discontinuous coefficients and singular sources, SIAM J. Numer. Anal., 31(4), 1019-1044, doi:10.1137/0731054.
30. LeVeque, R. J., and Z. Li (1997), Immersed interface methods for Stokes flow with elastic boundaries or surface tension, SIAM J. Sci. Comput., 18(3), 709-735, doi:10.1137/S1064827595282532.
31. Lin, S. C., and P. E. van Keken (2006), Dynamics of thermochemical plumes: 1. Plume formation and entrainment of a dense layer, Geochem. Geophys. Geosyst., 7, Q02006, doi:10.1029/2005GC001071.
32. Liu, X.-D., S. Osher, and T. Chan (1994), Weighted essentially non-oscillatory schemes, J. Comput. Phys., 115(1), 200-212, doi:10.1006/jcph. 1994. 1187.
33. Liu, X.-D., R. P. Fedkiw, and M. Kang (2000), A boundary condition capturing method for Poisson's equation on irregular domains, J. Comput. Phys., 160(1), 151-178, doi:10.1006/jcph.2000.6444.
34. Malladi, R., J. A. Sethian, and B. C. Vemuri (1995), Shape modeling with front propagation: A level set approach, IEEE Trans. Pattern Anal. Mach. Intell., 17(2), 158-175, doi:10.1109/34.368173.
35. Manga, M., and H. A. Stone (1993), Buoyancy-driven interactions between two deformable viscous drops, J. Fluid Mech., 256, 647-683, doi:10.1017/ S0022112093002915.
36. Manga, M., and H. A. Stone (1994), Interactions between bubbles in magmas and lavas: Effects of bubble deformation, J. Volcanol. Geotherm. Res., 63(3-4), 267-279, doi:10.1016/0377-0273(94)90079-5.
37. Manga, M., and H. A. Stone (1995), Collective hydrodynamics of deformable drops and bubbles in dilute low Reynolds number suspensions, J. Fluid Mech., 300, 231-263, doi:10.1017/S0022112095003673.
38. Manga, M., H. A. Stone, and R. J. O'Connell (1993), The interaction of plume heads with compositional discontinuities in the Earth's mantle, J. Geophys. Res., 98(B11), 19,979-19,990, doi:10.1029/93JB00441.
39. Marsh, B. D. (1988), Crystal capture, sorting, and retention in convecting magma, Geol. Soc. Am. Bull., 100(11), 1720-1737, doi:10.1130/0016- 7606(1988)100<1720:CCSARI>2.3.CO;2.
40. Mayo, A. (1984), The fast solution of Poisson's and the biharmonic equations on irregular regions, SIAM J. Numer. Anal., 21(2), 285-299, doi:10.1137/0721021.
41. Mayo, A. (1985), Fast high order accurate solution of Laplace's equation on irregular regions, SIAM J. Sci. Stat. Comput., 6(1), 144-157, doi:10.1137/ 0906012.
42. Moresi, L., F. Dufour, and H. B. Mühlhaus (2003), A Lagrangian integration point finite element method for large deformation modeling of viscoelastic geomaterials, J. Comput. Phys., 184(2), 476-497, doi:10.1016/ S0021-9991(02)00031-41
43. Osher, S., and J. A. Sethian (1988), Fronts propagating with curvaturedependent speed: Algorithms based on Hamilton-Jacobi formulations, J. Comput. Phys., 79(1), 12-49, doi:10.1016/0021-9991(88)90002-2.
44. Römer, M. M., and H. J. Neugebauer (1991), The salt dome problem: A multilayered approach, J. Geophys. Res., 96(B2), 2389-2396, doi:10.1029/ 90JB02193.
45. Royden, L. H., and L. Husson (2006), Trench motion, slab geometry and viscous stresses in subduction systems, Geophys. J. Int., 167(2), 881-905, doi:10.1111/j.1365-246X.2006.03079.x. (Pubitemid 44697446)
46. Schaeffer, N., and M. Manga (2001), Interaction of rising and sinking mantle plumes, Geophys. Res. Lett., 28(3), 455-458, doi:10.1029/ 2000GL012215.
47. Schenk, O., and K. Gärtner (2004), Solving unsymmetric sparse systems of linear equations with PARDISO, J. Future Gen. Comput. Syst., 20(3), 475-487, doi:10.1016/j.future.2003.07.011.
48. Schenk, O., and K. Gärtner (2006), On fast factorization pivoting methods for sparse symmetric indefinite systems, Electr. Trans. Numer. Anal., 23, 158-179.
49. Schmalzl, J., and A. Loddoch (2003), Using subdivision surfaces and adaptive surface simplification algorithms for modeling chemical heterogeneities in geophysical flows, Geochem. Geophys. Geosyst., 4(9), 8303, doi:10.1029/2003GC000578.
50. Schmeling, H. (1987), On the relation between initial conditions and late stages of Rayleigh-Taylor instabilities, Tectonophysics, 133(1-2), 65-80, doi:10.1016/0040-1951(87)90281-2.
51. Schmeling, H., et al. (2008), A benchmark comparison of spontaneous subduction models-Towards a free surface, Phys. Earth Planet. Inter., 171(1-4), 198-223, doi:10.1016/j.pepi.2008.06.028.
52. Schmid, D., and Y. Podladchikov (2003), Analytical solutions for deformable elliptical inclusions in general shear, Geophys. J. Int., 155(1), 269-288, doi:10.1046/j.1365-246X.2003.02042.x. (Pubitemid 37293471)
53. Schott, B., D. A. Yuen, and H. Schmeling (2000), The significance of shear heating in continental delamination, Phys. Earth Planet. Inter., 118(3-4), 273-290, doi:10.1016/S0031-9201(99)00159-4. (Pubitemid 30218494)
54. Sethian, J. A. (1996), A fast marching level set method for monotonically advancing fronts, Proc. Natl. Acad. Sci. U. S. A., 93(4), 1591-1595, doi:10.1073/pnas.93.4.1591. (Pubitemid 26075170)
55. Sethian, J. A. (1997), Level set methods: An act of violence, Am. Sci., 85(3), 1-23.
56. Sethian, J. A. (1999a), Level Set Methods and Fast Marching Methods, 2nd ed., Cambridge Univ. Press, Cambridge, U. K.
57. Sethian, J. A. (1999b), Fast marching methods, SIAM Rev., 41(2), 199-235, doi:10.1137/S0036144598347059.
58. Shu, C. W., and S. Osher (1988), Efficient implementation of essentially non-oscillatory shock capturing schemes, J. Comput. Phys., 77(2), 439-471, doi:10.1016/0021-9991(88)90177-5.
59. Sussman, M., and E. Fatemi (1999), An efficient, interface-preserving level set redistancing algorithm and its application to interfacial incompressible fluid flow, SIAM J. Sci. Comput., 20(4), 1165-1191, doi:10.1137/ S1064827596298245. (Pubitemid 32251212)
60. Sussman, M., P. Smereka, and S. Osher (1994), A level set approach for computing solutions to incompressible two-phase flow, J. Comput. Phys., 114(1), 146-159, doi:10.1006/jcph.1994.1155. (Pubitemid 124013478)
61. Sussman, M., E. Fatemi, P. Smereka, and S. Osher (1998), An improved level set method for incompressible two-phase flows, Comput. Fluids, 27(5-6), 663-680, doi:10.1016/S0045-7930(97)00053-4. (Pubitemid 28437604)
62. Tackley, P. J., and S. D. King (2003), Testing the tracer ratio method for modeling active compositional fields in mantle convection simulations, Geochem. Geophys. Geosyst., 4(4), 8302, doi:10.1029/2001GC000214.
63. Tornberg, A. K., and B. Engquist (2003), Regularization techniques for numerical approximation of PDEs with singularities, J. Sci. Comput., 19(1-3), 527-552, doi:10.1023/A:1025332815267.
64. van Keken, P. E. (1993), Numerical modeling of thermochemically driven fluid flow, applied to the Earth's lithosphere and mantle, Ph.D. thesis, Univ. of Utrecht, Utrecht, Netherlands.
65. van Keken, P. E., S. D. King, H. Schmeling, U. R. Christensen, D. Neumeister, and M.-P. Doin (1997), A comparison of methods for the modeling of thermochemical convection, J. Geophys. Res., 102(B10), 22,477-22,495, doi:10.1029/97JB01353.
66. Woidt, W. D. (1978), Finite element calculations applied to salt dome analysis, Tectonophysics, 50(2-3), 369-386, doi:10.1016/0040-1951(78) 90143-9.