Articles | Volume 13, issue 12
https://doi.org/10.5194/gmd-13-6265-2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-13-6265-2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Model description paper
| 
10 Dec 2020
Model description paper |
| 10 Dec 2020
| 10 Dec 2020
A fast and efficient MATLAB-based MPM solver: fMPMM-solver v1.1
Institute of Earth Sciences, University of Lausanne, 1015 Lausanne, Switzerland
Yury Alkhimenkov
Institute of Earth Sciences, University of Lausanne, 1015 Lausanne, Switzerland
Swiss Geocomputing Centre, University of Lausanne, 1015 Lausanne, Switzerland
Faculty of Mechanics and Mathematics, Lomonosov Moscow State University, Moscow, 119899, Russia
Michel Jaboyedoff
Institute of Earth Sciences, University of Lausanne, 1015 Lausanne, Switzerland
Swiss Geocomputing Centre, University of Lausanne, 1015 Lausanne, Switzerland
Yury Y. Podladchikov
Institute of Earth Sciences, University of Lausanne, 1015 Lausanne, Switzerland
Swiss Geocomputing Centre, University of Lausanne, 1015 Lausanne, Switzerland
Faculty of Mechanics and Mathematics, Lomonosov Moscow State University, Moscow, 119899, Russia
Related authors
Clément Hibert, François Noël, David Toe, Miloud Talib, Mathilde Desrues, Emmanuel Wyser, Ombeline Brenguier, Franck Bourrier, Renaud Toussaint, Jean-Philippe Malet, and Michel Jaboyedoff
Earth Surf. Dynam., 12, 641–656, https://doi.org/10.5194/esurf-12-641-2024, https://doi.org/10.5194/esurf-12-641-2024, 2024
Short summary
Short summary
Natural disasters such as landslides and rockfalls are mostly difficult to study because of the impossibility of making in situ measurements due to their destructive nature and spontaneous occurrence. Seismology is able to record the occurrence of such events from a distance and in real time. In this study, we show that, by using a machine learning approach, the mass and velocity of rockfalls can be estimated from the seismic signal they generate.
Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
Geosci. Model Dev., 14, 7749–7774, https://doi.org/10.5194/gmd-14-7749-2021, https://doi.org/10.5194/gmd-14-7749-2021, 2021
Short summary
Short summary
We propose an implementation of the material point method using graphical processing units (GPUs) to solve elastoplastic problems in three-dimensional configurations, such as the granular collapse or the slumping mechanics, i.e., landslide. The computational power of GPUs promotes fast code executions, compared to a traditional implementation using central processing units (CPUs). This allows us to study complex three-dimensional problems tackling high spatial resolution.
Xiong Tang, Wei Liu, Siming He, Lei Zhu, Michel Jaboyedoff, Huanhuan Zhang, Yuqing Sun, and Zenan Huo
Geosci. Model Dev., 18, 4743–4758, https://doi.org/10.5194/gmd-18-4743-2025, https://doi.org/10.5194/gmd-18-4743-2025, 2025
Short summary
Short summary
The paper presents an explicit stabilized two-phase material point method (MPM) based on the one-point two-phase MPM scheme. The novelty of the work lies in the employment of stabilized techniques, including the strain smoothing method and the multi-field variational principle. With its effective and easy-to-implement stabilized techniques, the proposed model offers an effective and reliable approach for simulating both static and dynamic processes in solid–fluid porous media.
Yury Alkhimenkov and Yury Y. Podladchikov
Geosci. Model Dev., 18, 563–583, https://doi.org/10.5194/gmd-18-563-2025, https://doi.org/10.5194/gmd-18-563-2025, 2025
Short summary
Short summary
The accelerated pseudo-transient (APT) method is an efficient way to solve partial differential equations, particularly well-suited for parallel computing. This paper explores the APT method's effectiveness in solving elastic, viscoelastic, and hydromechanical problems, focusing on quasi-static conditions in 1D, 2D, and 3D. The study examines the best numerical settings for fast and accurate solutions. The paper shows how the APT method can handle complex problems in high-resolution models.
Yury Alkhimenkov, Lyudmila Khakimova, and Yury Podladchikov
EGUsphere, https://doi.org/10.5194/egusphere-2024-3237, https://doi.org/10.5194/egusphere-2024-3237, 2024
Short summary
Short summary
This study examines stress drops and earthquake nucleation in elasto-plastic media using 2D simulations, highlighting the importance of high temporal and spatial resolutions in capturing stress evolution and strain fields. Stress drops reflect fault rupture mechanics and emulate earthquake behavior. The non-Gaussian distribution of stress drop amplitudes resembles "solid turbulence." Elasto-plastic models simulate key earthquake processes and could improve seismic hazard assessment.
Yury Alkhimenkov and Yury Y. Podladchikov
EGUsphere, https://doi.org/10.5194/egusphere-2024-3238, https://doi.org/10.5194/egusphere-2024-3238, 2024
Short summary
Short summary
This paper presents a rigorous derivation of Gassmann's equations, grounded in thermodynamic principles and conservation laws, addressing gaps and potential inconsistencies in the original formulation. It also explores Biot's poroelastic equations, demonstrating that Gassmann's equations are a specific case within Biot’s framework. The study affirms the robustness of Gassmann's equations when assumptions are met, and symbolic Maple routines are provided to ensure reproducibility of the results.
Charlotte Wolff, Marc-Henri Derron, Carlo Rivolta, and Michel Jaboyedoff
Geosci. Instrum. Method. Data Syst., 13, 225–248, https://doi.org/10.5194/gi-13-225-2024, https://doi.org/10.5194/gi-13-225-2024, 2024
Short summary
Short summary
The remote-sensing InSAR technique is vital for monitoring slope instabilities but requires understanding. This paper delves into differences between satellite and GB-InSAR. It offers a tool to determine the optimal GB-InSAR installation site, considering various technical, meteorological, and topographical factors. By generating detailed maps and simulating radar image characteristics, the tool eases the setup of monitoring campaigns for effective and accurate ground movement tracking.
Clément Hibert, François Noël, David Toe, Miloud Talib, Mathilde Desrues, Emmanuel Wyser, Ombeline Brenguier, Franck Bourrier, Renaud Toussaint, Jean-Philippe Malet, and Michel Jaboyedoff
Earth Surf. Dynam., 12, 641–656, https://doi.org/10.5194/esurf-12-641-2024, https://doi.org/10.5194/esurf-12-641-2024, 2024
Short summary
Short summary
Natural disasters such as landslides and rockfalls are mostly difficult to study because of the impossibility of making in situ measurements due to their destructive nature and spontaneous occurrence. Seismology is able to record the occurrence of such events from a distance and in real time. In this study, we show that, by using a machine learning approach, the mass and velocity of rockfalls can be estimated from the seismic signal they generate.
François Noël, Michel Jaboyedoff, Andrin Caviezel, Clément Hibert, Franck Bourrier, and Jean-Philippe Malet
Earth Surf. Dynam., 10, 1141–1164, https://doi.org/10.5194/esurf-10-1141-2022, https://doi.org/10.5194/esurf-10-1141-2022, 2022
Short summary
Short summary
Rockfall simulations are often performed to make sure infrastructure is safe. For that purpose, rockfall trajectory data are needed to calibrate the simulation models. In this paper, an affordable, flexible, and efficient trajectory reconstruction method is proposed. The method is tested by reconstructing trajectories from a full-scale rockfall experiment involving 2670 kg rocks and a flexible barrier. The results highlight improvements in precision and accuracy of the proposed method.
Ludovic Räss, Ivan Utkin, Thibault Duretz, Samuel Omlin, and Yuri Y. Podladchikov
Geosci. Model Dev., 15, 5757–5786, https://doi.org/10.5194/gmd-15-5757-2022, https://doi.org/10.5194/gmd-15-5757-2022, 2022
Short summary
Short summary
Continuum mechanics-based modelling of physical processes at large scale requires huge computational resources provided by massively parallel hardware such as graphical processing units. We present a suite of numerical algorithms, implemented using the Julia language, that efficiently leverages the parallelism. We demonstrate that our implementation is efficient, scalable and robust and showcase applications to various geophysical problems.
Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
Geosci. Model Dev., 14, 7749–7774, https://doi.org/10.5194/gmd-14-7749-2021, https://doi.org/10.5194/gmd-14-7749-2021, 2021
Short summary
Short summary
We propose an implementation of the material point method using graphical processing units (GPUs) to solve elastoplastic problems in three-dimensional configurations, such as the granular collapse or the slumping mechanics, i.e., landslide. The computational power of GPUs promotes fast code executions, compared to a traditional implementation using central processing units (CPUs). This allows us to study complex three-dimensional problems tackling high spatial resolution.
Martin Franz, Michel Jaboyedoff, Ryan P. Mulligan, Yury Podladchikov, and W. Andy Take
Nat. Hazards Earth Syst. Sci., 21, 1229–1245, https://doi.org/10.5194/nhess-21-1229-2021, https://doi.org/10.5194/nhess-21-1229-2021, 2021
Short summary
Short summary
A landslide-generated tsunami is a complex phenomenon that involves landslide dynamics, wave dynamics and their interaction. This phenomenon threatens numerous lives and infrastructures around the world. To assess this natural hazard, we developed an efficient numerical model able to simulate the landslide, the momentum transfer and the wave all at once. The good agreement between the numerical simulations and physical experiments validates our model and its novel momentum transfer approach.
Cited articles
Abe, K., Soga, K., and Bandara, S.: Material point method for coupled
hydromechanical problems, J. Geotechn. Geoenviron.
Eng., 140, 04013033, https://doi.org/10.1061/(ASCE)GT.1943-5606.0001011, 2014. a
Anderson Jr., C. E.: An overview of the theory of hydrocodes, Int.
J. Impact Eng., 5, 33–59, 1987. a
Bandara, S. and Soga, K.: Coupling of soil deformation and pore fluid flow
using material point method, Comput. Geotech., 63, 199–214, 2015. a
Bardenhagen, S., Brackbill, J., and Sulsky, D.: The material-point method for
granular materials, Comput. Method. Appl. M.,
187, 529–541, 2000. a
Baumgarten, A. S. and Kamrin, K.: A general fluid–sediment mixture model and
constitutive theory validated in many flow regimes, J. Fluid
Mechan., 861, 721–764, 2019. a
Beuth, L., Benz, T., Vermeer, P. A., and Więckowski, Z.: Large
deformation analysis using a quasi-static material point method, J.
Theor. Appl. Mechan., 38, 45–60, 2008. a
Coombs, W. M., Augarde, C. E., Brennan, A. J., Brown, M. J., Charlton, T. J.,
Knappett, J. A., Motlagh, Y. G., and Wang, L.: On Lagrangian mechanics and
the implicit material point method for large deformation elasto-plasticity,
Comput. Method. Appl. Mechan., 358, 112622, https://doi.org/10.1016/j.cma.2019.112622, 2020. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p
Cortis, M., Coombs, W., Augarde, C., Brown, M., Brennan, A., and Robinson, S.:
Imposition of essential boundary conditions in the material point method,
Int. J. Num. Method., 113, 130–152,
2018. a
Davis, T. A.: Suite Sparse,
available at: https://people.engr.tamu.edu/davis/research.html (last access: 6 October 2020), 2013. a
de Koster, P., Tielen, R., Wobbes, E., and Möller, M.: Extension of
B-spline Material Point Method for unstructured triangular grids using
Powell–Sabin splines, Comput. Part. Mechan., 1–16, https://doi.org/10.1007/s40571-020-00328-3, 2020. a
de Souza Neto, E. A., Peric, D., and Owen, D. R.: Computational methods for
plasticity: theory and applications, John Wiley & Sons, John Wiley & Sons Ltd, The Atrium, Southern Gate, Chichester, West Sussex, PO19 8SQ, United Kingdom, 2011. a
de Vaucorbeil, A., Nguyen, V., and Hutchinson, C.: A Total-Lagrangian
Material Point Method for solid mechanics problems involving large
deformations, Computer Methods in Applied Mechanics and Engineering, 360, https://doi.org/10.1016/j.cma.2019.112783, 2020. a
Dunatunga, S. and Kamrin, K.: Continuum modelling and simulation of granular
flows through their many phases, J. Fluid Mechan., 779, 483–513,
2015. a
Fern, J., Rohe, A., Soga, K., and Alonso, E.: The Material Point Method for Geotechnical Engineering. Boca Raton: CRC Press, https://doi.org/10.1201/9780429028090, 2019. a
Gan, Y., Sun, Z., Chen, Z., Zhang, X., and Liu, Y.: Enhancement of the material
point method using B-spline basis functions, Int. J.
Num. Method., 113, 411–431, 2018. a
Gaume, J., van Herwijnen, A., Gast, T., Teran, J., and Jiang, C.:
Investigating the release and flow of snow avalanches at the slope-scale
using a unified model based on the material point method, Cold Reg.
Sci. Technol., 168, 102847, https://doi.org/10.1016/j.coldregions.2019.102847, 2019. a
Gracia, F., Villard, P., and Richefeu, V.: Comparison of two numerical
approaches (DEM and MPM) applied to unsteady flow, Comput. Part.
Mechan., 6, 591–609, 2019. a
Homel, M. A., Brannon, R. M., and Guilkey, J.: Controlling the onset of
numerical fracture in parallelized implementations of the material point
method (MPM) with convective particle domain interpolation (CPDI) domain
scaling, Int. J. Num. Method., 107,
31–48, 2016. a
Iaconeta, I., Larese, A., Rossi, R., and Guo, Z.: Comparison of a material
point method and a galerkin meshfree method for the simulation of
cohesive-frictional materials, Materials, 10, 1150, 2017. a
Leavy, R., Guilkey, J., Phung, B., Spear, A., and Brannon, R.: A
convected-particle tetrahedron interpolation technique in the material-point
method for the mesoscale modeling of ceramics, Comput. Mechan., 64,
563–583, 2019. a
Moler, C.: MATLAB Incorporates LAPACK,
available at: https://ch.mathworks.com/de/company/newsletters/articles/matlab-incorporates-lapack.html?refresh=true (last access: 6 October 2020),
2000. a
Ni, R. and Zhang, X.: A precise critical time step formula for the explicit
material point method, Int. J. Num. Method., 121, 4989–5016, 2020. a
Steffen, M., Wallstedt, P., Guilkey, J., Kirby, R., and Berzins, M.:
Examination and analysis of implementation choices within the material point
method (MPM), Comput. Model. Eng. Sci., 31, 107–127,
2008b. a
Stomakhin, A., Schroeder, C., Chai, L., Teran, J., and Selle, A.: A material
point method for snow simulation, ACM Transactions on Graphics (TOG), 32,
1–10, 2013. a
Sulsky, D., Zhou, S.-J., and Schreyer, H. L.: Application of a particle-in-cell
method to solid mechanics, Comput. Phys. Commun., 87, 236–252,
1995. a
Vardon, P. J., Wang, B., and Hicks, M. A.: Slope failure simulations with MPM,
J. Hydrodynam., 29, 445–451, 2017. a
Vermeer, P. A. and De Borst, R.: Non-associated plasticity for soils, concrete
and rock, HERON, 29, 1984, 163–196, 1984. a
Wallstedt, P. C. and Guilkey, J.: An evaluation of explicit time integration
schemes for use with the generalized interpolation material point method,
J. Computat. Phys., 227, 9628–9642, 2008. a
Wang, B., Hicks, M., and Vardon, P.: Slope failure analysis using the random
material point method, Géotech. Lett. 6, 113–118,
2016a. a
Wang, B., Vardon, P., and Hicks, M.: Investigation of retrogressive and
progressive slope failure mechanisms using the material point method,
Comput. Geotech., 78, 88–98, 2016b. a
Wang, L., Coombs, W. M., Augarde, C., Cortis, E. M., Charlton, T. J., Brown, M. J., Knappett, J., Brennan, A., Davidson, C., Richards, and Blake, D. A.: On the use of
domain-based material point methods for problems involving large distortion,
Comput. Method. Appl. Mechan. Eng., 355, 1003–1025, 2019. a, b, c, d
Więckowski, Z.: The material point method in large strain engineering
problems, Comput. Method. Appl. Mechan. Eng., 193,
4417–4438, 2004. a
Wyser, E., Alkhimenkov, Y., Jayboyedoff, M., and Podladchikov, Y.:
fMPMM-solver, Zenodo, https://doi.org/10.5281/zenodo.4068585, 2020a. a
Wyser, E., Alkhimenkov, Y., Jayboyedoff, M., and Podladchikov, Y.:
fMPMM, available at: https://bitbucket.org/ewyser/fmpmm-solver/src/master/, last access: 6 October 2020. a
York, A. R., Sulsky, D., and Schreyer, H. L.: The material point method for
simulation of thin membranes, Int. J. Num. Method., 44, 1429–1456, 1999. a
Short summary
In this work, we present an efficient and fast material point method (MPM) implementation in MATLAB. We first discuss the vectorization strategies to adapt this numerical method to a MATLAB implementation. We report excellent agreement of the solver compared with classical analysis among the MPM community, such as the cantilever beam problem. The solver achieves a performance gain of 28 compared with a classical iterative implementation.
In this work, we present an efficient and fast material point method (MPM) implementation in...
