  We present results for dihadron fragmentation functions from the NJL-jet
model evolved from the model scale to a typical experimental scale of $4
\mathrm{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 \to \pi^+\pi^-$), as
well as a singly favored ($u \to \pi^+K^-$) process involving light quarks.
Finally, we explore the production of $K^+ K^-$ pairs from an initial $u$, $d$
or $s$ quark.
