Bone ingrowth into a porous implant is necessary for its long-term fixation. Although attempts have been made to quantify the peri-implant bone growth using finite element (FE) analysis integrated with mechanoregulatory algorithms, bone ingrowth into a porous cellular hip stem has scarcely been investigated. Using a three-dimensional (3D) FE model and mechanobiology-based numerical framework, the objective of this study was to predict the spatial distribution of evolutionary bone ingrowth into an uncemented novel porous hip stem proposed earlier by the authors. A CT-based FE macromodel of the implant-bone structure was developed. The bone material properties were assigned based on CT grey value. Peak musculoskeletal loading conditions, corresponding to level walking and stair climbing, were applied. The geometry of the implant-bone macromodel was divided into multiple submodels. A suitable mapping framework was used to transfer maximum nodal displacements from the FE macromodel to the cut boundaries of the FE submodels. CT grey value-based bone materials properties were assigned to the submodels. Thereafter, the submodels were solved and simulations of bone ingrowth were carried out using mechanoregulatory principle. A gradual increase in the average Young's modulus, from 1200 to 1500 MPa, of the bone tissue layer was observed considering all the submodels. The distal submodel exhibited 82% of bone ingrowth, whereas the proximal submodel experienced 65% bone ingrowth. Equilibrium in the bone ingrowth process was achieved in 7 weeks postoperatively, with a notable amount of bone ingrowth that should lead to biological fixation of the novel hip stem.