Variations in threshold voltage (VTH) and charge distributions in noncircular cells of three-dimensional (3D) NAND flash memory are investigated. The electric field difference in the circular and spike regions causes nonuniform trapped electron density in the charge trap layer (CTL) and influences the VTH variation. Such less-trapped electron (LT) regions in CTL between the circular and spike regions exhibit a lower electric field during the program, resulting in a higher current flow through that region. For the two-spike cells, the charge distribution and VTH variation are analyzed at different heights (HSpike) and angles between spikes. These irregular cells decrease VTH as the HSpike increases.