Extending Esh3D Code to Solve Interacting Eshelby's Inhomogeneity Problems

摘要

I present a new development of an open source code Esh3D in solving interacting Eshelby's inhomogeneity problems in whole, half, and finite spaces. The code resolves the interaction between multiple inhomogeneity bodies by assembling and solving a global linear system. For half or finite space problems, the code imposes Neumann and Dirichlet boundary conditions and resolves the surface‐body interaction by hybridizing Eshelby's analytical solutions with a numerical part. I demonstrate the impact of multibody interaction and surface‐body interaction by comparing the results of different problems in whole and finite spaces. This analytical‐numerical hybrid approach, compared to conventional finite element method, is inexpensive and flexible in modeling multiple and evolving inhomogeneity bodies. Plain Language Summary In this work I further develop an open source code, Esh3D, to solve interacting Eshelby's inhomogeneity problems in whole, half, and finite spaces. This new development improves the code's applicability to geophysical, as well as to many other, problems.