Mathematical modeling and numerical solutions for functionally dependent bone remodeling