This research examines diapycnal tracer mixing in Eulerian numerical simulations of stratified shear flow using techniques based on Lagrangian analysis. The first technique considers the use of quantitative density-tracer scatter plots, with their evolution providing a method of tracking the mixing of tracers across density ranges. The shapes of these scatter plots place constraints on the possible cross-isopycnal fluxes of a given tracer. For an initial layer of tracer in a typical stratified shear flow, it is observed that scatter plots tend to evolve toward a piecewise linear relationship. The second technique uses a variation of the Winters-D'Asaro-Nakamura effective density diffusivity formulation to predict tracer profile evolution. Individual fluid parcels are followed in the flow and cast to specific density classes to reconstruct a mean tracer profile. Geometric relationships between density and tracer gradients provide a forcing that affects the evolution of the mean tracer profile in smaller localized regions, while the bulk evolution of the profile is dominated by the effective density diffusivity contribution. The distribution of the final tracer profile is controlled primarily by the density diffusivity, the density range affected by mixing, and the tracer content within this density range.