Abstract: We propose a simulation method, thermonar dynamics (TD) simulation. Thermons in ideal gases can be simulated directly by TD method. A collision model between thermons is established and thermodynamic properties of thermon gas are obtained. By TD simulation method, equilibrium and non-equilibrium thermon gas of argon are simulated. Simulated pressures of equilibrium systems agree with the equation of state of thermon gas predicted by thermomass theory. Simulated thermal conductivities are consistent with experimental and theoretical data. Feasibility of the proposed thermonar dynamics method is shown.