The use of Neural Cellular Automata (NCAs) has been proven useful in simulating the morphogenetic process. NCAs has found its application primarily in the 2D domain, however, this project expands the use of NCAs to the 3D domain by incorporating 3D convolutions into the neural network architecture. This would enable the simulation of complex 3D structures and spatial dependencies in all three dimensions. The project create a CA model that defines global coordination out of local level interactions. The research findings highlight two significant contributions, namely, an expansion of NCA to 3D voxels and the development of a cellular automation technique for producing voxel structures with different levels of complexity. The project explores multiple aspects, including the structural decay of the system after running for additional iterations beyond the training phase. It also showcases the regeneration property of the system.