Three dimensional dielectric systems are simulated to obtain the relaxational polarization responses. The dielectric system consists of randomly distributed permanent dipoles that fluctuate due to thermal activation in double well potentials between two plane-parallel electrodes. The barrier height between the wells is assumed to be proportional to the dipole length. The electrodes' effect is considered with the method of images. The relaxational phenomena of two dielectric systems, one with fixed and the other with distributed dipole lengths having single and distributed relaxation times, respectively, are calculated with a Monte-Carlo simulation method. The local field at every dipole and consequently the transition probabilities are computed to determine the transient polarization response iteratively. Both dielectric systems are simulated over a wide range of temperatures for interacting and non-interacting cases. The mean relaxation times exhibit a Vogel-Fulcher law for interacting systems and an Arrhenius law for non-interacting systems. Index Terms-relaxational polarization, double well potentials, dipole-dipole interaction, temperature effect, method of images, computational materials science 2 SYSTEM MODEL Dielectric structures are modeled by assuming bidirectional dipoles that are randomly distributed in 3D space between the electrodes. These permanent dipoles fluctuate due to thermal