{"title":"Fully explicit numerical scheme for linearized wave propagation in nearly-incompressible soft hyperelastic solids","authors":"Giulia Merlini, Jean-Marc Allain, Sébastien Imperiale","doi":"10.1016/j.wavemoti.2025.103594","DOIUrl":null,"url":null,"abstract":"<div><div>The numerical approximation of wave propagation problems in nearly or pure incompressible solids faces several challenges such as locking and stability constraints. In this work we propose a stabilized Leapfrog scheme based on the use of Chebyshev polynomials to relax the stability condition, which is strongly limited by the enforcement of incompressibility. The scheme is fully explicit, second order accurate and energy-preserving. For the space discretization we use a mixed formulation with high-order spectral elements and mass-lumping. A strategy is proposed for an efficient and accurate computation of the pressure contribution with a new definition of the discrete Grad-div operator. Finally, we consider linear wave propagation problems in nearly-incompressible hyperelastic solids subject to static preload.</div></div>","PeriodicalId":49367,"journal":{"name":"Wave Motion","volume":"139 ","pages":"Article 103594"},"PeriodicalIF":2.1000,"publicationDate":"2025-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Wave Motion","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0165212525001052","RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"ACOUSTICS","Score":null,"Total":0}
引用次数: 0
Abstract
The numerical approximation of wave propagation problems in nearly or pure incompressible solids faces several challenges such as locking and stability constraints. In this work we propose a stabilized Leapfrog scheme based on the use of Chebyshev polynomials to relax the stability condition, which is strongly limited by the enforcement of incompressibility. The scheme is fully explicit, second order accurate and energy-preserving. For the space discretization we use a mixed formulation with high-order spectral elements and mass-lumping. A strategy is proposed for an efficient and accurate computation of the pressure contribution with a new definition of the discrete Grad-div operator. Finally, we consider linear wave propagation problems in nearly-incompressible hyperelastic solids subject to static preload.
期刊介绍:
Wave Motion is devoted to the cross fertilization of ideas, and to stimulating interaction between workers in various research areas in which wave propagation phenomena play a dominant role. The description and analysis of wave propagation phenomena provides a unifying thread connecting diverse areas of engineering and the physical sciences such as acoustics, optics, geophysics, seismology, electromagnetic theory, solid and fluid mechanics.
The journal publishes papers on analytical, numerical and experimental methods. Papers that address fundamentally new topics in wave phenomena or develop wave propagation methods for solving direct and inverse problems are of interest to the journal.