A simulation study on the Double-layer Beam Shaping Assembly (DBSA) system has been carried out. This study used fast neutron beam resulting from reactions of 30 MeV protons with beryllium target. The MCNPX code was utilized to design the DBSA and the phantom as well as to calculate neutron flux on the phantom. The distribution of epithermal neutron flux and gamma in the DBSA and phantom were computed using the PHITS code. The spectrum of radiation beams generated by the DBSA shows the characteristics that the typical epithermal neutron flux of 1.0 x109 n/(cm2.s), the ratio of epithermal to the thermal and fast neutron flux of 344 and 85, respectively and the ratio of gamma dose to the epithermal neutron flux of 1.82 x 10-13 Gy.cm2. The test of epithermal neutron beams irradiation on the water phantom shows that epithermal neutrons are thermalized and penetrate the phantom up to 12 cm in depth. The maximum value of neutron flux is 1.1 x 109 n/(cm2.s) at a depth of 2 cm in phantom.