Theoretical modeling of acoustofluidic systems faces extreme challenges as the thickness of the thermoviscous boundary layer is very small compared to the microscale fluid dimensions. The classical pressure acoustic model overcomes these difficulties and is extensively used in simulating three-dimensional (3D) or large two-dimensional (2D) acoustofluidic systems. However, this model cannot be applied to thermoviscous acoustofluidics, as it does not consider energy conservation. Modeling thermoviscous acoustofluidic systems is, therefore, difficult and restricted to small 2D systems only. Here, we have developed a thermo-pressure acoustic model that can effectively simulate thermoviscous acoustofluidic systems. The model has been validated with the full model by performing numerical simulations for a small 2D acoustofluidic system for which capturing the acoustic boundary layer effect is feasible using the full model. After successful validation, we demonstrate that the thermo-pressure acoustic model can also be applied to studying 3D acoustofluidic systems.