In this thesis the second-order QCD corrections to electroweak production of a Higgs boson in association with two jets through vector boson fusion are considered. This calculation is fully differential in the kinematics of the Higgs boson and of the final state jets. Infrared divergences are regulated using the antenna subtraction method. We detail the implementation of the process in the parton-level Monte Carlo integrator NNLOJET and present inclusive calculations as well as differential distributions for a wide range of observables at different center-of-mass energies.