We present results for dihadron fragmentation functions from the NJL-jet model evolved from the model scale to a typical experimental scale of 4 GeV 2 . The numerical method used in this evolution is discussed in detail. The effect of evolution on the shapes of the dihadron fragmentation functions is discussed for a doubly favored process (u → π + π − ), as well as a singly favored (u → π + K − ) process involving light quarks. Finally, we explore the production of K + K − pairs from an initial u, d or s quark.