Book contents
- Frontmatter
- Contents
- Acknowledgements
- Introduction
- 1 The continuity equation
- 2 Density and gravity
- 3 Numerical solutions of partial differential equations
- 4 Stress and strain
- 5 The momentum equation
- 6 Viscous rheology of rocks
- 7 Numerical solutions of the momentum and continuity equations
- 8 The advection equation and marker-in-cell method
- 9 The heat conservation equation
- 10 Numerical solution of the heat conservation equation
- 11 2D thermomechanical code structure
- 12 Elasticity and plasticity
- 13 2D implementation of visco-elasto-plastic rheology
- 14 The multigrid method
- 15 Programming of 3D problems
- 16 Numerical benchmarks
- 17 Design of 2D numerical geodynamic models
- Epilogue: outlook
- Appendix: MATLAB program examples
- References
- Index
Appendix: MATLAB program examples
Published online by Cambridge University Press: 05 June 2012
- Frontmatter
- Contents
- Acknowledgements
- Introduction
- 1 The continuity equation
- 2 Density and gravity
- 3 Numerical solutions of partial differential equations
- 4 Stress and strain
- 5 The momentum equation
- 6 Viscous rheology of rocks
- 7 Numerical solutions of the momentum and continuity equations
- 8 The advection equation and marker-in-cell method
- 9 The heat conservation equation
- 10 Numerical solution of the heat conservation equation
- 11 2D thermomechanical code structure
- 12 Elasticity and plasticity
- 13 2D implementation of visco-elasto-plastic rheology
- 14 The multigrid method
- 15 Programming of 3D problems
- 16 Numerical benchmarks
- 17 Design of 2D numerical geodynamic models
- Epilogue: outlook
- Appendix: MATLAB program examples
- References
- Index
Summary
Introduction
Program 1: Visualisation_is_important.m (Exercise Introduction.2) – visualisation of ‘sin’ and ‘cos’ functions with ‘plot’, ‘pcolor’, ‘contour’ and ‘surf’.
Chapter 1
Program 2: Divergence.m (Exercise 1.2) – computation and visualisation of velocity, divergence of velocity and time derivatives of density with ‘pcolor’ and ‘quiver’.
Chapter 2
Program 3: Periclase_EOS.m (Exercise 2.2) – computation and visualisation of density, thermal expansion and compressibility for periclase (MgO) using external Gibbs free energy function G_periclase.m.
Program 4: Density_map.m (Exercise 2.3) – loading from data files (m895_ro, morn_ro) and visualising density maps for pyrolite (m895_ro) and MORB (morn_ro) and density difference between pyrolite and MORB.
Chapter 3
Program 5: Poisson1D.m (Exercise 3.1) – solution of 1D Poisson equation with finite differences on a regular grid using direct solver ‘\’.
Program 6: Poisson2D_direct.m (Exercise 3.2) – solution of 2D Poisson equation with finite differences on a regular grid using direct solver ‘\’.
Program 7: Poisson2D_Gauss_Seidel.m (Exercise 3.3) – solution of 2D Poisson equation with finite differences on a regular grid using Gauss–Seidel iteration.
Program 8: Poisson2D_Jacobi.m (Exercise 3.4) – solution of 2D Poisson equation with finite differences on a regular grid using Jacobi iteration.
Chapter 4
Program 9: Strain_rate.m (Exercise 4.2) – computation and visualisation of velocity field, strain rate, deviatoric strain rate, and second strain rate invariant.
- Type
- Chapter
- Information
- Introduction to Numerical Geodynamic Modelling , pp. 319 - 325Publisher: Cambridge University PressPrint publication year: 2009