A finite element dual porosity approach to model deformation-induced fluid flow in cortical bone

Fluid flow through the osteocyte canaliculi network is widely believed to be a main factor that controls bone adaptation. The difficulty of in vivo measurement of this flow within cortical bone makes computational models an appealing alternative to estimate it. We present in this paper a finite element dual porosity macroscopic model that can contribute to evaluate the interstitial fluid flow induced by mechanical loads in large pieces of bone. This computational model allows us to predict the macroscopic fluid flow at both vascular and canalicular porosities in a whole loaded bone. Our results confirm that the general trend in the fluid flow field predicted is similar to the one obtained with previous microscopic models, and that in a whole bone model it is able to estimate the zones with higher bone remodeling.