Optical chirality has been recently suggested to complement the physically relevant conserved quantities of the well-known Maxwell's equations. This time-even pseudoscalar is expected to provide further insight in polarization phenomena of electrodynamics such as spectroscopy of chiral molecules. Previously, the corresponding continuity equation was stated for homogeneous lossless media only. We extend the underlying theory to arbitrary setups and analyse piecewise-constant material distributions in particular. Our implementation in a Finite Element Method framework is applied to illustrative examples in order to introduce this novel tool for the analysis of time-harmonic simulations of nano-optical devices.