Modeling the dynamical evolution of the canonical Eulalia family with the center a_c = 2.475 au. This solution used the bulk density ρ = 1.3 g cm⁻³, the surface thermal inertia Γ = 200 SI units, and the variable YORP timescale τ₀ = 1 My. The best match to the data (left panel) is obtained for T = 785 My and v₅ = 32 m s⁻¹ (recall that the characteristic initial ejection velocity is v_ej(D) = v₅ (5 km/D) in meters per second). The vertical dashed line in the left panel shows the C_⋆-value, though our data (black symbols with formal uncertainties) continue formally farther due to finite bin size. This solution provides a significant amount of D ≤ 3 km fragments directly injected to the J3:1 (right panel). Still, they represent ≃22% of the uncorrected population of fragments reaching the J3:1 in total. Both ordinates are calibrated to the adopted family identification without correcting for incompleteness.