Numerical Simulation of Electrical Properties in Cortical and Trabecular Bone: A Simplified Model
Abstract
The electrical properties of biological tissues depend on their composition and microstructure and determine their response to applied electric fields. In bone tissue, these properties are closely related to its composition, porosity, and microstructural organization. However, simplified models that link bone microstructure (trabecular and cortical) to its effective electrical properties while accounting for the electrical anisotropy of cortical bone remain scarce. In this work, a multiscale approach was developed, combining finite element simulations with Bruggeman-based homogenization models. Trabecular bone was modeled using microCT derived geometries, and its effective conductivity was analyzed as a function of BV/TV and free water content. Cortical bone was represented by a simplified three-dimensional network of Haversian canals (axial = z) and Volkmann canals (transverse), incorporating their directional organization through an anisotropic Bruggeman formulation. The models captured the conductivity range reported experimentally for trabecular bone (approximately 25-200 mS/m) and reproduced the anisotropic behavior observed in the cortical FEM simulations, with higher conductivity in the axial direction. In the generated cortical geometries, 65.5% of the total porosity was associated with axial Haversian canals. At the macroscale, the tibia model with lower trabecular and cortical BV/TV and reduced cortical thickness showed an approximately 38-39% reduction in impedance relative to the reference model. These results indicate that microstructural variations in bone tissue are reflected in its effective electrical properties and in the macroscopic impedance of the model, providing a physical and computational basis for future investigations into the sensitivity of electrical measurements to changes in bone microstructure.
Related
- Convergence properties of nonlinear RLC electrical networks
- Numerical Simulation of the Time-Fractional Nonlinear Schrödinger Equation and the Time-Fractional Schrödinger–Poisson System
- Numerical simulation of biological capsules based on interface-fitted polygonal meshes
- Space-Time Parallel Direct Numerical Simulation of Rayleigh-Bénard Convection
- Generating functions and efficient numerical simulation of aggregation-fragmentation kinetics
BibTeX
@misc{cervantes2026numerical,
title = {{Numerical Simulation of Electrical Properties in Cortical and Trabecular Bone: A Simplified Model}},
author = {María José Cervantes and Catalina A. Cely-Ortíz and C. Manuel Carlevaro and Ramiro M. Irastorza},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.31584},
eprint = {2609.31584},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.31584},
}
RIS
TY - UNPB
TI - Numerical Simulation of Electrical Properties in Cortical and Trabecular Bone: A Simplified Model
AU - María José Cervantes
AU - Catalina A. Cely-Ortíz
AU - C. Manuel Carlevaro
AU - Ramiro M. Irastorza
PY - 2026
DA - 2026/09/25
PB - arXiv
DO - 10.48550/arxiv.2609.31584
UR - https://arxiv.org/abs/2609.31584
AB - The electrical properties of biological tissues depend on their composition and microstructure and determine their response to applied electric fields. In bone tissue, these properties are closely related to its composition, porosity, and microstructural organization. However, simplified models that link bone microstructure (trabecular and cortical) to its effective electrical properties while accounting for the electrical anisotropy of cortical bone remain scarce. In this work, a multiscale approach was developed, combining finite element simulations with Bruggeman-based homogenization models. Trabecular bone was modeled using microCT derived geometries, and its effective conductivity was analyzed as a function of BV/TV and free water content. Cortical bone was represented by a simplified three-dimensional network of Haversian canals (axial = z) and Volkmann canals (transverse), incorporating their directional organization through an anisotropic Bruggeman formulation. The models captured the conductivity range reported experimentally for trabecular bone (approximately 25-200 mS/m) and reproduced the anisotropic behavior observed in the cortical FEM simulations, with higher conductivity in the axial direction. In the generated cortical geometries, 65.5% of the total porosity was associated with axial Haversian canals. At the macroscale, the tibia model with lower trabecular and cortical BV/TV and reduced cortical thickness showed an approximately 38-39% reduction in impedance relative to the reference model. These results indicate that microstructural variations in bone tissue are reflected in its effective electrical properties and in the macroscopic impedance of the model, providing a physical and computational basis for future investigations into the sensitivity of electrical measurements to changes in bone microstructure.
ER -