Neutral meson mixing and meson lifetimes are theory-side parametrised in terms four-quark operators which can be determined by calculating weak decay matrix elements using lattice Quantum Chromodynamics.While calculations of meson mixing matrix elements are standard, determinations of lifetimes typically suffer from complications in renormalisation procedures because dimension-6 four-quark operators can mix with operators of lower mass dimension and, moreover, quark-line disconnected diagrams contribute.We present work detailing the idea to use fermionic gradient flow to non-perturbatively renormalise matrix elements describing meson mixing or lifetimes, and combining it with a perturbative calculation to match to the MSbar scheme using the short flow time expansion.