Modern radiation treatment planning for photons includes full 3D modeling of the adsorbed dose distribution, accurate inclusion of the patient anatomy, and consideration of significant changes in material density and composition. Such efforts are founded in an accurate description of the radiation source and the beam delivery system. Modern fast neutron therapy facilities employ highly penetrating beams and isocentric beam delivery. Treatment planning is largely based on analytic models adapted from photon codes and interaction cross sections normalized to macroscopic attenuation. However, the recent PEREGRINE initiative at Lawrence Livermore Laboratory offers the possibility of fully stochastic modeling if the neutron source can be adequately described. In this article we report neutron source modeling of three high energy facilities. Neutron production is based on the intra-nuclear cascade model of the LAHET code while neutron transport through the beam delivery system is managed by MCNP using cross section libraries extended to 100 MeV neutron energy. PEREGRINE is then used to transport the neutron beam through typical phantoms. The resulting neutron sources are in excellent agreement with the limited experimental information and the measured phantom data are well described by the PEREGRINE transport using the LAHET/MCNP determined neutron sources.