A Chemo-poroelastic Analysis of Mechanically Induced Fluid and Solute Transport in an Osteonal Cortical Bone

Linguistic Analysis Identifies Emergent Biomaterial Fabrication Trends for Orthopaedic Applications

The expanse of publications in tissue engineering (TE) and orthopedic TE (OTE) over the past 20 years presents an opportunity to probe emergent trends in the field to better guide future technologies that can make an impact on musculoskeletal therapies. Leveraging this trove of knowledge, a hierarchical systematic search method and trend analysis using connected network mapping of key terms is developed. Within discrete time intervals, an accelerated publication rate for anatomic orthopedic tissue engineering (AOTE) of osteochondral defects, tendons, menisci, and entheses is identified. Within these growing fields, the top-listed key terms are extracted and stratified into evident categories, such as biomaterials, delivery method, or 3D printing and biofabrication. It is then identified which categories decreased, remained constant, increased, or emerged over time, identifying the specific emergent categories currently driving innovation in orthopedic repair technologies. Together, these data demonstrate a significant convergence of material types and descriptors used across tissue types. From this convergence, design criteria to support future research of anatomic constructs that mimic both the form and function of native tissues are formulated. In summary, this review identifies large-scale trends and predicts new directions in orthopedics that will define future materials and technologies.

A lumped model for long bone behavior based on poroelastic deformation and Darcy flow

The present paper provides a simplified model for compact bone behavior by accounting for bone fluid flow coupled to the elasticity of the porous structure. The lumped model considers the bone material as a layered poroelastic structure and predicts normal pressure versus displacement, i.e, a stress-strain curve. There is a parametric dependency on porosity and permeability but, in addition, on pressure history. Specifically, the pressure impulse (the integral of pressure versus time) plays a key role. This factor is alluded to in several past studies, but not highlighted in a simplified fashion. Based on a global flow balance, bone displacement depends on the fluid flow in a channel according to the classical Darcy model of 1856, and on the rate of change of fluid within the porous solid according to the 1941 classical model of Biot. The present results agree with those of Perrin et al. which, in turn, agree with results of a detailed numerical simulation.

Effects of Osteocyte Shape on Fluid Flow and Fluid Shear Stress of the Loaded Bone

This study was conducted to better understand the specific behavior of the intraosseous fluid flow. We calculated the number and distribution of bone canaliculi around the osteocytes based on the varying shapes of osteocytes. We then used these calculated parameters and other bone microstructure data to estimate the anisotropy permeability of the lacunar-canalicular network. Poroelastic finite element models of the osteon were established, and the influence of the osteocyte shape on the fluid flow properties of osteons under an axial displacement load was analyzed. Two types of boundary conditions (BC) that might occur in physiological environments were considered on the cement line of the osteon. BC1 allows free fluid passage from the outer elastic restraint boundary, and BC2 is impermeable and allows no free fluid passage from outer displacement constrained boundary. They both have the same inner boundary conditions that allow fluid to pass through. Changes in the osteocyte shape altered the maximum value of pressure gradient (PG), pore pressure (PP), fluid velocity (FV), and fluid shear stress (FSS) relative to the reference model (spherical osteocytes). The maximum PG, PP, FV, and FSS in BC2 were nearly 100% larger than those in BC1, respectively. It is found that the BC1 was closer to the real physiological environment. The fluid flow along different directions in the elongated osteocyte model was more evident than that in other models, which may have been due to the large difference in permeability along different directions. Changes in osteocyte shape significantly affect the degrees of anisotropy of fluid flow and porous media of the osteon. The model presented in this study can accurately quantify fluid flow in the lacunar-canalicular network.

Symmetry breaking and effects of nutrient walkway in time-dependent bone remodeling incorporating poroelasticity

Bone is an extraordinary biological material that continuously adapts its hierarchical microstructure to respond to static and dynamic loads for offering optimal mechanical features, in terms of stiffness and toughness, across different scales, from the sub-microscopic constituents within osteons—where the cyclic activity of osteoblasts, osteoclasts, and osteocytes redesigns shape and percentage of mineral crystals and collagen fibers—up to the macroscopic level, with growth and remodeling processes that modify the architecture of both compact and porous bone districts. Despite the intrinsic complexity of the bone mechanobiology, involving coupling phenomena of micro-damage, nutrients supply driven by fluid flowing throughout hierarchical networks, and cells turnover, successful models and numerical algorithms have been presented in the literature to predict, at the macroscale, how bone remodels under mechanical stimuli, a fundamental issue in many medical applications such as optimization of femur prostheses and diagnosis of the risk fracture. Within this framework, one of the most classical strategies employed in the studies is the so-called Stanford’s law, which allows uploading the effect of the time-dependent load-induced stress stimulus into a biomechanical model to guess the bone structure evolution. In the present work, we generalize this approach by introducing the bone poroelasticity, thus incorporating in the model the role of the fluid content that, by driving nutrients and contributing to the removal of wastes of bone tissue cells, synergistically interacts with the classical stress fields to change homeostasis states, local saturation conditions, and reorients the bone density rate, in this way affecting growth and remodeling. Through two paradigmatic example applications, i.e. a cylindrical slice with internal prescribed displacements idealizing a tract of femoral diaphysis pushed out by the pressure exerted by a femur prosthesis and a bone element in a form of a bent beam, it is highlighted that the present model is capable to catch more realistically both the transition between spongy and cortical regions and the expected non-symmetrical evolution of bone tissue density in the medium–long term, unpredictable with the standard approach. A real study case of a femur is also considered at the end in order to show the effectiveness of the proposed remodeling algorithm.