Bone is comprised of an enormously hierarchical construction that promotes transportation of necessary fluids and solids, guaranteeing accurate function and growth. Bone remodeling is a combined process of bone creation and destruction. A number of mathematical models have been developed for the balanced and imbalanced bone remodeling. A brief overview regarding mathematical modeling of bone remodeling is provided. The Lattice Boltzmann method (LBM) has widely been implemented in CFD simulations, and it is becoming more suitable in the application of image processing amongst several others. Mainly, the LBM simulates the communication between synthetic particles dispersed in a lattice. Canaliculi and tortuous channels that have more or less roughly circular structure link among oval bodies identified as lacunae, and are vital to the function of bone. As there is a lack of equipment to inspect flow in channels on the order of measure of canaliculi, so the use of computational methods are more advantageous to give perceptivities into the nature of the flows. In this article, the computational fluid dynamics analysis is descried, using the Lattice Boltzmann method, to examine the result of the microscopic surface roughness of the canalicular wall, which is formed by collagen fibrils, on the flow profiles in the pericellular space.