Numerical Simulations on Thermal Turbulent Flows with a Second-Order Closure Model