Skip to Main content Skip to Navigation
Journal articles

An unfitted radial basis function generated finite difference method applied to thoracic diaphragm simulations

Igor Tominec 1 Pierre-Frédéric Villard 2 Elisabeth Larsson 1 Víctor Bayona 3 Nicola Cacciani 4 
2 TANGRAM - Recalage visuel avec des modèles physiquement réalistes
Inria Nancy - Grand Est, UL - Université de Lorraine, LORIA - ALGO - Department of Algorithms, Computation, Image and Geometry, CNRS - Centre National de la Recherche Scientifique
Abstract : The thoracic diaphragm is the muscle that drives the respiratory cycle of a human being. Using a system of partial differential equations (PDEs) that models linear elasticity we compute displacements and stresses in a two-dimensional cross section of the diaphragm in its contracted state. The boundary data consists of a mix of displacement and traction conditions. If these are imposed as they are, and the conditions are not compatible, this leads to reduced smoothness of the solution. Therefore, the boundary data is first smoothed using the least-squares radial basis function generated finite difference (RBF-FD) framework. Then the boundary conditions are reformulated as a Robin boundary condition with smooth coefficients. The same framework is also used to approximate the boundary curve of the diaphragm cross section based on data obtained from a slice of a computed tomography (CT) scan. To solve the PDE we employ the unfitted least-squares RBF-FD method. This makes it easier to handle the geometry of the diaphragm, which is thin and non-convex. We show numerically that our solution converges with high-order towards a finite element solution evaluated on a fine grid. Through this simplified numerical model we also gain an insight into the challenges associated with the diaphragm geometry and the boundary conditions before approaching a more complex three-dimensional model.
Document type :
Journal articles
Complete list of metadata
Contributor : Pierre-Frédéric Villard Connect in order to contact the contributor
Submitted on : Tuesday, August 9, 2022 - 10:57:06 AM
Last modification on : Monday, August 22, 2022 - 2:06:44 PM


Files produced by the author(s)



Igor Tominec, Pierre-Frédéric Villard, Elisabeth Larsson, Víctor Bayona, Nicola Cacciani. An unfitted radial basis function generated finite difference method applied to thoracic diaphragm simulations. Journal of Computational Physics, Elsevier, 2022, pp.111496. ⟨10.1016/⟩. ⟨hal-03748102⟩



Record views


Files downloads