Nicolás A. Barnafi, Argyrios Petras, Luca Gerardo-Giorda
{"title":"Efficient reference configuration formulation in fully nonlinear poroelastic media","authors":"Nicolás A. Barnafi, Argyrios Petras, Luca Gerardo-Giorda","doi":"arxiv-2401.14536","DOIUrl":null,"url":null,"abstract":"Typical pipelines for model geometry generation in computational biomedicine\nstem from images, which are usually considered to be at rest, despite the\nobject being in mechanical equilibrium under several forces. We refer to the\nstress-free geometry computation as the reference configuration problem, and in\nthis work we extend such a formulation to the theory of fully nonlinear\nporoelastic media. The main steps are (i) writing the equations in terms of the\nreference porosity and (ii) defining a time dependent problem whose steady\nstate solution is the reference porosity. This problem can be computationally\nchallenging as it can require several hundreds of iterations to converge, so we\npropose the use of Anderson acceleration to speed up this procedure. Our\nevidence shows that this strategy can reduce the number of iterations up to\n80\\%. In addition, we note that a primal formulation of the nonlinear mass\nconservation equations is not consistent due to the presence of second order\nderivatives of the displacement, which we alleviate through adequate mixed\nformulations. All claims are validated through numerical simulations in both\nidealized and realistic scenarios.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":"181 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-01-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - CS - Numerical Analysis","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2401.14536","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
Typical pipelines for model geometry generation in computational biomedicine
stem from images, which are usually considered to be at rest, despite the
object being in mechanical equilibrium under several forces. We refer to the
stress-free geometry computation as the reference configuration problem, and in
this work we extend such a formulation to the theory of fully nonlinear
poroelastic media. The main steps are (i) writing the equations in terms of the
reference porosity and (ii) defining a time dependent problem whose steady
state solution is the reference porosity. This problem can be computationally
challenging as it can require several hundreds of iterations to converge, so we
propose the use of Anderson acceleration to speed up this procedure. Our
evidence shows that this strategy can reduce the number of iterations up to
80\%. In addition, we note that a primal formulation of the nonlinear mass
conservation equations is not consistent due to the presence of second order
derivatives of the displacement, which we alleviate through adequate mixed
formulations. All claims are validated through numerical simulations in both
idealized and realistic scenarios.