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

03/05/2021
by   Igor Tominec, et al.
0

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.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
03/06/2020

A least squares radial basis function finite difference method with improved stability properties

Localized collocation methods based on radial basis functions (RBFs) for...
research
11/11/2022

Local-basis Difference Potentials Method for elliptic PDEs in complex geometry

We develop efficient and high-order accurate finite difference methods f...
research
07/14/2022

Accurate Stabilization Techniques for RBF-FD Meshless Discretizations with Neumann Boundary Conditions

A major obstacle to the application of the standard Radial Basis Functio...
research
08/07/2023

An improved local radial basis function method for solving small-strain elasto-plasticity

Strong-form meshless methods received much attention in recent years and...
research
10/07/2021

Cantilevered, Rectangular Plate Dynamics by Finite Difference Methods

In this technical note, we consider a dynamic linear, cantilevered recta...
research
07/15/2020

An unfitted RBF-FD method in a least-squares setting for elliptic PDEs on complex geometries

Radial basis function generated finite difference (RBF-FD) methods for P...
research
02/23/2023

Walk on Stars: A Grid-Free Monte Carlo Method for PDEs with Neumann Boundary Conditions

Grid-free Monte Carlo methods based on the walk on spheres (WoS) algorit...

Please sign up or login with your details

Forgot password? Click here to reset