Open Access

ARTICLE

A Meshless LBIE/LRBF Method for Solving the Nonlinear Fisher Equation: Application to Bone Healing

by K. N. Grivas¹, M. G. Vavva¹, E. J. Sellountos², D. I. Fotiadis³, D. Polyzos¹

1 Department of Mechanical Engineering and Aeronautics, University of Patras, Patras, Greece.
2 Instituto Seperior Tecnico CEMAT, Lisbon, Portugal.
3 Department of Materials Science and Engineering, University of Ioannina.
4 Corresponding author. E-mail: polyzos@mech.upatras.gr; Tel: (0030) 2610969442

Computer Modeling in Engineering & Sciences 2015, 105(2), 87-122. https://doi.org/10.3970/cmes.2015.105.087

Abstract

A simple Local Boundary Integral Equation (LBIE) method for solving the Fisher nonlinear transient diffusion equation in two dimensions (2D) is reported. The method utilizes, for its meshless implementation, randomly distributed nodal points in the interior domain and nodal points corresponding to a Boundary Element Method (BEM) mesh, at the global boundary. The interpolation of the interior and boundary potentials is accomplished using a Local Radial Basis Functions (LRBF) scheme. At the nodes of global boundary the potentials and their fluxes are treated as independent variables. On the local boundaries, potential fluxes are avoided by using the Laplacian companion solution. Potential gradients are accurately evaluated without RBFs via a LBIE, valid for gradient of potentials. Nonlinearity is treated using the Newton-Raphson scheme. The accuracy of the proposed methodology is demonstrated through representative numerical examples. Fisher equation is solved here via the LBIE/LRBF method in order to predict cell proliferation during bone healing. Cell concentrations and their gradients are numerically evaluated in a 2D model of fractured bone. The results are demonstrated and discussed.

---

Keywords

---