{"title":"A pixel-based finite element implementation to estimate effective wave velocity in heterogeneous media","authors":"","doi":"10.1016/j.jappgeo.2024.105447","DOIUrl":null,"url":null,"abstract":"<div><p>In the present work, we present a 2D pixel-based Finite Element strategy to simulate the elastic wave propagation in heterogeneous media. An assembly-free approach is employed for the stiffness matrix, leveraging a pixel-based structured mesh to reduce the memory required to store computations. Additionally, a diagonal Lumped-Mass matrix technique is utilized to address challenges associated with the inversion and storage of the mass matrix. The Leap-frog integration method, known for its amalgamation of stability, precision, and efficiency, is adopted. The combination of these features is aimed at facilitating massive parallel computations for very large systems with 10<sup>8</sup> to 10<sup>9</sup> degrees of freedom. In that sense, the present work can be understood as a first step toward a very efficient massive parallel GPU-based voxel-based Finite Element implementation to treat very large digital images with personal computers. The implementation presented here has been validated against theoretical predictions and analytical results derived from classical wave propagation theory. Finally, the transmission test is simulated in two digital models, one representing a layered medium and another representing a medium with complex microstructue obtained via micro-tomography. For the first model, the results are compared with the so called Bakus average, while, for the second model, the results are compared with the corresponding outcomes acquired through an in-house developed static finite element homogenization implementation.</p></div>","PeriodicalId":54882,"journal":{"name":"Journal of Applied Geophysics","volume":null,"pages":null},"PeriodicalIF":2.2000,"publicationDate":"2024-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Applied Geophysics","FirstCategoryId":"89","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0926985124001630","RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"GEOSCIENCES, MULTIDISCIPLINARY","Score":null,"Total":0}
引用次数: 0
Abstract
In the present work, we present a 2D pixel-based Finite Element strategy to simulate the elastic wave propagation in heterogeneous media. An assembly-free approach is employed for the stiffness matrix, leveraging a pixel-based structured mesh to reduce the memory required to store computations. Additionally, a diagonal Lumped-Mass matrix technique is utilized to address challenges associated with the inversion and storage of the mass matrix. The Leap-frog integration method, known for its amalgamation of stability, precision, and efficiency, is adopted. The combination of these features is aimed at facilitating massive parallel computations for very large systems with 108 to 109 degrees of freedom. In that sense, the present work can be understood as a first step toward a very efficient massive parallel GPU-based voxel-based Finite Element implementation to treat very large digital images with personal computers. The implementation presented here has been validated against theoretical predictions and analytical results derived from classical wave propagation theory. Finally, the transmission test is simulated in two digital models, one representing a layered medium and another representing a medium with complex microstructue obtained via micro-tomography. For the first model, the results are compared with the so called Bakus average, while, for the second model, the results are compared with the corresponding outcomes acquired through an in-house developed static finite element homogenization implementation.
期刊介绍:
The Journal of Applied Geophysics with its key objective of responding to pertinent and timely needs, places particular emphasis on methodological developments and innovative applications of geophysical techniques for addressing environmental, engineering, and hydrological problems. Related topical research in exploration geophysics and in soil and rock physics is also covered by the Journal of Applied Geophysics.