Supercritical fluid extraction is powerful new technique in separation processes and environmentally is more appropriate than many other methods such as organic solvent extraction. Cinnamaldehyde and Eugenol are two active constituents of cinnamon and therapeutic for different disease especially diabetes and cancer. In this study, Cinnamaldehyde and Eugenol was extracted via (SC-CO 2 ) with methanol as an entrainer and subsequently the results were compared with Soxhlet method as the base case of 100% recovery. The extracted samples were analysed by gas chromatography (GC) and Response surface methodology was employed to optimize the extraction parameters. The effects of temperature, pressure, flow rate of CO 2 and extraction time on the recoveries of Cinnamaldehyde and Eugenol were investigated in the range of 40-80 , 10-30 MPa, 0.5-2.5 ml/min and40-120min, respectively. Cinnamaldehyde and Eugenol recovery results indicated that the data adequately fitted into a second-order polynomial model. R 2 and modified R 2 of the model are 95.26% and 91.11%, for Cinnamaldehyde and 97.93%, 96.12% for eugenol respectively. It was predicted that the optimum extraction conditions within the experimental ranges for cinnamaldehyde would be the pressure of 20.3 MPa, temperature of 68.2 o C, flow rate of1.8 ml/min, and extraction time of 95.7 min with maximum recovery 54.79 w/w%. and the eugenol recovery results indicated that the maximum recovery of 38.42 w/w% was obtained at optimal operating conditions of 42.04 o C, 20.7 MPa, 2.3 ml/min and 118 min (dynamic extraction time). Moreover, in the present study, mathematical modeling for (SC-CO 2 ) extraction of Cinnamaldehyde and Eugenol from cinnamon bark, was performed by two methods. The first model, broken and intact cell (BIC), is a general model on the base of mass transfer equations. first part of extraction was controlled by phase equilibrium and during the process, analysts diffuse on surface of particle and internal diffusion controlled by mass transfer coefficient in solid phase. Model parameters are such as volumetric mass transfer coefficient in fluid phase, volumetric mass transfer coefficient in solid phase and the analyst concentration ratio in intact ad broken cells. Two first parameters calculated by Empirical Correlations and third parameter achieved by experimental data's in order to minimize the absolute average deviations(AAD) between model results and experimental data's. In second model the effects of axial dispersion in fluid phase, film mass transfer resistance, pore diffusivity inside the particles has been investigated. The mass transfer parameters, i.e., effective pore diffusivity, film mass transfer coefficient and axial dispersion, along with the solubility parameter were chosen as the model parameters. The first three parameters were predicted using imperial correlations. Henry law was used to describe the equilibrium state of solid and pore fluid phases. Genetic algorithm was applied to determine the optimal value of Henry coefficient such that to minimize the absolute average deviations (AAD) between the experimental and predicted extraction recoveries. Finally, there was obvious good agreement between two methods of optimization (genetic algorithm and RSM).