Understanding plastic deformation in metallic materials necessitates modeling of dislocation multiplication. A well-known intragranular dislocation source, the Frank-Read (FR) source, plays an important role in size-dependent plasticity in metals. Despite numerous studies in this topic, a systematic investigation of multiple aspects of the FR source in refractory metals is lacking. In this work, we simulate dislocation multiplication from an FR source in six refractory metals with body-centered cubic structure –– Cr, Mo, Nb, Ta, V, and W –– via a phase-field dislocation dynamics (PFDD) method. Two main characteristics of the FR source, including the critical shear stress and critical dislocation configuration, are investigated. For the same metal and same initial dislocation segment, selected PFDD results are benchmarked against atomistic simulations. Our results shed light on developing more accurate continuum models for FR source activation.