Comparison of Eulerian and Lagrangian numerical techniques for the Stokes equations in the presence of strongly varying viscosity

Physics of the Earth and Planetary Interiors

Volumn 171, Issue 1-4, 2008, Pages 92-111

  Deubelbeiss, Y ^a   Kaus, B J P ^a,b  

^a ETH ZURICH   (Switzerland)

^b UNIVERSITY OF SOUTHERN CALIFORNIA   (United States)

Author keywords

Analytics; Finite difference methods; Finite element methods; Geodynamic modeling; Lagrangian and Eulerian grids; Numerical methods; Rotated staggered grid; Staggered grid; Stokes equations; Stream function

Indexed keywords

CONVERGENCE OF NUMERICAL METHODS; DIFFERENCE EQUATIONS; FINITE DIFFERENCE METHOD; FLUID DYNAMICS; GEOPHYSICS; HARMONIC ANALYSIS; HYDRAULIC CONDUCTIVITY; HYDRODYNAMICS; LAGRANGE MULTIPLIERS; NUMERICAL METHODS; PHOTOLITHOGRAPHY; PROBABILITY DENSITY FUNCTION; PROBLEM SOLVING; SOLUTIONS; STRESS CONCENTRATION; VISCOSITY;

ANALYTICS; FINITE DIFFERENCE METHODS; GEODYNAMIC MODELING; LAGRANGIAN AND EULERIAN GRIDS; ROTATED STAGGERED GRID; STAGGERED GRID; STOKES EQUATIONS; STREAM FUNCTION;

FINITE ELEMENT METHOD;

EULERIAN ANALYSIS; FINITE DIFFERENCE METHOD; FINITE ELEMENT METHOD; GEODYNAMICS; LAGRANGIAN ANALYSIS; NUMERICAL MODEL; STOKES FORMULA; VISCOSITY; VISCOUS FLOW;

EID: 55649110882     PISSN: 00319201     EISSN: None     Source Type: Journal    
DOI: 10.1016/j.pepi.2008.06.023     Document Type: Article

References (66)

