Numerical study of a cryogen-free vuilleumier type pulse tube cryocooler operating below 10 K