The particle formation by supercritical technology is an active research field. Rapid expansion of supercritical fluid (RESS) is a promising new technology for particle formation using supercritical fluids. In this study, numerical modeling of particle formation of TBTPP (5, 10, 15, 20-tetrakis (3, 5-bis (trifluoromethyl) phenyl) porphyrin) through RESS process is investigated. The RESS process contains two pre-expansion (nozzle part) and expansion vessel (supersonic free jet) parts. In the first section, mass, momentum and energy balances in addition to extended generalized Bender equation of state (egB-EoS) as accurate equation of state is used to calculate temperature, pressure, density and velocity of supercritical carbon dioxide (SC-CO2) based on mathematical modeling and computational fluid dynamics (CFD) simulation methods. Both models include friction in the nozzle and heat exchange in the expansion vessel and in the nozzle region. Using calculated pressure and temperature, the solubility of TBTPP in the supercritical CO2 is predicted by Peng-Robinson EoS with Kwak-Mansoori mixing rule and the results are compared with experimental data and a semi-empirical correlation. Supersaturation, homogenous nucleation and critical radius of TBTPP based on classical nucleation theory are calculated. The effect of pre-expansion temperature on supersaturation, nucleation and critical radius parameters is also studied. The results of hydrodynamic and particle formation modeling present the same trend as reported in the literature.