1. Albers M. A local mesh refinement multigrid method for 3D convection problems with strongly variable viscosity. Journal of Compuational Physics 160 1 (2000) 126-150
2. Auth C., and Harder H. Multigrid solution of convection problems with strongly variable viscosity. Geophysical Journal International 137 (1999) 793-804
3. Belytschko T., Moes N., Usui S., and Parimi C. Arbitrary discontinuitites in finite elements. International Journal for Numerical Methods in Engineering 50 4 (2001) 993-1013
4. Biot M. Theory of folding of stratified viscoelastic media and its implications in tectonics and orogenesis. Geological Society of America Bulletin 72 (1961) 1595-1620
5. Biot M., and Odé H. Theory of gravity instability with variable overburden and compaction. Geophysics 30 (1965) 153-166
6. Braun J., and Sambridge M. Dynamic Lagrangian remeshing (DLR): a new algorithm for solving large strain deformation problems and its application to fault-propagation folding. Earth and Planetary Science Letters 124 (1994) 211-220
7. Buiter, S.J.H., Babeyko, A., Ellis, S., Gerya, T., Kaus, B., Kellner, A., Schreurs, G., Yamada, Y., 2006. The numerical sandbox: comparison of model results for a shortening and an extension experiment. In: Buiter, S.J.H., Schreurs, G., (Eds.), Analogue and Numerical Modelling of Crustal-Scale Processes, vol. 253. Geological Society of London Special Publications, Geological Society of London, pp. 29-64.
8. Burg J.P., and Gerya T.V. The role of viscous heating in barrovian metamorphism of collisional orogens: thermomechanical models and application to the lepontine dome in the central alps. Journal of Metamorphic Geology 23 (2006) 75-95
9. Burg J.P., Kaus B., and Podladchikov Y. Dome structures in collision orogens mechanical investigation of the gravity/compression interplay. In: Whitney D., Teyssier C., and Siddoway C. (Eds). Gneiss Domes in Orogeny vol. 380 (2004), Geological Society of America, Boulder 47-66
10. Choblet G. Modelling thermal convection with large viscosity gradients in one block of the "cubed sphere". Journal of Computational Physics 205 1 (2005) 269-291
11. Connolly J.A.D., and Podladchikov Y. Decompaction weakening and channeling instability in ductile porous media: implications for asthenospheric melt segregation. Journal of Geophysical Research (2007)
12. Cuvelier C., Segal A., and van Steemhoven A. Finite element methods and Navier-Stokes equations. Mathematics and its Applications (1986), Reidel Publishing Company, Dordrecht
13. Dabrowski, M., Krotkiewski, M., Schmid, D.W., 2008. MILAMIN: MATLAB-based FEM solver for large problems. Geochem. Geophys. Geosyst., doi:10.1029/2007GC001719.
14. Elman H., Silvester D., and Wathen A. Finite elements and fast iterative solvers with applications to incompressible fluid dynamics. Numerical Mathematics and Scientific Computation (2005), Oxford University Press
15. Fletcher C. Computational Techniqus for Fluid Dynamics, vol. 2 (1991), Springer
16. Fullsack P. An arbitrary Lagrangian-Eulerian formulation for creeping flows and its application in tectonic models. Geophysical Journal International 120 1 (1995) 1-23
17. Gerya T.V., and Yuen D. Characteristics-based marker-in-cell method with conservative finite-differences schemes for modeling geological flows with strongly variable transport properties. Physics of the Earth and Planetary Interiors 140 (2003) 293-318
18. Gerya T.V., and Yuen D. Robust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems. Physics of the Earth and Planetary Interiors (2007)
19. Gorczyk W., Gerya T., Connolly J., and Yuen D. Growth and mixing dynamics of mantle wedge plumes. Geology 35 (2007) 587-590
20. Huismans R.S., and Beaumont C. Asymmetric lithospheric extension: the role of frictional plastic strain softening inferred from numerical experiments. Geology 30 3 (2002) 211-214
21. Jäger P., Steinmann P., and Kuhl E. On local tracking algorithms for the simulation of three-dimensional. Comput. Mech. 42 (2008) 395-406
22. Jamieson R., Beaumont C., Nguyen C.H., and Grujic D. Provenance of the greater himalayan sequence and associated rocks: predictions of channel flow models. In: Law R., Searle M., and Godin L. (Eds). Channel Flow, Ductile Extrusion and Exhumation of Lower-Mid Crust in Continental Collision Zones vol. 268 (2006), Geological Soc. London Spec. Publ
23. Kameyama M., Kageyama A., and Sato T. Multigrid iterative algorithm using pseudo-compressibility for three-dimensional mantle convection with strongly variable viscosity. Journal of Computational Physics 206 (2005) 162-181
24. Katz R.F., Knepley M.G., Smith B., Spiegelman M., and Coon E. Numerical simulation of geodynamic processes with the portable extensible toolkit for scientific Computation. Physics of th Earth and Planetary Interiors 163 1-4 (2007) 52-68
25. Katz R.F., Spiegelman M., and Holtzman B. The dynamics of melt and shear localization in partially molten aggregates. Nature 442 (2006) 676-679
26. Kaus B., Podladchikov Y., and Schmid D.W. Eulerian spectral/finite difference method for large deformation modelling of visco-elasto-plastic geomaterials. Bolletino Geofisico 45 (2004) 346-349
27. Kaus B.J.P., and Podladchikov Y. Initiation of localized shear zones in viscoelastoplastic rocks. Journal of Geophysical Research 111 B04412 (2006)
28. Landry, W., Hodkinson, L., Kientz, S., 2008. GALE User Manual, Version 1.2.2. http://www.geodynamics.org/cig/software/packages/long/gale.
29. Mancktelow N.S. How ductile are ductile shear zones?. Geology 34 5 (2006) 345-348
30. McKenzie D.P. The generation and compaction of partially molten rock. Journal of Petrology 25 3 (1984) 713-765
31. Moczo B., Kristek J., Vavrycuk V., Archuleta R., and Halada L. 3D heterogeneous staggered-grid finite difference modelling of seismic motion with volume harminic and arithmetic averaging of elastic moduli and densities. Bulletin of the Seismological Society of America 93 (2002) 3042-3066
32. Morency C., Huismans R., Beaumont C., and Fullsack P. A numerical model for coupled fluid flow and matrix deformation with applications to disequilibrium compaction and delta stability. Journal of Geophysical Research 112 (2007)
33. Moresi L., Dufour F., and Mühlaus H.B. Mantle convection modeling with viscoelastic/brittle lithosphere: numerical modeling and plate tectonic modeling. Pure and Applied Geophysics 159 (2002) 2335-2356
34. Moresi L., Dufour F., and Mühlhaus H.B. A Lagrangian integration point finite element method for large deformation modeling of viscoelastic geomaterials. Journal of Computational Physics 184 (2003) 476-497
35. Moresi L., Quenette V., Lemiale V., Meriaux C., Appelbe B., and Mühlhaus H.B. Computational approaches to studying non-linear dynamics of the crust and mantle. Physics of the Earth and Planetary Interiors 163 (2007)
36. Moresi L., and Solomatov V. Numerical investigations of 2D convection in a fluid with extremely large viscosity variations. Physics of Fluids 7 (1995) 2154-2162
37. Moresi L., Zhong S., and Gurnis M. The accuracy of finite element solutions of Stokes's flow with strongly varying viscosity. Physics of the Earth and Planetary Interiors 97 1-4 (1996) 83-94
38. Moskow S., Druskin V., Habashy T., Lee P., and Davydycheva S. A finite difference scheme for elliptic equations with rough coefficients using a Cartesian grid nonconforming to interfaces. SIAM Journal on Numerical Analysis 36 2 (1997) 442-464
39. Mühlhaus H., Bourgouin L., and Hale A. Free surfaces modeling based on Level sets. International Conference on Computational and Experimental Engineering and Sciences 3 4 (2007) 225-232
40. Petrini K., and Podladchikov Y. Lithospheric pressure-depth relationship in compressive regions of thickened crust. Journal of Metamorphic Geology 18 1 (2000) 67-77
41. Ranalli G. Rheology of the Earth. 2nd ed. (1995), Chapman and Hall, London
42. Saenger E., and Bohlen T. Anisotropic and viscoelastic finite-difference modeling using the rotated staggered grid. Geophysics 69 2 (2004) 583-591
43. Saenger E., Gold N., and Shapiro S. Modeling the propagation of elastic waves using a modified finite-difference grid. Wave Motion 31 (2000) 77-92
44. Schmalholz S.M., Podladchikov Y.Y., and Schmid D.W. A spectral/finite difference method for simulating large deformations of heterogeneous, viscoelastic materials. Geophysical Journal International 145 1 (2001) 199-208
45. Schmeling H., and Jacoby W.R. On modelling the lithosphere in mantle convection with nonlinear rheology. Journal of Geophysics 50 (1981) 89-100
46. Schmid, D.W., 2002. Finite and Infinite Heterogeneities Under Pure and Simple Shear. PhD. Swiss Federal Institute of Technology.
47. Schmid D.W., and Podladchikov Y. Analytical solutions for deformable elliptical inclusions in general shear. Geophysical Journal International 155 (2003) 269-288
48. Schott B., Yuen D.A., and Schmeling H. The diversity of tectonics from fluid-dynamical modeling of the lithosphere-mantle system. Tectonophysics 322 1-2 (2000) 35-51
49. Schubert G., Olson P., and Turcotte D.L. Mantle Convection in the Earth and Planets (2001), Cambridge University Press
50. Scott D. The competition between percolation and circulation in a deformable porous medium. Journal of Geophysical Research 93 B6 (1988) 6451-6462
51. Scott D., and Stevenson D.J. Magma ascent by porous flow. Journal of Geophysical Research 91 (1986) 9283-9296
52. Sharples W.K., Moresi L.N., May D.A., Giordani J., and Lemiale V. Gradient recovery methods for improved stresses and strain rates in geodynamic models. Eos Trans. AGU, Fall Meet. Suppl. 87 52 (2006) abstract T13C-530
53. Shaw H. Viscosity of basalt in the melting range. Journal of Petrology 10 (1969) 510-535
54. Spiegelman M. Linear analysis of melt band formation by simple shear. Geochemistry Geophysics Geosystems 4 9 (2003) 8615
55. Sussman M., Smereka P., and Osher S. A level set approach for computing solutions to incompressible two-phase flow. Journal of Computational Physics 114 1 (1994) 146-159
56. Tackley P. Self-consistent generation of tectonic plates in time-dependent, three-dimensional mantle convection simulations. 1. Pseudoplastic yielding. Geochem. Geophys. Geosyst. 01 23 (2000)
57. Tackley P., and King S.D. Testing the tracer ratio method for modeling active compositional fields in mantle convection simulations. Geochem. Geophys. Geosyst. 4 4 (2003)
58. Tannehill, J.C., Anderson, D.A., Pletcher, R.H., 1997. Computational fluid mechanics and heat transfer. Series in Computational and Physical Processes in Mechanics and Thermal Sciences.
59. Toussaint G., Burov E., and Avouac J.P. Tectonic evolution of a continental collision zone: a thermomechanical numerical model. Tectonics 23 (2004) TC6003
60. Trompert R., and Hansen U. The application of a finite volume multigrid method to 3D flow problems in a highly viscous fluid with variable viscosity. Geophysical and Astrophysical Fluid Dynamics 83 (1996) 261
61. Turcotte D.L., and Schubert G. Geodynamics (1982), Wiley
62. Vasilyev O.V., Podladchikov Y.Y., and Yuen D.A. Modeling of plume-lithosphere interaction using adaptive multilevel wavelet collocation algorithm. Geophysial Journal International 147 3 (2001) 579-589
63. Vermeer P.A. The orientation of shear bands in bi-axial tests. Geotechnique 40 2 (1990) 223-236
64. Weinberg R.F., and Schmeling H. Polydiapirs: multiwavelength gravity structures. Journal of Structural Geology 14 4 (1992) 425-436
65. Whitehead J. Fluid models of geological hotspots. Annual Reviews 20 (1988) 61-87
66. Zhong S., Yuen D., and Moresi L. Numerical methods in mantle convection. Treatise in Geophysics 7 (2007) 227-252