A method for implementing time-integral constitutive equations in commercial CFD packages