We employ the chirally rotated Schrodinger functional ($chi$SF) to study two-point fermion bilinear correlation functions used in the determination of $Z_{A,V,S,P,T}$ on a series of well-tuned ensembles. The gauge configurations, which span renormalisation scales from 4 to 70~GeV, are generated with $N_{rm f}=3$ massless flavors and Schrodinger Functional (SF) boundary conditions. Valence quarks are computed with $chi$SF boundary conditions. We show preliminary results on the tuning of the $chi$SF Symanzik coefficient $z_f$ and the scaling of the axial current normalization $Z_{rm A}$. Moreover we carry out a detailed comparison with the expectations from one-loop perturbation theory. Finally we outline how automatically $mathrm{O}(a)$-improved $B_{rm K}$ matrix elements, including BSM contributions, can be computed in a $chi$SF renormalization scheme.