We present an implementation of WZjj production via vector-boson fusion in the POWHEG BOX, a public tool for the matching of next-to-leading order QCD calculations with multi-purpose parton-shower generators. We provide phenomenological results for electroweak WZjj production with fully leptonic decays at the LHC in realistic setups and discuss theoretical uncertainties associated with the simulation. We find that beyond the leading-order approximation the dependence on the unphysical factorization and renormalization scales is mild. The two tagging jets are furthermore very stable against parton-shower effects. However, considerable sensitivities to the shower Monte-Carlo program used are observed for central-jet veto observables.