Designing scaffolds similar to the structure of trabecular bone requires specialised algorithms. Existing scaffold designs for bone tissue engineering have repeated patterns that do not replicate the random stochastic porous structure of the internal architecture of bones. In this research, the Voronoi tessellation method is applied to create random porous biomimetic structures. A volume mesh created from the shape of a Zygoma fracture acts as a boundary for the generation of random seed points by point spacing to create Voronoi cells and Voronoi diagrams. The Voronoi lattices were obtained by adding strut thickness to the Voronoi diagrams. Gradient Voronoi scaffolds of pore sizes (19.8 µm to 923 µm) similar to the structure of the trabecular bone were designed. A Finite Element Method-based computational fluid dynamics (CFD) simulation was performed on all designed Voronoi scaffolds to predict the pressure drops and permeability of non-Newtonian blood flow behaviour using the power law material model. The predicted permeability (0.33 × 10−9 m2 to 2.17 × 10−9 m2) values of the Voronoi scaffolds from the CFD simulation are comparable with the permeability of scaffolds and bone specimens from other research works.