This paper is devoted to studying the Bessel beam propagation in cylindrical coordinates using the Hankel transform beam propagation method (HT-BPM) and their behavior in different scenarios in the microscale and meter scale of propagation distances. The study compares the results obtained from the HT-BPM with another fast Fourier transform beam propagation method (FFT-BPM) to validate the accuracy and effectiveness of the HT-BPM in modeling Bessel beam propagation. The axial intensity of Bessel beam propagation is analyzed using the HT-BPM. The simulation results obtained from the HT-BPM are compared with those from the FFT-BPM to evaluate the agreement and consistency between the two methods in predicting the axial intensity of Bessel beam propagation. The results show that the HT-BPM is numerically faster than the FFT-BPM by ten times for different sampling points, furthermore, the FFT-BPM accuracy for evaluating the Bessel beam spot radius is 89.9% of the analytical value, while the HT-BPM is 99% relative to analytical value. The prediction of the axial intensity of the Bessel beam has been tested at different types of phase functions and different propagation distances: micrometer, centimeter, and meter scales. The results of the HT-BPM are matched with the analytical and experimental values. Finally, the HT-BPM is tested when the input light source takes different profiles.