Insulator-based DEP (iDEP) has been established as a powerful tool for manipulating particles in microfluidic devices. However, Joule heating may become an issue in iDEP microdevices due to the local amplification of electric field around the insulators. This results in an electrothermal force that can manifest itself in the flow field in the form of circulations, thus affecting the particle motion. We develop herein a transient, 3D, full-scale numerical model to study Joule heating and its effects on the coupled transport of charge, heat, and fluid in an iDEP device with a rectangular constriction microchannel. This model is validated by comparing the simulation results with the experimentally obtained fluid flow patterns and particle images that were reported in our recent works. It identifies a significant difference in the time scales of the electric, temperature, and flow fields in iDEP microdevices. It also predicts the locations of electrothermal flow circulations in different halves of the channel at the upstream and downstream of the constriction.