Source author record

Aaron L. Fogelson

Aaron L. Fogelson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis Quantitative Methods Cell Behavior Computational Engineering, Finance, and Science cond-mat.soft math.DG math.DS

Catalog footprint

What is connected

6works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

Immersed boundary simulations of cell-cell interactions in whole blood

We present a new method for the geometric reconstruction of elastic surfaces simulated by the immersed boundary method with the goal of simulating the motion and interactions of cells in whole blood. Our method uses parameter-free radial basis functions for high-order meshless parametric reconstruction of point clouds and the elastic force computations required by the immersed boundary method. This numerical framework allows us to consider the effect of endothelial geometry and red blood cell motion on the motion of platelets. We find red blood cells to be crucial for understanding the motion of platelets, to the point that the geometry of the vessel wall has a negligible effect in the presence of RBCs. We describe certain interactions that force the platelets to remain near the endothelium for extended periods, including a novel platelet motion that can be seen only in 3-dimensional simulations that we term "unicycling." We also observe red blood cell-mediated interactions between platelets and the endothelium for which the platelet has reduced speed. We suggest that these behaviors serve as mechanisms that allow platelets to better maintain vascular integrity.

preprint2020arXiv

A Mathematical Model of Platelet Aggregation in an Extravascular Injury Under Flow

We present the first mathematical model of flow-mediated primary hemostasis in an extravascular injury, which can track the process from initial deposition to occlusion. The model consists of a system of ordinary differential equations (ODE) that describe platelet aggregation (adhesion and cohesion), soluble-agonist-dependent platelet activation, and the flow of blood through the injury. The formation of platelet aggregates increases resistance to flow through the injury, which is modeled using the Stokes-Brinkman equations. Data from analogous experimental (microfluidic flow) and partial differential equation models informed parameter values used in the ODE model description of platelet adhesion, cohesion, and activation. This model predicts injury occlusion under a range of flow and platelet activation conditions. Simulations testing the effects of shear and activation rates resulted in delayed occlusion and aggregate heterogeneity. These results validate our hypothesis that flow-mediated dilution of activating chemical ADP hinders aggregate development. This novel modeling framework can be extended to include more mechanisms of platelet activation as well as the addition of the biochemical reactions of coagulation, resulting in a computationally efficient high throughput screening tool.

preprint2015arXiv

A Radial Basis Function (RBF)-Finite Difference Method for the Simulation of Reaction-Diffusion Equations on Stationary Platelets within the Augmented Forcing Method

We present a computational method for solving the coupled problem of chemical transport in a fluid (blood) with binding/unbinding of the chemical to/from cellular (platelet) surfaces in contact with the fluid, and with transport of the chemical on the cellular surfaces. The overall framework is the Augmented Forcing Point Method (AFM) (\emph{L. Yao and A.L. Fogelson, Simulations of chemical transport and reaction in a suspension of cells I: An augmented forcing point method for the stationary case, IJNMF (2012) 69, 1736-52.}) for solving fluid-phase transport in a region outside of a collection of cells suspended in the fluid. We introduce a novel Radial Basis Function-Finite Difference (RBF-FD) method to solve reaction-diffusion equations on the surface of each of a collection of 2D stationary platelets suspended in blood. Parametric RBFs are used to represent the geometry of the platelets and give accurate geometric information needed for the RBF-FD method. Symmetric Hermite-RBF interpolants are used for enforcing the boundary conditions on the fluid-phase chemical concentration, and their use removes a significant limitation of the original AFM. The efficacy of the new methods are shown through a series of numerical experiments; in particular, second order convergence for the coupled problem is demonstrated.

preprint2015arXiv

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.

preprint2014arXiv

A Radial Basis Function (RBF)-Finite Difference (FD) Method for Diffusion and Reaction-Diffusion Equations on Surfaces

In this paper, we present a method based on Radial Basis Function (RBF)-generated Finite Differences (FD) for numerically solving diffusion and reaction-diffusion equations (PDEs) on closed surfaces embedded in $\mathbb{R}^d$. Our method uses a method-of-lines formulation, in which surface derivatives that appear in the PDEs are approximated locally using RBF interpolation. The method requires only scattered nodes representing the surface and normal vectors at those scattered nodes. All computations use only extrinsic coordinates, thereby avoiding coordinate distortions and singularities. We also present an optimization procedure that allows for the stabilization of the discrete differential operators generated by our RBF-FD method by selecting shape parameters for each stencil that correspond to a global target condition number. We show the convergence of our method on two surfaces for different stencil sizes, and present applications to nonlinear PDEs simulated both on implicit/parametric surfaces and more general surfaces represented by point clouds.

preprint2012arXiv

A Study of Different Modeling Choices For Simulating Platelets Within the Immersed Boundary Method

The Immersed Boundary (IB) method is a widely-used numerical methodology for the simulation of fluid-structure interaction problems. The IB method utilizes an Eulerian discretization for the fluid equations of motion while maintaining a Lagrangian representation of structural objects. Operators are defined for transmitting information (forces and velocities) between these two representations. Most IB simulations represent their structures with piecewise-linear approximations and utilize Hookean spring models to approximate structural forces. Our specific motivation is the modeling of platelets in hemodynamic flows. In this paper, we study two alternative representations - radial basis functions (RBFs) and Fourier-based (trigonometric polynomials and spherical harmonics) representations - for the modeling of platelets in two and three dimensions within the IB framework, and compare our results with the traditional piecewise-linear approximation methodology. For different representative shapes, we examine the geometric modeling errors (position and normal vectors), force computation errors, and computational cost and provide an engineering trade-off strategy for when and why one might select to employ these different representations.