Computational Challenges in QCD Thermodynamics