Augmenting the Immersed Boundary Method with Radial Basis Functions (RBFs) for the Modeling of Platelets in Hemodynamic Flows

We present a new computational method by extending the Immersed Boundary (IB) method with a spectrally-accurate geometric model based on Radial Basis Function (RBF) interpolation of the Lagrangian structures. Our specific motivation is the modeling of platelets in hemodynamic flows, though we anticipate that our method will be useful in other applications as well. The efficacy of our new RBF-IB method is shown through a series of numerical experiments. Specifically, we compare our method with the traditional IB method in terms of convergence and accuracy, computational cost, maximum stable time-step size and volume loss. We conclude that the RBF-IB method has advantages over the traditional Immersed Boundary method, and is well-suited for modeling of platelets in hemodynamic flows.Comment: 25 pages, 4 figure

Doctor of Philosophy

dissertationPlatelet aggregation, an important part of the development of blood clots, is a complex process involving both mechanical interaction between platelets and blood, and chemical transport on and o the surfaces of those platelets. Radial Basis Function (RBF) interpolation is a meshfree method for the interpolation of multidimensional scattered data, and therefore well-suited for the development of meshfree numerical methods. This dissertation explores the use of RBF interpolation for the simulation of both the chemistry and mechanics of platelet aggregation. We rst develop a parametric RBF representation for closed platelet surfaces represented by scattered nodes in both two and three dimensions. We compare this new RBF model to Fourier models in terms of computational cost and errors in shape representation. We then augment the Immersed Boundary (IB) method, a method for uid-structure interaction, with our RBF geometric model. We apply the resultant method to a simulation of platelet aggregation, and present comparisons against the traditional IB method. We next consider a two-dimensional problem where platelets are suspended in a stationary fluid, with chemical diusion in the fluid and chemical reaction-diusion on platelet surfaces. To tackle the latter, we propose a new method based on RBF-generated nite dierences (RBF-FD) for solving partial dierential equations (PDEs) on surfaces embedded in 2D domains. To robustly tackle the former, we remove a limitation of the Augmented Forcing method (AFM), a method for solving PDEs on domains containing curved objects, using RBF-based symmetric Hermite interpolation. Next, we extend our RBF-FD method to the numerical solution of PDEs on surfaces embedded in 3D domains, proposing a new method of stabilizing RBF-FD discretizations on surfaces. We perform convergence studies and present applications motivated by biology. We conclude with a summary of the thesis research and present an overview of future research directions, including spectrally-accurate projection methods, an extension of the Regularized Stokeslet method, RBF-FD for variable-coecient diusion, and boundary conditions for RBF-FD

Hybrid finite difference/finite element immersed boundary method

The immersed boundary method is an approach to fluid-structure interaction that uses a Lagrangian description of the structural deformations, stresses, and forces along with an Eulerian description of the momentum, viscosity, and incompressibility of the fluid-structure system. The original immersed boundary methods described immersed elastic structures using systems of flexible fibers, and even now, most immersed boundary methods still require Lagrangian meshes that are finer than the Eulerian grid. This work introduces a coupling scheme for the immersed boundary method to link the Lagrangian and Eulerian variables that facilitates independent spatial discretizations for the structure and background grid. This approach employs a finite element discretization of the structure while retaining a finite difference scheme for the Eulerian variables. We apply this method to benchmark problems involving elastic, rigid, and actively contracting structures, including an idealized model of the left ventricle of the heart. Our tests include cases in which, for a fixed Eulerian grid spacing, coarser Lagrangian structural meshes yield discretization errors that are as much as several orders of magnitude smaller than errors obtained using finer structural meshes. The Lagrangian-Eulerian coupling approach developed in this work enables the effective use of these coarse structural meshes with the immersed boundary method. This work also contrasts two different weak forms of the equations, one of which is demonstrated to be more effective for the coarse structural discretizations facilitated by our coupling approach

Meshfree and Particle Methods in Biomechanics: Prospects and Challenges

The use of meshfree and particle methods in the field of bioengineering and biomechanics has significantly increased. This may be attributed to their unique abilities to overcome most of the inherent limitations of mesh-based methods in dealing with problems involving large deformation and complex geometry that are common in bioengineering and computational biomechanics in particular. This review article is intended to identify, highlight and summarize research works on topics that are of substantial interest in the field of computational biomechanics in which meshfree or particle methods have been employed for analysis, simulation or/and modeling of biological systems such as soft matters, cells, biological soft and hard tissues and organs. We also anticipate that this review will serve as a useful resource and guide to researchers who intend to extend their work into these research areas. This review article includes 333 references

Two Dimensional Compressible Flow Solver for Moving Geometries Using Immersed Boundary Method

RÉSUMÉ La méthode des frontières immergées (IB) a été mise en œuvre avec un certain succès pour différentes applications, y compris les géométries mobiles et stationnaires. Les travaux actuels portent sur l’extension de la méthode IB au traitement du déplacement des géométries pour les applications à écoulements compressible. En effet, pour les frontières mobiles, la mise en œuvre de l’approche IB comprend le changement du type de cellules, solide vers fluide, ce qui complexifie la méthodologie car il y a un manque dans l’historique de la solution numérique pour ces cellules. Afin de résoudre à cette problématique, deux approches sont disponibles dans la littérature. La première approche repose sur l’extrapolation des variables sur les cellules solides adjacentes (cellules fictives). La deuxième approche utilise la reconstruction de la solution toujours dans la région fluide. Un autre aspect de l’implémentation de la méthode IB est la représentation correcte d’une partie ou la totalité de la géométrie, en particulier lorsque l’échelle géométrique est inférieure à celle du maillage. Ce problème a été identifié et résolu pour les géométries stationnaires.----------ABSTRACT Immersed Boundary (IB) methods have been successfully implemented for different appli- cations including moving as well as stationary geometries. Present work focuses on the implementation of IB method for moving geometries for compressible flow cases. IB imple- mentation for moving boundaries includes the conversion of solid cells to fluid cells, which makes the problem a challenging one because of the abnormal values of the derivatives of pressure and velocity for the cells being converted. In order to solve this problem, mainly two different groups of methods are available. The first group of methods relies on extrapolating the variables on the adjacent solid cells(Ghost Cells) and the second group deals only with the interpolation in the fluid region. Another aspect of IB method implementation is the proper identification of the geometry or a part of the geometry, especially when the size of the geometry is smaller than the mesh size. This problem has been identified and solved for stationary geometries

Numerical modelling of additive manufacturing process for stainless steel tension testing samples

Nowadays additive manufacturing (AM) technologies including 3D printing grow rapidly and they are expected to replace conventional subtractive manufacturing technologies to some extents. During a selective laser melting (SLM) process as one of popular AM technologies for metals, large amount of heats is required to melt metal powders, and this leads to distortions and/or shrinkages of additively manufactured parts. It is useful to predict the 3D printed parts to control unwanted distortions and shrinkages before their 3D printing. This study develops a two-phase numerical modelling and simulation process of AM process for 17-4PH stainless steel and it considers the importance of post-processing and the need for calibration to achieve a high-quality printing at the end. By using this proposed AM modelling and simulation process, optimal process parameters, material properties, and topology can be obtained to ensure a part 3D printed successfully