Numerical algorithm for triphasic model of charged and hydrated soft tissues

This paper devises an efficient numerical algorithm for solving a two-dimensional triphasic model of charged and hydrated soft tissue by using the radial basis functions. The proposed numerical method is applied directly as a simple meshless collocation algorithm to approximate the solution of the governing system of continuity, momentum, and constitutive equations for the triphasic model. Since there is no requirement on meshing, the method can easily be applied to solve problems under complicated geometry. For verification, numerical simulations of stress, strain, and fluid flow patterns for a plane strain and an axisymmetric mechano-electrochemical coupling model with real synovial joint are given respectively. Classical domain decomposition technique is also combined successfully with the proposed method for solving large scale problems with numerical verification given in solving the axisymmetric case.